美联储模型策略
该宏观时机策略比较股票市场盈利收益率与10年期美国国债收益率。当股票收益率更高时持有股票ETF;当国债收益率更高时转入现金。每月对收益率差进行回归并预测下个月的值,从而减少噪声切换。
每个月末算法利用最近RegressionMonths个月的数据预测下一月的收益差。如果预测值为正则买入股票ETF,否则持有现金代理。只有当预测穿越零点时才调整仓位,降低换手率。
细节
- 入场条件:
- 月末对
(EarningsYield - BondYield)的RegressionMonths个观察值做回归并预测下一期。 - 如果预测为正且订单金额 ≥
MinTradeUsd,买入股票ETF。
- 月末对
- 多空方向:只做多股票或现金。
- 出场条件:当预测的收益差为负时卖出股票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()