Aan de slagGa gratis aan de slag

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

Cursus bekijken

Oefeninstructies

  • Het pakket astsa is vooraf geladen. x bevat de 200 AR(2)-waarnemingen.
  • Gebruik plot() om de gegenereerde data in x te plotten.
  • Plot het paar van de steekproef-ACF en -PACF met acf2() uit het astsa-pakket.
  • Gebruik sarima() om een AR(2) te fitten op de eerder gegenereerde data in x. 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

Code bewerken en uitvoeren