CommencezCommencez gratuitement

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

Voir le cours

Instructions de l’exercice

  • Utilisez diff() pour différencier les données (d_birth). Utilisez acf2() 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 dans dd_birth. Utilisez de nouveau acf2() 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

Modifier et exécuter le code