Recordemos que para que una serie de tiempo, \(y_t\), sea estacionaria (en covarianza), es necesario que:

• Su media y varianza sean constantes. Esto es, para todo \(t\in\mathbb{Z}\):

\[ \begin{align} \mathbb{E}(y_t) &= \mu;\\ \mathbb{V}ar(y_t)&= \gamma_0<\infty \end{align} \]

• Su función de autocovarianza dependa solo del desfase \(h\) y no del tiempo \(t\):

\[ \gamma(t,\tau) = \gamma(h); \ h=|t-\tau| \ \in \mathbb{Z} \]

Las siguientes series de tiempo podrían ser estacionarias:

¿Qué pasa si las series de tiempo presentan una tendencia? Por ejemplo:

Los modelos con tendencia violan directamente la condición de estacionariedad, esto es debido a que:

\[ \mathbb{E}(y_t) = \mu_t. \]

Esto es, la media es una función del tiempo.

Existen dos formas diferentes para incorporar la tendencia en el modelo:

• Tendencia determinista. Incorporamos explícitamente en el modelo una tendencia determinada por una función del tiempo.

• Tendencia estocástica. La tendencia proviene de una raíz unitaria.

Tendencia determinista vs. Tendencia estocástica

Aunque los dos tipos de tendencia son visualmente parecidos, las implicaciones teóricas y prácticas son diferentes:

• En un modelo con tendencia determinista, los shocks tienen efectos transitorios. En un modelo con tendencia estocástica, los shocks tienen efectos permanentes.

• El procedimiento para remover la tendencia y hacer la serie estacionaria es distinto para cada caso.

La tendencia determinista caracteriza el patrón gradual y de largo plazo de las variaciones de una serie de tiempo. Esta tendencia puede ser lineal o no lineal.

Si la tendencia es lineal, podemos modelarla como:

\[ Trend_t = \beta_1 + \beta_2 t \]

Por ejemplo:

Si la tendencia es no lineal, podemos intentar modelarla como una función cuadrática:

\[ Trend_t = \beta_1 + \beta_2 t +\beta_3 t^2 \]

Por ejemplo, para el volumen de transacciones en la bolsa de Nueva York:

También podemos modelar la tendencia de una serie de tiempo suponiendo que esta sigue un comportamiento exponencial o logarítmico. Por ejemplo, la siguiente figura muestra el volumen de transacciones en la bolsa de Nueva York en logaritmos:

Nota que el logaritmo del volumen de transacciones de la bolsa de Nueva York parece exhibir un comportamiento lineal. Esto es común para muchas series económicas. Si tenemos un comportamiento exponencial, entonces:

\[ Trend_t = \beta_1 \exp^{\beta_2 t} \]

Pero si aplicamos logaritmos:

\[ \ln Trend_t = \ln \beta_1 + \beta_2 t \]

Naturalmente el logaritmo de \(Trend_t\) es una función lineal del tiempo.

En ausencia de un componente cíclico, los coeficientes de los modelos de tendencia lineal y cuadrática (o de polinomios de mayor orden) pueden estimarse por MCO. Los otros tipos de ecuaciones no lineales también pueden estimarse por MCO siempre y cuando estos admitan una transformación que nos resulte en una ecuación lineal.

Nota

Para tendencias no lineales, también podemos estimar los coeficientes mediante el método de los Mínimos Cuadrados No Lineales (MCNL).

Obviamente, podríamos ajustar más de un tipo de tendencia a una serie de datos. Para elegir cuál modelo es el más adecuado, podemos utilizar los criterios de información Akaike (AIC) y Schwartz (BIC o SIC).

Como sabemos, tenemos que elegir el modelo con el AIC y el BIC más bajos posibles.

Ejercicio

Considera el modelo:

\[ y_t = \beta_1 + \beta_2 t + \varepsilon_t \]

Muestra que \(y_t\) no es estacionario.

Ejercicio: “detrending”

Considera el modelo:

\[ y_t = \beta_1 + \beta_2 t + \varepsilon_t \]

Claramente, la transformación:

\[ \begin{align} y_t^\ast &= y_t - \beta_1 - \beta_2 t \\ & = \varepsilon_t \end{align} \] es estacionaria. Esta transformación puede hacerse operacional reemplazando \(\beta_1\) y \(\beta_2\) por sus estimadores de MCO.

Nota

Una serie con tendencia estacionaria es aquella que se vuelve estacionaria una vez eliminada su tendencia determinista. En otras palabras, el componente estocástico de la serie es estacionario, y la no estacionariedad proviene únicamente de una tendencia predecible.

Ejemplo

Las técnicas habituales para analizar y pronosticar series de tiempo, aplican también para los modelos de tendencia (ver modelos ARMA). Por ejemplo, para la tasa de participación en el mercado laboral de Estados Unidos tenemos:

Como hemos visto, un proceso AR(\(p\)):

\[ \begin{align} y_t &= \phi_1 y_{t-1} + \phi_2 y_{t-2} + \dots + \phi_p y_{t-p} + \varepsilon_t\\ \varepsilon_t &\sim WN(0,\sigma^2), \end{align} \]

es estacionario si las raíces del polinomio \(\Phi(L)=1-\phi_1 L - \phi_2 L^2 - \dots - \phi_p L^p=0\) se encuentran fuera del círculo unitario.

Una realización de un AR(1) estacionario con \(\phi_1 = 0.5\) es:

Si al menos una de las raíces del polinomio \(\Phi(L)=0\) es unitaria (es decir, igual a \(1\)), tendremos que la serie de tiempo es no estacionaria.

Cuando esto ocurre tenemos que:

\[ \sum_{i=1}^p \phi_i = 1 \]

Ejercicio

Calcula la función de autocovarianzas para mostrar que un proceso AR(\(1\)):

\[ \begin{align} y_t &= \phi_1y_{t-1} + \varepsilon_t;\\ \varepsilon_t &\sim WN(0,\sigma^2) \end{align} \]

no es estacionario cuando \(\phi_1=1\). Este proceso se conoce como caminata aleatoria.

La siguiente es una gráfica de una realización de una caminata aleatoria:

Observa que, a diferencia de un proceso autorregresivo, una caminata aleatoria no presenta reversión a la media.

Ejercicio

Considera un proceso AR(\(1\)):

\[y_t = \phi y_{t-1} + \varepsilon\]

1. Muestra que, cuando \(|\phi|<1\) su función de impulso-respuesta (IRF) es \(\psi_h = \phi^h\).

2. Analiza qué ocurre con la IRF cuando \(\phi=1\). En particular: ¿en cuál caso el efecto de un shock es transitorio? ¿En cuál caso es permanente?

Sea \(y_t\) una serie de tiempo. El operador de diferencias \(\Delta\) se define como:

\[ \begin{align} \Delta y_t &= y_t - y_{t-1}\\ &=y_t - Ly_t\\ &=(1-L)y_t \end{align} \]

Este operador se puede aplicar varias veces. Por ejemplo:

\[ \begin{align} \Delta(\Delta y_t) &= \Delta^2 y_t \\ &= (1-L)^2 y_t \\ &=(1-2L+L^2)y_t\\ &= y_t - 2y_{t-1} + y_{t-2} \end{align} \]

Nota que:

\[ \Delta (y_t - \mu) = \Delta y_t \]

Y que:

\[ \Delta^2 (y_t-\mu-\beta t) = \Delta^2 y_t \]

En general, si aplicamos el operador de diferencias \(d\) veces, eliminamos toda la información sobre cualquier tendencia temporal de orden menor a \(d\).

Ejemplo

Si aplicamos primeras diferencias a una caminata aleatoria, obtenemos una serie estacionaria:

Lo anterior nos permite definir un proceso integrado. Se dice que una serie \(y_t\) es integrada de orden \(d\), \(y_t \sim I(d)\), si necesita ser diferenciada \(d\) veces para convertirse en un proceso estacionario.

Y entonces:

Una serie de tiempo no estacionaria \(y_t \sim I(d)\) sigue un proceso ARIMA\((p,d,q)\) si: \[ \Delta^d y_t \sim \text{ARMA}(p, q) \]

Ejercicio

Considera un proceso AR(\(1\)) con \(\phi_1 = 1\) (una caminata aleatoria):

\[ \begin{align} y_t &= y_{t-1} + \varepsilon_t;\\ \varepsilon_t &\sim WN(0,\sigma^2) \end{align} \]

Muestra que \(y_t \sim \text{ARIMA}(0,1,0)\).

Ejercicio

Considera una caminata aleatoria con deriva \(\mu\):

\[ \begin{align} y_t &= \mu + y_{t-1} + \varepsilon_t;\\ \varepsilon_t &\sim WN(0,\sigma^2) \end{align} \]

Muestra que \(y_t \sim \text{ARIMA}(0,1,0)\).

Para diferenciar entre los modelos con tendencia determinista y tendencia estocástica podemos aplicar la prueba de Dickey-Fuller. En su versión más simple, la prueba considera el siguiente proceso AR(1):

\[ \begin{align} y_t &= \phi_1 y_{t-1} + \varepsilon_t;\\ \varepsilon_t &\sim WN(0,\sigma^2) \end{align} \]

Si hacemos una regresión de \(y_t\) sobre \(y_{t-1}\), podemos aplicar una prueba de significancia para el coeficiente \(\phi_1\):

\[ \begin{align} H_0: \phi_1= 1\\ H_1: \phi_1\neq 1 \end{align} \]

El estadístico de la prueba es:

\[ \hat{\tau} = \frac{\hat{\phi_1}-1}{\text{s.e.}(\hat{\phi_1})} \]

Alternativamente:

\[ \begin{align} \Delta y_t &= (\phi_1 -1)y_{t-1} + \varepsilon_t \\ &=\gamma y_{t-1} + \varepsilon_t \\ \end{align} \]

Si hacemos una regresión de \(\Delta y_t\) sobre \(y_{t-1}\), podemos aplicar una prueba de significancia para el coeficiente \(\gamma\):

\[ \begin{align} H_0: \gamma= 0\\ H_1: \gamma \neq 0 \end{align} \]

El estadístico de la prueba es:

\[ \hat{\tau} = \frac{\hat{\gamma}}{\text{s.e.}(\hat{\gamma})} \]

Un problema encontrado por Dickey y Fuller es que el estadístico de la prueba (conocido como estadístico de Dickey-Fuller) no sigue una distribución \(t\), sino una distribución especial conocida en la literatura como distribución de Dickey-Fuller.

Los valores críticos de la prueba pueden consultarse en varios libros de texto, así como en varios paquetes de computadora.

El caso que hemos explicado es uno de tres posibles que pueden ser considerados por la prueba. Cada caso depende de la hipótesis alternativa:

• Caso 1: \(y_t\) sigue un proceso estacionario con media cero.

• Caso 2: \(y_t\) sigue un proceso estacionario con media diferente de cero. Esto es, un proceso estacionario con deriva: \[ \Delta y_t = \alpha_0 + \gamma y_{t-1} + \varepsilon_t \]

• Caso 3: \(y_t\) sigue un proceso con tendencia estacionaria. Esto es, un proceso estacionario con tendencia determinista. \[ \Delta y_t = \alpha_0 + \alpha_1 t + \gamma y_{t-1} + \varepsilon_t \]

La prueba de Dickey-Fuller no es válida si los errores no se comportan como un ruido blanco. La Prueba Aumentada de Dickey-Fuller (ADF) resuelve este problema considerando un proceso AR(\(p\)):

\[ y_t = \sum_{i=1}^p \phi_j y_{t-j} + \varepsilon_t \]

O, en primeras diferencias:

\[ \Delta y_t = \gamma y_{t-1} + \sum_{i=2}^p \beta_i \Delta y_{t-i+1} + \varepsilon_t \] Si tenemos una raíz unitaria, entonces \(\gamma = 0\) y \(y_t\) sigue un proceso AR(\(p-1\)) en primeras diferencias. El estadístico de la prueba Dickey-Fuller sigue la misma distribución asintótica que antes.

Igual que antes, podemos generalizar para incluir una media diferente que cero y una tendencia determinista. Entonces, la hipótesis alternativa considera:

• Caso 1, proceso estacionario con media cero.

• Caso 2, proceso estacionario con deriva: \[ \Delta y_t = \alpha_0 + \gamma y_{t-1} + \sum_{i=2}^p \beta_i \Delta y_{t-i+1} + \varepsilon_t \]

• Caso 3, proceso estacionario con tendencia determinista: \[ \Delta y_t = \alpha_0 + \alpha_1 t + \gamma y_{t-1} + \sum_{i=2}^p \beta_i \Delta y_{t-i+1} + \varepsilon_t \]

Nota

Podemos probar de manera conjunta la ausencia de deriva, tendencia temporal y raíz unitaria de forma conjunta utilizando un estadístico tipo F, conocido como la F de Dickey-Fuller.

En la práctica no es tan claro qué caso de la prueba ADF considerar. Una opción es seguir el siguiente procedimiento.

La principal diferencia entre las tendencias determinista y aleatoria es:

• Si la tendencia es determinista, los shocks externos solo afectan al componente estacionario, no a la tendencia. Entonces, en presencia de shocks externos, \(y_t\) se desvía temporalmente de su trayectoria, pero su naturaleza lo obliga a regresar a su trayectoria original. El efecto de los shocks es transitorio.

• Si la tendencia es estocástica, los shocks se incorporan permanentemente al nivel de la serie (a \(y_t\)). Entonces, la trayectoria futura de la serie se construye a partir de este nuevo nivel. El efecto de los shocks es permanente.

Aplicar la transformación incorrecta para hacer el modelo estacionario tiene consecuencias graves.

• Si una serie es estacionaria en tendencia y aplicamos diferencias, obtenemos: \[ \Delta y_t = \beta + \varepsilon_t−\varepsilon_{t-1} \] Hemos introducido artificialmente una estructura de media móvil en los errores que no es invertible. Esto se conoce como sobrediferenciación, lo cual complica la estimación y el pronóstico.

• Si una serie es una caminata aleatoria y le removemos una tendencia lineal, los residuos de esa regresión seguirán conteniendo la raíz unitaria. No hemos solucionado el problema de no estacionariedad, simplemente hemos ajustado una tendencia espuria.

Consideremos el caso de una caminata aleatoria:

Si hacemos una regresión de \(y_t\) contra el tiempo, encontramos un coeficiente de tendencia altamente significativo. Naturalmente, este es resultado completamente falso. Este es el fenómeno de la tendencia espuria.

Ejemplo

Naturalmente, una vez que la serie de tiempo es estacionaria, todo el análisis (ajuste, pronósticos, etc.) es idéntico al análisis de los modelos ARMA(\(p\), \(q\)). Por ejemplo, considera la tasa de cambio Yen-Dólar (en logaritmos) para el periodo enero 1973 a julio de 1996:

Ejemplo

Utilizamos una prueba de Dickey-Fuller aumentada para incluir seis retardos, constante y tendencia lineal. Esta prueba muestra que la serie de tiempo no es estacionaria:

## 
## ############################################################### 
## # Augmented Dickey-Fuller Test Unit Root / Cointegration Test # 
## ############################################################### 
## 
## The value of the test statistic is: -2.5779 2.9597 3.3887

Ejemplo

En cambio, las primeras diferencias de la serie de tiempo son estacionarias:

## 
## ############################################################### 
## # Augmented Dickey-Fuller Test Unit Root / Cointegration Test # 
## ############################################################### 
## 
## The value of the test statistic is: -5.9335 17.6073

Ejemplo

Un modelo ARIMA(0,1,1) resulta en: