Ajusta un modelo puramente estacional
Como con otros modelos, puedes ajustar modelos estacionales en R usando el comando sarima() del paquete astsa.
Para entender cómo funcionan los modelos puramente estacionales, lo mejor es usar datos simulados. Hemos generado 250 observaciones de un modelo puramente estacional dado por $$X_t = .9 X_{t-12} + W_t + .5 W_{t-12}\,,$$ que denotaríamos como un SARMA(P = 1, Q = 1)S = 12. Tienes representados tres años de datos y la ACF y PACF del modelo.
Vas a comparar los valores muestrales de ACF y PACF de los datos generados con los valores reales mostrados.
El paquete astsa ya está cargado y los datos generados están en x.
Este ejercicio forma parte del curso
Modelos ARIMA en R
Instrucciones del ejercicio
- Usa
acf2()para trazar la ACF y PACF muestrales de los datos generados hasta el rezago 60 y compáralas con los valores reales. Para estimar hasta el rezago 60, establece el argumentomax.lagigual a60. - Ajusta el modelo a los datos generados usando
sarima(). Además de los argumentosp,dyqen tu comandosarima(), especificaP,D,QyS(recuerda que R distingue mayúsculas y minúsculas).
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# Plot sample P/ACF to lag 60 and compare to the true values
acf2(___, max.lag = ___)
# Fit the seasonal model to x
sarima(x, p = 0, d = 0, q = 0, P = ___, D = 0, Q = ___, S = ___)