ARMA, kom d'r in
Je hebt inmiddels flink wat ervaring opgedaan met het fitten van ARMA-modellen op data, maar voordat je het viert, doe je nog één (soort van) zelfstandige oefening.
De data in oil zijn WTI-spotprijzen voor ruwe olie, FOB (in dollars per vat), wekelijkse gegevens van 2000 tot 2008. Gebruik je vaardigheden om een ARMA-model te fitten op de rendementen. De wekelijkse ruwe-olieprijzen (oil) zijn alvast voor je geplot. Werk in de hele oefening met de rendementen, die je zelf zult berekenen.
Zoals eerder is het astsa pakket voor je geladen. De data zijn als oil voor je ingeladen en geplot.
Deze oefening maakt deel uit van de cursus
ARIMA-modellen in R
Oefeninstructies
- Bereken de benaderde rendementen van de ruwe-olieprijs met
diff()enlog(). Zet de rendementen inoil_returns. - Plot
oil_returnsen let erop dat er een paar uitschieters zijn vóór 2004. Verzeker jezelf ervan dat de rendementen stationair zijn. - Plot de steekproef-ACF en -PACF van
oil_returnsmetacf2()uit hetastsa-pakket. - Uit het P/ACF-paar blijkt dat de correlaties klein zijn en de rendementen bijna ruis zijn. Maar het kan ook dat zowel de ACF als de PACF langzaam uitdoven. Als dat zo is, dan ligt een ARMA(1,1) voor de hand. Fit dit model op de olierendementen met
sarima(). Past het model goed? Zie je de uitschieters in de residu-plot?
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Calculate approximate oil returns
oil_returns <-
# Plot oil_returns. Notice the outliers.
# Plot the P/ACF pair for oil_returns
# Assuming both P/ACF are tailing, fit a model to oil_returns