IniziaInizia gratis

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

Visualizza il corso

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 in x.
  • Traccia i dati generati usando plot().
  • Traccia le ACF e PACF campionarie usando il comando acf2() del pacchetto astsa.
  • Usa sarima() di astsa per 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 è in x, per adattare un AR(1) ai dati usa sarima(x, p = 1, d = 0, q = 0) oppure semplicemente sarima(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

Modifica ed esegui il codice