CommencezCommencez gratuitement

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

Voir le cours

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 dans x.
  • Tracez les données générées avec plot().
  • Tracez l'ACF et la PACF empiriques à l'aide de la commande acf2() du paquet astsa.
  • Utilisez sarima() de astsa pour 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 dans x, pour ajuster un AR(1) aux données, utilisez sarima(x, p = 1, d = 0, q = 0) ou simplement sarima(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

Modifier et exécuter le code