Análisis de datos - tasa de natalidad
Ahora vas a aplicar tus nuevas habilidades para ajustar con cuidado un modelo SARIMA a la serie temporal birth de astsa. Los datos son nacimientos vivos mensuales (ajustados) en miles para Estados Unidos, 1948-1979, e incluyen el baby boom tras la Segunda Guerra Mundial.
Los datos de birth se han representado en tu consola de R. Fíjate en la tendencia a largo plazo (random walk) y en el componente estacional de los datos.
Este ejercicio forma parte del curso
Modelos ARIMA en R
Instrucciones del ejercicio
- Usa
diff()para diferenciar los datos (d_birth). Usaacf2()para ver la ACF y PACF muestrales de estos datos hasta el rezago 60. Observa la persistencia estacional. - Realiza otra llamada a
diff()para tomar la diferencia estacional de los datos. Guarda esto endd_birth. Usa otra llamada aacf2()para ver la ACF y la PACF de estos datos, de nuevo hasta el rezago 60. Concluye que un modelo SARIMA(0,1,1)x(0,1,1)12 parece razonable. - Ajusta el modelo SARIMA(0,1,1)x(0,1,1)12. ¿Qué ocurre?
- Añade un parámetro AR adicional (no estacional,
p = 1) para tener en cuenta correlación adicional. ¿Se ajusta bien el modelo?
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# Plot P/ACF to lag 60 of differenced data
d_birth <- diff(birth)
# Plot P/ACF to lag 60 of seasonal differenced data
dd_birth <- diff(d_birth, lag = 12)
# Fit SARIMA(0,1,1)x(0,1,1)_12. What happens?
# Add AR term and conclude