CommencezCommencez gratuitement

Réchauffement planétaire

Maintenant que vous avez un peu d'expérience avec l'ajustement d'un modèle ARIMA sur des données simulées, votre prochaine tâche consiste à appliquer vos compétences à des données réelles.

Les données globtemp (du paquet astsa) correspondent aux écarts annuels de température mondiale jusqu'en 2015. Dans cet exercice, vous utiliserez des techniques éprouvées pour ajuster un modèle ARIMA à ces données. Le graphique des données montre un comportement de marche aléatoire, ce qui suggère de travailler avec les données différenciées. Les données différenciées, diff(globtemp), sont aussi représentées.

Après avoir tracé la FAC et la FACP échantillonnales des données différenciées diff(globtemp), vous pouvez conclure que :

  1. La FAC et la FACP décroissent progressivement, ce qui suggère un modèle ARIMA(1,1,1).
  2. La FAC s'annule au décalage 2 et la FACP décroît progressivement, ce qui suggère un modèle ARIMA(0,1,2).
  3. La FAC décroît progressivement et la FACP s'annule au décalage 3, ce qui suggère un modèle ARIMA(3,1,0). Bien que ce modèle s'ajuste raisonnablement bien, c'est le pire des trois (vous pouvez le vérifier) parce qu'il utilise trop de paramètres pour des autocorrélations aussi faibles.

Après avoir ajusté les deux premiers modèles, comparez l'AIC et le BIC pour choisir le modèle privilégié.

Cette activité fait partie du cours

Modèles ARIMA en R

Voir le cours

Instructions de l’exercice

  • Tracez la FAC et la FACP échantillonnales des données différenciées diff(globtemp) afin de constater que deux modèles semblent raisonnables : un ARIMA(1,1,1) et un ARIMA(0,1,2).
  • Utilisez sarima() pour ajuster un modèle ARIMA(1,1,1) à globtemp. Tous les paramètres sont-ils significatifs?
  • Faites un autre appel à sarima() pour ajuster un modèle ARIMA(0,1,2) à globtemp. Tous les paramètres sont-ils significatifs? Quel modèle est le meilleur?

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

# 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?

Modifier et exécuter le code