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