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
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 dansxetacf2()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