Een AR(2)-model fitten
Voor deze oefening hebben we data gegenereerd uit het AR(2)-model, $$X_t = 1.5 X_{t-1} - .75 X_{t-2} + W_t,$$ met x <- arima.sim(model = list(order = c(2, 0, 0), ar = c(1.5, -.75)), n = 200). Bekijk de gesimuleerde data en het koppel van de steekproef-ACF en -PACF om de modelorde te bepalen. Fit daarna het model en vergelijk de geschatte parameters met de echte parameters.
Deze oefening maakt deel uit van de cursus
ARIMA-modellen in R
Oefeninstructies
- Het pakket astsa is vooraf geladen.
xbevat de 200 AR(2)-waarnemingen. - Gebruik
plot()om de gegenereerde data inxte plotten. - Plot het paar van de steekproef-ACF en -PACF met
acf2()uit hetastsa-pakket. - Gebruik
sarima()om een AR(2) te fitten op de eerder gegenereerde data inx. 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 AR(2) to the data and examine the t-table