ComeçarComece de graça

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

Ver curso

Instruções do exercício

  • Use diff() para diferenciar os dados (d_birth). Use acf2() 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 em dd_birth. Use outra chamada a acf2() 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

Editar e executar o código