Accrual Anomaly
Accrual Anomaly 策略实现了应计项异常因子。在每年五月的第一个交易日重新平衡,做多低应计股票并做空高应计股票。
测试显示年均收益约 12%,在美国股票市场表现最佳。
仓位每年调整一次,不使用日内信号。
详情
- 入场条件:参见实现中的应计项计算。
- 多空方向:双向。
- 退出条件:下一次计划的再平衡。
- 止损:无明确止损逻辑。
- 默认值:
Deciles = 10CandleType = TimeSpan.FromMinutes(5).TimeFrame()
- 过滤器:
- 分类: 基本面
- 方向: 双向
- 指标: 基本面
- 止损: 否
- 复杂度: 中等
- 时间框架: 日线
- 季节性: 是
- 神经网络: 否
- 背离: 否
- 风险等级: 中等
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>
/// Cross-sectional accrual anomaly strategy that trades the primary instrument against a synthetic accrual benchmark derived from the secondary instrument.
/// </summary>
public class AccrualAnomalyStrategy : Strategy
{
private readonly StrategyParam<string> _security2Id;
private readonly StrategyParam<int> _accrualLength;
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 _security2 = null!;
private ExponentialMovingAverage _primaryAccrualAverage = null!;
private ExponentialMovingAverage _secondaryAccrualAverage = null!;
private SimpleMovingAverage _spreadAverage = null!;
private StandardDeviation _spreadDeviation = null!;
private decimal _latestPrimaryAccrual;
private decimal _latestSecondaryAccrual;
private bool _primaryUpdated;
private bool _secondaryUpdated;
private decimal? _previousZScore;
private int _cooldownRemaining;
/// <summary>
/// Secondary security identifier.
/// </summary>
public string Security2Id
{
get => _security2Id.Value;
set => _security2Id.Value = value;
}
/// <summary>
/// Smoothing length for the synthetic accrual proxy.
/// </summary>
public int AccrualLength
{
get => _accrualLength.Value;
set => _accrualLength.Value = value;
}
/// <summary>
/// Lookback period for spread normalization.
/// </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>
/// Candle type used for both instruments.
/// </summary>
public DataType CandleType
{
get => _candleType.Value;
set => _candleType.Value = value;
}
/// <summary>
/// Initializes strategy parameters.
/// </summary>
public AccrualAnomalyStrategy()
{
_security2Id = Param(nameof(Security2Id), Paths.HistoryDefaultSecurity2)
.SetDisplay("Second Security Id", "Identifier of the secondary security", "General");
_accrualLength = Param(nameof(AccrualLength), 6)
.SetRange(2, 30)
.SetDisplay("Accrual Length", "Smoothing length for the synthetic accrual proxy", "Indicators");
_lookbackPeriod = Param(nameof(LookbackPeriod), 24)
.SetRange(10, 120)
.SetDisplay("Lookback Period", "Lookback period for spread normalization", "Indicators");
_entryThreshold = Param(nameof(EntryThreshold), 1.4m)
.SetRange(0.5m, 4m)
.SetDisplay("Entry Threshold", "Z-score threshold required to open a position", "Signals");
_exitThreshold = Param(nameof(ExitThreshold), 0.35m)
.SetRange(0m, 2m)
.SetDisplay("Exit Threshold", "Z-score threshold required to close a position", "Signals");
_cooldownBars = Param(nameof(CooldownBars), 12)
.SetRange(0, 100)
.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", "Candle series for both instruments", "General");
}
/// <inheritdoc />
public override IEnumerable<(Security sec, DataType dt)> 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();
_security2 = null!;
_primaryAccrualAverage = null!;
_secondaryAccrualAverage = null!;
_spreadAverage = null!;
_spreadDeviation = null!;
_latestPrimaryAccrual = 0m;
_latestSecondaryAccrual = 0m;
_primaryUpdated = false;
_secondaryUpdated = false;
_previousZScore = null;
_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("Secondary security identifier is not specified.");
_security2 = this.LookupById(Security2Id) ?? new Security { Id = Security2Id };
_primaryAccrualAverage = new ExponentialMovingAverage { Length = AccrualLength };
_secondaryAccrualAverage = new ExponentialMovingAverage { Length = AccrualLength };
_spreadAverage = new SimpleMovingAverage { Length = LookbackPeriod };
_spreadDeviation = new StandardDeviation { Length = LookbackPeriod };
_cooldownRemaining = 0;
var primarySubscription = SubscribeCandles(CandleType, security: Security);
var secondarySubscription = SubscribeCandles(CandleType, security: _security2);
primarySubscription
.Bind(ProcessPrimaryCandle)
.Start();
secondarySubscription
.Bind(ProcessSecondaryCandle)
.Start();
var area = CreateChartArea();
if (area != null)
{
DrawCandles(area, primarySubscription);
DrawCandles(area, secondarySubscription);
DrawOwnTrades(area);
}
StartProtection(
new Unit(2, UnitTypes.Percent),
new Unit(StopLoss, UnitTypes.Percent));
}
private void ProcessPrimaryCandle(ICandleMessage candle)
{
if (candle.State != CandleStates.Finished)
return;
_latestPrimaryAccrual = UpdateAccrualAverage(_primaryAccrualAverage, candle);
_primaryUpdated = true;
TryProcessSpread(candle.OpenTime);
}
private void ProcessSecondaryCandle(ICandleMessage candle)
{
if (candle.State != CandleStates.Finished)
return;
_latestSecondaryAccrual = UpdateAccrualAverage(_secondaryAccrualAverage, candle);
_secondaryUpdated = true;
TryProcessSpread(candle.OpenTime);
}
private decimal UpdateAccrualAverage(ExponentialMovingAverage average, ICandleMessage candle)
{
var accrualProxy = CalculateAccrualProxy(candle);
return average.Process(accrualProxy, candle.OpenTime, true).ToDecimal();
}
private decimal CalculateAccrualProxy(ICandleMessage candle)
{
var priceBase = Math.Max(candle.OpenPrice, 1m);
var range = Math.Max(candle.HighPrice - candle.LowPrice, Security?.PriceStep ?? 1m);
var body = candle.ClosePrice - candle.OpenPrice;
var bodyRatio = body / priceBase;
var closeLocation = ((candle.ClosePrice - candle.LowPrice) - (candle.HighPrice - candle.ClosePrice)) / range;
var volatilityRatio = range / priceBase;
return (bodyRatio * 18m) + (closeLocation * 0.8m) - (volatilityRatio * 6m);
}
private void TryProcessSpread(DateTime time)
{
if (!_primaryUpdated || !_secondaryUpdated)
return;
_primaryUpdated = false;
_secondaryUpdated = false;
if (!_primaryAccrualAverage.IsFormed || !_secondaryAccrualAverage.IsFormed)
return;
var spread = _latestPrimaryAccrual - _latestSecondaryAccrual;
var mean = _spreadAverage.Process(spread, time, true).ToDecimal();
var deviation = _spreadDeviation.Process(spread, time, true).ToDecimal();
if (!_spreadAverage.IsFormed || !_spreadDeviation.IsFormed || deviation <= 0)
return;
if (!IsFormedAndOnlineAndAllowTrading())
return;
if (_cooldownRemaining > 0)
_cooldownRemaining--;
var zScore = (spread - 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 accrual_anomaly_strategy(Strategy):
"""Cross-sectional accrual anomaly strategy using dual securities."""
def __init__(self):
super(accrual_anomaly_strategy, self).__init__()
self._security2_id = self.Param("Security2Id", "TONUSDT@BNBFT") \
.SetDisplay("Second Security Id", "Identifier of the secondary security", "General")
self._accrual_length = self.Param("AccrualLength", 6) \
.SetRange(2, 30) \
.SetDisplay("Accrual Length", "Smoothing length for the synthetic accrual proxy", "Indicators")
self._lookback_period = self.Param("LookbackPeriod", 24) \
.SetRange(10, 120) \
.SetDisplay("Lookback Period", "Lookback period for spread normalization", "Indicators")
self._entry_threshold = self.Param("EntryThreshold", 1.4) \
.SetRange(0.5, 4.0) \
.SetDisplay("Entry Threshold", "Z-score threshold required to open a position", "Signals")
self._exit_threshold = self.Param("ExitThreshold", 0.35) \
.SetRange(0.0, 2.0) \
.SetDisplay("Exit Threshold", "Z-score threshold required to close a position", "Signals")
self._cooldown_bars = self.Param("CooldownBars", 12) \
.SetRange(0, 100) \
.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", "Candle series for both instruments", "General")
self._security2 = None
self._primary_accrual_avg = None
self._secondary_accrual_avg = None
self._spread_average = None
self._spread_deviation = None
self._latest_primary_accrual = 0.0
self._latest_secondary_accrual = 0.0
self._primary_updated = False
self._secondary_updated = False
self._previous_z_score = None
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(accrual_anomaly_strategy, self).OnReseted()
self._security2 = None
self._primary_accrual_avg = None
self._secondary_accrual_avg = None
self._spread_average = None
self._spread_deviation = None
self._latest_primary_accrual = 0.0
self._latest_secondary_accrual = 0.0
self._primary_updated = False
self._secondary_updated = False
self._previous_z_score = None
self._cooldown_remaining = 0
def OnStarted2(self, time):
super(accrual_anomaly_strategy, self).OnStarted2(time)
sec2_id = str(self._security2_id.Value)
if not sec2_id:
raise Exception("Secondary security identifier is not specified.")
s = Security()
s.Id = sec2_id
self._security2 = s
accrual_len = int(self._accrual_length.Value)
lookback = int(self._lookback_period.Value)
self._primary_accrual_avg = ExponentialMovingAverage()
self._primary_accrual_avg.Length = accrual_len
self._secondary_accrual_avg = ExponentialMovingAverage()
self._secondary_accrual_avg.Length = accrual_len
self._spread_average = SimpleMovingAverage()
self._spread_average.Length = lookback
self._spread_deviation = StandardDeviation()
self._spread_deviation.Length = lookback
self._cooldown_remaining = 0
primary_subscription = self.SubscribeCandles(self.candle_type, True, self.Security)
secondary_subscription = self.SubscribeCandles(self.candle_type, True, self._security2)
primary_subscription.Bind(self.ProcessPrimaryCandle).Start()
secondary_subscription.Bind(self.ProcessSecondaryCandle).Start()
area = self.CreateChartArea()
if area is not None:
self.DrawCandles(area, primary_subscription)
self.DrawCandles(area, secondary_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_accrual = self.UpdateAccrualAverage(self._primary_accrual_avg, candle)
self._primary_updated = True
self.TryProcessSpread(candle.OpenTime)
def ProcessSecondaryCandle(self, candle):
if candle.State != CandleStates.Finished:
return
self._latest_secondary_accrual = self.UpdateAccrualAverage(self._secondary_accrual_avg, candle)
self._secondary_updated = True
self.TryProcessSpread(candle.OpenTime)
def UpdateAccrualAverage(self, average, candle):
accrual_proxy = self.CalculateAccrualProxy(candle)
result = process_float(average, accrual_proxy, candle.OpenTime, True)
return float(result)
def CalculateAccrualProxy(self, candle):
price_base = max(float(candle.OpenPrice), 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)
body = float(candle.ClosePrice) - float(candle.OpenPrice)
body_ratio = body / price_base
close_location = ((float(candle.ClosePrice) - float(candle.LowPrice)) - (float(candle.HighPrice) - float(candle.ClosePrice))) / range_val
volatility_ratio = range_val / price_base
return (body_ratio * 18.0) + (close_location * 0.8) - (volatility_ratio * 6.0)
def TryProcessSpread(self, time):
if not self._primary_updated or not self._secondary_updated:
return
self._primary_updated = False
self._secondary_updated = False
if not self._primary_accrual_avg.IsFormed or not self._secondary_accrual_avg.IsFormed:
return
spread = self._latest_primary_accrual - self._latest_secondary_accrual
mean_result = process_float(self._spread_average, spread, time, True)
mean = float(mean_result)
dev_result = process_float(self._spread_deviation, spread, time, True)
deviation = float(dev_result)
if not self._spread_average.IsFormed or not self._spread_deviation.IsFormed or deviation <= 0:
return
if not self.IsFormedAndOnlineAndAllowTrading():
return
if self._cooldown_remaining > 0:
self._cooldown_remaining -= 1
z_score = (spread - 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 accrual_anomaly_strategy()