Prognostisering med simulerad ARIMA
Nu när du behärskar anpassning av ARIMA-modeller kan du använda dina kunskaper för prognostisering. Börja med att arbeta med simulerade data.
Vi genererade 120 observationer från en ARIMA(1,1,0)-modell med AR-parametern 0,9. Data finns i y och de första 100 observationerna finns i x. Dessa observationer är redan plottade åt dig. Du ska anpassa en ARIMA(1,1,0)-modell till data i x och verifiera att modellen passar bra. Använd sedan sarima.for() från astsa för att prognostisera data 20 tidsperioder framåt. Därefter jämför du prognoserna med de faktiska data i y.
Grundsyntaxen för prognostisering är sarima.for(data, n.ahead, p, d, q), där n.ahead är ett positivt heltal som anger prognoshorisontens längd. De predikterade värdena och deras standardfel skrivs ut, data plottas i svart och prognoserna visas i rött tillsammans med 2 medelkvadratiska prediktionsfelsgränser som blå streckade linjer.
Paketet astsa är förladdat och data (x) samt differentierade data (diff(x)) är plottade.
Den här övningen är en del av kursen
ARIMA-modeller i R
Övningsinstruktioner
- Plotta stickprovets ACF och PACF för de differentierade data för att bestämma en lämplig modell.
- Använd
sarima()för att anpassa en ARIMA(1,1,0) till data. Granska utdata från dittsarima()-anrop för att utvärdera anpassningen och modelldiagnostiken. - Använd
sarima.for()för att prognostisera data 20 tidsperioder framåt. Jämför resultatet med de faktiska värdena.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# Plot P/ACF pair of differenced data
# Fit model - check t-table and diagnostics
# Forecast the data 20 time periods ahead
sarima.for(x, n.ahead = ___, p = ___, d = ___, q = ___)
lines(y)