IniziaInizia gratis

Gestire trend ed eteroschedasticità

Qui renderemo stazionari dati non stazionari calcolando il rendimento o tasso di crescita come segue.

Spesso le serie temporali sono generate come $$X_t = (1 + p_t) X_{t-1}$$, cioè il valore della serie osservato al tempo \(t\) è uguale al valore osservato al tempo \(t-1\) moltiplicato per una piccola variazione percentuale \(p_t\) al tempo \(t\).

Un semplice esempio deterministico è depositare denaro in banca con un tasso di interesse fisso \(p\). In questo caso, \(X_t\) è il valore del conto al tempo \(t\) con un deposito iniziale \(X_0\).

Di solito, \(p_t\) è chiamato rendimento o tasso di crescita di una serie temporale, e questo processo è spesso stabile.

Per ragioni che esulano dallo scopo di questo corso, si può dimostrare che il tasso di crescita \(p_t\) può essere approssimato da $$Y_t = \log X_t - \log X_{t-1} \approx p_t.$$

In R, \(p_t\) si calcola spesso come diff(log(x)) e lo si può rappresentare in un'unica riga con plot(diff(log(x))).

Questo esercizio fa parte del corso

Modelli ARIMA in R

Visualizza il corso

Istruzioni dell'esercizio

  • Come prima, i pacchetti astsa e xts sono già caricati.
  • Genera un grafico multifigure per (1) tracciare i dati trimestrali del GNP USA (gnp) e notare che non è stazionario, e (2) tracciare il tasso di crescita approssimato del GNP USA usando diff() e log().
  • Usa un grafico multifigure per (1) tracciare le chiusure giornaliere del DJIA (djia$Close) e notare che non è stazionario. I dati sono un oggetto xts. Poi (2) traccia i rendimenti approssimati del DJIA usando diff() e log(). Come si confrontano con il tasso di crescita del GNP?

Esercizio pratico interattivo

Prova a risolvere questo esercizio completando il codice di esempio.

# astsa and xts are preloaded 

# Plot GNP series (gnp) and its growth rate
par(mfrow = c(2,1))
plot(gnp)


# Plot DJIA closings (djia$Close) and its returns
par(mfrow = c(2,1))
plot(djia$Close)


 
Modifica ed esegui il codice