Quantstrat 不适用于 ADX

Quantstrat won't work with ADX

我是 运行 基于 R quantstrat 的脚本,摘自 Backtesting Strategies with R)。有用。直到我添加 ADX 作为指标和信号。如果是这样,那么我会收到以下错误:

Error in `colnames<-`(`*tmp*`, value = "ADXsig") : 
length of 'dimnames' [2] not equal to array extent

here in quantstrattrader and here in r.789695.n4.nabble.com 讨论的假定解决方案是将 ADX add.indicator 代码从 x=quote(Cl(mktdata)) 更改为 quote(Cl(mktdata)[,1]) .它不起作用,可能是因为 ADX 使用 HLC 而不是 Cl,并且 HLC 引用三列而不是一列。要清楚:将 HLC=quote(HLC(mktdata)) 更改为 quote(HLC(mktdata)[,1]) 无效。

下面是完整的工作代码,下面是单独的 ADX 代码:

# INSTALL PACKAGES
# install.packages("devtools")
# require(devtools)
# install_github("braverock/FinancialInstrument")
# install_github("joshuaulrich/xts") 
# install_github("braverock/blotter")
# install.packages("quantstrat", repos="http://R-Forge.R-project.org")
# install_github("braverock/PerformanceAnalytics")

# LIBRARIES
library(quantstrat)   

# INITIAL SETUP
Sys.setenv(TZ = "EST")
currency('USD') 
start_date <- "2015-01-01"
end_date <- "2016-12-31"
init_equity <- 1e4 # ,000
adjustment <- FALSE

# GET DATA
basic_symbols <- function() {symbols <- c("SPY")}
symbols <- basic_symbols()
getSymbols(Symbols = symbols, src = "google", index.class = "POSIXct", 
           from = start_date, to = end_date, adjust = adjustment)
stock(symbols,currency = "USD", multiplier = 1)

# DEFINE STRATEGY/PORTFOLIO/ACCOUNT NAMES
portfolio.st <- "Port.Luxor"
account.st <- "Acct.Luxor"
strategy.st <- "Strat.Luxor"

# REMOVE PRIOR STRATEGY/PORTFOLIO, INITIALIZE PORTFOLIO/ACCOUNT/STRATEGY, STORE STRATEGY
rm.strat(portfolio.st)
rm.strat(account.st)
initPortf(name = portfolio.st, symbols = symbols)
initAcct(account.st, portfolios = portfolio.st, initEq = init_equity)
initOrders(portfolio.st)
strategy(strategy.st, store = TRUE)

# INDICATORS & SIGNALS

add.indicator(strategy = strategy.st,
              name = "SMA",
              arguments = list(x = quote(Cl(mktdata)), 
                               n = 10),
              label = "nFast")

add.indicator(strategy = strategy.st, 
              name = "SMA", 
              arguments = list(x = quote(Cl(mktdata)), 
                               n = 30), 
              label = "nSlow")

add.signal(strategy = strategy.st,
           name="sigCrossover",
           arguments = list(columns = c("nFast", "nSlow"),
                            relationship = "gte"),
           label = "long")

add.signal(strategy = strategy.st,
           name="sigCrossover",
           arguments = list(columns = c("nFast", "nSlow"),
                            relationship = "lt"),
           label = "short")

# TRADING RULES
add.rule(strategy = strategy.st,
         name = "ruleSignal",
         arguments = list(sigcol = "long",
                          sigval = TRUE,
                          orderqty = 100,
                          ordertype = "stoplimit",
                          orderside = "long", 
                          threshold = 0.0005,
                          prefer = "High", 
                          TxnFees = -10, 
                          replace = FALSE),
         type = "enter",
         label = "EnterLONG")

add.rule(strategy.st, 
         name = "ruleSignal", 
         arguments = list(sigcol = "long", 
                          sigval = TRUE, 
                          orderside = "short", 
                          ordertype = "market", 
                          orderqty = "all", 
                          TxnFees = -10, 
                          replace = TRUE), 
         type = "exit", 
         label = "Exit2LONG")

# APPLY STRATEGY 
applyStrategy(strategy.st, portfolios = portfolio.st)

添加以下代码(用于 ADX 指标和信号)时会产生错误。即使策略未引用 - 它在 mktdata xts 对象中的存在也是产生错误的原因。 add.indicator 本身的存在不会导致错误,而是 add.signal 的存在会导致错误。认为错误可能是由 add.signal 引用 column = "ADX" 引起的,但不知道要引用哪个 ADX,因为 add.indicator 创建了三个 ADX 列。

add.indicator(strategy.st, name="ADX", 
              arguments=list(HLC=quote(HLC(mktdata)), n=14), 
              label="ADX")

add.signal(strategy.st, name = "sigThreshold",
           arguments = list(column = "ADX",
                            threshold = 30,
                            relationship = "gt",      
                            cross = TRUE),            
           label = "ADXsig")

找到解决方案:将 column = "ADX" 更改为 column = "ADX.ADX"。如下:

add.indicator(strategy.st, name="ADX", 
          arguments=list(HLC=quote(HLC(mktdata)), n=14), 
          label="ADX")

add.signal(strategy.st, name = "sigThreshold",
       arguments = list(column = "ADX.ADX",
                        threshold = 30,
                        relationship = "gt",      
                        cross = TRUE),            
       label = "ADXsig")