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 graduellement | S'annule après le retard \(q\) |
Décroît graduellement |
| PACF | S'annule après le retard \(p\) |
Décroît graduellement | Décroît graduellement |
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 l'ACF et la PACF empiriques pour déterminer l'ordre. Ensuite, vous ajusterez le modèle et comparerez les paramètres estimés aux paramètres réels.
Tout au long de ce cours, vous utiliserez sarima() du paquet astsa pour ajuster facilement des modèles aux données. La commande produit un graphique de diagnostic des résidus, que vous pouvez ignorer jusqu'à ce que les diagnostics soient abordés plus loin dans le chapitre.
Cette activité fait partie du cours
Modèles ARIMA en R
Instructions de l’exercice
- Le paquet astsa est préchargé.
- Utilisez la commande déjà é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 l'ACF et la PACF empiriques à l'aide de la commande
acf2()du paquetastsa. - Utilisez
sarima()deastsapour ajuster un AR(1) aux données générées. Examinez la table des t et comparez les estimations aux vraies valeurs. Par exemple, si la série 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 ce code d’exemple.
# 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