Dopasuj czysty model sezonowy
Podobnie jak inne modele, modele sezonowe można dopasowywać w R za pomocą polecenia sarima() z pakietu astsa.
Aby poczuć, jak działają czyste modele sezonowe, najlepiej zacząć od danych symulowanych. Wygenerowaliśmy 250 obserwacji z czystego modelu sezonowego opisanego wzorem $$X_t = .9 X_{t-12} + W_t + .5 W_{t-12}\,,$$ który oznaczamy jako SARMA(P = 1, Q = 1)S = 12. Na wykresie przedstawiono trzy lata danych oraz teoretyczne ACF i PACF modelu.
Porównasz wartości próbkowe ACF i PACF z wygenerowanych danych z wartościami teoretycznymi.
Pakiet astsa jest wstępnie załadowany, a wygenerowane dane znajdują się w zmiennej x.
To ćwiczenie jest częścią kursu
Modele ARIMA w R
Instrukcje do ćwiczenia
- Użyj funkcji
acf2(), aby wykreślić próbkowe ACF i PACF wygenerowanych danych do opóźnienia 60 i porównaj je z wartościami teoretycznymi. Aby wyznaczyć wartości do opóźnienia 60, ustaw argumentmax.lagna60. - Dopasuj model do wygenerowanych danych za pomocą
sarima(). Oprócz argumentówp,diqpodaj równieżP,D,QorazS(pamiętaj, że R rozróżnia wielkość liter).
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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 = ___)