Zacznij terazZacznij za darmo

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

Zobacz kurs

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 argument max.lag na 60.
  • Dopasuj model do wygenerowanych danych za pomocą sarima(). Oprócz argumentów p, d i q podaj również P, D, Q oraz S (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 = ___)
Edytuj i uruchom kod