Aan de slagGa gratis aan de slag

Opwarming van de aarde

Nu je wat ervaring hebt met het fitten van een ARIMA-model op gesimuleerde data, is je volgende taak om je vaardigheden toe te passen op echte data.

De data in globtemp (uit astsa) zijn de jaarlijkse wereldwijde temperatuurafwijkingen tot en met 2015. In deze oefening gebruik je beproefde technieken om een ARIMA-model op de data te fitten. Een grafiek van de data laat random walk-gedrag zien, wat suggereert dat je met de verschilde data moet werken. De verschilde data diff(globtemp) zijn ook geplot.

Na het plotten van de steekproef-ACF en -PACF van de verschilde data diff(globtemp), kun je stellen dat ofwel

  1. De ACF en de PACF beide langzaam uitdoven, wat wijst op een ARIMA(1,1,1)-model.
  2. De ACF afkapt bij lag 2 en de PACF langzaam uitdooft, wat wijst op een ARIMA(0,1,2)-model.
  3. De ACF langzaam uitdooft en de PACF afkapt bij lag 3, wat wijst op een ARIMA(3,1,0)-model. Hoewel dit model redelijk past, is het de slechtste van de drie (je kunt het controleren) omdat het te veel parameters gebruikt voor zulke kleine autocorrelaties.

Na het fitten van de eerste twee modellen, controleer je de AIC en BIC om het voorkeursmodel te kiezen.

Deze oefening maakt deel uit van de cursus

ARIMA-modellen in R

Cursus bekijken

Oefeninstructies

  • Plot de steekproef-ACF en -PACF van de verschilde data, diff(globtemp), en ontdek dat 2 modellen redelijk lijken: een ARIMA(1,1,1) en een ARIMA(0,1,2).
  • Gebruik sarima() om een ARIMA(1,1,1)-model op globtemp te fitten. Zijn alle parameters significant?
  • Gebruik nog een aanroep van sarima() om een ARIMA(0,1,2)-model op globtemp te fitten. Zijn alle parameters significant? Welk model is beter?

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

# Plot the sample P/ACF pair of the differenced data 


# Fit an ARIMA(1,1,1) model to globtemp


# Fit an ARIMA(0,1,2) model to globtemp. Which model is better?

Code bewerken en uitvoeren