Aan de slagGa gratis aan de slag

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

Cursus bekijken

Oefeninstructies

  • Het pakket astsa is voor je geladen. Er zijn 100 MA(1)-waarnemingen vooraf ingeladen als x.
  • Gebruik plot() om de gegenereerde data in x te plotten.
  • Plot de steekproef-ACF en -PACF met acf2() uit het astsa-pakket.
  • Gebruik sarima() uit astsa om 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

Code bewerken en uitvoeren