【发布时间】:2013-04-23 18:50:45
【问题描述】:
您好,我正在绘制一个 xts 对象:
我想将 xlim 延长到 2 月 30 日,因为我想添加来自 Arima 模型的预测:
par(mfrow=c(1,1))
pred <- predict(try2, n.ahead = 1,se.fit=T)
lim_sup <- pred$pred + 1.96 * pred$se
lim_inf <- pred$pred - 1.96 * pred$se
plot(heat["1996-02-01 00:00/1996-02-16 04:00"],type="l",main="Cross Validation")
points(exp(pred$pred),col="blue",lwd=2)
points(exp(lim_sup),lwd=2,col="green")
points(exp(lim_inf),lwd=2,col="green")
但我看不到情节中的任何线条。你能给我一些建议吗?
热时间序列是这样的
HC.f
1996-02-01 00:00:00 1437.000
1996-02-01 01:00:00 1441.600
1996-02-01 02:00:00 1489.300
1996-02-01 03:00:00 1501.300
1996-02-01 04:00:00 1568.400
1996-02-01 05:00:00 1629.400
而一小时步长的预测:
pred$pred
Time Series:
Start = 3002
End = 3002
Frequency = 1
[1] 7.480588
【问题讨论】:
-
能否提供
try2和heat让我们试一试?您对heat的索引可能有问题。 -
2 月 30 日不存在,恐怕。