Aan de slagGa gratis aan de slag

Pas een puur seizoensmodel

Net als bij andere modellen kun je in R seizoensmodellen fitten met het commando sarima() uit het astsa pakket.

Om gevoel te krijgen voor hoe pure seizoensmodellen werken, kun je het beste naar gesimuleerde data kijken. We hebben 250 observaties gegenereerd uit een puur seizoensmodel gegeven door $$X_t = .9 X_{t-12} + W_t + .5 W_{t-12}\,,$$ wat we zouden noteren als een SARMA(P = 1, Q = 1)S = 12. Drie jaar aan data en de model-ACF en -PACF zijn alvast voor je geplot.

Je vergelijkt de steekproef-ACF- en PACF-waarden van de gegenereerde data met de weergegeven ware waarden.

Het astsa-pakket is al voor je geladen en de gegenereerde data staan in x.

Deze oefening maakt deel uit van de cursus

ARIMA-modellen in R

Cursus bekijken

Oefeninstructies

  • Gebruik acf2() om de steekproef-ACF en -PACF van de gegenereerde data tot en met vertraging 60 te plotten en vergelijk met de werkelijke waarden. Om tot vertraging 60 te schatten, zet je het argument max.lag gelijk aan 60.
  • Fit het model op de gegenereerde data met sarima(). Geef naast de argumenten p, d en q in je sarima()-commando ook P, D, Q en S op (let op: R is hoofdlettergevoelig).

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

# Plot sample P/ACF to lag 60 and compare to the true values
acf2(___, max.lag = ___)

# Fit the seasonal model to x
sarima(x, p = 0, d = 0, q = 0, P = ___, D = 0, Q = ___, S = ___)
Code bewerken en uitvoeren