Data-analyse - werkloosheid II
Nu ga je verder met het fitten van een SARIMA-model op de maandelijkse Amerikaanse werkloosheidsreeks unemp door te kijken naar de steekproef-ACF en -PACF van de volledig gedifferentieerde reeks.
Let op: de lag-as in de steekproef-P/ACF-grafiek is uitgedrukt in jaren. Dus lag 1, 2, 3, … stelt 1 jaar (12 maanden), 2 jaar (24 maanden), 3 jaar (36 maanden), … voor.
Het pakket astsa is alweer voor je geladen.
Deze oefening maakt deel uit van de cursus
ARIMA-modellen in R
Oefeninstructies
- Differentieer de data volledig (zoals in de vorige oefening) en plot de steekproef-ACF en -PACF van de getransformeerde data tot lag 60 maanden (5 jaar). Houd rekening met het volgende:
- niet-seizoenscomponent: de PACF kapt af bij lag 2 en de ACF loopt langzaam uit.
- seizoenscomponent: de ACF kapt af bij lag 12 en de PACF loopt uit bij lag 12, 24, 36, …
- Stel een model voor en fit het met
sarima(). Controleer de residuen om te bevestigen dat het model goed past.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Plot P/ACF pair of fully differenced data to lag 60
dd_unemp <- diff(diff(unemp), lag = 12)
# Fit an appropriate model