ARIMA - prêt à l’emploi
Comme vous l’avez vu dans la vidéo, une série temporelle est dite ARIMA(\(p,d,q\)) si la série différenciée (d’ordre \(d\)) est ARMA(\(p,q\)).
Pour comprendre le fonctionnement du modèle, vous allez analyser des données simulées issues du modèle intégré $$ Y_t = .9 Y_{t-1} + W_t\, $$ où \(Y_t = \nabla X_t = X_t - X_{t-1}\). Ici, le modèle est un ARIMA(1,1,0) car les données différenciées suivent une autorégression d’ordre 1.
La série temporelle simulée est dans x et a été générée dans R comme suit
x <- arima.sim(model = list(order = c(1, 1, 0), ar = .9), n = 200).
Vous allez tracer les données générées ainsi que la FAC et la FACP empiriques correspondantes pour voir comment se comportent des données intégrées. Ensuite, vous appliquerez une différenciation pour rendre la série stationnaire. Vous tracerez les données différenciées et les FAC/FACP associées afin de constater l’effet de la différenciation.
Comme précédemment, le package astsa a été préchargé pour vous. Des données provenant d’un ARIMA(1,1,0) avec un paramètre AR de 0,9 sont enregistrées dans l’objet x.
Cet exercice fait partie du cours
Modèles ARIMA en R
Instructions
- Tracez les données générées.
- Utilisez
acf2()deastsapour tracer la paire FAC/FACP empirique des données générées. - Tracez les données différenciées.
- Faites un autre appel à
acf2()pour afficher la paire FAC/FACP empirique des données différenciées. Remarquez comment elles suggèrent un modèle AR(1) pour les données différenciées.
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# Plot x
# Plot the P/ACF pair of x
# Plot the differenced data
# Plot the P/ACF pair of the differenced data