Een ARMA-model fitten
Je bent nu klaar om het AR-model en het MA-model te combineren tot het ARMA-model. We hebben data gegenereerd uit het ARMA(2,1)-model, $$X_t = X_{t-1} - .9 X_{t-2} + W_t + .8 W_{t-1}, $$ x <- arima.sim(model = list(order = c(2, 0, 1), ar = c(1, -.9), ma = .8), n = 250). Bekijk de gesimuleerde data en het paar van de steekproef-ACF en -PACF om een mogelijk model te bepalen.
Onthoud dat voor ARMA(\(p, q\))-modellen zowel de theoretische ACF als PACF langzaam uitdoven. In dit geval zijn de ordes lastig af te leiden uit de data en is het misschien niet duidelijk of de steekproef-ACF of steekproef-PACF afkapt of juist uitdooft. In dit geval ken je de werkelijke modelordes, dus fit een ARMA(2,1) op de gegenereerde data. Algemene modelleerstrategieën komen later in de cursus aan bod.
Deze oefening maakt deel uit van de cursus
ARIMA-modellen in R
Oefeninstructies
- Het pakket astsa is al geladen. Er staan 250 ARMA(2,1)-waarnemingen in
x. - Gebruik, net als in de vorige oefeningen,
plot()om de gegenereerde data inxte plotten enacf2()om de steekproef-ACF en -PACF te bekijken. - Gebruik
sarima()om een ARMA(2,1) op de gegenereerde data te fitten. Bekijk de t-tabel en vergelijk de schattingen met de werkelijke 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 ARMA(2,1) to the data and examine the t-table