2013-04-23 88 views
3

你好我ploting一個XTS對象: enter image description here時間序列圖範圍

而且我想延長XLIM直到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 
+4

2月30日?? – Arun 2013-04-23 08:54:09

+0

您能否提供'try2'和'heat',以便我們試試看?你的'heat'索引可能有問題。 – plannapus 2013-04-23 09:14:24

+0

2月30日不存在,恐怕。 – Gx1sptDTDa 2013-04-23 10:51:38

回答

3

密謀合併之前,您的實際和預測的對象。那麼你將有一個時間序列在整個期望的時間段。

library(xts) 
set.seed(21) 
x <- xts(rnorm(10),Sys.Date()-10:1) 
y <- xts(rnorm(3),Sys.Date()+0:2) 
z <- merge(x,y) 
plot(z[,1], ylim=range(z,na.rm=TRUE)) 
lines(z[,2],col='red')