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
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’argumentmax.lagégal à60. - Ajustez le modèle aux données générées avec
sarima(). En plus des argumentsp,detqdans votre appel àsarima(), spécifiezP,D,QetS(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 = ___)