2015-08-14 77 views
2

我想创建一个信号状态时,我的价差计算比上/下布林带大于或小于差异(),但是我算了一笔账:R/quantmod - BBands之间runsd(EMA)的计算

pair <- c("qqq","iwm") 
start <- "2014-08-08" 
finish <- "2015-08-13" 
stckY <- suppressWarnings(getSymbols(pair[1], from = start, to = finish, auto.assign = FALSE)) 
stckX <- suppressWarnings(getSymbols(pair[2], from = start, to = finish, auto.assign = FALSE)) 

adY <- Ad(stckY) 
adX <- Ad(stckX) 

logY <- log(adY) 
logX <- log(adX) 

spread <- cbind(logY, logX) 
spread <- spread[complete.cases(spread),] #remove NAs 
spread$dailyDiff <- spread[,1] - spread[,2] 

ema <- EMA(spread[,1] - spread[,2], n=20) 
spread$UpBand <- (runSD(ema, n=20) * 2) + ema 
spread$LwBand <- ema - (runSD(ema, n=20) * 2) 

chartSeries(spread$dailyDiff, up.col = "white", 
      theme = chartTheme("black"), line.type = "l") 
addEMA(n = 20, col = "orange") 
addBBands(n = 20, sd = 2, maType = "EMA") 

与chartSeries显示屏上显示的band值不匹配,我无法弄清楚为什么?帮助文件指出,不使用SMA可能会导致“不一致”,所以这可能是问题的根源? chartSeries也使用EMA。

也许有更好的方法去做这件事?我不知道如何使用BBands(单独)来引用上/下带...

+1

这不是我清楚你比较什么,他们是如何不同。 [可重现的例子](http://stackoverflow.com/q/5963269/271616)比含糊的描述要好得多。 –

+0

您评论中的代码非常难以阅读。请将其编辑到您的问题中(使用正确的代码格式),然后删除您的评论。 –

+0

@JoshuaUlrich - 对不起 - 仍在学习本网站 - 我已更新原帖 – trock2000

回答

1

有2个问题:

  1. 你把移动平均的标准偏差时,你应该采取原始系列的标准偏差。
  2. BBands函数在runSD调用中使用sample = FALSE

此复制BBands函数的输出在图表中:

ema <- EMA(spread$dailyDiff, n=20) 
spread$UpBand <- runSD(spread$dailyDiff, n=20, sample=FALSE) * 2 + ema 
spread$LwBand <- ema - runSD(spread$dailyDiff, n=20, sample=FALSE) * 2 
tail(spread) 
#   QQQ.Adjusted IWM.Adjusted dailyDiff  UpBand  LwBand 
# 2015-08-06  4.704563  4.793060 -0.08849663 -0.06442705 -0.1381200 
# 2015-08-07  4.703295  4.786575 -0.08328008 -0.06687188 -0.1322478 
# 2015-08-10  4.714652  4.798267 -0.08361464 -0.06938022 -0.1267023 
# 2015-08-11  4.701752  4.789323 -0.08757113 -0.07421110 -0.1198771 
# 2015-08-12  4.705196  4.787408 -0.08221192 -0.07856667 -0.1126964 
# 2015-08-13  4.703566  4.784320 -0.08075361 -0.08283161 -0.1055975 
tail(BBands(spread$dailyDiff, n=20, maType="EMA")) 
#     dn  mavg   up  pctB 
# 2015-08-06 -0.1381200 -0.10127351 -0.06442705 0.6733800 
# 2015-08-07 -0.1322478 -0.09955985 -0.06687188 0.7490178 
# 2015-08-10 -0.1267023 -0.09804126 -0.06938022 0.7516764 
# 2015-08-11 -0.1198771 -0.09704411 -0.07421110 0.7074403 
# 2015-08-12 -0.1126964 -0.09563152 -0.07856667 0.8931942 
# 2015-08-13 -0.1055975 -0.09421457 -0.08283161 1.0912769 
+0

宾果,谢谢 - 你一直很有帮助,非常感谢 – trock2000