Een MA(1)-model passen
In deze oefening hebben we data gegenereerd uit een MA(1)-model, $$X_t = W_t - .8 W_{t-1} ,$$ x <- arima.sim(model = list(order = c(0, 0, 1), ma = -.8), n = 100). Bekijk de gesimuleerde data en de steekproef-ACF en -PACF om de orde te bepalen op basis van de tabel uit de eerste oefening. Pas daarna het model.
Onthoud dat bij pure MA(q)-modellen de theoretische ACF afkapt bij vertraging q, terwijl de PACF geleidelijk uitdooft.
Deze oefening maakt deel uit van de cursus
ARIMA-modellen in R
Oefeninstructies
- Het pakket astsa is voor je geladen. Er zijn 100 MA(1)-waarnemingen vooraf ingeladen als
x. - Gebruik
plot()om de gegenereerde data inxte plotten. - Plot de steekproef-ACF en -PACF met
acf2()uit hetastsa-pakket. - Gebruik
sarima()uitastsaom een MA(1) te passen op de eerder gegenereerde data. Bekijk de t-tabel en vergelijk de schattingen met de echte waarden.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# astsa is preloaded
# Plot x
# Plot the sample P/ACF of x
# Fit an MA(1) to the data and examine the t-table