Aan de slagGa gratis aan de slag

ARIMA - plug-and-play

Zoals je in de video zag, heet een tijdreeks ARIMA(\(p,d,q\)) als de gedifferentieerde reeks (van orde \(d\)) ARMA(\(p,q\)) is.

Om gevoel te krijgen voor hoe het model werkt, analyseer je gesimuleerde data uit het geïntegreerde model $$ Y_t = .9 Y_{t-1} + W_t\, $$ waarbij \(Y_t = \nabla X_t = X_t - X_{t-1}\). In dit geval is het model een ARIMA(1,1,0), omdat de gedifferentieerde data een autoregressie van orde 1 vormen.

De gesimuleerde tijdreeks staat in x en is in R gegenereerd als
x <- arima.sim(model = list(order = c(1, 1, 0), ar = .9), n = 200).

Je plot eerst de gegenereerde data en de steekproef-ACF en -PACF van deze data om te zien hoe geïntegreerde data zich gedragen. Daarna differentieer je de data om ze stationair te maken. Je plot de gedifferentieerde data en de bijbehorende steekproef-ACF en -PACF om te zien wat differentiëren oplevert.

Zoals eerder is het pakket astsa alvast voor je geladen. Data van een ARIMA(1,1,0) met AR-parameter .9 zijn opgeslagen in object x.

Deze oefening maakt deel uit van de cursus

ARIMA-modellen in R

Cursus bekijken

Oefeninstructies

  • Plot de gegenereerde data.
  • Gebruik acf2() uit astsa om het P/ACF-paar (steekproef) voor de gegenereerde data te plotten.
  • Plot de gedifferentieerde data.
  • Gebruik nog een aanroep van acf2() om het P/ACF-paar (steekproef) voor de gedifferentieerde data te bekijken. Let erop hoe ze een AR(1)-model voor de gedifferentieerde data suggereren.

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

# Plot x


# Plot the P/ACF pair of x


# Plot the differenced data


# Plot the P/ACF pair of the differenced data

Code bewerken en uitvoeren