IniziaInizia gratis

Analisi dei dati - tasso di natalità

Ora metterai in pratica le tue nuove competenze per stimare con cura un modello SARIMA per la serie temporale birth di astsa. I dati rappresentano le nascite vive mensili (corrette), in migliaia, negli Stati Uniti dal 1948 al 1979, e includono il baby boom post Seconda guerra mondiale.

I dati birth sono tracciati nella tua console R. Nota la tendenza di lungo periodo (random walk) e la componente stagionale dei dati.

Questo esercizio fa parte del corso

Modelli ARIMA in R

Visualizza il corso

Istruzioni dell'esercizio

  • Usa diff() per differenziare i dati (d_birth). Usa acf2() per visualizzare l'ACF e la PACF campionarie di questi dati fino al lag 60. Nota la persistenza stagionale.
  • Usa un'altra chiamata a diff() per effettuare la differenziazione stagionale dei dati. Salva il risultato in dd_birth. Usa di nuovo acf2() per visualizzare ACF e PACF di questi dati, ancora fino al lag 60. Concludi che un modello SARIMA(0,1,1)x(0,1,1)12 sembra ragionevole.
  • Stima il modello SARIMA(0,1,1)x(0,1,1)12. Cosa succede?
  • Aggiungi un parametro AR aggiuntivo (non stagionale, p = 1) per tenere conto di ulteriore correlazione. Il modello si adatta bene?

Esercizio pratico interattivo

Prova a risolvere questo esercizio completando il codice di esempio.

# 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

Modifica ed esegui il codice