ARIMA - plug and play
Come hai visto nel video, una serie temporale si chiama ARIMA(\(p,d,q\)) se la serie differenziata (di ordine \(d\)) è ARMA(\(p,q\)).
Per capire come funziona il modello, analizzerai dati simulati dal modello integrato $$ Y_t = .9 Y_{t-1} + W_t\, $$ dove \(Y_t = \nabla X_t = X_t - X_{t-1}\). In questo caso, il modello è un ARIMA(1,1,0) perché i dati differenziati sono un'autoregressione di ordine uno.
La serie temporale simulata è in x ed è stata generata in R come
x <- arima.sim(model = list(order = c(1, 1, 0), ar = .9), n = 200).
Tracerai i dati generati e le ACF e PACF campionarie dei dati generati per vedere come si comportano i dati integrati. Poi differenzierai i dati per renderli stazionari. Traccerai i dati differenziati e le corrispondenti ACF e PACF campionarie per vedere come la differenziazione fa la differenza.
Come prima, il pacchetto astsa è già caricato per te. Dati da un ARIMA(1,1,0) con parametro AR pari a .9 sono salvati nell'oggetto x.
Questo esercizio fa parte del corso
Modelli ARIMA in R
Istruzioni dell'esercizio
- Traccia i dati generati.
- Usa
acf2()daastsaper tracciare la coppia P/ACF campionaria per i dati generati. - Traccia i dati differenziati.
- Usa un'altra chiamata a
acf2()per visualizzare la coppia P/ACF campionaria per i dati differenziati. Nota come suggeriscano un modello AR(1) per i dati differenziati.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# Plot x
# Plot the P/ACF pair of x
# Plot the differenced data
# Plot the P/ACF pair of the differenced data