Series de Tiempo

Modelo GJR‑GARCH

Universidad Central del Ecuador Facultad de Ciencias Económicas Carrera de Estadística
Grupo 6  ·  Arboleda Melany · Hernández Angie · Mansilla Marley · Moncayo Bryan
Datos: Banco Central del Ecuador (2007–2026)
DESLIZA
01 · Marco teórico · Documento GARCH

Fundamento teórico: del ARCH al GARCH

Desarrollo íntegro del documento base GARCH, con sus cuatro apartados y ecuaciones.

1 · Introducción

Por qué la varianza no es constante

Los modelos econométricos clásicos asumen homocedasticidad: una varianza del error constante en el tiempo. Sin embargo, las series financieras muestran regularidades empíricas (Cont, 2001) que contradicen ese supuesto: agrupamiento de volatilidad (periodos de alta o baja volatilidad tienden a persistir), colas pesadas (más eventos extremos que en una normal) y heterocedasticidad condicional (la varianza, condicionada al pasado, cambia en el tiempo). Estas regularidades motivaron el desarrollo de los modelos ARCH y, posteriormente, GARCH.

2 · Modelo ARCH

Engle (1982)

El primer modelo formal de heterocedasticidad condicional autorregresiva ARCH fue propuesto por Robert F. Engle (1982). Engle (1982) introdujo una nueva clase de procesos estocásticos en los que la varianza condicional del término de error depende de los errores al cuadrado observados en periodos anteriores.

ecuación ARCH

Su principal limitación es que, para capturar bien la persistencia de la volatilidad, suele requerir muchos rezagos q, con el consecuente costo en grados de libertad (Engle, 1982). Esta limitación es precisamente lo que motiva la generalización propuesta por Bollerslev (1986).

3 · Modelo GARCH

Bollerslev (1986)

Con el fin de superar esta limitación, Bollerslev (1986) propuso una generalización del modelo ARCH, denominada GARCH, en la que la varianza condicional depende no solo de los errores al cuadrado pasados, sino también de sus propios valores rezagados. Esta extensión es conceptualmente análoga a la manera en que un proceso ARMA generaliza a un proceso autorregresivo puro (Bollerslev, 1986).

Un modelo GARCH(p,q) se define como:

ecuación GARCH(p,q)

Donde:

  • α₀ : constante o intercepto del modelo, con α₀ > 0.
  • αᵢ : coeficiente que pondera el efecto del shock rezagado i sobre la varianza actual, con αᵢ ≥ 0.
  • ε²t−i : error (innovación) al cuadrado observado i periodos atrás; captura el impacto de los shocks pasados.
  • σ²t−j : varianza condicional rezagada j periodos; es la memoria propia del modelo sobre su volatilidad pasada.
  • βj : coeficiente que pondera el efecto de la varianza rezagada j sobre la varianza actual, con βj ≥ 0.
  • p : número de rezagos de la propia varianza condicional incluidos en el modelo.

La especificación más utilizada en la práctica —y prácticamente un estándar en la modelación de series financieras— es el GARCH(1,1):

ecuación GARCH(1,1)

En esta versión, α₁ mide el efecto del último shock (ε²t−1) y β₁ mide cuánto persiste la volatilidad del periodo inmediatamente anterior (σ²t−1). La ventaja central del modelo GARCH frente al ARCH es su parsimonia: con un número reducido de parámetros —típicamente tres en la especificación (1,1)— es posible representar dinámicas de persistencia en la volatilidad que un modelo ARCH puro solo lograría replicar con un número mucho mayor de rezagos (Bollerslev, 1986).

4 · Limitación del modelo GARCH

La asimetría que no captura

El GARCH estándar responde solo a la magnitud del shock pasado (ε²t−1), sin distinguir su signo: una buena y una mala noticia de igual magnitud afectan igual a la volatilidad futura. Sin embargo, la evidencia empírica muestra una asimetría sistemática —el efecto apalancamiento (leverage effect), identificado por Black (1976)—: las malas noticias elevan la volatilidad futura más que las buenas de magnitud comparable, en parte porque una caída en el precio de las acciones incrementa el apalancamiento financiero de la empresa y, con ello, su riesgo percibido.

Como el GARCH simétrico no captura esta respuesta diferenciada, surgieron especificaciones alternativas como el modelo de Glosten, Jagannathan y Runkle (1993) —GJR-GARCH—, que añade un término para el impacto diferencial de los shocks negativos sobre la varianza condicional.

02 · Marco teórico · Documento GJR-GARCH

Modelo GJR-GARCH(1,1)

Introducción

Glosten, Jagannathan y Runkle (1993)

El modelo GJR-GARCH fue propuesto en 1993 por Lawrence Glosten, Ravi Jagannathan y David Runkle con el objetivo de mejorar el modelo GARCH tradicional. Su principal aporte fue corregir una limitación importante del modelo original, que no diferenciaba entre buenas y malas noticias en el mercado. Por eso, este modelo busca representar de forma más realista cómo cambia la volatilidad, especialmente en situaciones de caídas, lo que lo hace más útil para valorar derivados financieros y para la gestión del riesgo en inversiones.

Modelo matemático

Las dos ecuaciones

El modelo matemático debe presentarse en dos partes: la ecuación de la media (de dónde salen los errores) y la ecuación de la varianza condicional (el modelo GJR-GARCH).

La ecuación de la media condicional

El rendimiento del activo

Antes de calcular la volatilidad, se debe definir el comportamiento del rendimiento del activo rt:

media

Donde:

  • rt : Rendimiento en el tiempo t
  • μt : Media condicional
  • εt : El residuo, error o shock del mercado en el tiempo t

El residuo se define como:

residuo

donde Zt es una variable aleatoria independiente e idénticamente distribuida (i.i.d.) con media 0 y varianza 1:

znorm
La ecuación de la varianza condicional: GJR-GARCH(1,1)

El término asimétrico

El modelo GJR-GARCH(1,1) amplía la ecuación clásica de varianza añadiendo un término adicional diseñado específicamente para capturar choques asimétricos. La estructura matemática es la siguiente:

gjr
  • σt²: varianza o volatilidad en el periodo t
  • ω: Omega es la variación base a largo plazo.
  • α: Alfa es el impacto del choque previo (la parte simétrica).
  • β: Beta es la persistencia de la volatilidad (la memoria del modelo).
  • γ: Gamma es el parámetro de asimetría, el corazón del modelo.
  • ε²t-1: representa el error o shock del modelo, es decir, la diferencia entre el valor observado y el esperado.
  • It-1: variable indicadora
¿Cómo funciona la variable indicadora It-1?

Un interruptor para las malas noticias

It-1 es una función "dummy" o interruptor que toma los siguientes valores:

ind
  • It-1 = 1 si εt-1 < 0 (Mala noticia / Retorno negativo).
  • It-1 = 0 si εt-1 ≥ 0 (Buena noticia / Retorno positivo o cero).
Interpretación

Buenas vs. malas noticias

Cuando el shock anterior es positivo
pos

Es decir, cuando hubo una subida en el mercado, la volatilidad se calcula de manera normal con los términos omega, alfa y beta, que representan el nivel base, el impacto del shock pasado y la influencia de la volatilidad pasada.

Cuando el shock anterior es negativo
neg

El modelo GJR-GARCH incorpora una función indicadora que permite distinguir entre shocks positivos y negativos. Cuando el shock es negativo, se activa un término adicional representado por el parámetro γ, aumentando el impacto total sobre la volatilidad a α+γ.

Esto implica que las caídas del mercado generan un mayor incremento en la volatilidad que las subidas, fenómeno conocido como efecto de apalancamiento. En consecuencia, el modelo refleja que las malas noticias tienen un impacto más fuerte sobre el riesgo que las buenas.
Estimación de los parámetros

Máxima verosimilitud

Los parámetros ω, α, β y γ se estiman a partir de datos históricos mediante el método de máxima verosimilitud, el cual busca los valores que mejor explican los datos observados. En la práctica, este proceso se realiza con herramientas como R o Python, que calculan automáticamente los parámetros óptimos.

Ejemplo numérico sencillo

Igual choque, distinto efecto

Supongamos los siguientes valores: ω = 0.02, α = 0.05, β = 0.90 y γ = 0.08, con una varianza previa de 1.

Caso positivo
σ² = 0.97
0.02 + 0.05(1) + 0.90(1)
Caso negativo
σ² = 1.05
0.02 + (0.05 + 0.08)(1) + 0.90(1)

Este resultado muestra que, ante un choque de igual magnitud, la volatilidad es mayor cuando la noticia es negativa, lo que refleja el efecto de asimetría del modelo.

Estabilidad y estacionariedad

Para que el modelo funcione bien

Para que el modelo GJR-GARCH funcione bien y no dé errores, debe cumplir dos condiciones.

1 · Varianza positiva

La volatilidad nunca es negativa

La varianza siempre debe ser positiva, es decir, la volatilidad nunca puede ser negativa, porque no tendría sentido hablar de "inestabilidad negativa".

ω > 0    α ≥ 0    β ≥ 0    α + γ ≥ 0

Omega debe ser positivo, alfa y beta no pueden ser negativos, y la suma de alfa y gamma tampoco.

2 · Estacionariedad (persistencia)

Vuelve a un nivel normal

El modelo debe ser estable en el tiempo, lo que significa que la volatilidad no crece sin control, sino que con el tiempo vuelve a un nivel normal. Esto se cumple cuando la suma de los parámetros (alfa, beta y parte de gamma) es menor que uno.

pers

(Nota: Se multiplica gamma γ por ½ asumiendo que la distribución de los residuos es simétrica, por lo que las malas noticias ocurren el 50% de las veces).

El efecto apalancamiento

Por qué importa la asimetría

El efecto de apalancamiento explica que en los mercados financieros las malas noticias tienen un impacto más fuerte que las buenas. Es decir, cuando los precios bajan, la volatilidad aumenta más que cuando los precios suben en la misma proporción.

Qué lo causa · 1

Retroalimentación

Existe un efecto de retroalimentación donde la volatilidad genera más miedo, lo que hace que los inversionistas vendan y los precios sigan cayendo.

Qué lo causa · 2

Reacción a las pérdidas

Las personas reaccionan más fuertes ante las pérdidas que ante las ganancias, lo que provoca pánico en el mercado.

Qué lo causa · 3

Iliquidez

En momentos de crisis aparecen problemas como falta de liquidez o ventas obligadas, lo que aumenta aún más la inestabilidad.

Por qué es importante

Ignorar esta asimetría puede ser peligroso, ya que los modelos que no la consideran tienden a subestimar el riesgo cuando el mercado cae. Esto significa que se puede creer que el riesgo es menor de lo que realmente es, especialmente en momentos de crisis.

Relación con GJR-GARCH

El modelo GJR-GARCH es importante porque logra capturar esta asimetría de forma sencilla. A través de un parámetro adicional, permite reflejar que los choques negativos generan más volatilidad que los positivos.

Aplicación

¿Cuándo deberías usar GJR-GARCH?

El modelo GJR-GARCH se utiliza cuando se desea analizar la volatilidad en situaciones donde las malas noticias tienen un mayor impacto que las buenas. GJR-GARCH es especialmente valioso cuando:

Uso

Mercados bursátiles

Los efectos de apalancamiento son robustos en índices bursátiles y acciones individuales.

Uso

Previsión de crisis

Los shocks a la baja requieren respuestas de volatilidad diferentes.

Uso

Gestión de riesgos

Los cálculos VaR y ES se benefician de dinámicas asimétricas.

Uso

Valoración de opciones

Los sesgos de volatilidad implícita reflejan efectos de apalancamiento.

Limitaciones del modelo

A pesar de sus ventajas

A pesar de sus ventajas, el modelo GJR-GARCH presenta algunas limitaciones. En primer lugar, solo distingue entre choques positivos y negativos, sin considerar su magnitud relativa. En segundo lugar, puede ser sensible a valores extremos en los datos, lo que afecta la estimación de los parámetros. Finalmente, requiere una cantidad considerable de datos históricos para estimar correctamente el parámetro de asimetría.

Conclusiones del modelo

Síntesis teórica

El modelo GJR-GARCH es un modelo de volatilidad que incorpora asimetría, permitiendo capturar el mayor impacto de las malas noticias sobre el riesgo, siendo ampliamente utilizado en el análisis de mercados financieros.

El modelo GJR-GARCH captura eficazmente el efecto apalancamiento con un solo parámetro adicional, proporcionando una herramienta práctica e interpretable para gestores de riesgos e investigadores que necesitan modelar la volatilidad asimétrica en series de retornos financieros.

Además, destaca por su simplicidad e interpretabilidad, ya que con un solo parámetro adicional mejora significativamente la precisión del modelo, siendo ampliamente utilizado por gestores de riesgo e investigadores en finanzas.

02 · Ejercicio en Excel

Simulación Monte Carlo del GJR‑GARCH(1,1)

Se simula el modelo paso a paso para ver cómo la varianza condicional evoluciona y responde a los shocks. Descarga la hoja o experimenta con el simulador en vivo.

Parámetros en vivo
Persistencia α+β+γ/2 ⚠ ≥ 1
⬇ Descargar simulación (.xlsx)

Trayectoria de volatilidad simulada

Edita cualquier Zt en la tabla (celdas azules) o genera nuevos shocks, y ajusta los parámetros para ver cómo cambia la volatilidad.

ObsZtσ²tσtεtIPrecio
03 · Aplicación en Python

Ingresos petroleros del Ecuador

Datos mensuales del Banco Central del Ecuador (ene‑2007 a may‑2026). Se modela la volatilidad de los retornos logarítmicos con la librería arch.

0
Observaciones mensuales
0
Años de historia
0
Coef. de variación
0
Curtosis en exceso

Serie de ingresos por exportaciones petroleras públicas

Miles de USD · marcados los grandes choques externos.

Editor interactivo

El código completo del modelo ▶ Ejecución real · paso a paso

Este es el código original del Grupo 6 (38 celdas). Púlsale Ejecutar y correrá celda por celda en el servidor real (con arch, statsmodels y seaborn), mostrando cada salida y cada gráfica. Puedes editarlo antes de correr.

TRABAJO_GRUPO_6_MODELO_GJR_GARCH.py
▶ Pulsa Ejecutar para correr el análisis paso a paso; cada celda se ejecuta y muestra su salida.
Diagnóstico previo

¿Por qué un modelo de volatilidad?

Estacionariedad

ADF & KPSS

ADF = −23.24 (p<0.0001) y KPSS p>0.10: retornos estacionarios.

Efecto ARCH

ARCH‑LM

LM = 80.52 (p<0.0001): la volatilidad se agrupa. Justifica GARCH.

Normalidad

Colas pesadas

Curtosis 14.88 y JB p<0.0001 → se usa la t de Student.

Distribución de los retornos

dist

Densidad real (dorado) con colas más gruesas que la normal (rojo).

ACF de los retornos al cuadrado

acf

Autocorrelación = agrupamiento de volatilidad.

04 · Estimación y resultados

Comparación de modelos y efecto apalancamiento

Seis especificaciones (GARCH, GJR, EGARCH × normal / t) comparadas por AIC; luego se contrasta formalmente la asimetría.

Ranking de modelos por AIC (menor es mejor)

ModeloLog‑LikAICBICParámetros
La pregunta de investigación

¿Existe el efecto apalancamiento?

No basta el p‑valor de γ. Tres señales; la evidencia resulta débil y no concluyente.

Coeficiente γ
p = 0.121
No significativo
Razón de verosimilitud
p = 0.046
Marginalmente sig.
Test de Engle‑Ng
p = 0.184
Sin asimetría
Lectura honesta: solo la prueba LR respalda la asimetría, y apenas bajo 0.05. No hay evidencia robusta para rechazar H₀. LR contrasta  .
Parámetro (GJR‑GARCH‑t)Coeficientep‑valorSignificancia

Volatilidad condicional estimada

El riesgo se dispara en los grandes choques externos (%).

Pronóstico de volatilidad a 12 meses

Intervalo del 95% por simulación Monte Carlo (5.000 trayectorias).

Residuos estandarizados

residuos

QQ‑plot

qq

Índice de Riesgo Petrolero

semáforo
Modelo bien especificado: residuos sin autocorrelación (Ljung‑Box p=0.118) ni efecto ARCH remanente (p=0.984).
05 · Referencias

Bibliografía

Fuentes citadas en los documentos teóricos (GARCH y GJR-GARCH).

Documento GARCH

Black, F. (1976). Studies of stock price volatility changes. Proceedings of the 1976 Meeting of the Business and Economic Statistics Section, 177–181. American Statistical Association.

Bollerslev, T. (1986). Generalized autoregressive conditional heteroskedasticity. Journal of Econometrics, 31(3), 307–327.

Cont, R. (2001). Empirical properties of asset returns: Stylized facts and statistical issues. Quantitative Finance, 1(2), 223–236.

Engle, R. F. (1982). Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation. Econometrica, 50(4), 987–1007.

Glosten, L. R., Jagannathan, R., & Runkle, D. E. (1993). On the relation between the expected value and the volatility of the nominal excess return on stocks. The Journal of Finance, 48(5), 1779–1801.

Documento GJR-GARCH

Dube, F. (2026, abril). Building volatility models in MQL5 (part II): Implementing GJR-GARCH and TARCH in MQL5. MQL5 Community. https://www.mql5.com/en/articles/22258

Lakhloufi, H. (2025, mayo 8). GJR-GARCH asymmetric volatility calculator. MetricGate. https://metricgate.com/docs/gjr-garch-model/

Rao, M. (2025, mayo 7). GARCH vs. GJR-GARCH models in Python for volatility forecasting. QuantInsti. https://blog.quantinsti.com/garch-gjr-garch-volatility-forecasting-python/

V-Lab: GJR-GARCH Volatility documentation. (s/f). V-Lab. Recuperado el 10 de julio de 2026, de https://vlab.stern.nyu.edu/docs/volatility/GJR-GARCH