Adattare un modello AR(1)
Ricorda che usi la coppia ACF e PACF per aiutarti a identificare gli ordini \(p\) e \(q\) di un modello ARMA. La tabella seguente riassume i risultati:
| AR(\(p\)) | MA(\(q\)) | ARMA(\(p,q\)) | |
|---|---|---|---|
| ACF | Decresce gradualmente | Si azzera dopo il ritardo \(q\) |
Decresce gradualmente |
| PACF | Si azzera dopo il ritardo \(p\) |
Decresce gradualmente | Decresce gradualmente |
In questo esercizio genererai dati dal modello AR(1), $$X_t = .9 X_{t-1} + W_t,$$ osserverai i dati simulati e la coppia ACF e PACF campionarie per determinare l’ordine. Poi adatterai il modello e confronterai i parametri stimati con quelli veri.
Per tutto il corso userai sarima() dal pacchetto astsa per adattare facilmente i modelli ai dati. Il comando produce un grafico diagnostico dei residui che può essere ignorato finché non si parlerà di diagnostica più avanti nel capitolo.
Questo esercizio fa parte del corso
Modelli ARIMA in R
Istruzioni dell'esercizio
- Il pacchetto astsa è precaricato.
- Usa il comando già scritto
arima.sim()per generare 100 osservazioni da un modello AR(1) con parametro AR pari a .9. Salvale inx. - Traccia i dati generati usando
plot(). - Traccia le ACF e PACF campionarie usando il comando
acf2()del pacchettoastsa. - Usa
sarima()diastsaper adattare un AR(1) ai dati generati in precedenza. Esamina la t-table e confronta le stime con i valori veri. Ad esempio, se la serie temporale è inx, per adattare un AR(1) ai dati usasarima(x, p = 1, d = 0, q = 0)oppure semplicementesarima(x, 1, 0, 0).
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# Generate 100 observations from the AR(1) model
x <- arima.sim(model = list(order = c(1, 0, 0), ar = .9), n = 100)
# Plot the generated data
# Plot the sample P/ACF pair
# Fit an AR(1) to the data and examine the t-table