arima: Wie kann ich eine ARIMA-Zeitreihe bekommen?

Question

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) 

Es gibt 329 Beobachtungen in lill Objekt. Wie kann ich die letzte Beobachtung 328 statt der 330 Beobachtung vorhersagen? Vielen Dank.

Answer 1

Sie müssen nicht predict zur Vorhersage der beobachteten Daten aufrufen. Sie können:

fitted_values <- lill - m7$residuals 

Dies ist das angepasste ARIMA-Modell. Um den 328. Wert zu kontrollieren, tun

fitted_values[328] 

ich nicht Ihre Daten haben, so dass ich R-internen Datensatz LakeHuron als Spielzeug Demonstration.

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)