ARIMA - plug and play
Como viste en el vídeo, una serie temporal se denomina ARIMA(\(p,d,q\)) si la serie diferenciada (de orden \(d\)) es ARMA(\(p,q\)).
Para hacerte una idea de cómo funciona el modelo, vas a analizar datos simulados del modelo integrado $$ Y_t = .9 Y_{t-1} + W_t\, $$ donde \(Y_t = \nabla X_t = X_t - X_{t-1}\). En este caso, el modelo es un ARIMA(1,1,0) porque los datos diferenciados siguen una autorregresión de orden uno.
La serie temporal simulada está en x y se generó en R como
x <- arima.sim(model = list(order = c(1, 1, 0), ar = .9), n = 200).
Vas a representar los datos generados y la ACF y PACF muestrales de esos datos para ver cómo se comportan los datos integrados. Después, diferenciarás los datos para hacerlos estacionarios. Representarás los datos diferenciados y las ACF y PACF muestrales correspondientes para comprobar cómo la diferenciación marca la diferencia.
Como antes, el paquete astsa ya está cargado. Los datos de un ARIMA(1,1,0) con parámetro AR 0.9 están guardados en el objeto x.
Este ejercicio forma parte del curso
Modelos ARIMA en R
Instrucciones del ejercicio
- Representa los datos generados.
- Usa
acf2()deastsapara trazar el par ACF/PACF muestral de los datos generados. - Representa los datos diferenciados.
- Haz otra llamada a
acf2()para ver el par ACF/PACF muestral de los datos diferenciados. Observa cómo sugieren un modelo AR(1) para los datos diferenciados.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# Plot x
# Plot the P/ACF pair of x
# Plot the differenced data
# Plot the P/ACF pair of the differenced data