ComenzarEmpieza gratis

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

Ver curso

Instrucciones del ejercicio

  • Usa diff() para diferenciar los datos (d_birth). Usa acf2() 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 en dd_birth. Usa otra llamada a acf2() 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

Editar y ejecutar código