Aan de slagGa gratis aan de slag

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

Cursus bekijken

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 met lambda = 0. Sla dit op in fit.
  • 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 d en de seizoensdifferencing aan D.
  • Plot voorspellingen voor de komende 2 jaar met het gefitte model. Stel h hierop 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 %>% ___ %>% ___
Code bewerken en uitvoeren