Gesimuleerde ARIMA
Voordat je echte tijdreeksgegevens gaat analyseren, is het goed om met een iets ingewikkelder model te oefenen.
Hier hebben we 250 observaties gegenereerd uit het ARIMA(2,1,0)-model met drift, gegeven door $$Y_t = 1 + 1.5 Y_{t-1} - .75 Y_{t-2} + W_t\,$$ waarbij \(Y_t = \nabla X_t = X_{t} - X_{t-1}\).
Je gaat de bekende technieken gebruiken om een model op de gegevens te fitten.
Het pakket astsa is vooraf geladen en de gegenereerde gegevens staan in x. De reeks x en de gedifferentieerde reeks y <- diff(x) zijn al geplot.
Deze oefening maakt deel uit van de cursus
ARIMA-modellen in R
Oefeninstructies
- Plot de steekproef-ACF en -PACF met
acf2()van de gedifferentieerde gegevensdiff(x)om een model te bepalen. - Fit een ARIMA(2,1,0)-model met
sarima()op de gegenereerde gegevens. Bekijk de t-tabel en andere uitvoerinformatie om de model-fit te beoordelen.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Plot sample P/ACF of differenced data and determine model
# Estimate parameters and examine output