CommencerCommencer gratuitement

Analyse des données – taux de natalité

Vous allez maintenant mettre en pratique vos nouvelles compétences pour ajuster avec soin un modèle SARIMA à la série temporelle birth de astsa. Ces données correspondent aux naissances vivantes mensuelles (ajustées), en milliers, aux États‑Unis de 1948 à 1979, et incluent le baby‑boom d’après‑guerre.

Les données birth sont tracées dans votre console R. Remarquez la tendance de long terme (marche aléatoire) et la composante saisonnière des données.

Cet exercice fait partie du cours

Modèles ARIMA en R

Afficher le cours

Instructions

  • 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 retard 60. Notez la persistance saisonnière.
  • Utilisez un autre appel à diff() pour effectuer la différence saisonnière des données. Enregistrez cela dans dd_birth. Utilisez à nouveau acf2() pour afficher l’ACF et la PACF de ces données, encore jusqu’au retard 60. Concluez qu’un modèle SARIMA(0,1,1)x(0,1,1)12 semble approprié.
  • 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 cet exemple de code.

# 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