Réchauffement climatique
Maintenant que vous avez de l’expérience avec l’ajustement d’un modèle ARIMA sur des données simulées, votre prochaine mission est d’appliquer ces compétences à des données réelles.
Les données globtemp (du package astsa) correspondent aux écarts annuels de température globale jusqu’en 2015. Dans cet exercice, vous utiliserez des techniques éprouvées pour ajuster un modèle ARIMA aux données. Le graphique des données montre un comportement de marche aléatoire, ce qui suggère de travailler sur les données différenciées. Les données différenciées diff(globtemp) sont également tracées.
Après avoir tracé l’ACF et la PACF échantillonnales des données différenciées diff(globtemp), vous pouvez conclure que :
- L’ACF et la PACF décroissent progressivement, ce qui suggère un modèle ARIMA(1,1,1).
- L’ACF s’annule au retard 2 et la PACF décroît, ce qui suggère un modèle ARIMA(0,1,2).
- L’ACF décroît et la PACF s’annule au retard 3, ce qui suggère un modèle ARIMA(3,1,0). Bien que ce modèle s’ajuste raisonnablement bien, c’est le moins bon des trois (vous pouvez le vérifier), car il utilise trop de paramètres pour des autocorrélations aussi faibles.
Après avoir ajusté les deux premiers modèles, vérifiez l’AIC et le BIC pour choisir le modèle préféré.
Cet exercice fait partie du cours
Modèles ARIMA en R
Instructions
- Tracez l’ACF et la PACF échantillonnales des données différenciées,
diff(globtemp), pour constater que deux modèles sont plausibles : 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 ? - Utilisez 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 cet exemple de code.
# 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?