CommencerCommencer gratuitement

Ajuster un modèle ARMA

Vous êtes maintenant prêt à combiner le modèle AR et le modèle MA en un modèle ARMA. Nous avons généré des données à partir du modèle ARMA(2,1), $$X_t = X_{t-1} - .9 X_{t-2} + W_t + .8 W_{t-1}, $$ x <- arima.sim(model = list(order = c(2, 0, 1), ar = c(1, -.9), ma = .8), n = 250). Examinez les données simulées ainsi que la paire ACF et PACF empiriques pour déterminer un modèle possible.

Rappelez-vous que, pour les modèles ARMA(\(p, q\)), l’ACF et la PACF théoriques décroissent progressivement. Dans ce cas, les ordres sont difficiles à déduire des données et il peut ne pas être évident de voir si l’ACF ou la PACF échantillonnale se tronque ou décroît. Ici, vous connaissez les ordres réels du modèle ; ajustez donc un ARMA(2,1) aux données générées. Les stratégies générales de modélisation seront abordées plus loin dans le cours.

Cet exercice fait partie du cours

Modèles ARIMA en R

Afficher le cours

Instructions

  • Le package astsa est préchargé. 250 observations ARMA(2,1) sont dans x.
  • Comme dans les exercices précédents, utilisez plot() pour tracer les données générées dans x et acf2() pour afficher la paire ACF et PACF échantillonnales.
  • Utilisez sarima() pour ajuster un ARMA(2,1) aux données générées. Examinez la table des t et comparez les estimations aux vraies valeurs.

Exercice interactif pratique

Essayez cet exercice en complétant cet exemple de code.

# astsa is preloaded

# Plot x


# Plot the sample P/ACF of x


# Fit an ARMA(2,1) to the data and examine the t-table

Modifier et exécuter le code