Automatische ARIMA-modellen voor niet-seizoensgebonden tijdreeksen
In de video heb je geleerd dat de functie auto.arima() een passend autoregressive integrated moving average (ARIMA)-model kiest voor een tijdreeks, net zoals de functie ets() dat doet voor ETS-modellen. De functie summary() kan extra inzichten geven:
> # p = 2, d = 1, p = 2
> summary(fit)
Series: usnetelec
ARIMA(2,1,2) with drift
...
In deze oefening kies je automatisch een ARIMA-model voor de vooraf ingeladen austa-reeks, die het jaarlijkse aantal internationale bezoekers aan Australië van 1980–2015 bevat. Vervolgens controleer je de residuen (onthoud dat een p-waarde groter dan 0,05 aangeeft dat de data op white noise lijkt) en maak je enkele voorspellingen. Afgezien van de modelleringsfunctie is dit identiek aan wat je deed met ETS-voorspellen.
Deze oefening maakt deel uit van de cursus
Voorspellen in R
Oefeninstructies
- Pas een automatisch ARIMA-model toe op de
austa-reeks met de nieuw geïntroduceerde functie. Sla dit op alsfit. - Gebruik de juiste functie om te controleren of de residuen van het resulterende model op white noise lijken. Ken
TRUEtoe (als de residuen op white noise lijken) ofFALSE(als dat niet zo is) aanresidualsok. - Pas
summary()toe op het model om de geschatte coëfficiënten te zien. - Op basis van de resultaten met
summary(): wat is de AICc-waarde tot op twee decimalen? Hoeveel verschillen zijn er gebruikt? Ken deze respectievelijk toe aanAICcend. - Plot ten slotte, met de pipe-operator, voorspellingen voor de volgende 10 perioden vanuit het gekozen model.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Fit an automatic ARIMA model to the austa series
fit <- ___
# Check that the residuals look like white noise
___
residualsok <- ___
# Summarize the model
___
# Find the AICc value and the number of differences used
AICc <- ___
d <- ___
# Plot forecasts of fit
fit %>% forecast(h = ___) %>% ___()