CommencezCommencez gratuitement

Ajuster un modèle purement saisonnier

Comme pour les autres modèles, vous pouvez ajuster des modèles saisonniers en R avec la commande sarima() du paquet astsa.

Pour bien comprendre le fonctionnement des modèles purement saisonniers, il est préférable d'examiner des données simulées. Nous avons généré 250 observations à partir d'un modèle purement saisonnier défini par $$X_t = .9 X_{t-12} + W_t + .5 W_{t-12}\,,$$ que nous noterions SARMA(P = 1, Q = 1)S = 12. Trois années de données ainsi que la FAC et la FACP du modèle sont tracées pour vous.

Vous comparerez les valeurs de la FAC et de la FACP issues de l'échantillon généré aux valeurs exactes affichées.

Le paquet astsa est déjà chargé pour vous et les données générées sont dans x.

Cette activité fait partie du cours

Modèles ARIMA en R

Voir le cours

Instructions de l’exercice

  • Utilisez acf2() pour tracer la FAC et la FACP échantillonnales des données générées jusqu'au délai 60 et comparez-les aux valeurs réelles. Pour estimer jusqu'au délai 60, fixez l'argument max.lag à 60.
  • Ajustez le modèle aux données générées avec sarima(). En plus des arguments p, d et q dans votre commande sarima(), indiquez aussi P, D, Q et S (notez que R est sensible à la casse).

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

# Plot sample P/ACF to lag 60 and compare to the true values
acf2(___, max.lag = ___)

# Fit the seasonal model to x
sarima(x, p = 0, d = 0, q = 0, P = ___, D = 0, Q = ___, S = ___)
Modifier et exécuter le code