Dopasowywanie modelu MA(1)
W tym ćwiczeniu wygenerowaliśmy dane z modelu MA(1): $$X_t = W_t - .8 W_{t-1} ,$$ x <- arima.sim(model = list(order = c(0, 0, 1), ma = -.8), n = 100). Przyjrzyj się symulowanym danym oraz próbkowym wykresom ACF i PACF, aby określić rząd modelu na podstawie tabeli z pierwszego ćwiczenia. Następnie dopasuj model.
Pamiętaj, że dla czystych modeli MA(q) teoretyczna ACF zanika po opóźnieniu q, natomiast PACF wygasa stopniowo.
To ćwiczenie jest częścią kursu
Modele ARIMA w R
Instrukcje do ćwiczenia
- Pakiet astsa jest już załadowany. 100 obserwacji z modelu MA(1) zostało wczytanych jako
x. - Użyj funkcji
plot(), aby wykreślić wygenerowane dane zx. - Wykreśl pary próbkowych ACF i PACF za pomocą funkcji
acf2()z pakietuastsa. - Użyj funkcji
sarima()z pakietuastsa, aby dopasować model MA(1) do wygenerowanych danych. Przeanalizuj tablicę t i porównaj oszacowania z prawdziwymi wartościami.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# astsa is preloaded
# Plot x
# Plot the sample P/ACF of x
# Fit an MA(1) to the data and examine the t-table