如何在R中繪製預測模型？

Question

我最近將R從2.x升級到3.1，我感覺有些變化我不知道與繪圖模型有關。

airpass_model = arima(log(AirPassengers) , order=c(0,1,1), 
    seasonal=list(order=c(0,1,1), period=12)) 
plot(airpass_model, n1=c(1969, 1), n.ahead=24, pch=19, ylab='Log(Air Passengers)') 

執行第二行後出現以下錯誤。如果有人能解釋爲什麼會很棒？

Error in xy.coords(x, y, xlabel, ylabel, log) : 
    'x' is a list, but does not have components 'x' and 'y' 

Answer 1

最好使用Arima()功能在forecast庫，然後繪製之前調用forecast：

library(forecast) 
airpass_model = arima(log(AirPassengers) , order=c(0,1,1), 
    seasonal=list(order=c(0,1,1), period=12)) 
plot(forecast(airpass_model,h=24), pch=19, ylab='Log(Air Passengers)') 

畢竟，當你談論 「繪製模式」，你可能是在談論繪製ACF/PACF或許多其他信息。

最後，this may be enlightening - 非常推薦。