R plot.xts错误栏

时间:2016-09-16 17:33:05

标签: r xts timeserieschart

尝试使用plot.xts

在时间序列图上放置误差线
> myTS[1:20]
           [,1]
2013-07-01   29
2013-07-03   24
2013-07-03   16
2013-07-03   16
2013-07-03   12
2013-07-03   12
2013-07-03   16
2013-07-03   21
2013-07-03   21
2013-07-03   16
2013-07-05   12
2013-07-05   12
2013-07-05   12
2013-07-05   12
2013-07-08   16
2013-07-08   23
2013-07-08   16
2013-07-08   12
2013-07-09   16
2013-07-09   12

我已使用myTSquarterly = apply.quarterly(myTS,mean)

汇总了此内容
> myTSquarterly 
               [,1]
2013-09-30 24.50829
2013-12-31 23.79624
2014-03-31 24.15170
2014-06-30 24.57641
2014-09-30 23.71467
2014-12-31 22.99500
2015-03-31 24.50423
2015-06-30 25.19950
2015-09-30 24.76330
2015-12-31 24.65810
2016-03-31 25.35616
2016-06-30 22.71066
2016-07-27 20.63636

我可以使用plot.xts(myTSquarterly)轻松绘图:

enter image description here

我也可以使用apply.quarterly(myTS,sd)

轻松计算标准偏差

我想将这些标准偏差信息作为错误条添加到图中,但我找不到一个方法吗?

1 个答案:

答案 0 :(得分:2)

如果您打算使用xts进行绘图,您可能需要考虑使用xts的开发版本,该版本位于v0.10(我认为CRAN版本仍然较低),其中绘制了{的功能{1}}已得到改进。 。然后检查简单示例plot.xts

您想要做的事情如下:

?plot.xts

enter image description here

(如果您想继续使用较旧的# Want xts v0.10.0 library(devtools) install_github("joshuaulrich/xts") library(quantmod) getSymbols("GOOG") myTS <- GOOG[, 4] myTSquarterly_mean <- apply.quarterly(myTS,mean) myTSquarterly_sd <- apply.quarterly(myTS, sd) c <- 2 plot(myTSquarterly_mean) lines(myTSquarterly_mean - c * myTSquarterly_sd, col = "red") lines(myTSquarterly_mean + c * myTSquarterly_sd, col = "red", pch = 17) points(myTSquarterly_mean + c * myTSquarterly_sd, col = "red", pch = 17) 库在一个图上绘制多条曲线,请考虑将时间序列转换为动物园类型(xts),然后绘制)