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
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 usandodiff()elog(). - Usa un grafico multifigure per (1) tracciare le chiusure giornaliere del DJIA (
djia$Close) e notare che non è stazionario. I dati sono un oggettoxts. Poi (2) traccia i rendimenti approssimati del DJIA usandodiff()elog(). 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)