arima:如何安装 ARIMA 时间序列?

arima: How can I get fitted ARIMA time series?

m7 = arima(lill,order=c(0,0,1),
          seasonal=list(order=c(1,0,0),period=22),
          xreg=data.frame(lpGDP))

preds = predict(m7,n.ahead = 1, newxreg = 1)

lill 对象中有 329 个观测值。我如何预测最后一次观察 328,而不是 330 观察?谢谢。

您不需要调用 predict 来预测观测数据。你可以这样做:

fitted_values <- lill - m7$residuals

这是拟合的 ARIMA 模型。要检查第 328 个值,请执行

fitted_values[328]

我没有你的数据,所以我使用R的内置数据集LakeHuron作为玩具演示。

fit <- arima(LakeHuron, order = c(2,0,0), xreg = time(LakeHuron) - 1920)
fitted_values <- LakeHuron - fit$residuals
ts.plot(LakeHuron)  ## observed time series (black)
lines(fitted_values, col = 2)  ## fitted time series (red)