Automatische ARIMA-modellen voor seizoensreeksen
Zoals je in de video hebt geleerd, werkt de functie auto.arima() ook met seizoensdata. Let op: als je lambda = 0 instelt in de functie auto.arima() — wat neerkomt op een logtransformatie — dan wordt het model gefit op de getransformeerde data en worden de voorspellingen teruggetransformeerd naar de oorspronkelijke schaal.
Na het toepassen van summary() op dit soort gefitte modellen kun je een uitvoer zien zoals hieronder, wat overeenkomt met \((p,d,q)(P,D,Q)[m]\):
ARIMA(0,1,4)(0,1,1)[12]
In deze oefening ga je met deze functies het vooraf ingeladen h02-bestand modelleren en voorspellen. Deze gegevens bevatten de maandelijkse verkoop van corticosteroïde medicijnen in Australië.
Deze oefening maakt deel uit van de cursus
Voorspellen in R
Oefeninstructies
- Gebruik de standaard plotfunctie om de gelogde
h02-data te plotten en te controleren of de variantie stabiel is. - Fit een seizoens-ARIMA-model op de
h02-reeks metlambda = 0. Sla dit op infit. - Vat het gefitte model samen met de juiste methode.
- Welke niveaus van differencing zijn in het model gebruikt? Ken de hoeveelheid lag-1-differencing toe aan
den de seizoensdifferencing aanD. - Plot voorspellingen voor de komende 2 jaar met het gefitte model. Stel
hhierop in.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Check that the logged h02 data have stable variance
h02 %>% ___ %>% ___
# Fit a seasonal ARIMA model to h02 with lambda = 0
fit <- ___
# Summarize the fitted model
___
# Record the amount of lag-1 differencing and seasonal differencing used
d <- ___
D <- ___
# Plot 2-year forecasts
fit %>% ___ %>% ___