我们现在在学R软件来简单的分析

我先分享一下我的程序
我的模型是
fit=arima(IPI,order=c(1,1,1),seasonal=list(order=c(0,0,2),period=12))
所做的预测
IPI1=ts(IPI)
dbnp.prev.dyn=ts(start=1000,end=1120)
dbnp.se.dyn=ts(start=1000,end=1120)
for (i in 999:1120){
dbnp.mod.temp=arima(window(IPI1,start=1,end=i),order=c(1,1,1),include.mean=F)
     dbnp.prev.dyn[i-999]=predict(dbnp.mod.temp,n.ahead=1)$pred
    dbnp.se.dyn[i-999]=predict(dbnp.mod.temp, n.ahead=1)$se
}
IPI2=ts(dbnp.prev.dyn,start=c(2002,4),end=c(2012,4),frequency=12)
dbnp.se.dyn2=ts(dbnp.se.dyn,start=c(2002,4),end=c(2012,4),frequency=12)
par(mfrow=c(1,1))
plot(window(IPI,start=c(2002,4),end=c(2012,4)))
lines(IPI2,col="red")
lines(IPI2+1.00*dbnp.se.dyn2,col="blue")
lines(IPI2-1.00*dbnp.se.dyn2,col="blue")

上面方法是老师讲的
不知道有没有大神能用更加简洁的方法
对ARIMA模型进行样本内的预测和样本外的预测
想知道详细的过程
先在这里谢过了