Modelli ARIMA automatici per serie temporali stagionali
Come hai visto nel video, la funzione auto.arima() funziona anche con dati stagionali. Nota che impostare lambda = 0 nella funzione auto.arima() — applicando quindi una trasformazione logaritmica — significa che il modello verrà adattato ai dati trasformati e che le previsioni verranno riportate alla scala originale tramite retro-trasformazione.
Dopo aver applicato summary() a questo tipo di modello adattato, potresti vedere un output simile a quello seguente che corrisponde a \((p,d,q)(P,D,Q)[m]\):
ARIMA(0,1,4)(0,1,1)[12]
In questo esercizio userai queste funzioni per modellare e prevedere i dati h02 già caricati, che contengono le vendite mensili di farmaci corticosteroidi in Australia.
Questo esercizio fa parte del corso
Previsioni in R
Istruzioni dell'esercizio
- Usando la funzione di plot standard, traccia i dati
h02in scala logaritmica per verificare che la varianza sia stabile. - Adatta un modello ARIMA stagionale alla serie
h02conlambda = 0. Salvalo infit. - Riassumi il modello adattato usando il metodo appropriato.
- Quali livelli di differenziazione sono stati utilizzati nel modello? Assegna la quantità di differenziazione al ritardo 1 a
de la differenziazione stagionale aD. - Disegna le previsioni per i prossimi 2 anni usando il modello adattato. Imposta
hdi conseguenza.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# Check that the logged h02 data have stable variance
h02 %>% ___ %>% ___
# Fit a seasonal ARIMA model to h02 with lambda = 0
fit <- ___
# Summarize the fitted model
___
# Record the amount of lag-1 differencing and seasonal differencing used
d <- ___
D <- ___
# Plot 2-year forecasts
fit %>% ___ %>% ___