Fitování modelu MA(1)
V tomto cvičení jsme vygenerovali data 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). Prohlédni si simulovaná data a výběrové ACF a PACF, abys na základě tabulky z prvního cvičení určil/a řád modelu. Potom model nafituj.
Pamatuj, že u čistých modelů MA(q) teoretické ACF odřízne po zpoždění q, zatímco PACF postupně doznívá.
Toto cvičení je součástí kurzu
Modely ARIMA v R
Pokyny k cvičení
- Balíček astsa je předem načten. 100 pozorování z MA(1) je k dispozici jako
x. - Pomocí
plot()vykresli vygenerovaná data zx. - Vykresli výběrové ACF a PACF pomocí
acf2()z balíčkuastsa. - Pomocí
sarima()z balíčkuastsanafituj model MA(1) na dříve vygenerovaná data. Prohlédni si t-tabulku a porovnej odhadnuté hodnoty se skutečnými.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# astsa is preloaded
# Plot x
# Plot the sample P/ACF of x
# Fit an MA(1) to the data and examine the t-table