Fedモデル戦略
このマクロ・タイミング・システムは、株式市場の益回りと10年物米国債利回りを比較します。株式の方が高い利回りを提供している場合は株式ETFを保有し、債券の利回りが高い場合は現金に移行します。毎月の利回り格差に対する回帰分析が翌月の値を予測し、ノイズによる売買切替を減らします。
毎月末に、アルゴリズムは直近1年のデータを使って翌月の利回りスプレッドを予測します。予測が正であれば株式を買い、そうでなければ現金プロキシを保有します。予測がゼロをまたいだ時のみポジションが変わり、売買回転を最小化します。
詳細
- エントリー条件:
- 月末に、
(EarningsYield - BondYield)の直近RegressionMonths個の観測値で回帰を行い、次の値を予測する。 - 予測がゼロを上回り、かつ注文が
MinTradeUsdを満たす場合に株式ETFを買う。
- 月末に、
- ロング/ショート: 株式または現金のロングのみ。
- エグジット条件: 予測された利回りスプレッドがマイナスに転じたら株式ポジションを決済する。
- ストップ: なし。
- デフォルト値:
Universe– [株式ETF、オプションで現金ETF]。BondYieldSym– 10年物国債利回り系列。EarningsYieldSym– 株式市場の益回り。RegressionMonths= 12。CandleType= 1日。MinTradeUsd– 最低取引金額。
- フィルター:
- カテゴリ: マクロ。
- 方向: ロングのみ。
- 時間軸: 月次。
- リバランス: 毎月。
using System;
using System.Collections.Generic;
using Ecng.Common;
using StockSharp.Algo.Indicators;
using StockSharp.Algo.Strategies;
using StockSharp.BusinessEntities;
using StockSharp.Configuration;
using StockSharp.Messages;
namespace StockSharp.Samples.Strategies;
/// <summary>
/// Fed model strategy that trades the primary instrument when its synthetic earnings yield exceeds a synthetic bond yield benchmark.
/// </summary>
public class FedModelStrategy : Strategy
{
private readonly StrategyParam<string> _security2Id;
private readonly StrategyParam<int> _yieldLength;
private readonly StrategyParam<int> _lookbackPeriod;
private readonly StrategyParam<decimal> _entryThreshold;
private readonly StrategyParam<decimal> _exitThreshold;
private readonly StrategyParam<int> _cooldownBars;
private readonly StrategyParam<decimal> _stopLoss;
private readonly StrategyParam<DataType> _candleType;
private Security _benchmark = null!;
private ExponentialMovingAverage _earningsYield = null!;
private ExponentialMovingAverage _bondYield = null!;
private SimpleMovingAverage _gapAverage = null!;
private StandardDeviation _gapDeviation = null!;
private decimal _latestPrimaryGap;
private decimal _latestBenchmarkGap;
private decimal? _previousZScore;
private bool _primaryUpdated;
private bool _benchmarkUpdated;
private int _cooldownRemaining;
/// <summary>
/// Benchmark security identifier.
/// </summary>
public string Security2Id
{
get => _security2Id.Value;
set => _security2Id.Value = value;
}
/// <summary>
/// Smoothing length for synthetic yields.
/// </summary>
public int YieldLength
{
get => _yieldLength.Value;
set => _yieldLength.Value = value;
}
/// <summary>
/// Lookback period used to normalize the yield gap.
/// </summary>
public int LookbackPeriod
{
get => _lookbackPeriod.Value;
set => _lookbackPeriod.Value = value;
}
/// <summary>
/// Z-score threshold required to open a position.
/// </summary>
public decimal EntryThreshold
{
get => _entryThreshold.Value;
set => _entryThreshold.Value = value;
}
/// <summary>
/// Z-score threshold required to close a position.
/// </summary>
public decimal ExitThreshold
{
get => _exitThreshold.Value;
set => _exitThreshold.Value = value;
}
/// <summary>
/// Closed candles to wait before another position change.
/// </summary>
public int CooldownBars
{
get => _cooldownBars.Value;
set => _cooldownBars.Value = value;
}
/// <summary>
/// Stop loss percentage.
/// </summary>
public decimal StopLoss
{
get => _stopLoss.Value;
set => _stopLoss.Value = value;
}
/// <summary>
/// Type of candles used for processing.
/// </summary>
public DataType CandleType
{
get => _candleType.Value;
set => _candleType.Value = value;
}
/// <summary>
/// Initializes a new instance of the strategy.
/// </summary>
public FedModelStrategy()
{
_security2Id = Param(nameof(Security2Id), Paths.HistoryDefaultSecurity2)
.SetDisplay("Benchmark Security Id", "Identifier of the benchmark security", "General");
_yieldLength = Param(nameof(YieldLength), 12)
.SetRange(2, 80)
.SetDisplay("Yield Length", "Smoothing length for synthetic yields", "Indicators");
_lookbackPeriod = Param(nameof(LookbackPeriod), 24)
.SetRange(5, 120)
.SetDisplay("Lookback Period", "Lookback period used to normalize the yield gap", "Indicators");
_entryThreshold = Param(nameof(EntryThreshold), 1.1m)
.SetRange(0.2m, 5m)
.SetDisplay("Entry Threshold", "Z-score threshold required to open a position", "Signals");
_exitThreshold = Param(nameof(ExitThreshold), 0.25m)
.SetRange(0m, 2m)
.SetDisplay("Exit Threshold", "Z-score threshold required to close a position", "Signals");
_cooldownBars = Param(nameof(CooldownBars), 8)
.SetRange(0, 120)
.SetDisplay("Cooldown Bars", "Closed candles to wait before another position change", "Risk");
_stopLoss = Param(nameof(StopLoss), 2.5m)
.SetRange(0.5m, 10m)
.SetDisplay("Stop Loss %", "Stop loss percentage", "Risk");
_candleType = Param(nameof(CandleType), TimeSpan.FromHours(4).TimeFrame())
.SetDisplay("Candle Type", "Type of candles", "General");
}
/// <inheritdoc />
public override IEnumerable<(Security, DataType)> GetWorkingSecurities()
{
if (Security != null)
yield return (Security, CandleType);
if (!Security2Id.IsEmpty())
yield return (new Security { Id = Security2Id }, CandleType);
}
/// <inheritdoc />
protected override void OnReseted()
{
base.OnReseted();
_benchmark = null!;
_earningsYield = null!;
_bondYield = null!;
_gapAverage = null!;
_gapDeviation = null!;
_latestPrimaryGap = 0m;
_latestBenchmarkGap = 0m;
_previousZScore = null;
_primaryUpdated = false;
_benchmarkUpdated = false;
_cooldownRemaining = 0;
}
/// <inheritdoc />
protected override void OnStarted2(DateTime time)
{
base.OnStarted2(time);
if (Security == null)
throw new InvalidOperationException("Primary security is not specified.");
if (Security2Id.IsEmpty())
throw new InvalidOperationException("Benchmark security identifier is not specified.");
_benchmark = this.LookupById(Security2Id) ?? new Security { Id = Security2Id };
_earningsYield = new ExponentialMovingAverage { Length = YieldLength };
_bondYield = new ExponentialMovingAverage { Length = YieldLength };
_gapAverage = new SimpleMovingAverage { Length = LookbackPeriod };
_gapDeviation = new StandardDeviation { Length = LookbackPeriod };
var primarySubscription = SubscribeCandles(CandleType, security: Security);
var benchmarkSubscription = SubscribeCandles(CandleType, security: _benchmark);
primarySubscription
.Bind(ProcessPrimaryCandle)
.Start();
benchmarkSubscription
.Bind(ProcessBenchmarkCandle)
.Start();
var area = CreateChartArea();
if (area != null)
{
DrawCandles(area, primarySubscription);
DrawCandles(area, benchmarkSubscription);
DrawOwnTrades(area);
}
StartProtection(
new Unit(2, UnitTypes.Percent),
new Unit(StopLoss, UnitTypes.Percent));
}
private void ProcessPrimaryCandle(ICandleMessage candle)
{
if (candle.State != CandleStates.Finished)
return;
_latestPrimaryGap = UpdateYieldGap(_earningsYield, candle);
_primaryUpdated = true;
TryProcessGap(candle.OpenTime);
}
private void ProcessBenchmarkCandle(ICandleMessage candle)
{
if (candle.State != CandleStates.Finished)
return;
_latestBenchmarkGap = UpdateYieldGap(_bondYield, candle);
_benchmarkUpdated = true;
TryProcessGap(candle.OpenTime);
}
private decimal UpdateYieldGap(ExponentialMovingAverage average, ICandleMessage candle)
{
var yieldProxy = CalculateYieldProxy(candle);
return average.Process(yieldProxy, candle.OpenTime, true).ToDecimal();
}
private decimal CalculateYieldProxy(ICandleMessage candle)
{
var priceBase = Math.Max(candle.ClosePrice, 1m);
var range = Math.Max(candle.HighPrice - candle.LowPrice, Security?.PriceStep ?? 1m);
var normalizedRange = range / priceBase;
var closeLocation = (candle.ClosePrice - candle.LowPrice) / range;
return (1m / priceBase * 100m) + closeLocation - normalizedRange;
}
private void TryProcessGap(DateTime time)
{
if (!_primaryUpdated || !_benchmarkUpdated)
return;
_primaryUpdated = false;
_benchmarkUpdated = false;
var gap = _latestPrimaryGap - _latestBenchmarkGap;
var mean = _gapAverage.Process(gap, time, true).ToDecimal();
var deviation = _gapDeviation.Process(gap, time, true).ToDecimal();
if (!_gapAverage.IsFormed || !_gapDeviation.IsFormed || deviation <= 0m)
return;
if (!IsFormedAndOnlineAndAllowTrading())
return;
if (_cooldownRemaining > 0)
_cooldownRemaining--;
var zScore = (gap - mean) / deviation;
var bullishEntry = _previousZScore is decimal previousBullish && previousBullish < EntryThreshold && zScore >= EntryThreshold;
var bearishEntry = _previousZScore is decimal previousBearish && previousBearish > -EntryThreshold && zScore <= -EntryThreshold;
if (_cooldownRemaining == 0 && Position == 0)
{
if (bullishEntry)
{
BuyMarket();
_cooldownRemaining = CooldownBars;
}
else if (bearishEntry)
{
SellMarket();
_cooldownRemaining = CooldownBars;
}
}
else if (Position > 0 && zScore <= ExitThreshold)
{
SellMarket(Position);
_cooldownRemaining = CooldownBars;
}
else if (Position < 0 && zScore >= -ExitThreshold)
{
BuyMarket(Math.Abs(Position));
_cooldownRemaining = CooldownBars;
}
_previousZScore = zScore;
}
}
import clr
clr.AddReference("StockSharp.Messages")
clr.AddReference("StockSharp.Algo")
clr.AddReference("StockSharp.BusinessEntities")
clr.AddReference("StockSharp.Algo.Indicators")
clr.AddReference("StockSharp.Algo.Strategies")
from System import TimeSpan, Math
from StockSharp.Messages import DataType, CandleStates, Unit, UnitTypes
from StockSharp.Algo.Indicators import ExponentialMovingAverage, SimpleMovingAverage, StandardDeviation
from StockSharp.Algo.Strategies import Strategy
from StockSharp.BusinessEntities import Security
from indicator_extensions import *
class fed_model_strategy(Strategy):
"""Fed model strategy that trades the primary instrument when its synthetic earnings yield exceeds a synthetic bond yield benchmark."""
def __init__(self):
super(fed_model_strategy, self).__init__()
self._security2_id = self.Param("Security2Id", "TONUSDT@BNBFT") \
.SetDisplay("Benchmark Security Id", "Identifier of the benchmark security", "General")
self._yield_length = self.Param("YieldLength", 12) \
.SetRange(2, 80) \
.SetDisplay("Yield Length", "Smoothing length for synthetic yields", "Indicators")
self._lookback_period = self.Param("LookbackPeriod", 24) \
.SetRange(5, 120) \
.SetDisplay("Lookback Period", "Lookback period used to normalize the yield gap", "Indicators")
self._entry_threshold = self.Param("EntryThreshold", 1.1) \
.SetRange(0.2, 5.0) \
.SetDisplay("Entry Threshold", "Z-score threshold required to open a position", "Signals")
self._exit_threshold = self.Param("ExitThreshold", 0.25) \
.SetRange(0.0, 2.0) \
.SetDisplay("Exit Threshold", "Z-score threshold required to close a position", "Signals")
self._cooldown_bars = self.Param("CooldownBars", 8) \
.SetRange(0, 120) \
.SetDisplay("Cooldown Bars", "Closed candles to wait before another position change", "Risk")
self._stop_loss = self.Param("StopLoss", 2.5) \
.SetRange(0.5, 10.0) \
.SetDisplay("Stop Loss %", "Stop loss percentage", "Risk")
self._candle_type = self.Param("CandleType", DataType.TimeFrame(TimeSpan.FromHours(4))) \
.SetDisplay("Candle Type", "Type of candles", "General")
self._benchmark = None
self._earnings_yield = None
self._bond_yield = None
self._gap_average = None
self._gap_deviation = None
self._latest_primary_gap = 0.0
self._latest_benchmark_gap = 0.0
self._previous_z_score = None
self._primary_updated = False
self._benchmark_updated = False
self._cooldown_remaining = 0
@property
def candle_type(self):
return self._candle_type.Value
def GetWorkingSecurities(self):
result = []
if self.Security is not None:
result.append((self.Security, self.candle_type))
sec2_id = str(self._security2_id.Value)
if sec2_id:
s = Security()
s.Id = sec2_id
result.append((s, self.candle_type))
return result
def OnReseted(self):
super(fed_model_strategy, self).OnReseted()
self._benchmark = None
self._earnings_yield = None
self._bond_yield = None
self._gap_average = None
self._gap_deviation = None
self._latest_primary_gap = 0.0
self._latest_benchmark_gap = 0.0
self._previous_z_score = None
self._primary_updated = False
self._benchmark_updated = False
self._cooldown_remaining = 0
def OnStarted2(self, time):
super(fed_model_strategy, self).OnStarted2(time)
sec2_id = str(self._security2_id.Value)
if not sec2_id:
raise Exception("Benchmark security identifier is not specified.")
s = Security()
s.Id = sec2_id
self._benchmark = s
yield_len = int(self._yield_length.Value)
lookback = int(self._lookback_period.Value)
self._earnings_yield = ExponentialMovingAverage()
self._earnings_yield.Length = yield_len
self._bond_yield = ExponentialMovingAverage()
self._bond_yield.Length = yield_len
self._gap_average = SimpleMovingAverage()
self._gap_average.Length = lookback
self._gap_deviation = StandardDeviation()
self._gap_deviation.Length = lookback
primary_subscription = self.SubscribeCandles(self.candle_type, True, self.Security)
benchmark_subscription = self.SubscribeCandles(self.candle_type, True, self._benchmark)
primary_subscription.Bind(self.ProcessPrimaryCandle).Start()
benchmark_subscription.Bind(self.ProcessBenchmarkCandle).Start()
area = self.CreateChartArea()
if area is not None:
self.DrawCandles(area, primary_subscription)
self.DrawCandles(area, benchmark_subscription)
self.DrawOwnTrades(area)
self.StartProtection(
Unit(2, UnitTypes.Percent),
Unit(float(self._stop_loss.Value), UnitTypes.Percent)
)
def ProcessPrimaryCandle(self, candle):
if candle.State != CandleStates.Finished:
return
self._latest_primary_gap = self.UpdateYieldGap(self._earnings_yield, candle)
self._primary_updated = True
self.TryProcessGap(candle.OpenTime)
def ProcessBenchmarkCandle(self, candle):
if candle.State != CandleStates.Finished:
return
self._latest_benchmark_gap = self.UpdateYieldGap(self._bond_yield, candle)
self._benchmark_updated = True
self.TryProcessGap(candle.OpenTime)
def UpdateYieldGap(self, average, candle):
yield_proxy = self.CalculateYieldProxy(candle)
result = process_float(average, yield_proxy, candle.OpenTime, True)
return float(result)
def CalculateYieldProxy(self, candle):
price_base = max(float(candle.ClosePrice), 1.0)
price_step = float(self.Security.PriceStep) if self.Security is not None and self.Security.PriceStep is not None else 1.0
range_val = max(float(candle.HighPrice) - float(candle.LowPrice), price_step)
normalized_range = range_val / price_base
close_location = (float(candle.ClosePrice) - float(candle.LowPrice)) / range_val
return (1.0 / price_base * 100.0) + close_location - normalized_range
def TryProcessGap(self, time):
if not self._primary_updated or not self._benchmark_updated:
return
self._primary_updated = False
self._benchmark_updated = False
gap = self._latest_primary_gap - self._latest_benchmark_gap
mean_result = process_float(self._gap_average, gap, time, True)
mean = float(mean_result)
dev_result = process_float(self._gap_deviation, gap, time, True)
deviation = float(dev_result)
if not self._gap_average.IsFormed or not self._gap_deviation.IsFormed or deviation <= 0:
return
if not self.IsFormedAndOnlineAndAllowTrading():
return
if self._cooldown_remaining > 0:
self._cooldown_remaining -= 1
z_score = (gap - mean) / deviation
entry_thresh = float(self._entry_threshold.Value)
exit_thresh = float(self._exit_threshold.Value)
cooldown = int(self._cooldown_bars.Value)
bullish_entry = self._previous_z_score is not None and self._previous_z_score < entry_thresh and z_score >= entry_thresh
bearish_entry = self._previous_z_score is not None and self._previous_z_score > -entry_thresh and z_score <= -entry_thresh
if self._cooldown_remaining == 0 and self.Position == 0:
if bullish_entry:
self.BuyMarket()
self._cooldown_remaining = cooldown
elif bearish_entry:
self.SellMarket()
self._cooldown_remaining = cooldown
elif self.Position > 0 and z_score <= exit_thresh:
self.SellMarket(self.Position)
self._cooldown_remaining = cooldown
elif self.Position < 0 and z_score >= -exit_thresh:
self.BuyMarket(Math.Abs(self.Position))
self._cooldown_remaining = cooldown
self._previous_z_score = z_score
def CreateClone(self):
return fed_model_strategy()