Een gemengd seizoensmodel fitten
Zuiver seizoensafhankelijke patronen, zoals eerder in dit hoofdstuk besproken, komen relatief weinig voor. De meeste seizoensreeksen hebben een gemengde afhankelijkheid: slechts een deel van de variatie wordt verklaard door seizoenspatronen.
Onthoud dat het volledige seizoensmodel wordt genoteerd als SARIMA(p,d,q)x(P,D,Q)S, waarbij hoofdletters de seizoensordes aangeven.
Zoals eerder, vraagt deze oefening je om het steekproef-P/ACF-paar te vergelijken met de werkelijke waarden voor gesimuleerde seizoensdata en een model te fitten op de data met sarima(). Deze keer zijn de gesimuleerde data afkomstig van een gemengd seizoensmodel, SARIMA(0,0,1)x(0,0,1)12. De grafieken tonen drie jaar aan data, evenals de model-ACF en -PACF. Let op dat, in tegenstelling tot het zuiver seizoensmodel, er correlaties zijn op zowel de niet-seizoensmatige vertragingen als de seizoensmatige vertragingen.
Zoals altijd is het pakket astsa al geladen. De gegenereerde data staan in x.
Deze oefening maakt deel uit van de cursus
ARIMA-modellen in R
Oefeninstructies
- Plot de steekproef-ACF en -PACF van de gegenereerde data tot lag 60 (
max.lag = 60) en vergelijk ze met de echte waarden. - Fit het model op de gegenereerde data (
x) metsarima(). Zorg er, net als in de vorige oefening, voor dat je de extra seizoensargumenten meegeeft in jesarima()-opdracht.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Plot sample P/ACF pair to lag 60 and compare to actual
# Fit the seasonal model to x