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
Istruzioni dell'esercizio
- Usa
diff()per differenziare i dati (d_birth). Usaacf2()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 indd_birth. Usa di nuovoacf2()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