Extreme Value Theory
Extreme Value Theory
Extreme Value Theory
Seminario de Titulo
La medida de riesgo mas popular es Value at Risk (VaR), la cual ha sido criticada
fuertemente en el ultimo tiempo por no ser una medida real, debido a que supone
normalidad en los retornos. Dado este problema nosotros desarrollamos Extreme Value
Theory, para el calculo de VaR y Expected Shortfall (medida de riesgo para eventos
extremos), en Matrix Laboratory (Matlab), específicamente para el índice Nasdaq,
dándonos dos escenarios de distribución de los retornos, uno de normalidad y otro de DGP
(distribución general de pareto). Aquí se aplican los procedimientos planteados en el
trabajo realizado por Manfred Gilli y Evis Këllezi, para el calculo de VaR y ES. De lo
cual mostramos que el cálculo de VaR asumiendo normalidad en la distribución de los
retornos subestima la pérdida potencial para un activo financiero con respecto al VaR
obtenido bajo una DGP.
3
Introducción
El presente trabajo aplica una nueva teoría acerca de la medición de riesgo, que se
refiere al modelamiento de valores extremos. Hasta hoy, la medida de riesgo mas popular
es el VaR (Value at Risk), teoría que revolucionó el mercado financiero, pero que tiene
algunas fallas importantes que nos hacen pensar que su calculo no es exacto, lo que puede
llevar a grandes pérdidas en el mercado financiero.
Aquellos agentes que suponen normalidad son aquellos que creen que todas las
cosas que pasan en el mundo son normales, por lo tanto siguen la distribución de
probabilidad Normal.
Por el contrario los que no creen en la normalidad creen que cada fenómeno tiene
su propia distribución, y que suponer una distribución normal para todos los fenómenos, si
bien tiene a favor la simpleza de su cálculo, no entrega resultados acertados sino
aproximaciones que pueden sobreestimar o subestimar ganancias o pérdidas para los
inversionistas.
Por otra parte, durante la historia han existido diversos eventos extremos. Nos
referimos a ellos como episodios de la historia en que se han producido pérdidas extremas
en el mercado financiero. Ejemplos de eventos extremos existen por montones, entre ellos
destacamos las crisis financieras del año 1929, Long Term Capital Management, los
ataques terroristas del 11 de Septiembre del 2001 y diversos huracanes que han azotado las
costas de Estados Unidos, entre otros.
Con la información anterior, se nos plantea una duda, como poder calcular las
pérdidas máximas posibles en cierto periodo de tiempo, con un nivel de confianza dado y
1
Analizado en el trabajo “Implementación del Value at Risk Condicional” de Dr. Rafael Romero M y Dr.
Sigifredo Leangle S
4
aplicando la distribución de retornos que corresponda. Además de una teoría que sea capaz
de cuantificar las perdidas máximas ante eventos extremos.
Ante esta pregunta nos damos cuenta de la existencia de Extreme Value Theory
(EVT), teoría que estudia la inquietud planteada en el párrafo anterior y que será
desarrollada en profundidad en este trabajo.
El EVT tiene sentido de ser calculado cuando la distribución de retorno tiene colas
gruesas, es decir, que hay más eventos extremos que en la distribución normal.
2
Wikipedia
5
Explicación y cálculo del VaR (Value at Risk)
Como primer acercamiento al concepto de VaR, podemos decir que éste es “la
máxima pérdida esperada en un periodo de tiempo y con un nivel de confianza dados, en
condiciones normales de mercado”3
Una definición mas formal que encontramos del VaR es la de Sharpe (1995)4, que
dice, dada una cartera P, un periodo temporal T y un nivel de probabilidad Q, se estima un
nivel de pérdidas L*, tal que existe una probabilidad Q de que las pérdidas efectivas L,
sean iguales o menores que L* durante el periodo T. A éste nivel de perdidas se le
denomina el Valor en Riesgo de una cartera (VaR)5. Formalmente:
Prob [L* ≥ L] = Q
Para la estimación del VaR hay dos temas que son claves para su cálculo, tanto
para posiciones lineales como no lineales. Estas son la validez de la hipótesis de
normalidad para el valor o rendimiento de un activo o una cartera (que incluye el estudio
del tema de las colas gruesas como el tema de la volatilidad variable con tendencia a
agruparse a lo largo del tiempo –volatility clustering-) y la hipótesis de no linealidad de los
activos para los cuales se quiere estimar el VaR.
Las medidas como VaR son utilizadas en diversos campos, como por ejemplo, para
medir el grado de exposición al riesgo de una operación de una unidad de negocio o de
varios agregados, y para comparar el riesgo de éstas, para cargar a cada posición el coste
3
Definición de Phillipe Jorion, ver Jorion 2000
4
Definición citada en el libro “La Gestión de Riesgo Financiero y Crédito” de Juan Ignacio Peña.
5
La empresa XXX tuvo beneficios diarios de US$ 7,6 MM. Sin embargo, el 5% de los días se observaron
pérdidas de US$ 11, 4 MM o superiores. Por tanto, podríamos decir que, si la distribución de pérdidas y
ganancias se establece en el tiempo, la máxima pérdida que podríamos esperar en 19 de cada 20 días es de
US$ 11,4 MM. Esta pérdida máxima esperada es el VaR para el nivel de confianza (95%) elegido. Usando la
formula expuesta: P[-US$ 11,4 MM ≥ L] =0,95
6
de capital adecuado a su valor de mercado y a su riesgo, protege a la empresa de los costes
de insolvencia, entre otros.
7
Desventajas: Sobreestima el VaR para niveles bajos de confianza, y lo
subestima para niveles altos de confianza (inconveniente derivado de la
suposición de normalidad en los retornos de la cartera y no capturar el
fenómeno de las colas gruesas). Otra desventaja es que la hipótesis de
linealidad condiciona a este método a ser aplicable a carteras lineales, en un
mundo donde los activos no lineales toman cada día mas fuerza. Por último,
incluso ampliando la aproximación del valor de la cartera a una cuadrática,
no se logra una buena precisión en la estimación del VaR de carteras no
lineales. Teniendo en cuenta que la aplicación complica el cálculo de VaR,
disminuyendo de esa forma una de las ventajas más claras que tiene el
método.
8
permite capturar el fenómeno de las colas gruesas, eliminando al mismo
tiempo, la necesidad de trabajar con volatilidades y correlaciones.
La simulación de Monte Carlo ofrece una descripción más realista del
riesgo ya que la distribución de las variaciones en los precios refleja un
abanico completo de todas las realizaciones y sus posibilidades.
Ventajas: Las ventajas fueron mencionadas anteriormente, y ubican a este
método de simulación por sobre los otros dos métodos analizados
anteriormente.
Desventajas: La mayor desventaja de este método es su lentitud, pero este
inconveniente se va solucionando rápidamente a través del tiempo debido al
desarrollo informático que se esta teniendo.
9
sobre que método es mejor, e incluso si es mejor o no utilizar estimaciones
implícitas.
5. Instrumentos no lineales: Los derivados presentan una no linealidad que dificulta
el cálculo del VaR. No existe ningún tipo de convención que establezca de que
forma se debe calcular el VaR en carteras con estas posiciones.
Dado que el VaR se centra en el estudio de los casos extremos (las colas), es
evidente entonces, que el hecho de que existan las colas gruesas complica el estudio y el
cálculo del VaR.
10
financieros reales sabemos que no es así, por lo que hay que estudiar cual puede ser la
frecuencia con que estos ocurran.6
Antes de entrar en la Teoría del Valor Extremo (EVT), debemos recordar que el
VaR se definía como la pérdida bajo situaciones normales de mercado, y que por lo tanto
los acontecimientos extremos (crisis de 1987, 1992, 1994, 1998, etc.) deben ser analizados
mediante instrumentos nuevos, entre los cuales destacamos EVT, que tiene como fin
calcular las pérdidas posibles bajo situaciones extremas de mercado.
6
Existe evidencia de que es así en Chew (1994), página 64, donde movimientos de tres desviaciones
estándar se registraron una mayor cantidad de veces que las que determina la distribución normal (0,3% de
veces)
11
acudir a la EVT, según la cual, la forma precisa de la distribución de las colas de los
rendimientos es perfectamente conocida y modelizable, y lo que es muy importante,
independiente de la distribución de probabilidad de los rendimientos. Por lo que, según
esta teoría, una única ley de probabilidad describe el comportamiento de las colas de casi
todas las distribuciones de probabilidad.
x
e
( x) e ,x
0, x 0
(x) 0
x
e ,x 0
( x)
e ,x 0
(x) 0
1, x 0
7
Fisher-Tippett 1928
8
Estas distribuciones son de densidad.
12
En el año 1995, Jenkinson propone la formula general que engloba a las tres
distribuciones anteriores, creando la Distribución Generalizada de los Valores Extremos:
1/
(1 X )
e si 0
H (x)
x
e e si 0
1
Obtenemos la función de distribución Frechet.
1
Si Obtenemos la función de distribución Weibull.
0 Obtenemos la función de distribución Gumbel.
13
El noventa por ciento de los eventos causan el diez por ciento de las pérdidas, y
el resto del diez por ciento de los eventos causa el noventa por ciento de las pérdidas.
(Una variación de la regla del noventa por ciento utilizada en el marketing)
Expected Shortfall
Constituye una medida de riesgo, la cual estima la pérdida por sobre el VaR,
formalmente sería:
ES P E( X / X VARP )
14
Aplicación de la Teoría del Valor Extremo al índice de precios accionarios NASDAQ:
En esta parte del trabajo calcularemos y compararemos las dos formas que
anteriormente mencionamos para el cálculo del VaR, estas son suponiendo normalidad en
la distribución de los retornos, y suponiendo que se distribuyen de una forma distinta, para
la cual ocuparemos la función de Distribución General de Pareto, además del Expected
Shortfall calculado para el índice.
Los datos usados para los cálculos en ambos casos fueron los retornos diarios de
índice NASDAQ entre el 28 de Octubre de 1996 hasta el 26 de Octubre de 2006, lo que
corresponde a una muestra de 2.517 datos diarios.
En los siguientes gráficos podemos observar la evolución de los precios del índice en
cuestión, así como la evolución de las rentabilidades de éste.
NASDAQ
6000
Valor del índice
5000 NASDAQ
4000
3000
2000
1000
0
10/29/96
10/29/97
10/29/98
10/29/99
10/29/00
10/29/01
10/29/02
10/29/03
10/29/04
10/29/05
Fecha
En este gráfico observamos la evolución de los precios del índice NASDAQ, como se puede
observar, hay una caída fuerte entre los años 1999 y 2000, donde el índice vuelve a los valores que
tenía el año 1996.
15
Rentabilidad del NASDAQ
14%
Rentabilidad del NASDAQ
9%
Rentabilidad
4%
-1%
-6%
10
10
10
10
10
10
10
10
10
10
/2
/2
/2
/2
/2
/2
/2
/2
/2
/2
9/
9/
9/
9/
9/
9/
9/
9/
9/
9/
-11%
96
97
98
99
00
01
02
03
04
05
Fechas
En este grafico observamos las rentabilidades diarias del índice entre Octubre de 1996 y Octubre de
2006.
6
Exponential Quantiles
-10 -5 0 5 10 15
Ordered Data
16
Histograma
400
350
300 Frecuencia
Frecuencia
250
200
150
100
50
...
00
01
02
03
04
05
06
07
07
08
09
10
11
12
y 0.13
0
6
5
1
.1
.0
.0
.0
.0
.0
.0
.0
.0
.0
or
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
-0
-0
-0
-0
-0
-0
-0
-0
-0
-0
ay
m
Clase
En la figura anterior observamos el histograma con las frecuencias de las variaciones diarias del
índice, como se puede observar, existen observaciones extremas, las que, en conjunto con el gráfico
anterior nos llevar a afirmar la existencia de colas gruesas en la distribución de retornos del
NASDAQ.
Como fue mencionado anteriormente, el calculo del VaR fue hecho tanto para la
cola izquierda (posición larga en el activo), como para la como para la cola derecha
(posición corta en el activo)
17
o Primer Caso “Posición larga en el activo”
Media 0,0004
Volatilidad 0,0182
X rentabilid ad . promedio
X Z 99% * , donde Z 99% valor.normal 2.3263
volatilida d
Media -0,00043
Volatilidad 0,02
9
El ultimo dato de la serie es 2379,1 correspondiente al 26 de Octubre de 2006, por lo tanto los VaR
calculados son para el 27 de Octubre de 2006
18
Luego aplicamos el mismo intervalo de confianza anterior, el cual nos arrojo un
VaR al 99% de 4,288%, lo cual significa que con un 99% de confianza el NASDAQ
aumentara a lo más un 4,288% o sea de 2379,1 a 2481,12 puntos, bajo situaciones
normales de mercado.
Hay que tener en cuenta que cuando el inversionista se encuentra en una posición
corta un aumento en el valor del activo representa una pérdida para el.
Resumen de resultados
Valor del índice hoy corresponde al ultimo dato de la serie (25 octubre 2006), y el valor del índice
mañana (26 octubre 2006) corresponde a la peor situación esperada bajo situaciones normales de
mercado con un 99% de confianza
19
Distribución General de Pareto en los retornos:
El método que aquí aplicaremos consta de los siguientes pasos: elegir un umbral
“u”, ajustar la Distribución General de Pareto a los puntos que exceden el umbral y luego
calcular puntos e intervalos para el VaR y para el ES (expected shortfall).
Dado que no existe un criterio aprobado mejor que el resto, cada analista elige
aquel que desee. En el caso de éste análisis, el método utilizado es aquel que deja el 5% de
las observaciones como observaciones extremas.
Los siguientes gráficos muestran los excesos de la media11 para los datos del índice
NASDAQ. En este caso, “u” debe ser 2.93% para el caso de la cola izquierda, dejando 125
observaciones por sobre este umbral. Para el caso de la cola derecha, “u” debe ser 2.783%,
dejando 125 observaciones por sobre el umbral (lo que corresponde al 5% de la muestra.
10
An application of Extreme Value Theory for Measuring Financial Risk.
11
Funcion de exceso sobre la media, se encuentra planteada en el Anexo 1.
20
10
0
-14.17 2.93 8.56
u
El grafico anterior corresponde a la cola izquierda
0
2.78
u
El grafico anterior corresponde a la cola derecha
21
Estimación de Máxima Verosimilitud:
Con los resultados anteriores nosotros suponemos que las observaciones sobre el
umbral, es decir, las observaciones de la cola, se distribuyen conforme a la distribución
General de Pareto (GPD). Existen diferentes métodos para estimar los parámetros de la
distribución GPD, ahora explicaremos el Método de Máxima Verosimilitud (ML).
Par una muestra y = {y1,..., yn} la función log-likelihood L (ε, σ/y) para la GPD es
el logaritmo de la función de densidad de las n observaciones.
n
1
n log ( 1) log( 1 yi ) Si, 0
i 1
L( , / y )
n
1
n log yi Si, 0
i 1
n (VaR p u) VaR p
VaR p u (( p) 1 u
Nu ES VaR p ^ ^ ^
1 1 1
12
Demostración en Anexo 1
22
La siguiente figura muestra como la GPD se ajusta los excedentes por sobre el
umbral.
0.995
0.99
0.985
0.98
0 5 10 15 20 25
Con ambos parámetros estimados, podemos calcular tanto el VaR como el ES con
una probabilidad dada. Para una probabilidad de 0.01, el VaR estimado es de 4.731% y el
ES estimado es de 6.195% para el caso de la cola Izquierda. Para la cola Derecha, el VaR
es de 5.082%, mientras que el ES es de 6.985%.
Luego de haber calculado los valores estimados para VaR y ES, calculamos los
intervalos de confianza de éstos.
23
Estimación de los intervalos:
n
( p) 1
Nu 1 0
1 (1 y) ,
VaR p u
G ,VaRp ( y)
n y
1 p * exp( ) 0
Nu VaR p u ,
n n
( p) 1 ( p) 1 1
Nu Nu 1 0
(1 y) ,
(VaR p u) VaR p u
g ,VaRp ( y)
n y
p * exp
Nu VaR p u 0
,
(VaR p u)
24
Asimismo, usando la siguiente reparametrización para ε ≠ 0, tenemos:
n
(
p) 1
G ,ES p ( y) Nu 1 0
1 (1 y
) ,
( ES p u )(1 )
n n
( p) 1 p) ( 1 1
g Nu Nu 0
, ES p ( y )
1
(1 y
) ,
( ES p u )(1 ) ( ES p u )(1 )
Las siguientes figuras muestran las regiones de confiabilidad para la cola izquierda
obtenidas a través del proceso anteriormente descrito. Estas estimaciones están hechas
sobre el VaR y ES a través de la función reparametrizada.
-1.921
La figura anterior nos muestra el perfil de la función log-likelihood y los intervalos de confianza para el VaR al
0.01
25
ML
BCa
0.359
0.142
-0.009
Esta figura nos muestra los intervalos de confianza individuales y conjuntos al 95% para las estimaciones de ε
para el VaR al 0.01. Los puntos representan 1000 estimaciones hechas por el método Booststrap (BCa) para el
NASDAQ.
-1.921
Esta figura nos muestra el perfil de la función log-likelihood y los intervalos de confianza para el ES al 0.01
26
ML
BCa
0.361
0.142
-0.017
Esta figura nos muestra los intervalos de confianza individuales y conjuntos al 95% para las estimaciones de ε
para el ES al 0.01. Los puntos representan 1000 estimaciones hechas por el método Booststrap para el
NASDAQ.
Como se observa en los gráficos anteriores, de los 1000 puntos estimados por
booststrapping, aproximadamente el 95% se encuentra dentro del intervalo conjunto, no
así de los intervalos individuales.
13
Bootstraping.
14
Máxima Verosimilitud.
27
Los resultados de la tabla anterior muestran que, con una probabilidad de 0.01, las
pérdidas de mañana, con una posición larga, serán mayores de 4.731%, y que la
correspondiente pérdida esperada, que es el promedio de las pérdidas en las situaciones en
cuando se excede el 4.731%, es 6.195%.
En este caso, los resultados de la tabla anterior muestran que, con una probabilidad
de 0.01, las pérdidas de mañana, con una posición corta en el índice, serán mayores de
5.082%, y que la correspondiente pérdida esperada, que es el promedio de las pérdidas en
las situaciones en cuando se excede el 5.082%, es 6.985%.
Es importante destacar que el limite superior del intervalo de confianza para ε es tal
que el momento de primer orden es finito (1/0.34 = 2.94 >1). Lo anterior garantiza que la
estimación del ES, que es un primer momento condicional, existe para ambas colas.
28
Conclusión
El cálculo del VaR asumiendo normalidad en los retornos es válido cuando la serie
sigue una distribución normal y no presenta colas gruesas, de no ser así estaríamos
subestimando el VaR, lo cual se ve reflejado en los resultados de nuestro trabajo.
Por lo cual la Teoría del Valor Extremo, es aplicable para cualquier escenario, sin
embargo esta puede ser mas compleja en su cálculo, debido a esto es preferible testear la
distribución para notar la existencia de colas gruesas, de ser así hay que utilizar EVT, de lo
contrario, ajustamos la serie a una distribución normal.
29
Expected ShortFall como medida de exposición ante eventos extremos
Una de la utilidades del cálculo del ES, es analizar que tanta volatilidad tiene un
activo ante cambios extremos en el mercado, como por ejemplo el NASDAQ tiene un ES
de 6.195, lo cual nos señala que ante cambios extremos las perdidas promedio esperada
será de un 6.195%.
Por otra parte el ES se puede utilizar para hacer un análisis comparativo entre dos
activos o portafolios, por ejemplo podemos encontrar que un portafolio A tiene un VaR
mayor pero un ES menor que el portafolio B, lo cual nos señala que la volatilidad del
portafolio A es mayor en situaciones normales de mercado, pero en situaciones extremas
el portafolio B tiene una mayor volatilidad.
VaR ES
Coca-Cola (KO) 4,259 6,552
Mc Donalds (MCD) 4,609 6,454
Estos resultados, nos dicen que Coca-Cola tiene una menor volatilidad que
Mc.Donalds bajo situaciones normales de mercado, por lo cual un inversionista debería
estar mas seguro con una acción de Coca-Cola ante fluctuaciones nórmales de mercado,
pero por el contrario un inversionista debería estar mas seguro con una acción de
Mc.Donalds, bajo movimientos extremos de Mercado.
30
Esta información adicional es de mucha relevancia a la hora de tomar la decisión
de donde invertir, pero mas que eso es conocer la máxima perdida esperada bajo distintas
situaciones del mercado.
Otra utilidad del ES es su aplicación para coberturas, sin embargo esta no se utiliza
en la realidad para cubrirse ya que no es lógico estar cubierto todos los días del año para
eventos que suceden una vez por año o menos.
31
Anexo 1
Demostración
F (u y ) F (u ) F ( x) F (u )
Fu ( y )
1 F (u ) 1 F (u )
Despejando F(X),
F ( x) (1 F (u )) Fu ( y) F (u )
^
1
^ Nu ^ Nu
F ( x) (1 (1 ^
( x u )) (1 )
n n
Simplificando
^
1
^ Nu ^
F ( x) 1 (1 ^
( x u )) ,
n
^
^ n ^
VaR p u ^
(( p) 1)
Nu
32
Por definición ES,
^ ^ ^ ^
ES p VaR p E ( x VaR p / x VaR p )
z
e( z ) E(x z/x z)
1
^ ^ ^ ^ ^ ^
^ ^ (VaR p u ) VaR p u
ES P VaR p ^ ^ ^
1 1 1
33
Bibliografía
34
- Alexander J. McNeil; May 1999; Extreme Value Theory for Risk
Managers.
- Philip Hua and Paul Wilmontt; Vale-at-Risk and Market Crashes.
- Alexander J. McNeil and Rüdiger Frey; June 1999; Estimating the Tail of
Loss Severity Distribution using Extreme Value Theory.
- Danielsson, Jon, and Casper de Vries; 1997; “Value at Risk and Extreme
Returns”; Discussion Paper No. 273, pp. 1-33 (London : School of
Economics and Political Science).
- Koedijk,Kees , Ronald Huisman, and Rachel Pownall; 1998; “VaR-x : fat
tails in financial risk management”; The Journal of Risk 1, pp. 47-63.
- Hull, Jhon, and Alan White; 1998; “Value at Risk When Daily Changes in
Market Variables are Not Normally Distributed ”; Journal of Derivatives,
Vol.5, No.3 (Spring), pp. 9-19.
35
36
37