CommencerCommencer 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 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

Afficher le cours

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 dans x.
  • Tracez les données générées avec plot().
  • Tracez les ACF et PACF empiriques à l’aide de la commande acf2() du package astsa.
  • Utilisez sarima() de astsa pour 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 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 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

Modifier et exécuter le code