LoslegenKostenlos loslegen

Globale Erwärmung

Nachdem du erste Erfahrungen mit dem Anpassen eines ARIMA-Modells an simulierte Daten gesammelt hast, wendest du deine Fähigkeiten nun auf echte Daten an.

Die Daten in globtemp (aus astsa) sind die jährlichen globalen Temperaturabweichungen bis 2015. In dieser Übung nutzt du gängige Verfahren, um ein ARIMA-Modell an diese Daten anzupassen. Ein Plot der Daten zeigt Random-Walk-Verhalten, was darauf hindeutet, dass du mit den differenzierten Daten arbeiten solltest. Die differenzierten Daten diff(globtemp) sind ebenfalls geplottet.

Nach dem Plotten der Stichproben-ACF und -PACF der differenzierten Daten diff(globtemp) kannst du sagen, dass entweder

  1. ACF und PACF beide langsam abklingen, was auf ein ARIMA(1,1,1)-Modell hindeutet.
  2. die ACF bei Lag 2 abschneidet und die PACF langsam abklingt, was auf ein ARIMA(0,1,2)-Modell hindeutet.
  3. die ACF langsam abklingt und die PACF bei Lag 3 abschneidet, was auf ein ARIMA(3,1,0)-Modell hindeutet. Obwohl dieses Modell ziemlich gut passt, ist es das schlechteste der drei (du kannst es überprüfen), weil es für so kleine Autokorrelationen zu viele Parameter verwendet.

Nachdem du die ersten beiden Modelle geschätzt hast, prüfe AIC und BIC, um das bevorzugte Modell zu wählen.

Diese Übung ist Teil des Kurses

ARIMA-Modelle in R

Kurs anzeigen

Anleitung zur Übung

  • Plotte die Stichproben-ACF und -PACF der differenzierten Daten diff(globtemp), um festzustellen, dass 2 Modelle sinnvoll erscheinen: ein ARIMA(1,1,1) und ein ARIMA(0,1,2).
  • Verwende sarima(), um ein ARIMA(1,1,1)-Modell an globtemp anzupassen. Sind alle Parameter signifikant?
  • Verwende einen weiteren Aufruf von sarima(), um ein ARIMA(0,1,2)-Modell an globtemp anzupassen. Sind alle Parameter signifikant? Welches Modell ist besser?

Interaktive Übung

Vervollständige den Beispielcode, um diese Übung erfolgreich abzuschließen.

# 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 bearbeiten und ausführen