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
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 argumentmax.laggelijk aan60. - Fit het model op de gegenereerde data met
sarima(). Geef naast de argumentenp,denqin jesarima()-commando ookP,D,QenSop (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 = ___)