时间序列分析作图

x=arima.sim(list(order=c(1,0,1),ar=0.6721,ma=-0.1467),n=35)

x.fit=arima(x,order=c(1,0,1))

plot(x.fit)

上述程序运行结果为：错误于min(x, na.rm = na.rm): 'type'(list)参数不对

请问大神应该如何修改？或者加载什么packages呢？

tsdiag(x.fit),你是想画这个图吗？还是什么图？

曲散人终 发表于 2014-12-22 09:57
tsdiag(x.fit),你是想画这个图吗？还是什么图？

很抱歉这个不是哦，我想画的图类似下图下图的代码如下：
library(TSA)

data(color)#color包含35个数据

m1.color=arima(color,order=c(1,0,1))

plot(m1.color)#加载TSA包后，会自动预测前置12期，并画出95%的预测极限