Análise de dados - taxa de natalidade
Agora você vai usar suas novas habilidades para ajustar cuidadosamente um modelo SARIMA à série temporal birth do pacote astsa. Os dados são nascimentos mensais (ajustados) nos Estados Unidos, em milhares, de 1948 a 1979, incluindo o baby boom após a Segunda Guerra Mundial.
Os dados birth estão plotados no seu console do R. Observe a tendência de longo prazo (random walk) e o componente sazonal dos dados.
Este exercício faz parte do curso
Modelos ARIMA em R
Instruções do exercício
- Use
diff()para diferenciar os dados (d_birth). Useacf2()para visualizar as ACF e PACF amostrais desses dados até a defasagem 60. Note a persistência sazonal. - Use outra chamada a
diff()para fazer a diferença sazonal dos dados. Salve emdd_birth. Use outra chamada aacf2()para visualizar a ACF e a PACF desses dados, novamente até a defasagem 60. Conclua que um modelo SARIMA(0,1,1)x(0,1,1)12 parece razoável. - Ajuste o modelo SARIMA(0,1,1)x(0,1,1)12. O que acontece?
- Adicione um parâmetro AR adicional (não sazonal,
p = 1) para explicar correlação adicional. O modelo se ajusta bem?
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# 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