Ajuster un modèle ARMA
Vous êtes maintenant prêt à combiner le modèle AR et le modèle MA dans le modèle ARMA. Nous avons généré des données selon le 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 FAC et FACP empiriques pour déterminer un modèle plausible.
Rappelez-vous que, pour les modèles ARMA(\(p, q\)), la FAC et la FACP théoriques décroissent graduellement. Dans ce contexte, il est difficile de discerner les ordres à partir des données et il se peut qu'il ne soit pas clair si la FAC ou la FACP empirique se coupe net 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.
Cette activité fait partie du cours
Modèles ARIMA en R
Instructions de l’exercice
- Le paquet 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 dansxet utilisezacf2()pour afficher la paire FAC et FACP empiriques. - 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 ce code d’exemple.
# 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