CommencerCommencer gratuitement

Ajuster un modèle purement saisonnier

Comme pour les autres modèles, vous pouvez ajuster des modèles saisonniers dans R à l’aide de la fonction sarima() du package astsa.

Pour comprendre le fonctionnement des modèles purement saisonniers, il est préférable de partir de données simulées. Nous avons généré 250 observations à partir d’un modèle purement saisonnier donné 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 l’ACF et la PACF du modèle sont tracées pour vous.

Vous comparerez les valeurs d’ACF et de PACF empiriques issues des données générées aux vraies valeurs affichées.

Le package astsa est préchargé et les données générées sont dans x.

Cet exercice fait partie du cours

Modèles ARIMA en R

Afficher le cours

Instructions

  • Utilisez acf2() pour tracer l’ACF et la PACF empiriques des données générées jusqu’au lag 60 et comparez-les aux valeurs réelles. Pour estimer jusqu’au lag 60, définissez l’argument max.lag égal à 60.
  • Ajustez le modèle aux données générées avec sarima(). En plus des arguments p, d et q dans votre appel à sarima(), spécifiez P, D, Q et S (notez que R est sensible à la casse).

Exercice interactif pratique

Essayez cet exercice en complétant cet exemple de code.

# 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