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 vous familiariser avec le fonctionnement du modèle, vous allez analyser des données simulées à partir du modèle intégré $$ Y_t = .9 Y_{t-1} + W_t\, $$ où \(Y_t = \nabla X_t = X_t - X_{t-1}\). Dans ce cas, 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 se trouve 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 de ces données pour observer le comportement d'une série intégrée. Ensuite, vous différencierez la série pour la rendre stationnaire. Vous tracerez les données différenciées et les FAC/FACP correspondantes afin de voir l'effet de la différenciation.
Comme précédemment, le paquet 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.
Cette activité fait partie du cours
Modèles ARIMA en R
Instructions de l’exercice
- Tracez les données générées.
- Utilisez
acf2()du paquetastsapour 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 cela suggère un modèle AR(1) pour les données différenciées.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# Plot x
# Plot the P/ACF pair of x
# Plot the differenced data
# Plot the P/ACF pair of the differenced data