创建真/假信号和相关函数

时间:2016-01-19 21:05:48

标签: r quantstrat

以下量子系统会产生以下错误:

Error in `colnames<-`(`*tmp*`, value = "hardStop") : 
  attempt to set 'colnames' on an object with less than two dimensions

这是我的“hardStop”功能/指标的问题。它是一个if / else函数,它接受当前基于ATR的停止的时间戳(由函数/指示符“stopLimit”生成并存储为mktdata中的列)。

当我为“buyTrigger”列运行以下代码时,它会产生令人满意的TRUE / FALSE信号。

applyIndicatorSignals(strategy=strategy.st, portfolio=portfolio.st, mktdata, sigcol="buyTrigger")

2015-10-21          0
2015-10-22          1
2015-10-23          0
2015-10-26          0
2015-10-27          0
2015-10-28          1
2015-10-29          0
2015-10-30          0

所以这是故障的功能。

hardStop <- function(strategy, portfolio, mktdata, sigcol) {
  if (sigcol==TRUE) {
    stopHere <- mktdata$loss.stopLimit
  }
  else {
    stopHere <- NA
  }
  stopHere <- na.locf(stopHere)
  out <- stopHere
  colnames(out) <- "hardStop"
  return(out)

我知道为什么会出现这个错误?

下面的完整代码,这是一个冗长的帖子。遗憾。

require(quantstrat)

initDate="1990-01-01"
from="2011-07-07"
to="2015-12-12"
options(width=70)

options("getSymbols.warning4.0"=FALSE)
rm(list=ls(.blotter), envir=.blotter)

currency('USD')
Sys.setenv(TZ="UTC")

symbols <- "SPY"

suppressMessages(getSymbols(symbols, from=from, to=to, src="yahoo", adjust=TRUE))
stock(symbols, currency="USD", multiplier=1)

#trade sizing and initial equity settings
tradeSize <- 100000
initEq <- tradeSize*length(symbols)

strategy.st <- portfolio.st <- account.st <- "NewerTry"
rm.strat(portfolio.st)
rm.strat(strategy.st)
initPortf(portfolio.st, symbols=symbols, initDate=initDate, currency='USD')
initAcct(account.st, portfolios=portfolio.st, initDate=initDate, currency='USD',initEq=initEq)
initOrders(portfolio.st, initDate=initDate)
strategy(strategy.st, store=TRUE)

#parameters
daysHigh <- 20

pctATR <- .01
period <- 14

#create functions
periodHigh <- function(HLC, n) {
  high <- (runMax(Hi(HLC), n=n))
  out <- lag(high, 1)
  colnames(out) <- "periodHighest"
  return(out)
}

stopLimit <- function(HLC, n, maType, pctATR) {
  ATR <- ATR(HLC, n=n, maType=maType)
  ATR <- ATR$atr
  close <- Cl(HLC)
  atrStopProfit <- close+(ATR*300*pctATR)
  atrStopLoss <- close-(ATR*100*pctATR)
  atrStop <- cbind(atrStopProfit, atrStopLoss)
  colnames(atrStop) <- c("profit", "loss")
  return(atrStop)
}

hardStop <- function(strategy, portfolio, mktdata, sigcol) {
  if (sigcol==TRUE) {
    stopHere <- mktdata$loss.stopLimit
  }
  else {
    stopHere <- NA
  }
  stopHere <- na.locf(stopHere)
  out <- stopHere
  colnames(out) <- "hardStop"
  return(out)
}

#indicators and signals
add.indicator(strategy.st, name="periodHigh",
              arguments=list(HLC=quote(HLC(mktdata)), n=daysHigh),
              label="periodHighest")

add.indicator(strategy.st, name="stopLimit",
              arguments=list(HLC=quote(HLC(mktdata)), n=period, wilder=TRUE, pctATR=pctATR),
              label="stopLimit")

add.signal(strategy.st, name="sigCrossover",
           arguments=list(columns=c("High", "periodHighest"), relationship="gt", cross=TRUE),
           label="buyTrigger")

applyIndicatorSignals(strategy=strategy.st, portfolio=portfolio.st, mktdata, sigcol="buyTrigger")

add.indicator(strategy.st, name="hardStop",
              arguments=list(strategy=strategy.st, portfolio=portfolio.st, mktdata=quote(HLC(mktdata)),
                             sigcol="buyTrigger"),
              label="hardStop")

#rules
add.rule(strategy.st, name="ruleSignal", 
         arguments=list(sigcol="buyTrigger", sigval=TRUE, ordertype="market", 
                        orderside="long", replace=FALSE, prefer="Open", 
                        orderqty=10000, orderset="orders"), 
         type="enter", path.dep=TRUE,
         label="newEntry")

#run
t1 <- Sys.time()
out <- applyStrategy(strategy=strategy.st,portfolios=portfolio.st)
t2 <- Sys.time()
print(t2-t1)

#set up analytics
updatePortf(portfolio.st)

#performance analytics
chart.Posn(portfolio.st)

1 个答案:

答案 0 :(得分:1)

为此,我最终在函数内部运行了一个函数。这样我就不必处理将列数据带回另一个函数的问题。虽然我仍然觉得paste()或grep()应该以某种方式工作。

stopLimit <- function(HLC, n, maType, pctATR) {
  ATR <- ATR(HLC, n=n, maType=maType)
  ATR <- ATR$atr
  close <- Cl(HLC)
  atrStopProfit <- close+(ATR*300*pctATR)
  atrStopLoss <- close-(ATR*100*pctATR)
  atrStop <- cbind(atrStopProfit, atrStopLoss)
  colnames(atrStop) <- c("profit", "loss")
  return(atrStop)
}

hardStop <- function(HLC, m, n, maType, pctATR) {
  dayHigh <- Hi(HLC)[,1]
  runHigh <- periodHigh(HLC, m)
  runHigh <- replace(runHigh, is.na(runHigh), 0)
  potsLoss <- stopLimit(HLC, n, maType, pctATR)$loss
  potsLoss <- ifelse(dayHigh > runHigh, potsLoss, NA)
  potsLoss <- na.locf(potsLoss)
  out <- potsLoss
  potsWin <- stopLimit(HLC, n, maType, pctATR)$profit
  potsWin <- ifelse(dayHigh > runHigh, potsWin, NA)
  print(potsWin <- na.locf(potsWin))
  pots <- cbind(potsLoss, potsWin)
  colnames(pots) <- c("potsloss", "potswin")
  return(pots)
}

add.indicator(strategy.st, name="hardStop",
              arguments=list(HLC=quote(HLC(mktdata)), m=daysHigh,
                             n=period, wilder=TRUE, pctATR=pctATR),
              label="potsSpot")

add.signal(strategy.st, name="sigCrossover",
           arguments=list(columns=c("Close", "potswin.potsSpot"), relationship="lt", cross=FALSE),
           label="sellTriggerProfit")