Analyse de données - taux de natalité
Vous allez maintenant mettre en pratique vos nouvelles compétences pour ajuster soigneusement un modèle SARIMA à la série chronologique birth du paquet astsa. Les données correspondent au nombre mensuel de naissances vivantes (ajusté), en milliers, aux États-Unis de 1948 à 1979, incluant le baby-boom d'après la Deuxième Guerre mondiale.
Les données birth sont tracées dans votre console R. Notez la tendance à long terme (marche aléatoire) et la composante saisonnière des données.
Cette activité fait partie du cours
Modèles ARIMA en R
Instructions de l’exercice
- Utilisez
diff()pour différencier les données (d_birth). Utilisezacf2()pour afficher l'ACF et la PACF empiriques de ces données jusqu'au décalage 60. Remarquez la persistance saisonnière. - Utilisez un autre appel à
diff()pour effectuer la différence saisonnière des données. Enregistrez le résultat dansdd_birth. Utilisez de nouveauacf2()pour afficher l'ACF et la PACF de ces données, encore jusqu'au décalage 60. Concluez qu'un modèle SARIMA(0,1,1)x(0,1,1)12 semble raisonnable. - Ajustez le modèle SARIMA(0,1,1)x(0,1,1)12. Que se passe-t-il?
- Ajoutez un paramètre AR supplémentaire (non saisonnier,
p = 1) pour tenir compte de la corrélation additionnelle. Le modèle s'ajuste-t-il bien?
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# Plot P/ACF to lag 60 of differenced data
d_birth <- diff(birth)
# Plot P/ACF to lag 60 of seasonal differenced data
dd_birth <- diff(d_birth, lag = 12)
# Fit SARIMA(0,1,1)x(0,1,1)_12. What happens?
# Add AR term and conclude