Cómo tratar la tendencia y la heterocedasticidad
Aquí vamos a forzar datos no estacionarios a ser estacionarios calculando el retorno o la tasa de crecimiento de la siguiente forma.
A menudo, las series temporales se generan como $$X_t = (1 + p_t) X_{t-1}$$, lo que significa que el valor de la serie observado en el instante \(t\) es igual al valor observado en \(t-1\) más un pequeño cambio porcentual \(p_t\) en el instante \(t\).
Un ejemplo determinista sencillo es ingresar dinero en un banco con un interés fijo \(p\). En este caso, \(X_t\) es el valor de la cuenta en el periodo \(t\) con un depósito inicial de \(X_0\).
Normalmente, a \(p_t\) se le llama retorno o tasa de crecimiento de una serie temporal, y este proceso suele ser estable.
Por razones que quedan fuera del alcance de este curso, se puede demostrar que la tasa de crecimiento \(p_t\) puede aproximarse por $$Y_t = \log X_t - \log X_{t-1} \approx p_t.$$
En R, \(p_t\) se suele calcular como diff(log(x)) y se puede representar en una línea con plot(diff(log(x))).
Este ejercicio forma parte del curso
Modelos ARIMA en R
Instrucciones del ejercicio
- Como antes, los paquetes astsa y xts están precargados.
- Genera una figura con múltiples paneles para (1) graficar los datos trimestrales del PIB de EE. UU. (
gnp) y comprobar que no es estacionaria, y (2) graficar la tasa de crecimiento aproximada del PIB usandodiff()ylog(). - Usa otra figura con múltiples paneles para (1) graficar los cierres diarios del DJIA (
djia$Close) y comprobar que no es estacionaria. Los datos son un objetoxts. Luego (2) grafica los retornos aproximados del DJIA usandodiff()ylog(). ¿Cómo se compara esto con la tasa de crecimiento del PIB?
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# 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)