Ajuster un modèle AR(1)
Rappelez-vous que vous utilisez le duo ACF et PACF pour aider à identifier les ordres \(p\) et \(q\) d’un modèle ARMA. Le tableau suivant résume les résultats :
| AR(\(p\)) | MA(\(q\)) | ARMA(\(p,q\)) | |
|---|---|---|---|
| ACF | Décroît progressivement | S’annule après le retard \(q\) |
Décroît progressivement |
| PACF | S’annule après le retard \(p\) |
Décroît progressivement | Décroît progressivement |
Dans cet exercice, vous allez générer des données selon le modèle AR(1), $$X_t = .9 X_{t-1} + W_t,$$ examiner les données simulées ainsi que le duo ACF et PACF empiriques pour déterminer l’ordre. Ensuite, vous ajusterez le modèle et comparerez les paramètres estimés aux vrais paramètres.
Tout au long de ce cours, vous utiliserez sarima() du package astsa pour ajuster facilement des modèles aux données. Cette commande produit un graphique de diagnostic des résidus que vous pouvez ignorer jusqu’à la section sur les diagnostics plus loin dans le chapitre.
Cet exercice fait partie du cours
Modèles ARIMA en R
Instructions
- Le package astsa est préchargé.
- Utilisez la commande préécrite
arima.sim()pour générer 100 observations d’un modèle AR(1) avec un paramètre AR de 0,9. Enregistrez le résultat dansx. - Tracez les données générées avec
plot(). - Tracez les ACF et PACF empiriques à l’aide de la commande
acf2()du packageastsa. - Utilisez
sarima()deastsapour ajuster un AR(1) aux données précédemment générées. Examinez la table des t et comparez les estimations aux vraies valeurs. Par exemple, si la série temporelle est dansx, pour ajuster un AR(1) aux données, utilisezsarima(x, p = 1, d = 0, q = 0)ou simplementsarima(x, 1, 0, 0).
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# Generate 100 observations from the AR(1) model
x <- arima.sim(model = list(order = c(1, 0, 0), ar = .9), n = 100)
# Plot the generated data
# Plot the sample P/ACF pair
# Fit an AR(1) to the data and examine the t-table