Modèles ARIMA automatiques pour des séries chronologiques non saisonnières
Dans la vidéo, vous avez vu que la fonction auto.arima() choisit un modèle autorégressif intégré à moyenne mobile (ARIMA) approprié à partir d'une série chronologique, tout comme la fonction ets() le fait pour les modèles ETS. La fonction summary() peut fournir des précisions supplémentaires :
> # p = 2, d = 1, p = 2
> summary(fit)
Series: usnetelec
ARIMA(2,1,2) with drift
...
Dans cet exercice, vous allez choisir automatiquement un modèle ARIMA pour la série austa déjà chargée, qui contient le nombre annuel de visiteurs internationaux en Australie de 1980 à 2015. Vous vérifierez ensuite les résidus (rappelez‑vous qu'une valeur p supérieure à 0,05 indique que les données ressemblent à du bruit blanc) et produirez des prévisions. À part la fonction de modélisation, c'est identique à ce que vous avez fait avec les prévisions ETS.
Cette activité fait partie du cours
Prévision en R
Instructions de l’exercice
- Ajustez automatiquement un modèle ARIMA à la série
austaà l'aide de la nouvelle fonction. Enregistrez le résultat dansfit. - Utilisez la fonction appropriée pour vérifier que les résidus du modèle obtenu ressemblent à du bruit blanc. Assignez
TRUE(si les résidus ressemblent à du bruit blanc) ouFALSE(si ce n'est pas le cas) àresidualsok. - Appliquez
summary()au modèle pour voir les coefficients ajustés. - D'après les résultats de
summary(), quelle est la valeur AICc à deux décimales? Combien de différenciations ont été utilisées? Assignez ces valeurs àAICcetd, respectivement. - Enfin, au moyen de l'opérateur pipe, tracez les prévisions des 10 prochaines périodes à partir du modèle choisi.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# Fit an automatic ARIMA model to the austa series
fit <- ___
# Check that the residuals look like white noise
___
residualsok <- ___
# Summarize the model
___
# Find the AICc value and the number of differences used
AICc <- ___
d <- ___
# Plot forecasts of fit
fit %>% forecast(h = ___) %>% ___()