Aan de slagGa gratis aan de slag

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

Cursus bekijken

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 in x te plotten en acf2() 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

Code bewerken en uitvoeren