SciELO - Scientific Electronic Library Online

 
vol.22 número53Comercio intra-industrial y variaciones del empleo en la industria colombiana, 1991-2015Factores de en uso de los medios de pago móviles millennials y centennials índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • En proceso de indezaciónCitado por Google
  • No hay articulos similaresSimilares en SciELO
  • En proceso de indezaciónSimilares en Google

Compartir


Semestre Económico

versión impresa ISSN 0120-6346versión On-line ISSN 2248-4345

Semest. Econ. vol.22 no.53 Medellín oct./dic. 2019

https://doi.org/10.22395//seec.v22n53a3 

Artículos

Valoración de riesgo mediante modelos GARCH y simulación Montecarlo: evidencia del mercado accionario colombiano*

Value at Risk Through GARCH Models and Montecarlo Simulation: Evidence from the Colombian Stock Market

Avaliação do risco mediante modelos GARCH e simulação Montecarlo: evidência do mercado acionista colombiano

María Ines Barbosa Camargo** 

Alejandra Salazar Sarmiento*** 

Kelly Jhohana Peñaloza Gómez**** 

** Economista y magíster en Ciencias Económicas, Universidad Nacional de Colombia, Bogotá, Colombia. Doctora en Análisis Económico Aplicado e Historia Económica, Universidad de Sevilla, Sevilla, España. Docente investigadora, Universidad de la Salle, Bogotá, Colombia. Integrante del Grupo de Investigación en Estudios Sociales, Financieros e Internacionales (ESFI), Bogotá, Colombia. Correo electrónico: mibarbosa@unisalle.edu.co. Orcid: https://orcid.org/0000-0002-7705-7983

*** Profesional en Finanzas y Comercio Internacional, Universidad de la Salle, Bogotá, Colombia. Correo electrónico: salejandra64@unisalle.edu.co. Orcid: https://orcid.org/0000-0001-8445-5202

**** Profesional en Finanzas y Comercio Internacional, Universidad de la Salle, Bogotá, Colombia. Correo electrónico: kpenaloza01@unisalle.edu.co. Orcid: http://orcid.org/0000-0001-6835-8920


RESUMEN

Este documento evalúa el comportamiento de varios modelos de volatilidad en estimaciones de un día del valor en riesgo (VaR) de veinticuatro series de retornos de acciones en Colombia con diferentes distribuciones. Al considerar que todas las series de retornos presentan clúster de volatilidad y memoria de largo plazo, se utilizan modelos tipo GARCH que incluyen diferentes distribuciones: normal, T-Student y GED. Los hallazgos corroboran la dificultad de elegir un único modelo para el cálculo del VaR, pero validan el uso de modelos paramétricos con distribución normal y simulación Montecarlo en mercados financieros emergentes.

CLASIFICACIÓN JEL: G17, G32, C58

CONTENIDO: Introducción; 1. Breve revisión de la literatura: valor en riesgo; 2. Metodología; 3. Resultados; Conclusiones; Bibliografía.

PALABRAS CLAVE: Valor en riesgo; modelos GARCH; simulación; Colombia

ABSTRACT

This paper evaluates the performance of several volatility models for estimating one-day-ahead Value-at-Risk (VaR) of twenty-four stocks return series in Colombia, using a number of distributional assumptions. Because all return series exhibit volatility clustering and long-range memory GARCH-type models including models under normal, T-Student and generalized error distribution are examined. The findings corroborate the difficulty of choosing a single model for calculating VaR, but validate the use of parametric models with normal distribution and Montecarlo simulation in emerging financial markets.

JEL CLASSIFICATION: G17, G32, C58

CONTENT: Introduction; 1. Brief literature review; Value-at-Risk; 2. Methodology; 3. Results; Conclusions; Bibliography.

KEYWORDS: Value-at-Risk; GARCH-type models; simulation; Colombia

RESUMO

Este documento avalia o comportamento de vários modelos de volatilidade em estimativas de um dia do Value at Risk (VaR) de 24 séries de retornos das ações na Colômbia com diferentes distribuições. Ao considerar que todas as séries de retornos apresentam cluster de volatilidade e memória de longo prazo, são utilizados modelos tipo GARCH que incluem diferentes distribuições: normal, t-student e GED. Os achados corroboram a dificuldade de escolher um único modelo para calcular o VaR, mas validam o uso de modelos paramétricos com distribuição normal e simulação Montecarlo em mercados financeiros emergentes.

CLASSIFICAÇÃO JEL: G17, G32, C58

CONTEÚDO: Introdução; 1. Breve revisão da literatura: value at risk; 2. Metodologia; 3. Resultados; Conclusões; Bibliografia.

PALAVRAS-CHAVE: Value at risk; modelos GARCH; simulação; Colômbia

INTRODUCCIÓN

El mercado de capitales es un tipo de mercado financiero en el que se realiza compra y venta de títulos valor. Dicha compra y venta de títulos valor se realiza en contextos de incertidumbre bajo la conocida relación rendimiento-riesgo extendida por Markovitz (1952), establece que la inversión en activos más riesgosos se ve compensada con mayores rentabilidades esperadas. A partir de esta relación, se da lugar a diversas aproximaciones de la medición del riesgo de los activos, que incluyen medir la volatilidad mediante el uso de técnicas, tanto estadísticas como econométricas, que buscan realizar una mejor gestión del riesgo. Una de las medidas más ampliamente utilizadas es el Valor en Riesgo (VaR por sus siglas en inglés).

Según Fan et al., (2008), el Valor en Riesgo (VaR) se convierte en una herramienta esencial para los mercados financieros al ofrecer una medición del riesgo mediante el uso de una distribución de probabilidad de variable aleatoria. Incluso, el Comité de Supervisión Bancaria de Basilea (BCBS por sus siglas en inglés) (2019) propone dicha medida dentro del conjunto de medidas acordadas internacionalmente como respuesta a la crisis financiera de 2008, con el objetivo de reforzar la regulación, supervisión y gestión del riesgo de los bancos. Su cálculo se puede realizar mediante la elaboración de modelos econométricos que resultan muy útiles en la medición de riesgo. Al respecto, Su y Hung, (2011) y Su, Lee, y Chiu (2014) señalan que esto es deseable dado que permite incorporar las características de las series de tiempo financieras, tales como colas anchas, asimetría, clúster de volatilidad y memoria de largo plazo. De hecho, varios autores coinciden en que capturar el comportamiento asimétrico es esencial para un pronóstico más preciso del VaR. Por ejemplo, Wong, Chin y Tan (2016) llegan a este resultado a partir de la comparación de volatilidad de tres índices accionarios. Abounoori, Elmi y Nademi (2016) hallan que los modelos con mejor desempeño con un horizonte de tiempo de un día son AR(2)-MRSGARCH-GED. Por su parte, Cheong (2008) encuentra para Malasia que la distribución de Pareto es la que mejor recoge el comportamiento asimétrico de ambas colas.

En efecto, existen diversos trabajos que abordan de manera empírica la elaboración de modelos en la gestión de riesgo. Algunos de estos se centran en la incorporación de volatilidad a través de diferentes tipos de modelos GARCH, por ejemplo, Mabrouk (2016) utiliza modelos de largo plazo en siete índices accionarios y tres tipos de cambio; Schaeffer et al., (2012) realizan una medición de la volatilidad en acciones del sector petrolero; Laporta, Merlo y Petrella (2018) utilizan diferentes modelos GARCH y cuantiles para el pronóstico del VaR para siete commodities (materias primas o bienes primarios) de energía.

Por su parte, trabajos como el de Youssef, Belkacem y Mokni (2015) confirman para el mercado del crudo y la gasolina que tener en cuenta elementos de memoria de largo plazo, asimetría y colas pesadas es importante para el manejo del riesgo y las estrategias de cobertura; Ewing y Malik (2017) resaltan la importancia de distinguir entre buenas y malas noticias e incorporar cambios estructurales en los modelos GARCH asimétricos para reducir la persistencia en volatilidad de los precios del crudo; Cheong (2009) utiliza modelos ARCH para estimar hechos estilizados de la volatilidad de los dos mayores mercados de crudo (WTI y Brent).

De igual manera, se encuentran estudios que utilizan modelos con aproximaciones diversas. Entre estos, aquellos implementados por Chiu y Chuang (2016) a las propiedades de memoria larga o persistencia en volatilidad de acciones; Bentes (2014) a los índices bursátiles; Tang y Shieh (2006) en tipos de cambio y mercados de futuros; Aloui y Mabrouk (2010) en commodities de crudo y gas; Degiannakis y Potamia (2017) para acciones, commodities y tipos de cambio. Asimismo, existen otros que diferencian entre posición corta y posición larga, por ejemplo, Giot y Laurent (2003) estiman un modelo de memoria larga basado en distribución T; So y Yu (2006) estiman siete modelos GARCH para doce índices de mercado y cuatro tasas de cambio, resaltan la importancia de las colas anchas en la estimación del VaR. Todo esto evidencia la amplia variedad de trabajos que incorporan en sus aplicaciones empíricas el uso de modelos econométricos para el cálculo del VaR.

Sin embargo, la mayoría de estas aproximaciones se realizaron en mercados financieros desarrollados, son pocos los estudios que se centran en la gestión del riesgo en mercados emergentes. Se resalta el trabajo de Yang y Chang (2008) que señala la importancia de evaluar el valor y la estabilidad de la moneda local como parte del análisis de decisión de inversión en el mercado accionario; el de Reboredo, Rivera-Castro y Ugolini (2016) que estructura la dependencia de la relación entre los retornos de acciones y los tipos de cambio como cobertura de un portafolio; aquellos que analizan los índices accionarios y los efectos en los mercados locales en economías emergentes como Brasil, Chile, Colombia, India, México, Rusia, Suráfrica y Turquía (Assaf, 2015; Dimitrakopoulos et al., 2010; Su, 2015; Su y Knowles, 2006; Diamandis et al., 2011; Sadorsky, 2014); así como el de Slim, Koubaa y BenSaida (2017) que encuentra que en el caso de mercados emergentes, el modelo GARCH estándar es la mejor especificación para los pronósticos de estimaciones VaR.

Particularmente para el caso colombiano, son pocos los trabajos que buscan presentar el comportamiento del mercado accionario, y aún más escasos los que hacen uso de modelos econométricos para la gestión del riesgo. De manera puntual, destacan los trabajos de Uribe (2007), quien realiza una comparación con índices de la región y ciertos mercados desarrollados, con el ánimo de fijar el grado de desarrollo local; Rivera (2009) analiza los patrones estacionales de los antiguos índices de las Bolsas de Bogotá (IBB), Medellín (Ibomed) y Colombia (IGBC) y encuentra evidencia que contradice la hipótesis de mercados eficientes; Pino, Uribe y Jiménez (2017) proponen una descripción del mercado de acciones, realizan un análisis de topología de red con el fin de reconocer la importancia institucional; y Delgado y Durango (2018) evalúan la estructuración de portafolios a partir de comparaciones con la volatilidad del índice Colcap que tiene en cuenta que este es el principal indicador bursátil del país, está compuesto por las acciones más transadas del mercado.

No obstante, a pesar de la escasa literatura, se hallaron tres artículos similares a nivel nacional: el primero lo realizaron Ospina y Giraldo (2009), quienes calcularon el VaR para las acciones de mayor bursatilidad del mercado colombiano mediante el método paramétrico, como resultado, encontraron que los modelos más significativos fueron los GARCH asimétricos. El segundo, desarrollado por Kuri, Ojeda y Ovalle (2015), estimó el VaR para la acción de Ecopetrol mediante modelos de heterocedastidad condicionada, encontraron que el tipo de estimación que mejor modela los retornos de esta acción es el GARCH (1,1). El tercero -elaborado por Pérez y Fernández (2006)- considera la volatilidad del Índice General de la Bolsa de Valores de Colombia (IGBC).

En este contexto, esta investigación contribuye a la literatura de análisis de riesgo en mercados emergentes puesto que calcula el valor en riesgo (VaR) para cada una de las veinticuatro acciones de la canasta del índice Colcap que rige de mayo a julio de 2017, todo esto mediante la estimación de modelos econométricos que incorporan la volatilidad, asimetría y demás propiedades de las series financieras para el periodo comprendido entre agosto de 2013 y diciembre de 2017.

Asimismo, en mercados financieros emergentes como el colombiano, resulta de especial relevancia el análisis del riesgo de sus activos, teniendo en cuenta la vulnerabilidad que se presenta en pequeños mercados ante eventos extremos como las crisis financieras -de manera particular la crisis sub-prime de 2008 o la crisis hipotecaria de 1999-, así como las asimetrías de información que conllevan que los inversionistas realicen una toma de decisiones inadecuada.

En consecuencia, efectuar el cálculo del VaR para los activos propuestos y a partir de diversas metodologías, no solo permite identificar las principales características de las acciones del mercado de valores colombiano -y con ello aportar a la modelación de series financieras en mercados emergentes-, sino que también contribuye a situar la discusión de la gestión de riesgo frente a la importancia de contar con diversos métodos de análisis y formas de cálculo, en pos de una toma de decisiones más informada que conduzca a mejores rentabilidades y reducción del riesgo.

Este documento se estructura en cinco secciones que incluyen esta introducción. Posteriormente, se realiza una breve revisión de la literatura existente acerca del valor en riesgo, más adelante se describe la metodología utilizada para la elaboración de los modelos, el cálculo del VaR y el backtesting, seguido por una sección de resultados. La última sección presenta las conclusiones del artículo.

1. BREVE REVISIÓN DE LITERATURA: VALOR EN RIESGO

A menudo, la medición del riesgo de mercado se realiza a partir del VaR (valor en riesgo), el cual Jorion (2000, p. 13) define como "la peor pérdida esperada en un intervalo de tiempo determinado bajo condiciones normales del mercado ante un nivel de confianza dado". Al respecto, Sadegui y Shavvalpour (2006) describen el VaR como un método cuantitativo ampliamente adoptado por instituciones tales como los bancos, gestores de cartera y reguladores, al ser reconocida -según Dowd y Blake (2006)- como la medida más ampliamente utilizada en la práctica y cuyo cómputo resulta fundamental para el cálculo de otras medidas de riesgo basadas en cuantiles, como por ejemplo, el expected shortfall.

Existen diferentes formas de calcular el VaR, estas se pueden clasificar en paramétrica, no paramétrica y simulación Montecarlo. Según Hull y White (1998), el método paramétrico -también conocido como método analítico de la matriz de varianzas/ covarianzas- se basa en la distribución normal de los rendimientos de la cartera. El método no paramétrico se muestra independiente de cualquier distribución de probabilidad y por tal motivo, no se necesita estimar volatilidades ni correlaciones. El método de simulación Montecarlo permite la estimación de manera paramétrica y no paramétrica, ofrece una descripción más realista del riesgo a partir del ajuste de una distribución empírica de los retornos.

Según Cheong (2008), las estimaciones no paramétricas tienen las ventajas de la simplicidad y -lo que es más importante- la no necesidad de un supuesto de distribución específico para la serie de rendimientos, mientras que el cálculo del VaR basado en un modelo paramétrico requiere estimaciones precisas de la volatilidad y un cuantil correspondiente a la distribución empírica. En esta línea, Sadegui y Shavvalpour (2006) afirman que los estudios previos sobre la estimación del VaR basada en modelos se han centrado principalmente en la caracterización de los hechos estilizados que se encuentran con frecuencia en los datos financieros, especifican la estructura variable de la volatilidad y la forma de la distribución del rendimiento de los activos.

Dentro de la estimación paramétrica, el modelo más utilizado es el de RiskMetrics propuesto por J. P. Morgan, el cual asume que los rendimientos de los activos siguen una distribución normal con media cero y varianza condicional en un promedio móvil ponderado exponencialmente (EWMA) de los rendimientos históricos al cuadrado. Sin embargo, para González-Rivera, Lee y Mishra (2004), a pesar de su uso extendido, este modelo parece ser el que peor desempeño presenta. Lo anterior es consistente con el trabajo de Mabrouk y Saadi (2012) en cuanto que argumenta que dicha estimación del VaR es la menos precisa al no tener en cuenta las características más destacadas de las series de tiempo financieras: colas gruesas, asimetría, agrupamiento de volatilidad y memoria larga.

Además, es común que el primer supuesto (el de normalidad) no se presente en un contexto empírico, al tener en cuenta que las características de la distribución del retorno de activos involucran fenómenos de colas pesadas, leptocurtosis y asimetría. Por este motivo, el uso de una distribución condicional que no incluya dichas características puede llevar a un sesgo en la estimación del verdadero riesgo de mercado.

Por ello, estudios previos han hecho un gran esfuerzo para superar estos problemas al emplear conjuntos de distribuciones condicionales. Algunos de los trabajos pioneros son: Bollerslev (1987) con la T-Student, y Angelidis et al., (2004) con la distribución de errores generalizada (GED) de cola pesada, que permite capturar eventos extremos. Al respecto, Mabrouk (2016) añadió que ambas distribuciones son las más ampliamente utilizadas al presentar mejor desempeño al momento de capturar el exceso de curtosis. Sin embargo, las distribuciones GED y T-Student (que son simétricas) no pueden incluir el componente de asimetría presente en las innovaciones.

El segundo supuesto del promedio móvil ponderado exponencialmente para la varianza tampoco suele incorporar suficiente información con respecto a las propiedades empíricas de las dinámicas de volatilidad, como la memoria larga y el efecto de apalancamiento. Si bien para Hartz, Miitnik y Paolella (2006) existen diversos métodos que buscan estar en capacidad de incorporar las dos principales características dentro de los métodos que capturan estas dos características de una manera razonable -y que han sido ampliamente extendidas- se encuentran los modelos de heterocedasticidad condicionada, propuestos por Engle (1982) y generalizados por Bollerslev (1986).

Dichos modelos han proporcionado una buena primera aproximación a estos hechos estilizados, al modelar la estructura autorregresiva en la varianza condicional. En efecto, trabajos como el de Yang y Chang (2008) evidencian que la estimación de la volatilidad en el rendimiento de los activos resulta crucial para la fijación de precios de las opciones, en la previsión de umbrales de valor en riesgo (VaR) y en la selección óptima de la cartera, al encontrar que los modelos GARCH son adecuados para describir datos como clúster de volatilidad y comportamiento de cola gruesa.

Sin embargo, de acuerdo con autores como Andersen et al., (2002) y Gallant et al., (1997), los modelos tipo GARCH presentan alta persistencia y suelen no ser los más adecuados frente a eventos raros en retornos de activos. Estos saltos pueden llevar a una sobreestimación del riesgo de mercado real. De allí que se hayan desarrollado modelos que incorporan propiedades adicionales, como la asimetría, la cual se refiere a los impactos diferenciales de choques positivos y negativos de igual magnitud en la volatilidad. Según Cheong (2008), esta propiedad es muy importante en el análisis de riesgo en el cual las inversiones en posiciones largas y cortas durante un período de tiempo determinado se basan en gran medida en los comportamientos de las colas inferiores y superiores.

Para incorporar esta propiedad, Nelson (1991) desarrolló el modelo exponencial GARCH, EGARCH, que analiza los efectos positivos y negativos ocasionados por las noticias del mercado de valores, con choques negativos que afectan la volatilidad de manera diferente a los choques positivos de igual magnitud. Además, la posibilidad de apalancamiento que permite que los choques negativos aumenten la volatilidad -mientras que los choques positivos disminuyen la volatilidad- está presente en la formulación de los modelos EGARCH. Bajo la misma idea, se han desarrollado modelos asimétricos alternativos de volatilidad, como el umbral GARCH, TGARCH, de Zakoian (1994). Al respecto, Yang y Chang (2008) y Mabrouk (2016) señalan que tales modelos típicamente distinguen entre la asimetría de información, un umbral, y umbrales de valores críticos.

Consecuentemente, validaciones empíricas como la desarrollada por Su y Hung (2011) que resaltan, tanto las colas pesadas en el cálculo del VaR al 99 % en posiciones largas, como la asimetría y choques dinámicos, así como la de Bayer (2018), que diferencia entre ventanas y el uso de técnicas que combinan varios métodos, han mostrado mejoras en la eficiencia.

Por otra parte, la evaluación del pronóstico mayormente generalizada es aquella realizada a partir de backtesting, que el BCBS (2019) define como el proceso utilizado para comparar las ganancias y pérdidas diarias reales e hipotéticas, con las medidas de VaR generadas por el modelo para evaluar el conservadurismo de los sistemas de medición de riesgos.

La aproximación más utilizada es la propuesta por Kupiec (1995), también conocida como la versión incondicional de la prueba de cobertura, la cual se centra en si el número real de excepciones de VaR es igual al número esperado de excepciones, es decir, examina si la tasa de cobertura realizada (ex-post) es igual a la tasa de cobertura teórica. De acuerdo con Dimitrakopoulos et al., (2010) y Laporta et al., (2018), dicha prueba sigue un proceso de Bernoulli con parámetro p en la cual p es igual a la tasa teórica de cobertura del VaR, de esta manera, la variable toma el valor de 1 si la pérdida excede el VaR estimado y 0 en los otros casos.

Por definición, la tasa de fallas supone el número de veces que el retorno excede -en valor absoluto- el pronóstico del VaR. Si el modelo VaR se especifica correctamente, la tasa de fallas debería ser igual al nivel VaR especificado. La tasa de fallas empíricas f se define como una secuencia de observaciones sí/no, de tal manera que es posible probar si el modelo se especifica de manera adecuada. A un nivel del 95 %, un intervalo de confianza para la tasa de fallas se puede construir por:

El estadístico LR es

donde N es el número de violaciones del VaR, T es el número total de observaciones y a es la tasa de fallas teóricas. Bajo la hipótesis nula de que la tasa de fallas es correcta, el estadístico LR se distribuye como un chi-cuadrado con un grado de libertad.

2. METODOLOGÍA

2.1. Datos

El análisis de riesgo realizado para el mercado accionario colombiano utiliza para el cálculo del VaR las principales acciones que componen el índice Colcap, en este índice bursátil participan las acciones más liquidas que cotizan en la Bolsa de Valores de Colombia, el cual debe estar compuesto por mínimo veinte acciones de emisores diferentes.

El periodo de análisis comprende desde el 1° de agosto de 2013 hasta el 28 de diciembre de 2017, con un total de 1.074 observaciones para cada una de las 24 acciones que conforman la canasta del índice Colcap de mayo a julio de 2017. El periodo de estimación elegido busca tener la misma cantidad de observaciones por acción, de allí que considere el hecho de que la acción de Avianca Holdings cambió su nemotécnico el 19 de julio de 2013 de PFAVTA a PFAVH.

Los precios de cierre diarios de las acciones se obtienen de la Bolsa de Valores de Colombia (BVC). Los retornos diarios se calculan como la diferencia logarítmica en el precio de cierre diario, es decir, R t = log(Pr) - log(Pr-1) donde P t es el precio de cierre diario. La tabla 1 describe las empresas que componen el índice Colcap para el periodo de análisis.

Tabla 1 Empresas 

Fuente: elaboración propia a partir de las páginas oficiales de las empresas.

2.2. Estimación de modelos

Como se describió en la revisión de literatura, una mejor aproximación al riesgo debe incluir la estimación de modelos que incluyan características de las series financieras como clúster de volatilidad, asimetría y apalancamiento. Debido a la naturaleza de las series financieras presente en los retornos diarios, la estimación de los modelos contempló dos fases. En primer lugar, la definición de una ecuación para la media, habitualmente estimada mediante modelos autorregresivos de media móvil; en segundo lugar, la estimación de la varianza condicionada, la cual se puede dar mediante modelos de heterocedasticidad condicionada GARCH, E-GARCH y TARCH.

Para la estimación de dichos modelos se tuvo en cuenta que las series de retornos no se distribuyen de manera normal, por lo que se hizo uso de las distribuciones condicionales mayormente utilizadas, de acuerdo con lo propuesto en la literatura: la distribución T-Student y la distribución del error generalizado (GED). La primera es leptocúrtica, lo que significa que la parte central de la distribución es más delgada y tiene colas más anchas que la distribución normal, mientras que la segunda se extiende exponencialmente en la cola de la distribución (Castaño, Gómez y Gallón, 2008; Mendoza y Galvanovskis, 2014).

2.3. Cálculo del VaR

El cálculo del VaR se realizó bajo tres metodologías: no paramétrico, paramétrico y simulación Montecarlo. Las estimaciones se efectuadron con los softwares E-views®, y Risk Simulator®, respectivamente. Como su nombre lo indica, la simulación Montecarlo crea diversos modelos de posibilidades al calcular una y otra vez resultados (10.000 simulaciones) que comparan cada vez un grupo diferente de valores aleatorios de las funciones de probabilidad utilizadas y selecciona el mejor modelo de acuerdo con medidas de bondad de ajuste.

Una vez se obtuvieron los modelos para cada una de las series de retornos de los activos seleccionados, se continuó con el cálculo del VaR paramétrico a partir de la fórmula:

Donde µ i,t es el retorno condicional esperado del activo i, Z m,α corresponde al cuantil izquierdo α de la distribución usada en la serie de los residuales del modelo estimado en cada caso, y h it es la varianza condicional del activo i. Todos los cálculos se hicieron con un nivel de confianza del 95 % lo que según Jorion (2000), indica que el VaR se encuentra en el 5 % del lado izquierdo de la cola del histograma de cada una de las series.

Para determinar las acciones más riesgosas, cada medida de riesgo VaR se compara de acuerdo con la metodología empleada, a partir de la prueba de backtesting realizada, con el fin de validar la eficacia de los modelos implementados para la gestión de riesgo.

3. RESULTADOS

La tabla 2 presenta las ecuaciones de los modelos estimados para cada uno de los retornos de los precios de cierre de las 24 acciones que componen el índice Colcap. Algunos de estos retornos no requirieron ecuación para la media (Rt=0), pero todos incorporaron componentes de volatilidad bajo la especificación de modelos GARCH con distribuciones no-normales o con propiedades de asimetría.

Tabla 2 Modelos estimados 

Fuente: elaboración propia

La acción CLH presenta uno de los modelos más sencillos de los veinticuatro estimados, el cual es el ARCH. Por otra parte, veintiún modelos estimados corresponden al GARCH que es el modelo más usado para modelar series financieras. Asimismo, se puede evidenciar que la mayoría de los modelos son GARCH (1,1), lo que concuerda con los hallazgos de Slim, Koubaa y BenSaida (2017), quienes señalan para el caso de mercados emergentes que, el modelo GARCH estándar es la mejor especificación para los pronósticos de estimaciones VaR.

Las dos acciones restantes, que corresponden a Conconcreto y al Banco de Bogotá, requieren la estimación de modelos E-GARCH y TARCH respectivamente, los cuales corrigen los problemas que presentan los GARCH al incorporar el componente de asimetría. Lo anterior coincide con la importancia de estimar modelos que incluyan características propias de las series financieras y está ampliamente expuesto en la literatura (So y Yu, 2006; Yang y Chang, 2008; Youssef et al., 2015).

De otra parte, los resultados del cálculo del VaR se presentan en la tabla 3. Es de anotar que el método paramétrico implementado utiliza tanto las distribuciones de los modelos de heterocedasticidad condicionada estimados, como la normal y la simulación Montecarlo.

Tabla 3 Comparación del VaR 

Fuente: elaboración propia.

Del cálculo del VaR tradicional que hace uso de la distribución normal se resalta que las acciones que presentan mayor riesgo son: CNEC (5,41 %), PFAV (3,66 %) y ECOPETROL (3,38 %). En contraste, con el VaR paramétrico de los modelos de heteros-cedasticidad condicionada estimados -los cuales asumen distribuciones T-Student o GED- las acciones más riesgosas son: CLH (8,42 %), ETB (6,71 %) y NUTRESA (5,99 %). Por su parte, la tercera metodología paramétrica o simulación Montecarlo que asumió otras distribuciones como: Laplace, Gumbel mínimo, logística, parabólico y normal, evidenció como acciones más riesgosas a: CNEC (5,19 %), ECOPETROL (3,37 %) y CONCONCRET (3,08 %). Por último, según el VaR no paramétrico las acciones más riesgosas son: CNEC (4,62 %), ECOPETROL (3,22 %) y EXITO (3,05 %).

Si bien, hay diferencias en los riesgos evidenciados a partir del cálculo del VaR con las metodologías utilizadas, acciones como ECOPETROL y CNEC resultan las más riesgosas en tres de las cuatro metodologías utilizadas, lo que sugiere que habría que tratarse con cuidado su incorporación en decisiones de inversión. Por el contrario, acciones como NUTRESA o EEB se consideran menos riesgosas que las demás y por tanto, más atractivas para la conformación de carteras.

Para verificar cual metodología se ajusta mejor, se emplea la probabilidad arrojada por el backtesting, la cual nace a partir de: P en días, que es el número de fallas que soporta el modelo en días; M no P, que son las fallas que presentan cada una de las metodologías; y PROB, que es la probabilidad, la cual no puede ser inferior al 5 %, debido a que se tomó el 95 % como nivel de confianza para realizar esta prueba (tabla 4).

Tabla 4 Backtesting: modelos que no se rechazan 

Fuente: elaboración propia

Como resultado del backtesting tenemos que en el VaR no paramétrico el único valor que se rechaza es el de CLH. Por el contrario, en el VaR paramétrico -tanto de distribución normal como de simulación Montecarlo- no se rechaza ningún modelo. Finalmente, de la totalidad de modelos econométricos estimados se rechazan 15 debido a que no logran ajustar la volatilidad, lo que implica que dichos modelos subestiman el VaR. Este hallazgo refuerza la importancia de realizar análisis de gestión de riesgo mediante la incorporación de la mayor cantidad de propiedades de las series financieras recogidas, por ejemplo con la posibilidad de combinar técnicas y modelos. En este sentido, el uso de simulación Montecarlo resulta ser el método paramétrico más adecuado en cuanto su implementación considera el mejor ajuste entre varias distribuciones y simulaciones.

CONCLUSIONES

Este documento investigó los comportamientos de cola de las distribuciones de innovación para las veinticuatro principales acciones en Colombia. Para ello, se estimaron modelos ARMA que buscaban explicar los componentes de la media de los retornos, así como modelos de heterocedasticidad condicionada que incorporaran elementos propios de las series financieras. A pesar de que el modelo T-Student tipo GARCH es capaz de capturar la propiedad de cola pesada en la mayoría de los activos, no tuvo en cuenta el comportamiento asimétrico en dos de los casos, en los que resulta necesario la inclusión de modelos E-GARCH y TARCH.

Si bien en la gestión de riesgo el VaR sigue siendo la medida más utilizada por los profesionales, ya que al implementarla periódicamente y tener una metodología estructurada se pueden evitar desastres financieros, aún es difícil responder a la pregunta que ha suscitado un gran interés en la literatura financiera y que se refiere a cuál es el modelo más apropiado para pronosticar la volatilidad de los rendimientos de los activos, principalmente en mercados financieros emergentes.

Los resultados hallados para las principales acciones del mercado colombiano corroboran la dificultad de elegir un único modelo para el cálculo del VaR, pero dan luces acerca de la necesidad de incorporar en las estimaciones características propias de las series financieras, tales como volatilidad, asimetría, apalancamiento, entre otras. Esto solo se puede obtener por medio del uso de modelos econométricos o simulaciones particulares que diferencien entre activos y también entre distribuciones. Añádase a esto la validez del uso del método paramétrico de distribución normal, ampliamente utilizado de manera empírica.

Una vez calculado el VaR mediante las diferentes metodologías se evidencia que la simulación Montecarlo es la más acertada para valorar el riesgo accionario colombiano, debido a que sus resultados provienen de la incorporación de pruebas de bondad de ajuste que permiten seleccionar entre diversas distribuciones conocidas y frente a una amplia cantidad de repeticiones. El uso de econometría es igualmente válido en la medida en que se realice una comparación entre modelos que conlleve elegir el mejor, es decir aquel que efectivamente recoja las características de las series.

Finalmente, a partir de los resultados encontrados se llega a la conclusión de que al incorporar más información, el cálculo del VaR es más acertado, por lo que resulta fundamental realizar su cálculo con las diferentes metodologías de análisis existentes para obtener un panorama más completo y de esta manera mejorar la toma de decisiones de inversión.

REFERENCIAS

Abounoori, Esmaiel, Elmi, Zahra, y Nademi, Younes. (2016). Forecasting Tehran stock Exchange volatility; Markow switching GARCH approach, En: Physica A, N°. 445, p. 264-282. Doi: 10.1016/j.physa.2015.10.024 [ Links ]

Aloui, Chaker y Mabrouk, Samir. (2010). Value-at-risk estimations of energy commodities via long-memory, asymmetry and fat-tailed GARCH models. En: Energy Policy, N°. 38, p. 2326-2339. Doi:10.1016/j.enpol.2009.12.020 [ Links ]

Andersen, Torben, Benzoni, Luca y Lund, Jesper. (2002). An empirical investigation of continuous-time equity return models. En: Journal of Finance, Vol. 57, N°. 3, p. 1239-1284. Doi: 10.1111/1540-6261.00460 [ Links ]

Angelidis, Timotheos, Benos, Alexandros y Degiannakis, Stavros. (2004). The use of GARCH models in VaR estimation. En: Statistical Methodology, N°. 1, p. 105-128. [ Links ]

Assaf, Ata. (2015). Value-at-Risk analysis in the MENA equity markets: Fat tails and conditional asymmetries in return distributions. En: Journal of Multinational Financial Management, N°. 29, p. 30-45. Doi: 10.1016/j.mulfin.20T14.11.002 [ Links ]

Bayer, Sebastian. (2018). Combining Value-at-Risk forecasts using penalized quantile Regressions, En: Econometrics and Statistics, N°. 8, p. 56-77. [ Links ]

Bentes, Sónia. (2014). Measuring persistence in stock market volatility using the FIGARCH approach. En: Physica A, N°. 408, p. 190-197. Doi: 10.1016/j.physa.2014.04.032 [ Links ]

Bollerslev, Tim. (1986). Generalized autoregressive conditional heteroskedasticity. En: Journal of Econometrics, N°. 31, p. 307-327. Doi: 10.1016/0304-4076(86)90063-1 [ Links ]

Bollerslev, T. (1987). A conditionally heteroscedastic time series model for speculative prices and rates of return. En: Review of Economics and Statistics, N°. 69, p. 542-547. [ Links ]

Castaño, Elkin,. Gómez, Karoll y Gallón, Santiago (2008). Pronóstico y estructuras de volatilidad multiperíodo de la tasa de cambio del peso colombiano. En: Cuadernos de Economía, Vol. XXVII, No. 48, p. 241-266 [ Links ]

Cheong, Chin Wen. (2008). Heavy-tailed value-at-risk analysis for Malaysian stock exchange. En: Physica A, N°. 387, p. 4285-4298. Doi: 10.1016/j.physa.2008.01.075 [ Links ]

Cheong, Chin Wen. (2009). Modeling and forecasting crude oil markets using ARCH-type models. En: Energy Policy, N°. 37, p. 2346-2355. Doi: 10.1016/j.enpol.2009.02.026 [ Links ]

Chiu, Yen-Chen y Chuang, I-Yuan. (2016). The performance of the switching forecast model of value-at-risk in the Asian stock markets. En: Finance Research Letters, N°. 18, p. 43-51. Doi: 10.1016/j.frl.2016.03.019 [ Links ]

Comité de Supervisión Bancaria de Basilea (BCBS) (2019). Minimum capital requirements for market risk. Suiza: BCBS. [ Links ]

Degiannakis, Stavros y Potamia, Artemis. (2017). Multiple-days-ahead value-at-risk and expected shortfall forecasting for stock indices, commodities and exchange rates: Inter-day versus intra-day data. En: International Review of Financial Analysis, N°: 49, p. 176-190. Doi: 10.1016/j.irfa.201610.008 [ Links ]

Delgado, Luis David y Durango, María Patricia (2018). Estructuración de un portafolio de inversiones con acciones colombianas, En: Semestre Económico, Vol. 21, N°. 46, p. 167-183. Doi: 10.22395/seec.v21n46a7. [ Links ]

Diamandis, Panayiotis, Drakos, Anastassios, Kouretas, Georgios y Zarangas, Leonidas. (2011). Value-at-risk for long and short trading positions: Evidence from developed and emerging equity markets. En: International Review of Financial Analysis, N°. 20, p. 165-176. Doi: 10.1016/j.irfa.201L02.009 [ Links ]

Dimitrakopoulos, Dimitris, Kavussanos, Manolis y Spyrou, Spyros. (2010). Value at risk models for volatile emerging markets equity portfolios. En: The Quarterly Review of Economics and Finance, N°. 50, p. 515-526. Doi: 10.1016/j.qref.2010.06.006 [ Links ]

Dowd, Kevin y Blake, David. (2006). After Var: the theory, estimation and insurance applications of quantile-based risk measures. En: The Journal of Risk and Insurance, Vol. 73, N°. 2, p. 193-229. Doi: 10.1111/j.1539-6975.2006.00171.x [ Links ]

Engle, Robert. (1982). Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation. En: Econométrica, Vol. 50, N°. 4, p. 987-1007. [ Links ]

Ewing, Bradley y Malik, Farroq. (2017). Modelling asymmetric volatility in oil prices under structural breaks. En: Energy Economics, N°. 63, p. 227-233. Doi: 10.1016/j.eneco.2017.03.001 [ Links ]

Fan, Ying, Zhang, Yue-Jung, Tsai, Hsien-Tang, y Wei, Yi-Ming. (2008). Estimating 'Value at Risk' of crude oil price and its spillover effect using the GED-GARCH approach. En: Energy Economics, N°. 30, p. 3156-3171 Doi: 10.1016/j.eneco.2008.04.002 [ Links ]

Gallant, Ronald, Hsieh, David y Tauchen, George. (1997). Estimation of stochastic volatility models with diagnostics. En: Journal of Econometrics, N°. 81, p. 159-192. [ Links ]

Giot, Pierre y Laurent, Sébastien. (2004). Modelling daily Value-at-Risk using realized volatility and ARCH type models. En: Journal of Empirical Finance, N°. 11, p. 379-398. Doi: 10.1016/j.jempfin.2003.04.003 [ Links ]

Gonzaléz-Rivera, Gloria, Lee, Tae-Hwy, y Mishra, Santosh. (2004). Forecasting volatility: A reality check based on option pricing, utility function, value-at-risk, and predictive likelihood. En: International Journal of Forecasting, N°. 20, p.629--645. Doi: 10.1016/j.ijforecast.2003.10.003 [ Links ]

Hartz, Christoph, Mittnik, Stefan, y Paolella, Marc. (2006). Accurate value-at-risk forecasting based on the normal-GARCH model, En: Computational Statistics and Data Analysis, N°. 51, p. 2295-2312. [ Links ]

Hull, John, White, Alan. (1998). Incorporating volatility updating into the historical simulation method for VaR. En: The Journal of Risk, N°. 1, p. 5-19. [ Links ]

Jorion, P. (2000). Valor en Riesgo: El nuevo paradigma para el control de riesgos con derivados. México, D.F.: Limusa. [ Links ]

Kupiec, Paul. (1995). Techniques for verifying the accuracy of risk measurement models. En: The Journal of Derivatives, Vol. 3, N°. 2, p. 73-84. Doi: 10.3905/jod.1995.407942. [ Links ]

Kuri, K., Ojeda, C., y Ovalle, D. (2015). Modelación de la serie de retornos diarios de la acción de ecopetrol en el periodo: 27/noviembre/2007-25/noviembre/2013. En: Heurística, N°. 17, p. 51-57. [ Links ]

Laporta, Alessandro, Merlo, Luca y Petrella, Lea. (2018). Selection of Value at Risk Models for Energy Commodities, En: Energy Economics, N°. 74, p.628-643. Doi: 10.1016/j.eneco.2018.07.009 [ Links ]

Mabrouk, Samir. (2016). Forecasting daily conditional volatility and h-step-ahead short and long Value-at-Risk accuracy: Evidence from financial data. En: The Journal of Finance and Data Science, N°. 2, p. 136-151. Doi: 10.1016/j.jfds.2016.06.00T [ Links ]

Mabrouk, Samir y Saadi, Samir. (2012). Parametric Value-at-Risk analysis: Evidence from stock indices. En: The Quarterly Review of Economics and Finance, N°. 52, p. 305-321. Doi: 10.1016/j.qref.2012.04.006 [ Links ]

Markovitz, Harry. (1952). Portfolio selection. En: The Journal of Finance, Vol. 7, N°. 1, p. 77-91. [ Links ]

Mendoza, Alfonso y Galvanovskis, Evalds (2014). La cópula GED bivariada. Una aplicación en entornos de crisis. En: El Trimestre Económico, Vol. 81, No. 323, p.721-746. [ Links ]

Nelson, Daniel. (1991). Conditional heteroscedasticity in asset returns: a new approach, En: Econometrica, N°. 59, p. 347-370. [ Links ]

Ospina D'Aleman, Federico y Giraldo Sánchez, David. (2009). Aplicación de los modelos GARCH a la estimacion del VaR de acciones colombianas. En: Soluciones de Postgrado EIA, N° 3. p. 11-24. [ Links ]

Pérez, Fredy y Fernández, Horacio. (2006). Análisis de la volatilidad del índice general de la bolsa de valores de Colombia utilizando modelos ARCH, En: Revista Ingenierías Universidad de Medellín, Vol. 5, N°. 8, p. 13-33. [ Links ]

Pino, Anderson, Uribe, Jorge Mario, y Jiménez, Diana. (2017). Relevancia de los inversionistas institucionales en el mercado accionario colombiano. Semestre Económico, Vol. 20, n°. 44, p. 45-65. Doi: 10.22395/seec.v20n44a3. [ Links ]

Reboredo, Juan, Rivera-Castro, Miguel y Ugolini, Andrea. (2016). Downside and upside risk spillovers between exchange rates and stock prices. En: Journal of Banking & Finance, N°. 62, p.76-96. [ Links ]

Rivera, David. (2009). Modelación del efecto del día de la semana para los índices accionarios de Colombia mediante un modelo STAR GARCH. En: Revista de Economía del Rosario, Vol. 12 N°. 1, p. 1-24. [ Links ]

Sadegui, Mehdi y Shavvalpour, Saeed. (2006). Energy risk management and value at risk modeling, En: Energy Policy, N°. 34, p. 3367-3373. [ Links ]

Sadorsky, Perry. (2014). Modeling volatility and correlations between emerging market stock prices and the prices of copper, oil and wheat. En: Energy Economics, N°. 43, p. 72-81. Doi: 10.1016/j.eneco.2014.02.014 [ Links ]

Schaeffer, Roberto, Borba, Bruno, Rathmann, Régis, Szklo, Alexandre y Castelo-Branco, David. (2012). Dow Jones sustainability index transmission to oil stock market returns: A GARCH approach, En: Energy, N°. 45, p.933-943. [ Links ]

Slim, Skander, Koubaa, Yosra y BenSaida, Ahmed. (2017). Value-at-Risk under Lévy GARCH models: Evidence from global stock markets. En: Journal of International Financial Markets, Institutions & Money, N°. 46, p. 30-53. Doi: 10.1016/j.intfin.2016.08.008 [ Links ]

So, Mike y Yu, Philip. (2006). Empirical analysis of GARCH models in value at risk estimation. International. En: Financial Markets, Institutions & Money, N°. 16, p. 180-197. [ Links ]

Su, Ender y Knowles, Thomas. (2006). Asian Pacific stock market volatility modeling and value at risk analysis. En: Emerging Markets Finance and Trade, N°. 42, p. 18-62. [ Links ]

Su, Jung Bin. (2015). Value-at-risk estimates of the stock indices in developed and emerging markets including the spillover effects of currency market. En: Economic Modelling, N°. 46, p. 204-224. Doi: 10.1016/j.econmod.2014.12.022 [ Links ]

Su, Jung-Bin y Hung, Jui-Cheng. (2011). Empirical analysis of jump dynamics, heavy-tails and skewness on value-at-risk estimation. En: Economic Modelling, N°. 28, p. 1117-1130. [ Links ]

Su, Jung-Bin, Lee, Ming-Chih, y Chiu, Chien-Liang (2014). Why does skewness and the fat-tail effect influence value-at-risk estimates? Evidence from alternative capital markets. En: International Review of Economics and Finance, N°.31, 59-85. Doi: 10.1016/j.iref.2013.12.001 [ Links ]

Tang, Ta Lun y Shieh, Shwu-Jane. (2006). Long memory in stock index futures markets: A value--at-risk approach. En: Physica A, N°: 366, p. 437-448. Doi: 10.1016/j.physa.2005.10.017 [ Links ]

Uribe, Jorge Mario. (2007). Caracterización del mercado accionario colombiano, 2001-2006: un análisis comparativo. En: Borradores de Economía, N°. 456, p. 1-35. [ Links ]

Wong, Zhen Yao, Chin, Wen Cheon y Tan, Siow Hooi. (2016). Daily value-at-risk modeling and forecast evaluation: The realized volatility approach. En: The Journal of Finance and Data Science, Vol. 2, N°. 3, p.171-187. [ Links ]

Yang, Yung-Lieh y Chang, Chia-Lin. (2008). A double-threshold GARCH model of stock market and currency shocks on stock returns, En: Mathematics and Computers in Simulation, N°. 79, p. 458-474. [ Links ]

Youssef, Manel, Belkacem, Lotfi y Mokni, Khaled. (2015). Value-at-Risk estimation of energy commodities: A long-memory GARCH-EVT approach. En: Energy Economics, N°. p. 51, 99-110. Doi: 10.1016/j.eneco.2015.06.010 [ Links ]

Zakoian, Jean Michel. (1994). Threshold heteroskedastic models, En: Journal of Economic Dynamics and Control, Vol. 18, N°. 5, p. 931-955. Doi: 10.1016/0165-1889(94)90039-6 [ Links ]

* El presente artículo es resultado de la investigación desarrollada en el marco del proyecto Modelación de los mercados financieros de Colombia, adscrito al Grupo de Investigación en Estudios Sociales, Financieros e Internacionales (ESFI), clasificado en categoría C en Colciencias (2018) y financiado por la Universidad de la Salle. Se realizó entre abril de 2018 y noviembre de 2019. El artículo se basado en el trabajo de grado obligatorio para optar al título de Profesional en Finanzas y Comercio Internacional de la Universidad de La Salle.

Recibido: 30 de Noviembre de 2019; Aprobado: 12 de Diciembre de 2019

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons