Introducción
- La mayoría de las variables macroeconómicas y financieras son no estacionarias
- Se desplazan hacia arriba con el tiempo y a menudo exhiben características de una tendencia estocástica
- Engle & Granger propuso que el orden de integración combinado de varias variables puede ser menor que su orden de integración individual
- Permite el análisis de datos no estacionarios, donde podemos retener información relacionada con el nivel
- Facilita un análisis de equilibrio, donde las variables exhiben desviaciones de la tendencia no estacionaria (estado estacionario)
Cointegración definida
En la mayoría de los casos, cuando las variables
son no estacionarios \(y_{1,t}\) y \(y_{2,t}\) \(I(1)\) variables, una combinación lineal de estas variables también será no estacionariaSin embargo, en algunos casos la combinación lineal de las variables puede ser estacionaria.
Esto ocurre cuando las variables comparten la misma tendencia estocástica
El efecto de la tendencia estocástica común puede estar contenido cuando se combinan las variables
En estos casos, decimos que las variables están cointegradas
Cointegración definida
En la mayoría de los casos, cuando las variables \(y_{1,t}\) y \(y_{2,t}\) son no estacionarias \(I(1)\) variables, una combinación lineal de estas variables también será no estacionaria. Sin embargo, en algunos casos la combinación lineal de estas variables puede ser estacionaria. Esto sucede cuando las variables comparten las mismas tendencias estocásticas, que se cancelan cuando se combinan. En estos casos, decimos que las variables están cointegradas.
Considere dos variables, \(y_{1,t}\) y \(y_{2,t}\), que están integrados de primer orden, \(I(1)\)
Al retroceder estas variables entre sí, podríamos reorganizar el modelo de regresión lineal, de modo que
\[u_t=y_{1,t}-\beta_1 y_{2,t}\]
Ahora si el término de error,\(u_t\) es estacionario \(I(0)\), entonces, por definición, la combinación Lineal \(y_{1,t}-\beta_1 y_{2,t}\) también debe ser estacionario, ya que las propiedades del lado izquierdo deben ser iguales a las propiedades del lado derecho. Por lo tanto, mientras ambos \(y_{1,t} , y_{2,t}\) tienen tendencias estocásticas, decimos que las variables \(y_{1,t}\) y \(y_{2,t}\) están cointegrados, como la combinación lineal \(y_{1,t}-\beta_1 y_{2,t}\) y tienen las mismas propiedades estadísticas que una \(I(0)\) variable.
Cointegración definida
Usando álgebra matricial, si dejamos \(\mathbf{y_t}=(y_{1,t} , y_{2,t} )'\) denona un vector de \(2\times 1\) de variables\(I(1)\) y \(\beta=(1,-\beta_1)'\)
Se deduce que la relación \(\beta'y_t= y_{1,t} -\beta_1 y_{2,t}\)
Entonces el sistema se cointegra cuando \(\beta'\mathbf{y_t}\) es \(I(0)\)
En este caso las variables en el vector \(\mathbf{y_t}\) comparten una tendencia estocástica común y se desplazarán juntos en el equilibrio a largo plazo
El vector \(\beta\) se denomina vector cointegrante
Resumen
- El modelado de ecuaciones simples de variables no estacionarias da resultados espurios: a menos que las variables cointegraran
- La cointegración implica que las variables comparten una tendencia común
- La cointegración describe la relación a largo plazo entre variables
- A menudo, tales relaciones están dadas por la teoría económica.
- Para determinar si las variables comparten una o más relaciones de cointegración, podemos hacer lo siguiente:
- Determine si las series temporales son estacionarias o no estacionarias (análisis gráfico, prueba de raíz unitaria)
- Si las series son estacionarias, proceda a estimar el modelo dinámico en los niveles de los datos.
- Si las series no son estacionarias, proceda a probar la cointegración
Resumen
- Para probar la cointegración, se pueden usar dos enfoques:
- Enfoque de ecuación única: procedimiento de dos pasos de Engel-Granger
- Enfoque multivariante: prueba de Johansen
- Como la cointegración es inherentemente una propiedad del sistema, generalmente se prefiere el enfoque multivariante
- Si somos capaces de rechazar la hipótesis nula de no cointegración, proceda a estimar un modelo de corrección de equilibrio ya sea por:
- Modelo de ecuación única: ECM o ARDL
- VECM multivariante
- Si no se puede rechazar el nulo de no cointegración, proceda a estimar el modelo en primeras diferencias
Prueba de Cointegración
Johansen propuso 2 tipos de pruebas para \(r\) :
La prueba eigenvalor máximo:
Esta prueba está basada en la razón de máxima verosimilitud \(\ln \left[L_{M V}(r) / L_{M V}(r+1)\right],\) y se efectua secuencialmente para \[ r=0,1, \ldots, n-1 \] - Esta prueba corrobora la hipótesis nula de que el rango de cointegración es \(r\) versus la alterna de que el rango de cointegración es \(r+1\)
- El estadístico de prueba es \[ \ell_{r+1}^{*}-\ell_{r}^{*}=-\frac{T}{2} \ln \left(1-\widehat{\lambda}_{r+1}\right) \]
La prueba de la traza:
Esta prueba se basa en la razón de máxima verosimilitud \(\ln \left[L_{M V}(r) / L_{M V}(n)\right]\) y es efectuada secuencialmente para \(r=n-1, \ldots, 1,0\) Esta prueba comprueba la hipótesis nula de que el rango de cointegración es \(r\) frente a la alternativa que el rango de cointegración es \(n\) El estadístico de prueba es \[ \ell_{A}^{*}-\ell_{0}^{*}=-\frac{T}{2} \sum_{i=r+1}^{n} \log \left(1-\widehat{\lambda}_{i}\right) \]
library(urca)data(denmark)
sjd <- denmark[, c("LRM", "LRY", "IBO", "IDE")]
sjd.vecm <- ca.jo(sjd, ecdet = "const", type="eigen", K=2, spec="longrun",
season=4)
summary(sjd.vecm)##
## ######################
## # Johansen-Procedure #
## ######################
##
## Test type: maximal eigenvalue statistic (lambda max) , without linear trend and constant in cointegration
##
## Eigenvalues (lambda):
## [1] 4.331654e-01 1.775836e-01 1.127905e-01 4.341130e-02 6.927550e-16
##
## Values of teststatistic and critical values of test:
##
## test 10pct 5pct 1pct
## r <= 3 | 2.35 7.52 9.24 12.97
## r <= 2 | 6.34 13.75 15.67 20.20
## r <= 1 | 10.36 19.77 22.00 26.81
## r = 0 | 30.09 25.56 28.14 33.24
##
## Eigenvectors, normalised to first column:
## (These are the cointegration relations)
##
## LRM.l2 LRY.l2 IBO.l2 IDE.l2 constant
## LRM.l2 1.000000 1.0000000 1.0000000 1.000000 1.0000000
## LRY.l2 -1.032949 -1.3681031 -3.2266580 -1.883625 -0.6336946
## IBO.l2 5.206919 0.2429825 0.5382847 24.399487 1.6965828
## IDE.l2 -4.215879 6.8411103 -5.6473903 -14.298037 -1.8951589
## constant -6.059932 -4.2708474 7.8963696 -2.263224 -8.0330127
##
## Weights W:
## (This is the loading matrix)
##
## LRM.l2 LRY.l2 IBO.l2 IDE.l2 constant
## LRM.d -0.21295494 -0.00481498 0.035011128 2.028908e-03 -1.726523e-13
## LRY.d 0.11502204 0.01975028 0.049938460 1.108654e-03 9.428195e-14
## IBO.d 0.02317724 -0.01059605 0.003480357 -1.573742e-03 4.143714e-14
## IDE.d 0.02941109 -0.03022917 -0.002811506 -4.767627e-05 7.781415e-14
#data(finland)
sjf <- finland
sjf.vecm <- ca.jo(sjf, ecdet = "none", type="trace", K=2,
spec="longrun", season=4)
summary(sjf.vecm)##
## ######################
## # Johansen-Procedure #
## ######################
##
## Test type: trace statistic , with linear trend
##
## Eigenvalues (lambda):
## [1] 0.30932660 0.22599561 0.07308056 0.02946699
##
## Values of teststatistic and critical values of test:
##
## test 10pct 5pct 1pct
## r <= 3 | 3.11 6.50 8.18 11.65
## r <= 2 | 11.00 15.66 17.95 23.52
## r <= 1 | 37.65 28.71 31.52 37.22
## r = 0 | 76.13 45.23 48.28 55.43
##
## Eigenvectors, normalised to first column:
## (These are the cointegration relations)
##
## lrm1.l2 lny.l2 lnmr.l2 difp.l2
## lrm1.l2 1.0000000 1.000000 1.0000000 1.000000
## lny.l2 -0.9763252 -1.323191 -0.9199865 1.608739
## lnmr.l2 -7.0910749 -2.016033 0.2691516 -1.375342
## difp.l2 -7.0191097 22.740851 -1.8223931 -15.686927
##
## Weights W:
## (This is the loading matrix)
##
## lrm1.l2 lny.l2 lnmr.l2 difp.l2
## lrm1.d 0.033342108 -0.020280528 -0.129947614 -0.002561906
## lny.d 0.022544782 -0.005717446 0.012949130 -0.006265406
## lnmr.d 0.053505000 0.046876449 -0.007367715 0.002173242
## difp.d 0.005554849 -0.017353903 0.014561151 0.001531004