CommencezCommencez gratuitement

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

Voir le cours

Instructions de l’exercice

  • Ajustez automatiquement un modèle ARIMA à la série austa à l'aide de la nouvelle fonction. Enregistrez le résultat dans fit.
  • 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) ou FALSE (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 à AICc et d, 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 = ___) %>% ___()
Modifier et exécuter le code