DE LA FUNDACIN C / 219
Teora de cpulas.
Introduccin
y aplicaciones
a Solvencia II
Teora de cpulas
Introduccin y aplicaciones
a Solvencia II
10%
www.fundacionmapfre.org
ISBN: 978-84-9844-647-0
Depsito Legal: M-18220-2017
Maquetacin y produccin editorial: Cyan, Proyectos Editoriales, S.A.
A M Carmen, Mnica y Noemi
Luis Latorre Llorns es doctor en Ciencias Econmicas, Actuario de seguros y licen-
ciado en Derecho. En 1970 ingres en el cuerpo de Inspeccin de Seguros del
Ministerio de Hacienda y en 1980 en el de Inspectores de Finanzas del Estado. Ha
desempeado los cargos de Subdirector General de Inspeccin y Director del
Gabinete de Estudios y Relaciones Internacionales de la Direccin General de
Seguros y Fondos de Pensiones. Ha sido tambin ponente decano en el Tribunal
Econmico Administrativo Central, en donde prest servicios hasta su jubilacin, en
el ao 2007.
Constituye para m una gran satisfaccin tener la oportunidad de prologar esta nue-
va obra de Luis Latorre. Lo es por la personalidad del autor y lo es tambin por la
posibilidad de conocer el resultado de un anlisis en profundidad en una materia del
mayor inters en el anlisis de riesgos.
Luis Latorre, inspector de Seguros del Estado e inspector de Hacienda del Estado,
ocup cargos de relevancia en la Direccin General de Seguros en una poca difcil
de cambios para el sector asegurador. Aunque profesionalmente orient su activi-
dad hacia otros mbitos de la funcin pblica, hemos tenido siempre la oportunidad
de contar con su anlisis riguroso y reflexivo -caractersticas ambas tan definidoras
de su personalidad- en una gran diversidad de asuntos en el mbito actuarial y de
los seguros en general.
Son varias las personas que me han ayudado en la redaccin de este trabajo, ya sea
dndome tiles orientaciones, ya manifestando su buena disposicin a colaborar en
mi propsito de una forma u otra. ste ha sido el caso de los profesores Jos Antonio
Gil Fana, Antonio Heras Martnez (ambos de la Facultad de Ciencias Econmicas y
Empresariales de la Universidad Complutense de Madrid) y Mara Victoria Rivas
Lpez (profesora del Centro Universitario Villanueva, adscrito a dicha universidad) as
como tambin de Lorenzo Esteban Jdar, Carlos Montalvo Rebuelta y Daniel Prez
Ortiz, los tres inspectores de la Direccin General de Seguros y Fondos de Pensiones
que en aquellas fechas prestaban servicios en EIOPA (European Insurance and
Occupational Pensions Authority). Agradezco profundamente a todos ellos su actitud
generosa.
He recibido otras ayudas que, de forma concreta, sirvieron para allanar obstculos,
corregir errores en el texto o, en general, mejorarlo. Aqu debo citar a Guillermo
Ayala Gallego (profesor del Departamento de Estadstica e Investigacin Operativa de
la Universidad de Valencia) que tuvo la amabilidad de leer el borrador y cuyas obser-
vaciones sobre el mismo fueron muy tiles para m. Enrique Llorns Fuster (profesor
del Departamento de Anlisis Matemtico de la misma universidad) no slo me asis-
ti en la resolucin de ciertas cuestiones matemticas sino que, adems, me puso en
contacto con el profesor Ayala. Por su parte, Jos Luis Llorns Fuster (profesor de la
Universidad Politcnica de Valencia) me hizo valiosas sugerencias sobre la presenta-
cin del texto, lo que me permiti suplir, al menos en gran parte, mi torpeza al res-
pecto. Jos Luis Vilar Zann (profesor de la Facultad de Ciencias Econmicas y
Empresariales de la Universidad Complutense de Madrid) tuvo asimismo la paciencia
de leer el trabajo y sus valiosas observaciones contribuirn a mejorar su claridad.
INTRODUCCIN 23
1. Dependencia y correlacin 43
2. Distribuciones esfricas 51
3. Distribuciones elpticas 53
II. CPULAS 65
3. Cpulas y dependencia 82
4. Cpulas arquimedianas 93
La cpula de Gumbel 96
La cpula de Clayton 97
La cpula de Frank 98
Cpula t 156
Cpula t 227
Figura 1.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
Figura 1.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
Figura 1.3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
Figura 1.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
Figura 1.5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
Figura 1.6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
Figura 1.7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
Figura 2.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
Figura 2.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
Figura 2.3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
Figura 2.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
Figura 2.5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
Figura 2.6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
Figura 2.7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
Figura 2.8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
Figura 2.9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
Figura 2.10 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
Figura 2.11 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
Figura 2.12 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
Figura 2.13 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
Figura 2.14 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
Figura 2.15 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
Figura 2.16 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
Figura 2.17 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
Figura 2.18 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
Figura 2.19 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
Figura 2.20 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
Figura 2.21 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
Figura 2.22 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 117
Figura 2.23 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
Figura 2.24 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
Figura 2.25 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 121
Figura 2.26 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124
Figura 2.27 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 126
Figura 2.28 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 127
Figura 2.29 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 128
Figura 2.30 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 131
Figura 2.31 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
Figura 2.32 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
Figura 3.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 145
Figura 3.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 150
Figura 3.3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 156
Figura 3.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 158
Figura 4.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 165
Figura 4.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 165
Figura 4.3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 166
Figura 4.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 168
Figura 4.5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 171
Figura 4.6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 172
Figura 4.7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 179
Figura 4.8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 192
Figura 4.9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 194
Figura 4.10 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 195
Figura 4.11 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 198
Figura 5.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 213
Figura 5.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 213
Figura 5.3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 216
Figura 5.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 218
Figura 5.5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 219
Figura 5.6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 220
Figura 5.7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 220
Figura 5.8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 222
Figura 5.9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 223
Figura 5.10 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 225
Figura 5.11 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 227
Figura 5.12 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 228
Figura 5.13 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 230
Figura 5.14 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 233
Figura 5.15 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 236
Figura 5.16 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 236
Figura 5.17 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 242
Figura 5.18 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 243
Figura 5.19 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 244
NDICE DE CUADROS
Cuadro 1.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
Cuadro 2.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
Cuadro 3.1a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 145
Cuadro 3.1b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 146
Cuadro 3.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 149
Cuadro 3.3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 162
Cuadro 4.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 165
Cuadro 4.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 166
Cuadro 4.3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 167
Cuadro 4.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 167
Cuadro 4.5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 168
Cuadro 4.6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 181
Cuadro 4.7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 182
Cuadro 4.8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 183
Cuadro 4.9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 183
Cuadro 4.10 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 185
Cuadro 4.11 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 186
Cuadro 4.12 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 187
Cuadro 4.13 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 187
Cuadro 4.14 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 191
Cuadro 4.15 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 193
Cuadro 4.16 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 194
Cuadro 4.17 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 195
Cuadro 4.18 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 197
Cuadro 4.19 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 197
Cuadro 5.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 214
Cuadro 5.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 217
Cuadro 5.3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 219
Cuadro 5.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 221
Cuadro 5.5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 226
Cuadro 5.6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 229
Cuadro 5.7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 231
Cuadro 5.8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 232
Cuadro 5.9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 233
Cuadro 5.10 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 233
Cuadro 5.11 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 234
Cuadro 5.12 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 238
Cuadro 5.13 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 239
Cuadro 5.14 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 243
Cuadro 5.15 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 244
Cuadro 5.16 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 249
Cuadro 5.17 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 250
NDICE DE EXPRESIONES
(1.1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
(1.2) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
(1.3) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
(1.4) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
(1.5) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
(1.6) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
(2.1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
(2.2a) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
(2.2b) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
(2.3) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
(2.4) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
(2.5) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
(2.6) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
(2.7) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
(2.8) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
(2.9) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
(2.10) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
(2.11) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
(2.12a) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
(2.12b) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
(2.13) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
(2.14) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
(2.15) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
(2.16) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
(2.17) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
(2.18) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
(2.19) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
(2.20) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
(2.21) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
(2.22) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
(2.23) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
(2.24) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
(2.25) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
(2.26) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
(2.27) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
(2.28) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
(2.28a) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
(2.28b) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
(2.29) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
(2.30) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
(2.31) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
(2.32) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
(2.33) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
(2.34) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 137
(3.1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
(3.2) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
(3.3) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 148
(4.1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 170
(4.2) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 170
(4.3) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 173
(4.4a) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 174
(4.4a) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 174
(4.5a) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 177
(4.5b) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 177
(4.6) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 192
(5.1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 203
(5.2) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 204
(5.3) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 205
(5.4) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 208
(5.4a) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 216
(5.4b) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 231
(5.5) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 245
(5.6) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 247
(5.7) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 248
INTRODUCCIN
David X. Li naci en una zona rural de China1 en donde vivi hasta iniciar sus estu-
dios de Economa en la prestigiosa Universidad de Nankai. En 1987 y tras apren-
der francs en cuatro meses viaj a Canad y obtuvo un MBA en la Universidad de
Quebec, un mster en matemtica actuarial y un doctorado en estadstica en la de
Waterloo (Ontario). En 1997 empez a trabajar en el Canadian Imperial Bank of
Commerce y tres aos ms tarde se traslad a Nueva York. All prest sus servi-
cios en la unidad de Riskmetrics de J.P. Morgan y public el famoso artculo ya
mencionado. En este trabajo aplicaba la cpula gaussiana a la valoracin de los
CDOs (Collaterized Debt Obligations); en primer lugar, basaba el clculo de las
probabilidades de insolvencia de los diversos activos en el valor de mercado de su
correspondiente spread o prima de riesgo, teniendo en cuenta tambin la pro-
porcin de recobro sobre la inversin total, en caso de insolvencia. Obtenida la
probabilidad de insolvencia durante un ao, deduca el tanto instantneo corres-
pondiente y con l la funcin de densidad de la variable expresiva del tiempo
23
hasta la insolvencia del activo en cuestin (supervivencia). Aceptaba la distribu-
cin exponencial como modelo para dicha variable. Pero, como dice Embrechts
(2009), al contrario de lo que ocurre con la normal o la t multivariante, no existe una
funcin de distribucin estndar exponencial d-dimensional. Para obtener una fun-
cin de distribucin multivariante cuyos mrgenes fueran exponenciales, utiliz la
cpula Gauss2. Este modelo meta-gaussiano contina diciendo P. Embrechts
se hizo enormemente popular y caus problemas porque el mercado crey dema-
siado en l. Se lleg a decir que la cpula gaussiana es el peor invento para la
gestin del riesgo de crdito que ha habido nunca.
En agosto de 2004 Moodys incorpor la cpula gaussiana a sus mtodos para va-
lorar CDOs3. Dice Sam Jones que los bancos ya no necesitaban mantener los ries-
gos, sino que, al ser capaces de evaluarlos mediante complejas matemticas, po-
dan empaquetarlos y comercializarlos como cualquier otra inversin. En lugar
de ir cobrando los intereses y las amortizaciones de un prstamo hipotecario, los
bancos los agrupaban y troceaban, clasificando cada tramo segn su calificacin
crediticia y transfiriendo los productos a entidades especiales que, a su vez, emi-
tan bonos para obtener efectivo. La facilidad de liquidar estos activos pudo incen-
tivar la realizacin de ms operaciones cuyos crditos subyacentes tuvieran una
insuficiente calidad. As, en el ao 2000 el total de CDOs emitidos era de unos
70.000 millones de dlares, mientras que hacia 2006-2007 la cifra alcanz los
500.000 millones por ao. En 2010 descenda a slo 10.000 millones.4
http://cyrusfarivar.com/docs/li.defaultcorrelation.pdf
2
E. Forslund y D. Johansson: Gaussian Copulas. What happens when models fail?. Noviembre 2012.
4
24
por Li fue utilizado en muchos casos por quienes no lo entendan del todo y que
probablemente no slo se us, sino que se abus de l5 al considerarlo como una
frmula para medir los riesgos y su dependencia, sin preocuparse de si el modelo
y los datos introducidos eran o no adecuados a la realidad circundante.
Probablemente la mayor crtica que cabe hacer se refiere a la utilizacin del mo-
delo gaussiano y, con l, de la correlacin lineal para medir la dependencia
entre riesgos, sin que ello fuera correcto en este caso por las caractersticas y
comportamiento de tales riesgos. En agosto de 1999 P. Embrechts, A. McNeil y D.
Straumann haban ya publicado un importante artculo significativamente titulado
Correlation and Dependence in Risk Management: Properties and Pitfalls; gran
parte de este trabajo se dedicaba a estudiar en qu condiciones la correlacin es
una medida de la dependencia y a advertir del peligro de extender el uso de aquel
concepto ms all del mbito de las distribuciones elpticas. Ms recientemente,
en enero de 2010, C. Donnelly y P. Embrechts han publicado el artculo The Devil
is in the tails: actuarial mathematics and the subprime mortgage crisis. Se trata de
un trabajo en que analizan la cuestin, describen el modelo de Li y concluyen que
la cpula gaussiana, al carecer de la dependencia en colas6 no es capaz de re-
flejar el fenmeno de la mayor dependencia entre insolvencias en situaciones de
crisis, como lo son otros modelos de cpula.
5
Confessions of a Risk Manager, The Economist, 7 de agosto de 2008.
6
Nos ocuparemos de este concepto en el Captulo II, 3.
25
construirse una funcin de distribucin multivariante. Es un concepto muy til en
muchas situaciones en que interese modelizar la dependencia entre variables
aleatorias y construir una gran variedad de funciones de distribucin multivarian-
tes. No tiene nada de extrao que encuentre muchas aplicaciones en el mbito de
las finanzas, la gerencia de riesgos y el seguro.
26
unos y otros. Las cpulas ofrecen la posibilidad de tener en cuenta las posibles
dependencias entre tales riesgos.
En los seguros distintos al de vida, es clara la posible aplicacin de las cpulas; pen-
semos, por ejemplo, en contratos que cubran riesgos dependientes entre s, o en la
posibilidad de acontecimientos que afecten simultneamente a varias plizas.
27
dependencia entre estos riesgos. Si fueran perfectamente dependientes, de ma-
nera que variasen todos al unsono, entonces se ve intuitivamente que el capital
total tendra que ser la suma 200 + 50 + 100 unidades monetarias. Pero, fuera de
este supuesto tan poco realista, la respuesta a aquella pregunta ya no es tan in-
mediata y dependera de la forma en que esos tres riesgos se comportan conjun-
tamente. Y aqu es donde las cpulas podran constituir un interesante instrumen-
to de anlisis.
A principios de los aos 2000 Suiza implant el sistema conocido como SST (Swiss
Solvency Test), tambin encaminado a exigir de cada asegurador el nivel de capi-
tales de solvencia que exija la medida de los riesgos a que est sometido. Dicho
sistema de solvencia, as como el de Basilea establecido para la Banca, influy
fuertemente en la normativa de la Comunidad Europea conocida como Solvencia
II, para diferenciarlo del sistema que lo ha venido precediendo durante muchos
aos. Este primitivo sistema no pretenda fundarse en la totalidad de los riesgos
que afectan al asegurador, sino slo en los cubiertos por las plizas; no se plan-
teaba problema alguno de dependencia entre varios riesgos (separaba los riesgos
de los seguros de vida y los distintos al de vida). El nuevo sistema de Solvencia II,
por el contrario, trata de cubrir varios tipos de riesgos de los que afectan a las
empresas de seguros (operacional, suscripcin de seguros distintos del seguro de
vida, suscripcin del seguro de vida, suscripcin del seguro de enfermedad,
28
riesgo de mercado y riesgo de incumplimiento de la contraparte o riesgo de cr-
dito). A su vez, y a travs de la estructura modular a la que nos referiremos en el
Captulo V, cada uno de stos se divide en otros, todo lo cual plantea el problema
de cmo agregarlos hasta llegar al capital total o Capital de Solvencia Obligatorio,
conocido tambin por las siglas CSO o SCR, estas ltimas del ingls Solvency
Capital Requirement. Cmo se hace esta agregacin? En el modelo estndar
que establece Solvencia II, se logra mediante los coeficientes de correlacin de
cada par de riesgos del mismo nivel y mdulo (o sub-mdulo); en la mayora de los
casos se supone que las variables aleatorias representativas de cada riesgo se
distribuyen normalmente con media nula y la varianza que el propio sistema apor-
ta; otras veces se recurre igualmente a dichas correlaciones, las cuales permiten
determinar la desviacin tpica de la variable suma de varios riesgos, y se adopta
el modelo logartmico-normal, lo que conduce al capital correspondiente a dicha
suma. Es decir, que Solvencia II no es excepcin a la regla de que los sistemas de
solvencia basados en varios riesgos exigen analizar la dependencia entre ellos. El
sistema resuelve la cuestin con los coeficientes de correlacin lineal como ins-
trumento para medir la dependencia. Como veremos en el Captulo I, tal proceder
es correcto slo en el mbito de las distribuciones elpticas, lo que apunta a la
teora de cpulas como eficaz alternativa para modelizarlos. En efecto, el coefi-
ciente de correlacin mide el grado de dependencia lineal de una distribucin a lo
largo de todo su rango; pero esta dependencia puede ser mucho mayor en una
zona de la distribucin que en el resto. En la prctica, encontramos en ocasiones
variables o riesgos cuya dependencia es dbil o nula en condiciones normales,
mientras que, bajo situaciones extremas (por ejemplo, una crisis econmica o sa-
nitaria, una catstrofe natural, etctera) la correlacin es considerable. Un ejem-
plo es la grfica siguiente, en la que cabe distinguir, a simple vista, un ncleo in-
tegrado por la mayora de los valores, que no aparentan correlacin alguna; pero,
en la zona positiva extrema de ambas variables (para valores de cada una superio-
res a 4, aproximadamente), el comportamiento conjunto es claramente distinto, al
mostrar una correlacin positiva considerable. Si calculamos el coeficiente de co-
rrelacin separadamente para ambas zonas, tendramos que, en la primera, ste
es prcticamente nulo, mientras que en la segunda es, aproximadamente, 0,9. La
correlacin del conjunto de todos estos puntos sera 0,08. Observemos que la co-
rrelacin alta se da precisamente en la cola de la distribucin, es decir, en la
zona en que tanto X como Y toman valores elevados, que en muchas aplicaciones
29
es la ms relevante. Este ejemplo manifiesta dos cosas: primera, el peligro de
utilizar el coeficiente de correlacin en ciertos casos y la segunda, que una posible
solucin para evitarlo es ajustar dicho coeficiente, es decir, elevarlo para que
refleje la dependencia lineal en la zona extrema, ms que en el conjunto. Conviene
que nos detengamos brevemente en esta cuestin.
7
A Global Framework for Insurer Solvency Assessment. Research Report of the Insurer Solvency
Assessment Working Party.
Diversification. Technical Paper, 31 de octubre de 2005.
8
30
Actuarial Internacional proponga su uso como el mtodo tericamente correcto
para calcular los efectos de la diversificacin; (dichos efectos aadimos aqu
estn en funcin de la dependencia existente entre los riesgos en juego: a mayor
dependencia, menor ser la atenuacin que sobre el riesgo total genere la agru-
pacin de todos y recprocamente; de ah que medir los efectos de la diversifica-
cin exija conocer la dependencia entre las variables o riesgos intervinientes).
Ahora bien, en dicho informe se reconocen las dificultades de clculo que las c-
pulas pueden plantear y se ofrece como alternativa a las mismas la utilizacin de
una matriz de correlaciones ajustada, es decir, modificada para que refleje las
correlaciones entre las variables en la cola de la distribucin. El problema que
plantea esta alternativa es cmo cuantificar estos coeficientes de correlacin li-
neal en la cola. Los datos de que habitualmente se dispone en esa zona de las
distribuciones suelen ser escasos, insuficientes para estimar los coeficientes y la
nica solucin viable consiste en determinarlos acudiendo al buen juicio y expe-
riencia. Es frecuente graduar las correlaciones en la cola en valores que sean
mltiplos de un cierto valor, por ejemplo 0,25, para llegar as a la gama de valores
0, 0,25, 0,50, 0,75 y 1, segn se califique la correlacin de nula, escasa, apreciable-
moderada, alta y total. sta es la solucin adoptada en la normativa de Solvencia
II referente al llamado modelo estndar; no se recurre en l a las cpulas, sino
que se opta quiz por el principio de proporcionalidad9 por la solucin de apli-
car coeficientes de correlacin en la cola de las distribuciones.
9
Segn el prembulo de la Directiva 2009/138/CE del Parlamento Europeo y del Consejo, de 25 de
Noviembre de 2009 (Directiva de Solvencia II), el principio de proporcionalidad consiste en no impo-
ner requisitos que puedan resultar demasiado gravosos para las pequeas o medianas empresas
de seguros.
31
Quiz uno de los trabajos ms frecuentemente citados en materia de teora de
cpulas sea el artculo de Frees y Valdez Understanding Relationships Using
Copulas. Los autores muestran cmo el planteamiento de diversas situaciones
propias del seguro de vida conduce a diversos tipos de cpulas. Desarrollan un
ejemplo numrico para estudiar la dependencia que en el seguro de responsabi-
lidad civil se da entre el importe de la indemnizacin y el de los gastos generados
por la liquidacin de las prestaciones (Allocated loss adjustment expense o ALAE).
Se trata de un ejercicio basado en datos reales, obtenidos del Insurance Services
Office, Inc. Una vez elegida la cpula y distribuciones marginales que mejor se
ajustan a los datos utilizados, se aportan dos ejemplos de cmo puede aprove-
charse el conocimiento de la funcin de distribucin bivariante que resulta de
aqullas. Uno, es la tarificacin de un reaseguro de exceso de prdida con deter-
minadas caractersticas. Y el otro consiste en elaborar varias funciones de regre-
sin cuantlica, que permiten obtener informacin acerca del coste del ALAE para
cada valor de la indemnizacin.
Otro caso es Staudt, (2010), en cuyo artculo Tail Risk, Systemic Risk and Copulas
desarrolla un ejemplo de cmo determinar el tail conditional expectation o expec-
ted shortfall10 al 95 % de una cartera compuesta por 4 lneas o tipos de seguros:
responsabilidad civil general, responsabilidad civil de vehculos comerciales, multi-
rriesgos sobre los bienes comerciales y responsabilidad civil comercial. Aplica dos
tipos de cpulas (gaussiana y t) a una serie de datos consistentes en ratios de si-
niestralidad observados entre 1986 y 2008, recopilados por el Departamento de
Seguros del Estado de Texas.
10
Segn se examinar en el Captulo V, ambas denominaciones inglesas significan el valor medio
condicionado a que la variable supere el cuantil prefijado o bien dficit medio. Se refiere al valor
medio de la variable X, condicionado a que X>VaR. Por ejemplo, el VaR al 99 % es el valor que tiene
la propiedad de que X slo lo supera con una probabilidad del 1 %, pero nada informa acerca de la
cuanta aquellos valores de X superiores a dicho VaR. Este inconveniente trata de salvarlo dicho
dficit medio.
32
Exposures Using Copulas) se dedica a estudiar una aplicacin de la cpula t a la
cobertura del riesgo de catstrofes naturales, que pueden originar dependencia
entre los resultados de tres lneas de seguro: daos a las viviendas, daos a ins-
talaciones comerciales y automviles.
33
correlacin por rangos es la misma que la de dicha distribucin bivariante. En
cuarto lugar se aborda el estudio de las cpulas, reconocindose as la utilidad de
estos modelos en el anlisis financiero dinmico.
34
Solvencia II con una base de datos similar a la del trabajo anterior, aunque exten-
dida a 2010, bajo diversas hiptesis en cuanto a las matrices de correlacin, mos-
trndose la alta sensibilidad del resultado ante cambios en stas. Definido un
modelo interno basado en la prediccin del resultado de cada lnea de seguro, se
aplican las cpulas de Gauss y de Student con diversos parmetros. Este trabajo
contiene una frmula de credibilidad para los coeficientes de correlacin basada
en la Estadstica bayesiana.
J. De Mey, (2007).
14
35
Warranty16. Los autores describen estas operaciones, examinan sus clases y sus
efectos. En los bonos de catstrofes (o CAT bonds), el asegurador o reasegurador
que desea proteccin frente al riesgo de catstrofes promueve la emisin de tales
bonos a travs de una entidad especfica (special purpose vehicle); estos ttulos tie-
nen dos ventajas principales para el suscriptor (que, habitualmente, es un fondo de
riesgo o hedge fund, un fondo de pensiones, una entidad aseguradora, etctera): por
una parte, le ofrecen, a cambio del riesgo que incorporan, una rentabilidad mayor
que los bonos convencionales y, por otra, las fluctuaciones en el valor de mercado
de estos bonos de catstrofes guardan muy baja correlacin con la de los restantes
de renta fija, lo que es un incentivo para el inversor que desee diversificar su carte-
ra. En caso de que ocurra la catstrofe a la que el bono est asociado, el inversor
soporta el riesgo en los trminos preestablecidos, bien dejando de percibir el im-
porte de los cupones, bien, adems, con la prdida del derecho a recuperar todo o
parte del principal. En los Industry Loss Warranties, el comprador de esta forma de
transferencia alternativa de riesgos (empresa de seguros o reaseguros que desea
cobertura para el caso de siniestralidad extraordinaria de los tipos de seguro a que
la operacin se refiere) paga una prima al vendedor (que puede ser un hedge
fund o fondo de riesgo, o bien un reasegurador). A cambio, obtiene el derecho a per-
cibir una cantidad de dinero en caso de que ocurra una catstrofe del tipo preesta-
blecido (por ejemplo, huracn, terremoto o cualquier tipo de catstrofe natural, et-
ctera) en el territorio a que la operacin se refiere y dentro del perodo tambin
fijado, que haya ocasionado unos daos totales al sector asegurador superiores a la
cifra establecida (por ejemplo, 20.000 millones de dlares); estos daos se determi-
nan por medio de unos ndices de prdidas (por ejemplo, los elaborados por Property
Claims Services para diversas zonas de los Estados Unidos). Si la catstrofe rene
los requisitos as fijados, el comprador de la garanta tendr derecho a percibir la
suma de dinero establecida en el contrato (por ejemplo, 100 millones de dlares).
En el trabajo sobre los Industry Loss Warranties, los autores se proponen determinar
el importe de la prima de una de estas operaciones, referentes a cobertura de da-
os extraordinarios por catstrofe natural que afecte a dos lneas de seguro (daos
y responsabilidad civil). Ante la dificultad de obtener datos reales, se ha llevado a
cabo una simulacin de las cuantas de los daos. El clculo se hace por medio de
simulaciones, adoptando diversas hiptesis acerca de las distribuciones marginales
36
(log-normal, Pareto, logstica, paralogstica, exponencial) y de las estructuras de
dependencia. El trabajo contiene una exposicin de la teora de cpulas.
Las cpulas, por tanto, irrumpen en el mundo actuarial como un instrumento que
puede aportar soluciones en muchos casos. Pero ello no significa que carezcan de
inconvenientes. Ya hemos hecho referencia a las dificultades de clculo que pue-
den plantear, pero hay otros. Quiz el principal sea que suele resultar problem-
tico determinar el modelo de cpula ms adecuado para una situacin dada. Los
datos requeridos para el ajuste de cpulas han de ser muy numerosos y no suele
contarse con informacin suficiente. Es raro que un solo asegurador disponga de
experiencia estadstica en forma de varios cientos de observaciones del vector
que interese para el caso (por ejemplo, pares de valores (x, y) significativos de la
siniestralidad observada en dos tipos de seguro) y eso representa una grave difi-
cultad a la hora, por ejemplo, de elaborar un modelo interno de solvencia. Y an si
dispusiera de una larga experiencia, habra que evitar utilizar datos que, por su
antigedad, carecieran de representatividad sobre los fenmenos que pretendie-
ran modelizarse. De ah que sea importante organizar la obtencin de estadsticas
procedentes de un gran nmero de empresas, de forma que sean aptas como ba-
ses de datos para este tipo de trabajos.
Los cinco Captulos que siguen contienen una introduccin bsica a las cpulas, re-
dactada con la intencin de lograr claridad en la exposicin y carcter prctico, con
multitud de ejemplos. Va destinada a los actuarios, estudiantes de materias actuaria-
les y, en general a quienes deseen tener una idea de qu son y para qu sirven las
cpulas. Se pretende aportar los conceptos necesarios para que, quien est interesa-
do, pueda profundizar en el estudio con mayor facilidad. Fieles a dicho carcter prc-
tico, incluimos en este trabajo una aplicacin de la teora de cpulas a la determina-
cin del capital de solvencia en un supuesto de seguros de responsabilidad civil,
segn criterios de la normativa de Solvencia II y valindonos de datos reales.
37
lineal puede conducir a conclusiones errneas. Se determina en qu mbitos pue-
de utilizarse este coeficiente como medida de la dependencia estadstica, lo que
exige estudiar las distribuciones esfricas y elpticas. Ello se hace mediante un
ejemplo y varias grficas. El Captulo concluye con el examen de otras medidas de
la dependencia, como los coeficientes de Kendall y de Spearman.
El siguiente Epgrafe de este Captulo II describe cmo unir o vincular una cpula
con las funciones de distribucin marginal, para obtener as una funcin de distri-
bucin multivariante. Otra cuestin de inters que examinamos en este Captulo
es la cpula de supervivencia; se explica el concepto, se aportan ejemplos y se
presenta la llamada de cpula HRT, til en aplicaciones actuariales.
Expresin (2.23)
17
38
Al principio de este Captulo II se hace referencia a las cpulas elpticas como
ejemplos de las implcitas; en el Epgrafe 9 se aportan ms casos de cmo obtener
cpulas a partir de distribuciones bivariantes que responden a problemas prcti-
cos. Uno de estos casos es la cpula de Marshall-Olkin, que se examina con cier-
to detenimiento. Finaliza este Captulo con una breve referencia a las cpulas de
valores extremos.
39
El ltimo Captulo, tambin de ndole prctica, se destina a la aplicacin de las
cpulas al clculo del capital de solvencia en el caso concreto de dos lneas de
seguro (responsabilidad civil general y responsabilidad civil de vehculos de mo-
tor), siguiendo los criterios establecidos por la normativa de Solvencia II. Los da-
tos utilizados en este caso han sido obtenidos de unas estadsticas publicadas en
forma de grfico en un documento de CEIOPS;18 no se trata, por tanto, de datos
simulados, sino reales. Tras una exposicin resumida de las normas de Solvencia
II ms relevantes, se desarrolla esta aplicacin prctica paso a paso, siguiendo el
proceso que qued esbozado en Captulos anteriores, es decir, se obtiene la cpu-
la emprica y las distribuciones marginales, a las que se busca un modelo terico.
Se comparan varios modelos de cpula por los mtodos ya enumerados y se de-
termina que la que mejor se ajusta a los datos disponibles es la cpula de Frank.
A partir de ella y de dichos modelos para las distribuciones marginales, se llega a
la distribucin de probabilidades y se determina el VaR y el TailVaR al 99,5% (que
es la probabilidad exigida por la normativa europea). La conclusin obtenida es
una cifra de capital de solvencia inferior a la que resultara de dicha normativa, lo
que pone de manifiesto que en el caso contemplado no queda justificada la eleva-
cin de coeficientes de correlacin para atender a la existente en la cola de las
distribuciones.
40
el texto. Un ejemplo de esto lo encontramos en la Nota a pie de pgina nmero 24.
Con los caracteres insertados en las grficas ocurrir lo mismo.
Los Captulos quedan divididos en Epgrafes numerados. Para citar uno de ellos
se escribe el nmero del Captulo y el del Epgrafe, separados ambos por un pun-
to. As, III.2 es el Epgrafe nmero 2 del Captulo III.
Anlogamente, para identificar los Cuadros se acude tambin a dos nmeros se-
parados por un punto: el primero significa el Captulo en que se encuentra y el
segundo, el nmero de orden dentro del Captulo. Y lo mismo ocurre con las
Figuras. As, la Figura 4.2 es la segunda dentro del Captulo IV.
Se ha procurado que los ejemplos queden destacados del texto general y para lo-
grarlo, cada uno de ellos va encabezado por un ttulo en negrita con el nmero,
segn el mismo sistema que el empleado para las Figuras y la expresin del ob-
jetivo o finalidad del Ejemplo propuesto. Adems, el texto del ejemplo cuenta con
unos mrgenes laterales algo superiores a los generales, para que el lector pueda
determinar fcilmente dnde empieza y dnde termina. Se ha seguido este crite-
rio a lo largo de todo texto salvo en los Epgrafes dedicados a contener ejemplos,
por haberse considerado innecesario; an en estos casos dichos ejemplos quedan
destacados del texto. Tal es el caso del Captulo IV, Epgrafe 4.
***
Con ello se intenta que salte a la vista la extensin de esas demostraciones, cuyo
seguimiento puede no convenir al lector ms interesado en el resultado que en el
procedimiento para llegar a l.
A las frmulas o expresiones que lo requieran, se asocian dos nmeros entre pa-
rntesis, separados por un punto. Por ejemplo, la expresin (2.5) est en el Captulo
II y es la nmero 5.
41
I. LA DEPENDENCIA ENTRE RIESGOS
1. DEPENDENCIA Y CORRELACIN
En efecto, la funcin de densidad conjunta del modelo normal bivariante (0,1) es:
Sobre las distribuciones elpticas, vase el Epgrafe 3 de este Captulo. Las distribuciones elpti-
20
cas bivariantes se caracterizan porque los puntos en que la densidad es constante, forman una
elipse. A esta clase pertenecen la distribucin normal y la t de Student con 0.
43
Si hacemos =0 podemos expresar f(x,y) como producto de las densidades margi-
nales de X y de Y, lo que demuestra su independencia.
Pues bien, si =0 resulta una funcin que no puede expresarse como producto de
las densidades marginales, faltando as la condicin necesaria y suficiente para la
independencia entre X e Y en este caso.
En este Epgrafe examinaremos cuatro ejemplos que ilustran el error que se co-
mete al utilizar la correlacin para medir la independencia en situaciones inade-
cuadas. En los dos primeros el error proviene de que las distribuciones no son
elpticas, mientras que en los otros dos, las distribuciones son elpticas pero el
coeficiente se utiliza de forma incorrecta.
44
La funcin de distribucin de Y es:
Por tanto,
Si, por ejemplo, s=6 obtenemos que Prob(X+Y 6) = 0,983 y Prob(X+Y > 6) =
0,017.
Este caso est sugerido por el que proponen Embrechts, McNeil y Straumann
(2002). Dos distribuciones bivariantes con el mismo coeficiente de correlacin
(en ambos casos, = 0,7) e iguales distribuciones marginales (gamma con
parmetros 3 y 1), presentan distintos grados de dependencia; para mostrar-
lo, representan grficamente 1000 puntos simulados a partir de cada una de
45
ambas distribuciones de probabilidad; en la primera parte de la grfica slo
unos pocos puntos estn ms all de los lmites x=8, y=8 mientras que en la
segunda los puntos que superan ambos valores que podramos considerar
como umbrales de peligrosidad son muchos ms.
21
La distribucin simulada en la Figura 1.1 tiene como estructura de dependencia la cpula gaus-
siana con =0,666; la de la Figura 1.2, la cpula es Gumbel con parmetro 1,543. En ambas las dis-
tribuciones marginales son exponenciales con media 5.
46
Figura 1.2: Simulacin de 2000 puntos; cpula Gumbel y mrgenes exponenciales.
En cada una de las Figuras 1.3 y 1.4 se representan 2000 puntos simulados de
una distribucin bivariante con correlacin =0,6 y mrgenes con media cero
47
y varianza igual a 2 en ambos casos. La diferencia est en que, en la Figura 1.3
la distribucin bivariante es la normal, mientras que en la Figura 1.4 los pun-
tos representados proceden de una t de Student con n=4 grados de libertad.
Dentro de la clase de las distribuciones elpticas a que ambas pertenecen, con
un mismo coeficiente de correlacin se originan distribuciones con muy dis-
tinta peligrosidad, como ilustra el distinto nmero de puntos (x, y) tales que
x+y>7 en uno y otro caso.
Figura 1.3: Simulacin de 2000 puntos de una normal bivariante (0, ) con =0,6.
Figura 1.4: Simulacin de 2000 puntos de una t bivariante con =0,6 y n=4 grados de libertad.
48
Ejemplo 1.4: Otro caso de incorrecta aplicacin de la correlacin lineal
49
unidades es de 1,16 milsimas y por tanto muy superior24. Concluimos una vez
ms que el coeficiente no sirve para determinar si existe o no independencia
entre variables aleatorias cuando la distribucin conjunta de ellas no sea el
modelo gaussiano, si bien desempea, en las dems distribuciones elpticas,
como es la t de Student, el papel de medir la dependencia, aunque de otra
manera (sin que incorrelacin equivalga a independencia). Lo comprobare-
mos empricamente siguiendo con este Ejemplo 1.4: en la distribucin biva-
riante de Student que acabamos de utilizar, hagamos que =0,5. En tal caso,
la dependencia es mayor que antes y la probabilidad de que ambas siniestra-
lidades (la de A y la de B) superen simultneamente el umbral de 25 unidades
es de 3,5 milsimas, es decir, aproximadamente el triple que si =025.
Adems de las limitaciones del coeficiente de correlacin que hemos venido ilus-
trando mediante ejemplos, an pueden citarse otros inconvenientes de esta me-
dida de dependencia:
50
2) Requiere que las varianzas sean finitas; de lo contrario, el coeficiente de corre-
lacin no est definido. Esto impide utilizar el coeficiente en casos como el se-
alado por Moscadelli (2004) referente al riesgo operacional, en donde apare-
cen distribuciones con media y varianza infinita.
(X,Y)[LN(X), LN(Y)].
2. DISTRIBUCIONES ESFRICAS
El valor de la densidad ser el mismo para cualquier par (x1,x2) tal que
x12 + x22 = 2
Si aplicamos a (x1, x2) una transformacin que traslade dicho punto a otro dentro
de la misma circunferencia obtenemos igual densidad. Esta transformacin ha
de ser ortogonal y se logra multiplicando el vector (X1,X2) por una matriz
51
ortogonal A (es ortogonal si A = A1). Por tanto, si A es una matriz ortogonal n x
n y X un vector de orden n,
A.X =d X 26
y tambin:
X =d A.X. 27
26
El smbolo =d significa que las variables que figuran a ambos lados del mismo siguen la misma
ley de distribucin estadstica.
27
En el caso de la normal bivariante, cualquier vector (X1,X2) tendr una distribucin uniforme si sus
componentes determinan puntos de la misma circunferencia de radio, por ejemplo, igual a 1; si hace-
mos ahora X1 = cos(), X2 = sen() (siendo el ngulo que forma el radio que pasa por (X1,X2) con el eje
de abscisas, resulta que es una variable uniforme (0, 2). Por otra parte, cualquier punto aleatorio
(N1,N2) del espacio bidimensional puede expresarse como producto de una variable R0 por (X1,X2), es
decir, que N1=R.X1, N2=R.X2. Este hecho puede aprovecharse para simular nmeros aleatorios distribui-
dos normalmente a partir de nmeros aleatorios (o, ms exactamente, pseudoaleatorios) uniformes
(0,1). Para obtener N1 simulamos el radio de la circunferencia, R; como R2 es la suma de los cuadrados
de dos variables independientes y normales (0, 1), su funcin de distribucin es F(x) = Prob(R2 x) = 1-
e-x/2, es decir, la de una variable 22. Pero F(X) es uniforme (0,1) y por tanto, si simulamos un nmero
aleatorio U1, uniforme (0,1) podemos hacer que U1 =1- e-x/2 y, por tanto, x= 2.LN(1 U1) es un valor
posible de R2; la simulacin aleatoria del radio R ser x=(2.LN(1 U1))0,5. Por otra parte, al ser uni-
forme (0, 2), resulta que 2U2 ser una simulacin de , si U2 es asimismo uniforme e independiente de
U1. Por tanto, el producto R.X =(-2 LN(U1))0,5 . cos (2U2) ser una simulacin del nmero aleatorio N1 y
R.Y =(-2.LN(U1))0,5 . sen (2U2) lo ser de N2 (es irrelevante sustituir 1U1 por U1 en estas expresiones).
stas son las conocidas frmulas utilizadas para convertir un par de nmeros aleatorios uniformes (0,1)
e independientes, en otro par N1 y N2, ambos normales (0,1) y tambin independientes entre s.
52
Pero si X tiene distribucin esfrica, X =d A.X (siendo A ortogonal) y por tanto,
3. DISTRIBUCIONES ELPTICAS
Una distribucin elptica resulta de hacer una transformacin afn sobre una esf-
rica. Dicha transformacin lleva a cabo, sobre la variable esfrica, un cambio en
las unidades, un giro y una traslacin. La ecuacin de la transformacin afn de un
vector X en n es
Y= C.X + M
53
Si aplicamos esta transformacin a la circunferencia de radio 2, x12 + x22 =4 ob-
tenemos la elipse 0,25 y12+ 1,5625 y22 = 4, representada en la Figura 1.5.
En tercer lugar, para trasladar dicha elipse, por ejemplo, 3 unidades horizon-
talmente y 1 unidad en sentido vertical, en ambos casos en sentido positivo
(es decir, hacia la derecha y hacia arriba), sumamos el vector:
54
Figura 1.6: Giro de la elipse por transformacin ortogonal.
W = B . A. X + M = C . X + M (1.2)
siendo
55
Figura 1.7: Elipse desplazada por suma de un vector.
es una matriz simtrica, porque equivale al producto de una matriz (C) por su
traspuesta. Por tanto, queda XX= (WM). -1. (WM). Observemos que -1 es
tambin simtrica. En efecto, si lo es, entonces = y en consecuencia -1 =
()-1 = -1).
56
Hay que dividir entre el jacobiano de la transformacin, que es la raz cuadrada del
determinante de .
Pero E(ei TCX )es la funcin caracterstica de la variable esfrica X con vector de
parmetros C.T, que es (ver expresin (1.1)):
hx(TCCT) = hx(T T)
57
Mediante unas transformaciones algebraicas algo tediosas expresaremos
esta funcin de densidad en una forma ms reconocible como propia de la t
bivariante modificada, para que tenga los valores medios y varianzas elegidos
(ver expresin (4.6)28); en este caso la media de W1 es 3 y la de W2 es 1. Las
varianzas, en el caso general de la distribucin t, son iguales a /(-2) pero en
nuestro ejemplo ello no tiene por qu ser as; en efecto, recordemos que la
presente distribucin bivariante resulta de una esfrica (t con 4 grados de li-
bertad y correlacin nula) deformada y desplazada caprichosamente, lo que
determina que las dispersiones marginales tengan valores diferentes de 4/
(4-2)=2.
f(W) = j(WW) =
Sera errneo concluir que, de la misma forma que 3 y 1 son, en nuestro ejem-
plo, los valores medios de W1 y W2, 3,16 y 1,48 son las varianzas respectivas; ni
lo son, ni tampoco es la matriz de covarianzas. Eso s, sta resulta de multi-
plicar a por una constante k de manera que la transformacin equivalente sea
Pgina 192.
28
58
(i = 1,2) en donde fi es la respectiva funcin de densidad marginal y mi la media.
Los resultados son:
Adems, calculamos que Cov(W1, W2) = E[(W13).(W2 1)] 2,91 29; por tanto,
valor que coincide con el que figura en (1.5). Por tanto, la matriz de covarian-
zas 0 y su inversa, son, respectivamente:
lo que nos dice, por comparacin con (1.3), que la constante k = 2 (igual a 4/
(4-2)):
0 = 2
29
Obtenemos este resultando integrando el producto (w1 3)(w2 1) por la funcin de densidad
conjunta, (expresin (1.5)).
59
Esta funcin expresa la densidad de una t bivariante con matriz de covarian-
zas 0, 4 grados de libertad y media (3, 1); su coeficiente de correlacin es
0,673 y las desviaciones tpicas de los mrgenes, 2,514 y 1,720. Es la misma
funcin (1.5), aunque con apariencia algo diferente, porque ahora figura ex-
presamente la matriz de covarianzas.
Por ejemplo, el coeficiente de correlacin mide la dependencia lineal; los coeficientes S (de
31
Spearman) y (de Kendall) que veremos ms adelante, miden la dependencia, entendida como con-
cordancia.
60
distribuciones marginales. Uno de estos tipos de dependencia es la concordancia.
Dos variables X e Y son concordantes si los valores grandes (pequeos) de X tien-
den a presentarse con valores grandes (pequeos) de Y. Este tipo de dependencia
no est en funcin de la forma en que se distribuyan, por separado, X e Y.
Sean 4 pares (X, Y) como los que figuran en las dos primeras columnas del
Cuadro 1.1.
X Y N(X) N(Y) E(X) E(Y)
0,1 0,3 0,54 0,62 0,095 0,26
0,2 0,5 0,58 0,69 0,18 0,39
0,3 0,4 0,62 0,65 0,26 0,33
1,0 1,2 0,84 0,88 0,63 0,70
Cuadro 1.1: La concordancia es invariante a la escala (ejemplo).
61
consecuencia, concordantes), al menor valor posible de X correspondera el me-
nor posible de Y, al siguiente de X el siguiente de Y y as sucesivamente hasta los
mayores de ambas, que iran igualmente emparejados. Y lo mismo ocurrira con
los valores u = FX(x) y v = FY(y). Pues bien, el coeficiente de correlacin de rangos
de Spearman S(X,Y) es igual al coeficiente de correlacin lineal entre las varia-
bles U y V (siendo stas las transformadas de X e Y segn sus respectivas funcio-
nes de distribucin marginal, es decir: U = FX(X), V= FY(Y)) y, por tanto, S(X,Y) =
(U, V) = 12cov(U, V). El valor de este coeficiente no depende de dichas distribu-
ciones marginales, porque slo mide el grado de coincidencia en el orden de X y
de Y. As, si tenemos slo 3 pares (X,Y) el coeficiente ser S = 1 tanto en el caso
(1, 2), (2, 3), (3, 4) como si los pares son (1, 2), (20, 3), (21, 4).
El valor de , referido a una muestra de n pares, es, si tenemos en cuenta que para
variables continuas el nmero de pares posibles es la suma de concordantes (c)
ms discordantes (d):
en donde (X1, Y1) y (X2,Y2) son dos pares independientes de la variable bivariante
(X, Y).
62
Dos variables aleatorias (X,Y) que estn asociadas de tal manera que para cuales-
quiera (X1,Y1) y (X2,Y2), (X1<X2) (Y1<Y2), son comonotnicas. Si, por el contrario
(X1<X2) (Y1>Y2), son contra-comonotnicas. La comonotonicidad implica, por
tanto, que = S = 1, mientras que el caso opuesto significa que ambos coeficientes
valen 1. Es decir, si X e Y son concordantes, son comonotnicas y ambas medi-
das son iguales a 1. Y lo contrario ocurre si son discordantes.
Segunda: Segn acabamos de ver, ambas medidas estn comprendidas entre 1 y +1.
63
b) Para T(X) montona decreciente, (T(X),Y) = (X,Y).
= 2.Prob[(X1X2)(Y1Y2)>0] 1 = 4.Prob[(X1<X2)(Y1<Y2)] 1
Por tanto,
64
II. CPULAS
Sea F(x1 x2 ... xd) la funcin de distribucin conjunta de d variables continuas, X1,
X2, ..., Xd. Esta funcin describe el comportamiento conjunto de dichas variables y
as, nos permite obtener la probabilidad de que una variable no supere cierto va-
lor, con independencia del que tomen las otras (distribucin marginal) o cono-
ciendo el valor asumido por stas (distribucin condicionada) o, por ltimo, que
las otras tomen valores inferiores o superiores a los que se especifiquen. Es decir,
que dentro de la funcin F(x1 x2 ... xd) hay dos tipos de informacin:
65
b) Sobre el comportamiento conjunto de las variables. Este comportamiento radi-
ca en la estructura de dependencia entre las variables, la cual viene modeliza-
da por la cpula: la naturaleza no paramtrica de la dependencia entre dos varia-
bles aleatorias es expresada por la cpula.34
ui = Fi(xi) (i = 1, 2, ... d)
34 Nelsen, (2006), pgina 52. Es decir, la cpula es un concepto ajeno o distinto a la distribucin
marginal y contiene la ley de dependencia entre las variables intervinientes en la distribucin mul-
tivariante.
35 Observemos que F -1(u ) se introduce en la funcin de distribucin conjunta para anular o desha-
i i
cer el contenido de sta en cuanto a la distribucin marginal de xi.
66
Las funciones de distribucin marginal de X e Y son, respectivamente:
u = FX(x) = 1 exp(1.x)
v = FY(y) = 1 exp(2.y)
x = FX-1(u) = LN(1u)/1
y = FY-1(v) =LN(1v)/2
C(u,v) = u.v
F(x1, x2, ... , xd) = F(F1-1(u1), F2-1(u2), ... , Fd-1(ud)) = C(u1,u2 , ... , ud)
La funcin C(u1,u2, ... , ud) o cpula se llama as porque conecta o acopla las distribu-
ciones marginales ui = Fi(xi) a la estructura de dependencia, en que la cpula consiste,
para llegar a la funcin de distribucin conjunta. La cpula es una funcin de distribu-
cin de variables uniformes (0,1), U1, ... , Ud y convierte cualquier punto del cubo de
lado unitario en el espacio de d dimensiones, en un nmero real del intervalo [0, 1]. Es
decir, [0, 1]d [0, 1]. En efecto, tengamos en cuenta que las Ui = Fi(Xi) son variables
uniformes (0,1) y que la transformacin ui = Fi(xi) es montona creciente:
67
Veamos ahora dos ejemplos especialmente importantes.
x = -1(0,6) = 0,253
y = -1(0,3) = 0,524
68
Figura 2.1: Funcin de densidad de la cpula gaussiana con parmetro = 0,5.
69
y su funcin de densidad,
t3-1(0,6) = 0,277
t3-1(0,3) =0,584
70
Figura 2.2: Funcin de densidad de la cpula t con parmetros n = 3 y = 0,5.
71
Esta nueva funcin de distribucin tiene la misma cpula que la bivariante original
F(x,y) pero nuevos mrgenes, X(x) y Y(y). La posibilidad de combinar una estruc-
tura de dependencia con distintos mrgenes o viceversa, hace de las cpulas un
instrumento generador de gran nmero de modelos multivariantes, para adaptar-
se a la realidad lo mejor posible. En el Epgrafe 7 de este mismo Captulo desarro-
llaremos este punto, pero ahora propondremos unos ejemplos en d=2 dimensio-
nes, de cmo la separacin entre estructura de dependencia o cpula, por un lado
y distribuciones marginales, por otro, permite unir diversos modelos de una y
otras para dar lugar a una amplia variedad de distribuciones multivariantes.
(a): Cpula gaussiana con =0,4; mrgenes (b): Cpula gaussiana con =0,4; mrgenes
lognormales con media y desviacin tpica 2 exponenciales con media y desviacin tpica 2
72
En resumen:
Una cpula C(u) = C(u1 u2 ... un) es una funcin de distribucin de n variables alea-
torias uniformes (0,1). Por tanto, a cada vector en [0,1]n hace corresponder un
nmero en [0,1]. Para que C(u1 u2 ... un) sea, en estas condiciones, una funcin de
distribucin, ha de reunir los requisitos siguientes:
73
1. Ser no-decreciente en cada uno de sus componentes ui36.
2. La distribucin marginal de Ui. que es C(1, 1, ... ui ... 1), equivale a ui. Por otra
parte C(1,1, ... 1) = 1 como no poda menos de suceder (hacemos ui=1).
3. La funcin C es n-creciente. En general, se dice que una funcin F(x1, x2 , ... , xn)
es n-creciente si, para cualquier hiper-rectngulo [x1(1), x1(2)]. [x2(1), x2(2)]. . . [xn(1),
xn(2)] (siendo xi(1) xi(2) para i =1,2, ...,n) cuyos vrtices se encuentran en el domi-
nio de definicin de F se verifica que:
F(x1(2), x2(2), x3(2), ... , xn-1(2), xn(2)) F(x1(1), x2(2), x3(2), ... , xn-1(2), xn(2)) F(x1(2), x2(1),
x3(2), ... , xn-1(2), xn(2)) F(x1(2), x2(2), x3(1), ... , xn-1(2), xn(2)) F(x1(2), x2(2), x3(2),
... , xn-1(2), xn(1)) + F(x1(1), v2(1), x3(2), ... , xn-1(2), xn(2)) + F(x1(2), x2(1), x3(1), ... , xn-1(2),
xn(2)) + . + F(x1(2), x2(2), x3(2), ... , xn-1(1), xn(1)) F(x1(1), x2(1), x3(1), ... , xn-1(1),
xn(1)) 0.
Veamos lo anterior en dos dimensiones, para lo que nos valdremos de la Figura 2.4.
36 Puede sustituirse esta primera propiedad por la de ser C(0,...,0, u , 0,...,0)=0. El resultado es
i
equivalente. Para el caso bidimensional, tenemos que, al ser C(1,v)=v (segn exige la Condicin 2),
C(1,0) = C(0,1) = 0 y por ser C no decreciente en u y en v resulta que C(u,0) y C(0,v) han de ser cero
necesariamente (al ser u,v 1). Y, recprocamente, para el rectngulo (0, u) x (v, v+k) (k > 0), resulta
que C(u, v+k) C(0, v+k) C(u, v) + C(0,v) 0 o bien C(u, v+k) C(u, v) 0 lo que demuestra que
C(u,v) es no-decreciente en v y lo mismo para u.
74
Prob(A) = C(b,c) C(a,c).
Figura 2.4.
75
continuas cuya cpula es C(u1, u2 , . . . , ud). Si aplicamos a cada variable Xi una
transformacin Ti montona creciente, las variables Yi = Ti(Xi) tienen la misma
cpula C(u1 , u2 , . . . , ud).
***
C[F1(T1-1(x1)) . . . Fd(Td-1(xd))].
***
Los argumentos Gi(xi) pueden tomar cualquier valor comprendido en [0. 1] al ser
continuas las variables Xi.
76
efecto en el seno de la distribucin conjunta consiste en transformar a las varia-
bles uniformes U1, U2, ... Ud segn la respectiva funcin montona creciente Ui =
Fi(Xi). El coeficiente de correlacin depende de la cpula, pero tambin de las
distribuciones marginales. En cambio, los coeficientes de Kendall y de
Spearman estn en funcin slo de la cpula.
Sean ahora dos variables X1 y X2 cuya cpula es gaussiana de parmetro ; sus dis-
tribuciones marginales son F1(x1) y F2(x2). El coeficiente de correlacin entre ambas,
es decir, (X,Y) no es (salvo que dichas distribuciones marginales se identificaran
con la normal). En estas condiciones, existe una frmula que relaciona el coeficien-
te entre X e Y, y el parmetro de su cpula; dicha frmula es:
39 La demostracin de estas frmulas puede verse en McNeil, Frey y Embrechts, (2005), pgina 215.
77
Ejemplo 2.6: Variable normal bivariante a la que se aplica transformacin
Y = (X). Significado del parmetro
78
Para obtener el anterior resultado tenemos en cuenta que la densidad de U es
constante e igual a 1 en [0, 1], al ser variable uniforme (0,1).
En general,
Por otra parte, el lmite del denominador dividido tambin entre hn-1, es la funcin
de densidad conjunta de las variables U2,...,Un; si llamamos C1(u2, u3, . . . , un) a la
funcin de distribucin marginal de dichas variables, la correspondiente funcin
de densidad conjunta ser
79
siempre el mismo valor, es decir, U = V. En tal caso Prob [U u; V v]= Prob = [U
u; U v] = min(u,v), que es el valor mximo posible de la funcin de distribucin
para el punto (u,v). Si V = 1 U tenemos el caso opuesto:
El lmite inferior es, por tanto, u + v 1, sin que pueda ser negativo. De ah que las
cotas sean:
max [(u + v 1), 0] C(u,v) min(u,v).
80
Figura 2.6(a): Cpula mxima. Figura 2.6(b): Cpula mnima.
Figura 2.7(a): Cpula mxima, curvas de nivel Figura 2.7(b): Cpula mnima, curvas de nivel
81
La cpula M(u,v) significa la dependencia total entre dos variables X e Y, entendida
dicha dependencia en el sentido de comonotonicidad41, mientras que W(u,v) re-
presenta la contramonotonicidad: Cuando X es grande, Y es pequea y viceversa.
3. CPULAS Y DEPENDENCIA
41 Ver Captulo I.
42 McNeil, Frey, Embrechts, (2005), pgina 200.
82
***
Para demostrar esta equivalencia observamos que, dados dos puntos x, z cua-
lesquiera en la recta, (Figura 2.8) la distancia entre ambos es
Figura 2.8.
(2.8)
Primer trmino: A + B + C + D
Segundo trmino: C + D
Tercer trmino: B + C
Cuarto trmino: C
83
Consideremos ahora un par de variables (X1, X2) con funcin de distribucin
conjunta F(x1,x2) y mrgenes F1 y F2; el coeficiente de correlacin lineal entre
X1 y X2 es ; hacemos tambin E(X1) = 1; E(X2) = 2.
Si (Z1 Z2) es un par independiente del anterior, con la misma distribucin con-
junta, tendremos:
Observemos ahora que, por ejemplo E[1(s>Z1)] es el valor medio de una varia-
ble dicotmica Z1 que puede ser mayor o menor que s; por tanto, dicho valor
medio equivale a F1(s). De aqu que, en virtud de la independencia entre los
pares (X1, X2) por un lado y (Z1, Z2), por otro,
Figura 2.9.
84
Por otra parte,
puesto que los valores medios de los productos (X1 1) (Z2 2) y (X2 2)
(Z1 1) se anulan y los valores medios de (X1 1)(X2 2) y de (Z1 1)(Z2
2) coinciden, al tener (X1, X2) la misma distribucin conjunta que (Z1 , Z2).
Por tanto,
***
85
Nelsen (2006) propone un interesante ejemplo de aplicacin de esta medida,43 re-
ferida a dos variables (X,Y) perfectamente dependientes entre s de forma que
ambas son las coordenadas de puntos de una circunferencia de radio 1, siendo
constante la densidad en dichos puntos. Pues bien, en este caso el resultado es
= 0,25. La perfecta dependencia no conduce en este caso a un coeficiente de de-
pendencia igual a la unidad porque X e Y no son comonotnicas.
Como observa Nelsen, esta integral puede interpretarse como el valor medio de
la variable C(U,V), en donde C es la funcin de distribucin conjunta de dos varia-
bles uniformes (0, 1) 44.
43 Pgina 209.
44 Nelsen (2006), pgina 162.
45 Figura 2.10(a), cpula Frank con = 5,74; Figura 2.10(b), cpula Gumbel con = 2 .
86
la distribucin (superior o inferior, segn los casos; en este ejemplo es de mayor
inters la dependencia en el extremo superior). sta es la dependencia en colas,
cuyo coeficiente se define (para las variables X, Y con distribucin conjunta F(x,y)
y mrgenes FX y FY como:
(a) (b)
Figura 2.10: Simulacin de 500 puntos de dos cpulas diferentes, ambas con coeficiente =0,5
siendo .
87
igual manera a lo que ocurre con los coeficientes y S. Y tambin significa que el
valor que tome radicar en el comportamiento asinttico de dicha cpula.
y, por tanto,
***
Adems, ; e igualmente:
88
Recordemos ahora expresin (2.5) y razonamiento que nos condujo a ella
que si F(x,y) es una funcin de distribucin derivable,
(siendo f(x) la funcin de densidad de X). Si aplicamos esto a una cpula C(u,v)
obtenemos = Prob(VvU=u).1; por tanto,
y en consecuencia,
47Comparemos este resultado con la expresin (2.11), en donde el condicionante era V>, pero el
lmite no se duplicaba, como ocurre aqu.
89
Por tanto, la variable (YX=x) es es decir, una variable nor-
mal con media .x y desviacin tpica y la probabilidad buscada ser,
llamando a la funcin de distribucin normal (0, 1),
(salvo que =1). Por tanto, en la cpula de Gauss no hay dependencia en las
colas si <1.
90
y de la que se resta .x. Es decir:
es una variable A(x,y) que se distribuye como una t con +1 grados de li-
bertad.
***
48 Si dicha probabilidad es menor que 1, ello significa que la expresin (2.13) dar un valor del co-
eficiente mayor que cero. Hacer que y tienda a infinito equivale a que, en dicha expresin, 1.
91
El Cuadro 2.1 contiene valores de para algunos pares (, ).
49 Los ceros de esta columna no significan ausencia total de dependencia, sino valores inferiores a
una milsima.
92
4. CPULAS ARQUIMEDIANAS
Las arquimedianas son una clase de cpulas cuya gran importancia radica en dos
circunstancias: En primer lugar, sus interesantes propiedades, que hacen de ellas un
instrumento de fcil aplicacin, al menos en el caso bivariante; y en segundo lugar, su
notable variedad, que posibilita construir modelos adaptables a muchas situaciones.
Si hacemos (t) =LN(t) resultar que -1(t) = e-(t). Por tanto u.v= e-(u)(v) y esto,
expresado en funcin de (.) y de su inversa da lugar a:
Es decir que -1(t) = exp(t1/) y por tanto la cpula resultante es, en virtud de (2.14):
93
Esta funcin resulta ser una cpula (llamada cpula de Gumbel) que veremos ms
adelante.
3. Que (1) = 0.
En la Figura 2.12 se representan dos tipos de funcin [designada como (t)] que
renen estas caractersticas. La funcin ( en la grfica) es generador estric-
to si (0) = ; en este caso -1 est definida para cualquier valor de la ordenada
positiva. En otro caso (Figura 2.12(b)) -1 slo est definida para valores de la or-
denada comprendidos en [0,(0)] y de ah que haya que definirla convencional-
mente para ordenadas superiores a (0); ello se logra haciendo -1(t) = 0 para
t(0). A esta inversa con definicin parcialmente convencional se denomina
pseudo-inversa y ha de designarse por [-1], porque la notacin -1(t) es espec-
fica del generador estricto. Por tanto, si el generador es estricto, [-1] = -1.
Una vez caracterizada la funcin generador (t) falta demostrar que las funciones
de la forma (2.14) obtenidas con tales generadores renen los requisitos propios
de las cpulas. El de ser creciente52 puede demostrarse que se cumple por ser el
generador convexo (Demostracin en Nelsen (2006), pgina 111).
Adems:
94
Figura 2.12(a): Generador estricto. Figura 2.12(b): Generador no estricto.
y, por tanto,
95
Otra caracterstica de las cpulas arquimedianas es que,
el signo de desigualdad est justificado por ser (y, por tanto, tambin [-1]) decre-
ciente.
La cpula de Gumbel
Su generador es
(2.18).
53 Slo examinaremos tres casos de cpulas arquimedianas. Nelsen (2006), pginas 116 a 118,
incluye una tabla con 22 tipos de estas cpulas.
54 Epgrafe 6.
96
Figura 2.13(a): Densidad Gumbel, =2. Figura 2.13(b): 1000 puntos Gumbel, =2.
La cpula de Clayton
Con la cpula de Clayton ocurre lo contrario que con la anterior, de Gumbel: Las
mayores densidades estn en el ngulo inferior izquierdo y de ah que se utilice
cuando la dependencia entre las variables sea mayor en el caso de que valores
bajos de una tiendan a coincidir con valores bajos de la otra variable (Figura 2.14).
97
Figura 2.14(a): Densidad Clayton, =2. Figura 2.14(b): 1000 puntos Clayton, =2.
La cpula de Frank
Esta cpula es simtrica radialmente, igual que ocurre con las cpulas elpticas;
por tanto, las densidades estn igualmente repartidas en zonas prximas a (0,0) y
a (1,1) del cuadrado [0, 1]2. En las Figuras 2.15(a) y (b) se observan estas caracte-
rsticas.
y la funcin C(u,v):
98
Figura 2.15(a): Densidad Frank, =8. Figura 2.15(b): 1000 puntos Frank, =8.
55 Pgina 114.
99
Nelsen se ocupa de esta cuestin56 y explica que el origen est en el axioma de
Arqumedes: Dados dos nmeros reales positivos a, b, existe siempre un entero n
tal que n.a>b. Y las cpulas arquimedianas tienen una propiedad que recuerda
este axioma. De similar forma que la multiplicacin de n por a genera un nuevo
nmero, la introduccin de u y v en la funcin C(u,v) da lugar a un tercer nmero.
Si reiteramos esta operacin a partir de un nmero inicial u [0,1] obtenemos,
aplicando (2.14) y (2.15):
Ahora bien, lo sorprendente es que no slo las arquimedianas, sino otras muchas
cpulas cumplen la condicin que atribuye su denominacin a las primeras. Sea
C(u,v) una cpula tal que para todo punto (u,v) perteneciente al cuadrado (0,1]2
cumpla la condicin de ser menor que la cpula mxima, es decir:
C(u,v)< min(u,v)
100
Sobre esta base podemos elaborar la sucesin siguiente:
C1 = C(u, u)< u
C2 = C(u, C1)< C1
C3 = C(u, C2)< C2
..........
Cn+1 = C(u,Cn)< Cn
Tenemos, por tanto, u> C1 > C2> . . . >Cn > . . . que es una sucesin de nmeros
positivos, infinita, decreciente e inferiormente acotada en 0 (puesto que toda c-
pula ha de cumplir C(u, 0) = 0. Por tanto, la sucesin Cn0 y siempre ser posible
encontrar un nmero entero n tal que Cn < v, por pequeo que v sea.
Abordamos ahora el estudio de una frmula de gran utilidad que es, quiz, una de
las principales razones de que el uso de las cpulas arquimedianas est tan ex-
tendido. Esta frmula es:
(2.23)
101
en donde KC(t) = Prob[C(U,V) t] y +(t) es la derivada por la derecha de la funcin
generador.
***
57La cpula emprica ser objeto de estudio en el Captulo III. De momento podemos identificar
cpula emprica con la versin muestral de la cpula.
58 La demostracin que sigue se basa en la de Nelsen, (2006) pginas 127 y siguientes.
102
Figura 2.16.
B = [0,5<u<0,6] x [0,2<v<0,3].
La M-medida de B ser:
VM(B) =
M(0,6; 0,3) M(0,6; 0,2) M(0,5; 0,3) + M(0,5; 0,2) =0,3 0,2 0,3 + 0,2 = 0.
Ello nos indica que en B no hay probabilidad, porque ningn punto de la recta
que la concentra est contenido en B. La cpula M presenta, encima de B, la
forma de una superficie plana descendente hacia el eje u.
103
Figura 2.17.
59 En esta grfica se representa como (kw/n)^[-1] lo que en el texto designamos como (kw/n)[-1]
104
Ordenada Abscisa correspondiente
v(n) =0 u(n) = 1
v(n-1) =w/n u(n1)
v(n-2) =2w/n u(n2)
............. ............
v(k+1) = (nk1)w/n u(k+1) = u[n(nk1)]
v(k) = (nk)w/n u(k) = u[n(nk)]
v(k-1) = (nk+1)w/n u(k-1) = u[n (nk+1)]
.............. ..............
v(1) = (n1)w/n u(1) = u[n(n1)]
v(0) = nw/n u(0) = u[nn] = t(0)
Figura 2.18.
105
o bien, despejando v:
Si u = 1,
Y si v = 1,
1 = -1[(t(0))(u)]
Figura 2.19.
106
u(k) = -1 ; la ordenada de los dos vrtices superiores es, en virtud de
(2.24) y teniendo en cuenta que [t(0)] = w (ver Figura 2.18),
El rectngulo R(k)= Rk es, por tanto: Rk = [u(k1), u(k)] x [0, v(k1)]; y por con-
siguiente, VC(Rk) ser:
107
La derivada obtenida es la de la funcin -1(w) y es la derivada por la izquier-
da. Equivale a la C-medida de la superficie rayada (Figura 2.20) y tambin
puede expresarse as:
KC(t) = t (t)/(t)+.
Figura 2.20.
***
108
Por ejemplo, en una cpula de Gumbel con =2, la expresin anterior ser:
= 4.E[C(U,V] 1
(2.25)
Por otra parte, (2.25) permite obtener fcilmente el valor de en las cpulas ar-
quimedianas, como vamos a ver en seguida.
(2.26)
(2.27)
c) En el caso de la cpula Frank el generador viene dado por (2.21). Su derivada es:
109
Demostraremos a continuacin que, en este modelo de cpula, la expresin de
en funcin de viene dada por (2.28).
***
Definimos ahora las funciones A(x) y B(x) y la integral C, de manera que sea
110
Por tanto, C = D1().
(2.28)
***
Obtenidas las frmulas para en estos tres modelos de cpulas, finalizamos este
Epgrafe con la expresin de los coeficientes de dependencia en colas. Al tratarse
de cpulas explcitas, no hay ms que aplicar (2.12a) y (2.12b) para sup y inf. Para
llegar al coeficiente respectivo obtenemos la funcin f() = C(, ) y su derivada y
aplicamos a continuacin la regla de LHpital en el clculo del lmite en las cita-
das expresiones (2.12a y 2.12b).
a) Cpula de Gumbel.
111
(2.28a)
b) Cpula de Clayton.
(2.28b)
c) Cpula de Frank.
112
La simetra de esta cpula es congruente con la igualdad de ambos coeficientes
de dependencia que, adems, son nulos como ocurra en la cpula gaussiana.
Ya hemos visto que el Teorema de Sklar61 nos dice que a partir de una cpula y de
unas distribuciones marginales puede obtenerse una distribucin multivariante.
ste es el origen de la aplicacin ms interesante de las cpulas. Con ellas pode-
mos formar funciones de distribucin. Al separar las piezas que integran la dis-
tribucin multivariante (cpula por una parte y mrgenes por otra) y ser stas
susceptibles de combinarse entre s, podemos formar una gran variedad de fun-
ciones de distribucin.
Sea C(u1,u2, . . . , ud) una cpula de d dimensiones y las Fi (i = 1, 2, ... , d) sus res-
pectivas distribuciones marginales continuas; si hacemos ui = Fi(xi) obtenemos
C(u1,u2, . . . , ud) = C(F1, F2, ..., Fd) = F(x1, x2, ..., xd)
(2.29)
en donde las fi(xi) son las densidades marginales. Es decir, que la funcin de densi-
dad conjunta es igual a la funcin de densidad de la cpula, multiplicada por el pro-
ducto de las densidades marginales. La expresin (2.29) ha de manejarse con cui-
dado porque, tal como acabamos de escribirla, en el primer miembro aparecen x1,
x2, ..., xd como variables independientes y en el tercero, adems de ellas, u1, u2, ...,ud.
61 Epgrafe 2.
113
univariante); las funciones r(u) y s(v) son las transformadas inversas segn dicha
normal univariante, es decir que:
r(u) = -1(u), s(v) = -1(v),
siendo (.) la funcin de distribucin normal univariante. Si nos ceimos al caso
bivariante y deseamos aplicar (2.29) cuando se trata de dos variables X, Y con c-
pula gaussiana y mrgenes con funciones de distribucin FX(x) y FY(y) y densida-
des fX(x), fY(y) resultar:
La funcin f(x, y) representada en la Figura 2.3(b) por medio de sus curvas de nivel
es la densidad de una distribucin con cpula Gauss con = 0,4 y mrgenes expo-
nenciales62 con parmetro a=2.
114
Las dos ltimas fracciones equivalen, cada una, a la correspondiente densi-
dad exponencial, puesto que u = 1ex/2; como la densidad es ex/2/2, sta se
identifica con la primera de ambas fracciones y otro tanto cabe decir respecto
de Y en cuanto a la ltima fraccin.
Camino parecido hay que seguir para obtener la funcin de densidad de la otra
cpula implcita que hemos visto, la cpula t.
63 En Venter (2002) hay varios grficos que permiten comparar las densidades de diversas dis-
115
de la densidad de la cpula, c(u,v), por el de las densidades marginales. En este
caso:
Figura 2.21: Curvas de nivel de la densidad de la distribucin con cpula de Gumbel de parmetro
=1,355 y mrgenes exponenciales de media y desviacin tpica iguales a 2.
Si, por ejemplo, =1,355 (con lo que =0,26264) resulta, para a=2, la densidad f(x,y)
representada por medio de curvas de nivel que aparecen en la Figura 2.21.
64 Hemos elegido = 0,262 por ser ste el de la cpula gaussiana de las distribuciones representa-
116
8. CPULAS DE SUPERVIVENCIA
Figura 2.22.
117
Figura 2.23: Funcin de supervivencia conjunta de la distribucin con cpula Gumbel con = 2 y
mrgenes exponenciales de parmetro 1.
Su grfica, con los mismos datos para la estructura de dependencia que los de la
Figura 2.23, est en la Figura 2.24.
118
Figura 2.24: Funcin (u,v) = Prob(U>u; V>v) correspondiente a la cpula de Gumbel de
parmetro = 2.
Que (u,v) es una cpula puede confirmarse al comprobar que cumple las propie-
dades de estas funciones:
En efecto,
119
ya que: 1a > 1 b y 1 c > 1d y que C es una cpula.
En primer lugar, que la funcin (u,v) = Prob(U>u; V>v) no es una cpula, segn
se ve claramente en la Figura 2.24. Y en segundo trmino, que en esta misma fun-
cin (u,v)=1 u v + C(u,v) hemos introducido 1u = y 1v = en el lugar de
u y v respectivamente, con lo cual (1u, 1v) = u + v 1 + C(1u, 1v). Es ste
un paso fundamental en la gnesis de porque significa pasar a contemplar o
medir a desde el ngulo superior derecho del cuadrado [0, 1]2 y no, como origi-
nariamente, desde el inferior izquierdo. Es decir, que en aparecen u y v que,
evidentemente, son variables cuyo punto cero se encuentra en dicho ngulo infe-
rior izquierdo. Al introducir y en lugar de u y de v trasladamos este origen
desde (0,0) a (1,1) y con ello logramos una funcin que guarda, con respecto a C,
la simetra que veremos ms adelante.
120
(a) (b)
(c) (d)
Figura 2.25: Densidad de la cpula de Clayton (a) y (c) y de su cpula de supervivencia (b) y
(d)para =3.
121
Ya hemos visto que (u,v) = (F,G) = ( , ) y que, por tanto,
( , ) = (F,G) = (x,y).
La densidad en el punto (u,v) segn (u,v) coincide con la densidad en (1u, 1v)
segn C(u,v). Por eso algunos autores denominan a la cpula de supervivencia
flipped copula65. Una forma alternativa de obtener (u,v), seguira los siguientes
pasos:
122
Queda comprobado, por tanto, que la cpula correspondiente a la densidad sim-
trica (es decir, a c(1 u, 1v)) es precisamente la cpula de supervivencia.
Para hacer tal comparacin, partimos de la hiptesis de que ambas cpulas cu-
yos parmetros designamos como G y C segn se trate de la cpula de Gumbel
o de la de Clayton tienen el mismo coeficiente de Kendall. Si ponemos el res-
pectivo parmetro en funcin del coeficiente [expresiones (2.26) y (2.27)] y susti-
tuimos cada expresin resultante en (2.28a) y (2.28b), resulta, designando a los
respectivos coeficientes de dependencia en la cola superior, de Gumbel y de
Clayton (cpula de supervivencia) como G y C:
G = 22 1-
C = 2 (1-)/(2)
123
Figura 2.26: Dependencia en la cola derecha de la cpula de Gumbel y de la cpula de
supervivencia de Clayton, en funcin del coeficiente .
Hemos visto ya algunos ejemplos de cpulas implcitas (la gaussiana, la cpula t).
Expondremos a continuacin unos casos ms en que, por un camino u otro, se
llega a la cpula a partir de la distribucin conjunta.
124
funcin de distribucin, G(.). Realizamos un experimento que genera un valor para
cada variable. Ambos resultados son independientes entre s. Pero si conocira-
mos el menor de estos dos valores y el importe del mismo, ganaramos informa-
cin acerca del otro, que necesariamente no sera inferior al primero. Por tanto,
las variables N = min(X, Y) y M = max(X, Y) no son independientes y nos propone-
mos obtener su funcin de distribucin conjunta y su cpula. A dicha funcin de
distribucin conjunta la designaremos F(x,y) = Prob(Nx; My). Pues bien, F(x,y) es:
x = G-1(u)
y = G-1(v)
125
y sustituyendo x e y en F(x,y) obtenemos la cpula. Si admitimos que G(.) es una
funcin de distribucin uniforme (0,1), entonces u =2xx2; v =y2. En consecuencia,
(a) (b)
Figura 2.27: Cpula correspondiente al primer ejemplo examinado.
126
Figura 2.28: Densidad de la cpula correspondiente al primer ejemplo examinado.
Para valores (u,v) tales que la densidad es nula, pues esta des-
igualdad significa que y<x; si N toma un valor en el entorno de x, M no puede to-
mar un valor prximo a y.
Por ltimo, si x es bajo e y alto (caso que se corresponde con el ngulo inferior
izquierdo de la grfica, cerca del punto 0,1) la densidad es baja porque estamos
ante un supuesto en que N est prximo a 0 y los valores posibles de M pueden
oscilar en casi todo su rango.
127
Figura 2.29: 2000 puntos simulados de la cpula representada en las Figuras 2.27 y 2.28.
En la Figura 2.29 tenemos 2000 simulaciones de esta cpula, que permiten obser-
var su congruencia con la forma de la densidad segn la Figura 2.28. Curva de
discontinuidad y carcter no simtrico de esta cpula.
Prob(Zi t) =1 eit.
128
Llamamos X e Y a las variables representativas del tiempo o duracin del funcio-
namiento de los elementos 1 y 2. Por tanto, X=min(Z1, Z12); Y=min(Z2 , Z12). La pro-
babilidad de que (X>x, Y>y) ser:
Hacemos: = u, = v:
(x,y) = (1 u, 1 v) = (u,v) =
129
Hacemos ahora
12/(1 + 12) = ; 12/(2 + 12) = .
Si ambos parmetros son bajos, ocurre lo contrario y se produce una gran disper-
sin entre X e Y.
130
***
Si u v C(u, v) = u.v1-
Si u v C(u, v) = v.u1-
Figura 2.30.
Pues bien, el valor de la integral para cualquier par (u,v) tal que v<u/, es decir,
131
ser la parte continua de C(u,v) y equivaldr a la suma de las probabilidades
que, segn la cpula C(u,v) corresponden a los recintos A1 + A2 + B. En otras
palabras, equivaldr a la suma de la C-medida68 de A1, ms la de A2, ms la
de B, que designamos respectivamente por C(A1), C(A2), C(B). La curva de la
Figura 2.30 tiene como ecuacin v = u/; u(v) es igual a v/, siendo v la ordena-
da del recinto de integracin. Es decir, u(v) es la abscisa que corresponde,
segn la curva, a dicha ordenada mxima. La probabilidad del recinto de inte-
gracin, es decir, C(u, v) (para v<u/) est, pues, compuesta por los tres su-
mandos C(A1), C(A2) y C(B) siguientes:
Por tanto,
(2.32)
132
Observemos ahora que v.u1- es la expresin de la cpula para v u. Por tan-
to, el segundo trmino, con signo negativo, es el importe de la probabilidad
concentrada en el segmento de curva v = u/ comprendido en el recinto de
integracin. Y observemos tambin que, al ser los mrgenes continuos, en
dicha curva no hay ningn punto que tenga masa de probabilidad. El segundo
trmino de A puede tambin escribirse:
(2.33)
Es decir, que para (u,v) tal que u > v el valor de C(u,v) = v.u1- est integrado por:
De manera similar obtendramos, para (u,v) tal que u < v, que el valor de
C(u,v) = u.v1- se compone de:
Parte singular =
Parte continua, el resto.
***
133
La parte singular es
La Figura 2.31 muestra cmo cambia el aspecto de los puntos simulados de esta
cpula al variar los parmetros. Si stos estn ms prximos a la unidad, los pun-
tos aparecen ms concentrados en la curva, lo que indica una ms alta dependen-
cia. Y si > , como ocurre en la parte (a) de la Figura, la curva es convexa desde
abajo.
134
Hagamos que dicho parmetro sea, a su vez, aleatorio y distribuido segn la fun-
cin gamma (, ). X e Y son independientes entre s, si sus distribuciones respec-
tivas se condicionan a un valor de . Pero si vara, X e Y son dependientes, al
estar su distribucin en funcin del parmetro comn .
135
nos permiten concluir que el anterior resultado conduce a la cpula de supervi-
vencia de Clayton72. En efecto, si , entonces la cpula
de supervivencia, denominada en este caso cpula HRT73, es:
Otra forma de construir una cpula consiste en combinar otras para llegar a una
nueva. ste es el caso de la cpula de Cuadras-Aug, que consiste en la media
geomtrica de la cpula mxima y la cpula de independencia, es decir:
136
ya que el hecho de que M y N no superen a x y a y requiere que cada una de las n
variables (independientes) de cada serie sean menores o iguales que x e y. Adems,
la probabilidad de que M x y N y ser:
u = FM(x)
v= GN(y).
es decir, la cpula de las variables (Xi, Yi) y la de las variables (M, N) son la misma.
En tal caso de dice que la cpula C(u,v) es max-estable. Por ejemplo, la cpula de
Gumbel es max-estable:
C(u,v) = (u + v 1)1/
(C(u1/n,v1/n))n = (u/n + v /n1)n/ C(u,v)
137
En La Figura 2.32 se representa la densidad de (C(u1/n,v1/n))n (n = 5), siendo C la
cpula de Clayton con parmetro = 3, para que sea comparable con la densidad
representada en la Figura 2.25(a), a la misma escala. La comparacin sugiere que
al aumentar n, la cpula resultante se ira aproximando a la cpula de indepen-
dencia, puesto que el efecto es equivalente a si hacemos que disminuya el par-
metro en el modelo de Clayton.
138
C(u,v) = limn (C(u1/n,v1/n))n = limm (C(u1/(mr),v1/(mr)))mr =
= (C(u1/r,v1/r))r.
139
III. NMEROS Y CPULAS
1. LA CPULA EMPRICA
El problema que plantea la cpula emprica es que las variables estadsticas a las
que se refiere no son directamente observables. La muestra consiste en n pares
(xi, yi) de manera que slo a travs de la distribucin marginal de cada variable
podemos obtener las frecuencias acumuladas ui, vj es decir, las variables estads-
ticas a las que se refiere la distribucin de frecuencias o cpula emprica. De ah
que la definicin usual de la cpula emprica se haga a partir de la muestra de n
pares (x1, y1), ... , (xn, yn) como la funcin :
emprica respectiva. Para lograr que esta funcin marginal sea continua, hace-
mos la hiptesis de que, entre dos puntos sucesivos cualesquiera observados en
la muestra, xi, xi+1, la densidad es constante.77
Captulo V, Epgrafe 5.
76
Captulo V, Epgrafe 4.
77
141
Segn la anterior definicin, la cpula emprica de una muestra de n pares (x,y)
ser una matriz n x n. Pero puede tambin representarse con una matriz m x m
(m<n), en cuyo caso el elemento (r,s) significar la frecuencia con que se ha pre-
sentado el resultado xa, yb (a y b son los mayores enteros que cumplen la con-
dicin a n.r/m; b n.s/m. As, para los 10 pares (x, y):
obtenemos una cpula emprica que adopta la forma de la siguiente matriz 10 x 10:
Por ejemplo, el elemento (6, 7), igual a 0,6, significa la frecuencia de pares (x,y)
cuyos dos miembros cumplen la condicin de ser x x6, y y y7. En efecto, x6 e y7
142
son los valores que ocupan el lugar 6 y 7 del respectivo estadstico de orden; as,
x6 = 2,898...; y7 = 3,382...78.
Cada columna o fila supone ahora una frecuencia acumulada de 0,25 de manera
que, por ejemplo, el elemento (2, 3), es decir, 0,5 significa que el 50 por 100 de los
pares (x,y) cumplen la condicin de ser no superiores al estadstico de orden n-
mero 5 y 7, respectivamente, es decir, a 2,506 y 3,38279.
Segn dijimos anteriormente, transformar los pares observados (x,y) en sus co-
rrespondientes (ui,vj) requerira conocer las funciones de distribucin marginal;
como ello no es posible en la prctica, es necesario hacer alguna hiptesis sobre
ellas. Cuando adoptamos como tales a las distribuciones empricas, la matriz
definida anteriormente y aquella cuyo elemento (i, j) significa el nmero
de pares (ui, vj) tales que (ui i/n; vjj/n), coinciden.
143
expresara frecuencias. Con las cifras del ejemplo anterior, la matriz sera:
2. AJUSTE DE CPULAS
144
50x50 elementos. Para obtener una idea del coeficiente de dependencia en la
cola superior, tratamos de determinar valores, segn la cpula emprica, de
la frmula (2.11) aunque sin paso al lmite (es decir, valores del cociente
para varios valores de prximos a 1). Designaremos a estos valores por .
En el Cuadro 3.1a se recogen los datos obtenidos con referencia a dicha cola
superior, a partir de = 0,74. En la Figura 3.1 se representan los puntos de
abscisa y ordenada coincidente con la correspondiente cifra de la tercera fila
del Cuadro 3.1a, es decir, con los valores .
0,74 0,76 0,78 0,80 0,82 0,84 0,86 0,88 0,90 0,92 0,94 0,96 0,98
C(,) 0,661 0,68 0,70 0,73 0,76 0,78 0,80 0,83 0,86 0,89 0,92 0,95 0,97
0,70 0,69 0,65 0,64 0,64 0,63 0,60 0,58 0,64 0,65 0,68 0,65 0,55
145
efecto, recordemos que inf = C(,)/(); si a la vista de la matriz de la cpula
emprica 50x50, elaboramos un cuadro similar al 3.1a, pero referido a valores
de cercanos a 0, obtenemos el resultado que resume el Cuadro 3.1b. A la
vista de l, es evidente que para valores de ms cercanos a 0, el valor de
tambin se aproxima a cero. El coeficiente de dependencia en la cola inferior
es nulo en la cpula de Gumbel.
0,06 0,04 0,02
C(, ), 0,013 0,005 0,003
0,014 0,005 0,003
Cuadro 3.1b: Aproximaciones del coeficiente de dependencia en la cola inferior. Cpula
emprica.
En segundo lugar, podemos valernos del mtodo grfico de Genest y Rivest para
seleccionar el modelo de cpula arquimediana ms adecuado para nuestros datos;
este mtodo utiliza las expresiones (2.23) y (2.25) y consiste en comparar la proba-
bilidad K(w) = Prob[C(U,V;)w] segn un modelo concreto de cpula, C(u, v; ), con
la frecuencia del suceso {cpula emprica w}. La idea es la siguiente: una vez cal-
culada la de Kendall a partir de la muestra, la expresin (2.25) permite determinar
el parmetro del modelo de cpula arquimediana que deseemos considerar; ob-
tenemos as una estimacin de dicho parmetro por un mtodo que recuerda el
conocido mtodo de los momentos. A continuacin, sustituimos esta estimacin
en el modelo de cpula y aplicamos (2.23), para obtener la funcin
Captulo 5, Epgrafe 7.
81
146
P(t, ) que se adapte a dichos n puntos (por ejemplo, un polinomio ajustado por
mnimos cuadrados) y aplicar la medida de Frees y Valdez (1998), que viene dada por
(3.1)
Esta integral da una medida de la distancia entre las funciones K(t, ) y P(t, ); en
el Captulo 5, Epgrafe 8, se estudia una aplicacin prctica de la misma, que sirve
de criterio para elegir la cpula ms adecuada.
=sen(./2) (3.2)
Captulo 5, Epgrafe 6.
82
147
(medida con el coeficiente ) al variar el tamao del cuadrado [0, z]2 segn el mo-
delo de cpula y el parmetro considerados. Si hacemos lo mismo con un modelo
distinto, obtendremos otra curva diferente y podremos observar sus diferencias
con la anterior.
Pues bien, si calculamos J(z) para la cpula emprica, resultar una funcin a la
que llamaremos y que compararemos con la de otros modelos para as determi-
nar cul es ms adecuado para nuestros datos.
(3.3).
148
Estos datos nos permiten calcular el parcial para unos cuantos valores de z,
por ejemplo, z=0,1; 0,2; ...; 0,9; 1. El resultado est en el Cuadro 3.2.
z 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1.0
-0,035 0,008 0,125 0,186 0,205 0,240 0,260 0,320 0.414 0,505
Cuadro 3.2: Varios valores del tau parcial.
(el valor muestral de es, por tanto, 0,505).
A continuacin, aplicamos (3.3) para obtener la funcin J(z) para tres modelos
de cpula: Gumbel, Clayton y Frank.
La Figura 3.2, adems de incluir las funciones J(z) de los tres modelos de c-
pula mencionados (todos ellos con parmetro elegido para que sea 0,5, es
decir, =2 en las cpulas de Gumbel y Clayton y 5,74 en la de Frank), se repre-
senta, mediante una lnea quebrada, la funcin (z). Puede verse cmo la c-
pula de Gumbel es, como era de esperar, la que mejor se ajusta a los datos,
sobre todo en valores de z superiores a 0,5. Para valores inferiores el mtodo
no resulta muy fiable por el escaso nmero de pares que entran en el clculo.
La funcin J(z) de la cpula de Clayton es constante e igual a 0,5, con los datos
introducidos.
T. Mikosch (2006) da una visin crtica del empleo de las cpulas y de los procedimientos para
85
149
Figura 3.2: Funcin J(z) para tres modelos de cpula y su comparacin con los valores empricos,
que vienen representados por la lnea quebrada.
Expresin (2.29)
86
150
b) La situacin habitual es la ignorancia de las distribuciones marginales, en
cuyo caso se presentan, a su vez, dos soluciones alternativas:87
87
Vase, por ejemplo, R. Shaw y G. Spivak (2009). A. Staudt, (2010) examina con ejemplos los erro-
res a que puede conducir una mala eleccin de las distribuciones marginales, al determinar la es-
tructura de dependencia.
88
A. Staudt, (2010).
151
ajustar un modelo adecuado a dicho extremo de la distribucin, tal como el de
Pareto o el logartmico-normal,89
Es obvio que estos tres enfoques que hemos examinado sobre el ajuste de cpulas
no son excluyentes entre s, sino complementarios.
Para comparar el ajuste de varios modelos a una base de datos suele utilizarse el
criterio de Akaike (Akaike Information Criterion, AIC), que consiste en calcular
y elegir como mejor el modelo al que corresponda un menor valor de AIC. Cuando
el nmero de parmetros estimados en los distintos modelos que se comparan es
el mismo, el criterio consiste, sencillamente, en elegir aqul que tenga una vero-
similitud ms alta.
3. SIMULACIN DE CPULAS
Es muy til simular cpulas para resolver ciertos problemas prcticos que pre-
sentan dificultades de clculo. La realizacin de un gran nmero de simulaciones
permite aproximarse a la probabilidad del fenmeno simulado.
152
sabemos que sta vincula a la distribucin conjunta con las marginales, de mane-
ra que un punto simulado de la cpula, (u1, u2, ..., ud), permite obtener una realiza-
cin de la variable conjunta (x1, x2, ..., xd) al aplicar la transformacin xi = Fi1(ui)
(siendo Fi la distribucin marginal de Xi).
En este Epgrafe examinaremos cmo simular cpulas, cuestin que puede divi-
dirse en otras dos, segn que la cpula que deseemos simular sea implcita o
explcita. En el primer caso, la forma de proceder es la siguiente:
1. Simular la variable (X1, X2,...,Xd) distribuida segn la ley cuya cpula implcita
se trata de simular (por ejemplo, la variable t multidimensional). Obtenemos
as d valores (x1, x2,...,xd).
ste es el caso de las cpulas gaussiana y t, que veremos con ms detalle a con-
tinuacin.
El primer paso es simular el vector Y = (Y1, Y2, ..., Yd) en donde las Yi son normales
(0,1) con matriz de covarianzas 90. En virtud de (1.2) podemos escribir:
Y=C.X
Al ser las variables Yi normales (0,1), la matriz de covarianzas coincide con la matriz de correlaciones.
90
Captulo I, Epgrafe 3.
91
153
Damos por resuelta la cuestin de cmo simular variables normales independien-
tes entre s;92 por tanto, si conociramos C obtendramos Y aplicando la anterior
expresin, pero no nos es dada dicha matriz, sino la de covarianzas, . Por (1.3)
sabemos que C.C = y el problema de obtener el vector Y puede resolverse con
la factorizacin de Cholesky. A partir de una matriz simtrica y definida positiva,
dicha factorizacin consiste en encontrar una matriz triangular inferior T tal que
T.T = . En efecto, si conocemos T, entonces Y = T.X, ya que
Observemos que T. T = .
Y2 = 0,8 X2 + 0,6 X3
Y3 = X3
154
Si, por ejemplo, obtuviramos 5 simulaciones del vector X (normales 0,1 e inde-
pendientes entre s) como las siguientes (cada fila es una simulacin de X)
En resumen, el algoritmo para simular un vector Y normal (0,1) con matriz de co-
varianzas y su correspondiente cpula, es el siguiente:
3. Calcular Y = T.X.
155
4. Aplicar la transformacin u =(y), es decir, determinar la probabilidad acumu-
lada que corresponde a cada valor simulado y segn el modelo normal (0,1).
Cpula t
156
simulacin de la cpula gaussiana. El resultado es la siguiente matriz de 3
columnas y 5 filas (cada fila es una simulacin del vector)
1. Obtener una simulacin del vector Y normal (0,1) con matriz de covarianzas o
correlaciones segn el proceso descrito para la cpula gaussiana.
5. Repetir el proceso.
157
Finalmente, simulamos 1000 puntos de la cpula t con parmetros = 0,7
y = 6.
En las proximidades de los vrtices (0, 1) y (1, 0) aparecen ms puntos que en los
mismos lugares de la grfica para la cpula gaussiana. Esto se debe a que la distri-
bucin t tiene unas colas ms pesadas que la distribucin normal. Si la correlacin
es moderada, pueden coincidir puntos de la cola derecha de X (lo que corresponde
a valores altos de U) con puntos de la cola izquierda de Y (lo que corresponde a va-
lores de V prximos a cero).
Procedimiento general
El mtodo para simular cpulas que esbozaremos ahora se basa en la idea de que
puede obtenerse un valor simulado de la variable X a partir de un nmero pseudo-
aleatorio uniforme (0,1), al que designamos Ri, mediante la inversin de la funcin
de distribucin u = F(x) de X, es decir, resolviendo la ecuacin xi = F-1(Ri); ello suele
exigir el recurso a mtodos numricos.
En el caso de las cpulas bivariantes C(u,v), aplicar dicho mtodo exige llegar a
una funcin de distribucin de una sola variable, lo que puede lograrse con la
158
funcin de distribucin condicionada C(uv); en las cpulas multivariantes se re-
curre a las distribuciones marginales
En el caso de una cpula multivariante C(u1, u2, ... , ud) ha de seguirse un mtodo
recurrente94 en el que intervienen las funciones de distribucin marginal de Ci(u1, ...,
ui, 1, ... , 1) (para i < d) y las funciones de distribucin condicionada, buscando siem-
pre en este proceso que aparezca una ecuacin con una sola incgnita, resultante
de igualar una de dichas funciones a un nmero pseudo-aleatorio previamente ge-
nerado y sustituir uno de los valores condicionantes por otro de estos nmeros.
En la cpula Frank puede despejarse u algebraicamente y por tanto no hay necesidad de aplicar
93
mtodos numricos:
Embrechts, McNeil y Straumann, (2002). Por otra parte, T. Schmidt, (2006) describe un mtodo
94
de simulacin aplicable a las cpulas arquimedianas que, como las examinadas en el presente tra-
bajo, pertenecen a la clase llamada cpulas arquimedianas de transformada de Laplace. En estas
cpulas, la inversa de la funcin generador puede representarse como la transformada de Laplace
de una cierta funcin, lo que puede aprovecharse para simular valores.
159
El proceso es el siguiente:
4. Calcular la raz u2 de dicha ecuacin. Resultan as los dos primeros valores si-
mulados del vector (U1, U2, ..., Ui), es decir, R1 y la raz u2, a los que llamaremos
1 y 2 respectivamente.
160
Al ser tres las variables de que depende la funcin C, necesitaremos dos fun-
ciones de distribucin marginal y sus correspondientes funciones de distribu-
cin condicionada.
C2(u,v)/v=
[2C(u, v, w)/uv]/[2C2(u,v)/(uv)] =
Hagamos, por ejemplo, = 10. Obtenidas las funciones que vamos a necesitar, el
primer paso es generar dos nmeros aleatorios95, que resultan R1= 0,3398 y R2
= 0,5127. A continuacin, en C(vu) sustituimos u por R1 e igualamos la funcin a
R2. Por procedimientos numricos determinamos la raz v, que es v= 0,3422.
161
Resolvemos la anterior ecuacin en w con el resultado w = 0,3130; por tanto,
la primera simulacin da (0,3398; 0,3422; 0,3130). El Cuadro 3.3 recoge las
cifras de las cinco simulaciones practicadas.
Simulacin n R1 R2 R3
1 0,3398 0,5127 0,2311 0,3398 0,3422 0,3130
2 0,9311 0,5745 0,9085 0,9311 0,9322 0,9452
3 0,8353 0,7978 0,2206 0,8353 0,8535 0,8288
4 0,6458 0,7176 0,4321 0,6458 0,6606 0,6501
5 0,1992 0,3890 0,6288 0,1992 0,1881 0,2153
Cuadro 3.3: Cinco simulaciones de la cpula de Gumbel trivariante (=10).
Cpula de Marshall-Olkin
El algoritmo para simularla sigue los pasos que, para llegar a ella, describimos
en el Captulo II, Epgrafe 9, es decir:
1 = 1; 2 = (1)/((1)); 12 = /(1).
162
IV. OTRAS CUESTIONES SOBRE CPULAS. EJEMPLOS
Llamamos cpula mixta C(u,v) a la que combina a otras dos: C1(u,v; 1) y C2(u,v; 2)
de la siguiente forma:
Como punto de partida para el ejemplo que vamos a estudiar, tomaremos 500
pares de nmeros (x,y) representativos de siniestros o prdidas observados en
relacin con dos riesgos, que denominaremos riesgo 1 y riesgo 297. Estos
Ver II,2.
96
163
datos quedan representados en la Figura 4.1 y algunas de sus caractersticas
se reflejan en el Cuadro 4.1.
Para centrarnos en la cuestin que ahora nos interesa, admitiremos que las
distribuciones marginales son del tipo Pareto, con parmetros 1 = 5,2; 1
=3,6; 2 = 2,6; 2 = 3,6. Ello determina que los valores medios son 2 y 1, respec-
tivamente para el riesgo 1 y 2, mientras que las desviaciones tpicas son 3 y
1,5. Llamamos F1 y F2 a dichas distribuciones marginales98 y aplicando las
transformaciones u=F1(x); v=F2(y) obtenemos los puntos (ui,vi) representados
en la Figura 4.2 y que corresponden a la cpula emprica. Nos proponemos
obtener un modelo que se ajuste a los datos que suponemos observados y
representados en dicha Figura 4.1.
164
Figura 4.1: 500 puntos representativos, cada uno de ellos, de las prdidas o siniestralidades de
dos riesgos, Riesgo 1 y Riesgo 2.
165
Log. mx.
Modelo en funcin de mximo-verosmil
verosimilitud
Frank 4,70 4,70 115,754
Gumbel 1,779 1,730 128,273
Cuadro 4.2: Parmetros estimados de dos modelos de cpula para los datos representados en
la Figura 4.2.
(a) (b)
Figura 4.3: Funcin K(w) para dos modelos de cpula y para los datos observados.
166
Tenemos, por tanto, tres parmetros, que procedemos a estimar por mxima
verosimilitud99, con el resultado que refleja el Cuadro 4.3.
167
observadas (para w = 0; 0,1; , ... , 0,9; 1) con los valores de K(w) de la cpula
mixta (obtenidos mediante simulacin) para dichas abscisas. Se observa en
4.4b la cercana de los puntos qq a la lnea recta.
(a) (b)
Figura 4.4: Funcin K(w) y grfica cuantil-cuantil para la cpula mixta.
Para confirmar la preferencia de la cpula mixta sobre los otros dos modelos
contemplados, aplicamos a las diferentes cpulas la medida (3.1) propuesta
por Frees y Valdez y calculamos tambin la de Akaike, con el siguiente resul-
tado:
168
2. REGRESIN MEDIANA
Si F(x,y) = Prob(X < x;Yy) la regresin de Y sobre X para el cuantil p ser el valor y
que satisfaga la ecuacin
F(yx) = p
169
y, en definitiva,
170
lo que nos da v 0,6744. Para este valor de v, la funcin de distribucin mar-
ginal de Y hace corresponder el importe y = 2,244 mientras que para u = 0,7 el
valor de x es x=1,20.
Figura 4.5: Curvas de regresin de Y sobre X para los cuantiles 0,50 y 0,90.
171
cobertura de cierto tipo de daos, pero as como el Riesgo 1 es la siniestrali-
dad del seguro directo, la del Riesgo 2 es la del reaseguro aceptado, cuya in-
formacin le llega con algn retraso. El asegurador tiene inters, por tanto,
en hacer predicciones sobre el coste que para l supondr el reaseguro acep-
tado, conociendo el incurrido en el seguro directo.
En la Figura 4.6 estn las curvas de regresin del Riesgo 2 (siniestralidad del
reaseguro aceptado) sobre el Riesgo 1 (del seguro directo) para los cuantiles
0,5, 0,95 y 0,99. Si, por ejemplo, en el seguro directo la siniestralidad ha sido
de 30 unidades, su importe en el reaseguro aceptado puede preverse que ser
de 10,5, aproximadamente, sin que supere las 18,7 unidades con una probabi-
lidad del 95%, ni las 22,6 con probabilidad del 99%.
Figura 4.6: Curvas de regresin de x=Riesgo 2 sobre y= Riesgo 1 para los cuantiles 0,50, 0,95 y
0,99.
Captulo II, 4.
100
172
consistira en aplicar reiteradamente la expresin (2.14) para ir incrementando el
nmero de dimensiones, es decir:
........
Cn = C(u1, u2, . . . , un) = -1[(Cn-1) + (un)] = -1[(u1) + (u2) + ... + (un)] (4.3).
Pero este procedimiento puede desembocar en funciones que no renen los re-
quisitos constitutivos de las cpulas. Un ejemplo que propone Nelsen (2006)101 es
el siguiente:
d = W(u2, v2, z2) W(u2, v2, z1) W(u2, v1, z2) W(u1, v2, z2) + W(u2, v1, z1) + W(u1, v2, z1)
+ W(u1, v1, z2) W(u1, v1, z1) (para u1<u2; v1<v2; z1<z2).
Si, por ejemplo (u2, v2, z2) = (1, 1, 1) y (u1, v1, z1) = (1/2, 1/2, 1/2), d = 0,5. La pregun-
ta que surge es, por tanto, qu condiciones tiene que reunir la funcin (t) para
que (4.3) sea una cpula103; la respuesta tiene que ver con el concepto de funcin
completamente montona, que es aqulla que, adems de ser continua, cumple la
condicin
Pgina 151.
101
Captulo II, 2.
102
173
Pues bien, hay un teorema, debido a Kimberling,104 que viene a afirmar que la fun-
cin Cn de (4.3) ser una cpula cuando (t) es un generador estricto y su inversa
-1(t) es completamente montona para t0.
C(u1, u2, ... , un) = exp{ [ [LN(u1)] + [LN(u2)] + ... +[LN(un)] ]1/}
es una cpula. Pero esta funcin C tiene el inconveniente de que sus n varia-
bles son intercambiables entre s en el sentido de que permutar el lugar que
ocupan en la funcin resulta irrelevante, al ser nico el generador empleado y,
con ello, tambin nico el parmetro ; como consecuencia, el grado de de-
pendencia entre cualquier par de variables (ui, uj) (i j) ser el mismo, lo que
limita la utilidad de esta cpula multivariante105. Ello nos lleva a plantear si es
posible construir una cpula arquimediana multivariante con distintos genera-
dores. Una posibilidad consiste en seguir el proceso siguiente106, en donde h
j significa h[j(t)]:
104
Nelsen, (2006), pgina 152.
La cpula arquimediana trivariante simulada en III,3 tena esta caracterstica, es decir, un solo
105
parmetro.
106
Seguimos aqu el trabajo citado en la anterior Nota 103.
107
Para llegar a esta expresin hemos hecho C3=2-1[23-1(3(u) + 3(v))+ 2(w)]
174
y as sucesivamente. Pero surge nuevamente la pregunta de si, en esta nueva si-
tuacin (es decir, diversidad de generadores), este mtodo para construir cpulas
arquimedianas con distintos generadores conduce siempre a una funcin Cn que
es una cpula o, por el contrario, no es suficiente la condicin de que las funciones
i(t) sean generadores estrictos y sus inversas completamente montonas. La
respuesta108 es que se necesitan ahora ms requisitos; as, cuando se trata de
la cpula de Gumbel, ha de ser 1 2 ... n109. Observemos ahora lo siguiente:
Si, a partir de (4.4a) obtenemos la cpula marginal C(u,v), resulta, teniendo en
cuenta que, por definicin de funcin generador, ha de ser 1(w) = 0, para w = 1,
Si no se cumple este requisito puede resultar una funcin que no es n-creciente. Por ejemplo,
109
en la cpula de Gumbel trivariante con 3-1=2 generadores distintos, cuya expresin es:
si hacemos 1=10 y 2 =1,1 (contraviniendo as el requisito de ser 1 2) tenemos que, para (u1, v1, w1)
= (0,5; 0,04; 0,05) y (u2, v2, w2) = (0,95; 0,95; 0,1) el volumen segn C(u, v, w) del recinto definido por
los anteriores puntos es negativo:
C(u2, v2, w2) C(u2, v2, w1) C(u2, v1, w2) C(u1, v2, w2) + C(u2, v1, w1) + C(u1, v2, w1) + C(u1, v1, w2) C(u1,
v1, w1) = 0,0032.
175
Llegamos, por tanto, al siguiente esquema para representar las relaciones de
dependencia entre U, V y W:
Si aplicamos esto mismo a (4.4b) llegamos a que en C(u, v, w, z) puede haber tres
relaciones de dependencia diferentes, de las cuatro existentes en total.
A) Dos riesgos
176
casos su media es 1, con una varianza de 2,778. Ello determina las siguientes fun-
ciones de distribucin de X e Y:
siendo f(y) la funcin de densidad del ratio Y. Por otra parte, la probabilidad con-
dicionada que figura en el integrando equivale a
y la funcin de densidad,
177
Con los tres pasos anteriores obtenemos la funcin de densidad f (x,y).
4. Sustituir x por sy, de forma que dicha funcin pasar a ser f (sy, y).
El VaR al 99,5% segn la funcin f (s) es 18,2, inferior, por tanto, al valor segn
la aproximacin lognormal. Si en lugar de esta aproximacin adoptramos la
normal (como sera propio si atendemos a la utilizacin de correlaciones como
medida de la dependencia) el VaR al 99,5%, con los anteriores datos sera 10,2.
178
Figura 4.7. Funcin de densidad de la suma de los ratios de siniestralidad X e Y.
B) Tres riesgos
Primer ejemplo
a. El seguro de enfermedad que se identifica por las siglas SLT (Similar to Life
insurance Technique, es decir, seguro de enfermedad planteado actuarialmen-
te segn tcnica similar a la del seguro de vida).
179
El capital obligatorio para el mdulo seguro de enfermedad, ha de calcularse
mediante la frmula siguiente, que lleva implcita la hiptesis de normalidad en la
distribucin de los riesgos de suscripcin de estos tres sub-mdulos111:
en donde SCR(.) son los capitales correspondientes a cada uno de los tres sub-
mdulos, calculados segn las reglas que se detallan en la normativa. Los coefi-
cientes de correlacin CorrH(i,j) quedan fijados por dichas normas y son los si-
guientes:
Entre los sub-mdulos SLT y NSLT, por una parte, y el submdulo de riesgo de
catstrofe, por otra, 13 = 23 = 0,25.
Pues bien, supongamos que el capital SCR o VaR al 99,5% para cada uno de los
tres sub-mdulos fuese igual a 10 unidades monetarias y nos preguntamos cul
sera el VaR o capital exigido para el mdulo riesgo de suscripcin del seguro de
salud en caso de que la estructura de dependencia entre los tres riesgos que lo
integran se apartara del modelo gaussiano. Concretamente, calcularemos dicho
capital en los dos casos siguientes de estructura de dependencia entre las tres
variables representativas de dichos riesgos: a) Cpula Gumbel con distintos valo-
res del parmetro de dependencia, y b) Cpula t, con varios grados de libertad.
Los mrgenes son, en todo caso, normales, con media cero y desviacin tpica
igual a 10/2,58 3,876.112 Es decir que referimos ambas distribuciones marginales
al ratio de prdida experimentada en el ejercicio econmico sobre el volumen de
negocio.
gaussiana con mrgenes normales. Que bajo el signo radical aparezcan capitales en lugar de des-
viaciones tpicas implica admitir que las medias marginales son nulas; el VaR al 99,5% coincide con
los capitales de solvencia exigidos.
112
Segn veremos en el Captulo V,1, la normativa de Solvencia II adopta como probabilidad mni-
ma de solvencia el 99,5%. En la distribucin normal (0,1), el cuantil correspondiente es, aproxima-
damente, 2,58. Por tanto, si el SCR calculado para cada sub-mdulo es 10, entonces 2,58 x desvia-
cin tpica de ese sub-mdulo = 10.
180
a) Cpula Gumbel
Calculamos el VaR o SCR del mdulo por medio de simulaciones (50.000 en cada
caso) para dos grupos de valores de los parmetros de la cpula trivariante (1,2),
elegidos de forma que las correlaciones 12, 13 y 23 sean, o bien similares a las
que establece la normativa (es decir, 0,5, 0,25 y 0,25 respectivamente) o bien ms
altas (0,7, 0,5 y 0,5).
Por tanto, si las correlaciones (segunda fila, columna 2 del Cuadro 4.6) son simi-
lares a las que prescriben las normas de Solvencia II, la aplicacin de stas infra-
valora el capital necesario para que el mdulo alcance una solvencia del 99,5%,
puesto que dicho capital necesario 25,15 representa el 112,5% del requerido
por la normativa, que es de 22,36. En el supuesto de que las correlaciones entre
los tres sub-mdulos fueran ms altas, la infravaloracin sera an ms impor-
tante (el 123,7 por 100).
181
b) Cpula t
Segundo ejemplo
182
Los parmetros de estas funciones estn en el Cuadro 4.8. Los coeficientes de
Kendall para cada par de riesgos figuran en el Cuadro 4.9.
Riesgos X Y Z
X - 0,5 0,3
Y 0,5 - 0,3
Z 0,3 0,3 -
Cuadro 4.9: Coeficientes entre X, Y, Z.
Cpula gaussiana
Cpula de Gumbel
Por tanto, en este ejemplo no nos referimos a la variable suma de los ratios de siniestralidad
113
183
En el caso de las cpulas gaussiana y t trivariantes determinaremos los parme-
tros XY, XZ y YZ mediante la expresin (3.2), que relaciona el coeficiente con los
coeficientes de correlacin en la cpula gaussiana; y as, tenemos que, con los
datos del Cuadro 4.9, resulta:
Por medio de simulaciones de las variables X, Y, Z114 obtenemos otros tantos va-
lores de la suma y computamos las veces en que sta no supera el importe desea-
do. El resultado de las 50.000 simulaciones realizadas para cada caso se recoge
en el Cuadro 4.10.
Se observa que, para los datos de este ejemplo, son muy similares los resultados
de la cpula gaussiana y los de la cpula t con 10 grados de libertad, mientras que
la de 3 grados de libertad se acerca a la cpula de Gumbel. Tambin queda mani-
fiesta la gran diferencia de todos estos resultados con los de independencia.
Ver III,3.
114
184
VaR segn modelo de cpula
t con 10 g. t con 3
Gaussiana Gumbel Independencia
de l. de g. de l.
99% 71,5 75 71,5 75 57
99,5% 88,7 95 89 96,5 71,2
99,9% 145 170 150 170 123
Cuadro 4.10: Importes del VaR para distintos niveles de probabilidad y estructuras de
dependencia entre los tres riesgos considerados (clculo por simulacin).
Por ltimo, compararemos el importe del VaR al 99,5 por 100 calculado ante-
riormente para cada hiptesis, con el VaR obtenido segn las normas de
Solvencia II. Aqu hemos de abrir un breve parntesis para aclarar que en tales
normas subyace la hiptesis de que las variables que representan los riesgos
de suscripcin del seguro novida se distribuyen segn la ley logartmico-nor-
mal; el procedimiento para calcular el capital de solvencia obligatorio o SCR en
este seguro es: Primero, determinar la desviacin tpica del ratio de siniestra-
lidad y gastos sobre el volumen de negocio, para el conjunto de lneas de segu-
ro mediante una frmula como la (5.4), es decir, la raz cuadrada de la suma de
los productos de las desviaciones tpicas por las correlaciones. A continuacin,
se ajusta una distribucin log-normal cuya media sea la unidad y desviacin
tpica la que acabamos de obtener. En tercer lugar, se determina el VaR al
99,5% segn esta funcin de distribucin y se le resta la unidad (la media), para
as llegar al SCR por unidad de volumen de negocio.115. Ya dijimos que en este
ejemplo las variables X, Y y Z no significan ratios como ocurre en el procedi-
miento de Solvencia II que acabamos de resumir sino cantidades de sinies-
tralidad medidas en unidades monetarias; por tanto, hemos de abandonar nece-
sariamente la hiptesis de que la media es igual a la unidad. El procedimiento
que vamos a seguir para llegar al VaR segn los criterios de Solvencia II adap-
tados a nuestro ejemplo, ser:
En el Captulo V, Epgrafe 1 hay una breve descripcin de las normas de Solvencia II ms rele-
115
vantes a nuestros efectos. Aqu nos interesan especialmente las expresiones (5.2), (5.3) y (5.4). En el
Epgrafe siguiente trataremos de la aproximacin logartmico-normal preconizada por Solvencia II
como modelo del ratio de siniestralidad global no-vida.
185
Probabilidad, para cada modelo de cpula, de superar el VaR
calculado para la hiptesis de independencia
VaR Gaussiana Gumbel t con 10 g. de l. con 3 de de l. Independencia
99% 0,021 0,022 0,018 0,02 0,01
99,5% 0,009 0,012 0,01 0,012 0,005
99,9% 0,002 0,002 0,002 0,002 0,001
Cuadro 4.11.
186
Hiptesis sobre el
t - 10 g de
modelo de Gaussiana Gumbel t - 3 g de lib. indep
lib.
dependencia
0,59/0,259/ 0,688/0,49/ 0,597/0,36/ 0,598/0,38/
XY / XZ / YZ; 0
0,22 0,53 0,36 0,51
VaR 99,5 % (segn
el modelo log- 90,30 101,9 95,23 98,60 72,90
normal)
Porcentaje del VaR
al 99,5% segn el
Cuadro 4.10, en 98,2 93,2 93.5 97,9 97,7
relacin con el de
la fila anterior.
Cuadro 4.12: Comparacin del VaR segn varias estructuras de dependencia entre los riesgos X, Y,
Z con el determinado bajo la aproximacin lognormal.
187
Para obtener las cifras de la tercera fila del Cuadro 4.13 hemos introducido en un
modelo normal la desviacin tpica resultante de combinar las desviaciones mar-
ginales con las correlaciones, as como la media igual a 18 y calcular el VaR al
99,5%. Por ejemplo, en el caso del modelo gaussiano (segunda columna) resulta
que los coeficientes de correlacin que para l contiene el Cuadro 4.13 y las des-
viaciones tpicas de las distribuciones marginales (recordemos que son 5, 5 y 10
respectivamente) conducen a una desviacin tpica global de 15,08 unidades,
aproximadamente. El cuantil al 99,5% de la distribucin normal (0,1) es, tambin
aproximadamente, 2,58. Por tanto, en la distribucin del riesgo global ese mismo
cuantil o VaR cuya media es 10+3+5=18 ser 15,08 x 2,58 + 18 56,91. El im-
porte de 88,7 que figura en el Cuadro 10 es el 155,86% de 56,91.
Segundo, clculo del ratio de la siniestralidad conjunta (es decir, de todas las
lneas de seguro novida), mediante la conocida frmula en que intervienen
188
los coeficientes de correlacin entre cada par de riesgos o lneas de seguro,
coeficientes que tambin aportan dichas normas.
Para facilitar las cosas supondremos, en los ejemplos que examinaremos inme-
diatamente, que la cartera novida cuyo riesgo de primas trataremos de medir, se
compone de dos lneas de seguro solamente, X e Y. Tambin admitiremos, para
mayor sencillez, que las distribuciones marginales de ambos riesgos son las mis-
mas. En tercer lugar, en todos los ejemplos que vamos a estudiar trataremos de
determinar el VaR al 99,5% de la variable suma de los dos riesgos X+Y. Dichos
riesgos, X e Y, consistirn en el ratio de la siniestralidad ms gastos sobre las
primas116 a excepcin del primero de estos ejemplos, en que las variables X e Y
representan la siniestralidad ms gastos medidos en unidades monetarias y no
como proporcin del volumen de negocio.
116
Prescindimos de la casustica contemplada en la normativa de Solvencia II acerca de cmo de-
terminar el volumen de primas que entra en el clculo.
189
Sobre esta base, y una vez elegidos los parmetros (coeficiente de correlacin)
y (nmero de grados de libertad) de la funcin de densidad conjunta t bivariante,
la expresin (4.5a) nos permitir calcular probabilidades Prob(X+Y<s) y, en defini-
tiva, llegar al importe de s que se identifique con el VaR al 99.5%. As, la integral
= 0,5
= 10.
Veamos ahora cul sera el VaR resultante para esta misma cartera si segui-
mos las hiptesis subyacentes en las normas de Solvencia II. Para su clculo
hemos de determinar los parmetros de una distribucin log-normal de me-
dia 40 y desviacin tpica igual a la de X+Y. La varianza de X (y de Y) es /(2)
por tratarse de una t de Student. Por tanto, la desviacin tpica de la suma,
teniendo en cuenta que la correlacin elegida es 0,5, ser:
El VaR segn los criterios implcitos en las reglas de Solvencia II, se determi-
nar mediante la integral de la densidad logartmico-normal (, ), es decir:
190
cuyos parmetros han de ser tales que la media sea 40 y la desviacin estn-
dar 1,94, o sea:
pues bien, la anterior integral vale muy aproximadamente 0,995 y, por tanto,
el VaR al 99,5% segn este segundo procedimiento es 45,27.
El Cuadro 4.14 resume los resultados de clculos como los anteriores, para
distintos valores del parmetro . El coeficiente de correlacin se ha fijado
en 0,5.
(3) (4)
(2) (5)
(1) Desviacin VaR 99,5%
VaR 99,5% 100 x (2)
(gr de lib) estndar Aprox
distrib exacta /(4)
cartera lognormal
3 50,30 3,00 48,38 104,0
6 46,42 2,12 45,78 101,4
10 45,49 1,94 45,27 100,5
20 44,93 1,82 44,93 100,0
30 44,765 1,79 44,85 99,8
Cuadro 4.14: Comparacin del VaR al 99,5% segn el modelo t bivariante y segn la
aproximacin lognormal (Ejemplo 1).
191
Figura 4.8: Porcentaje que representa el VaR al 99,5% sobre el calculado utilizando la
aproximacin lognormal.
(4.6)
192
en donde s representa la desviacin tpica elegida para las variables X e Y, el
coeficiente de correlacin y el nmero de grados de libertad. El procedi-
miento para obtener los resultados que se contienen en los Cuadros 4.15 y
4.16 es similar al del Ejemplo 1: Valindonos de la expresin (4.5a), en la que
introducimos la funcin f(x,y) anterior, calculamos el VaR al 99,5%, que com-
paramos con el que corresponde segn la aproximacin log-normal. sta ha
de tener media 2 (puesto que es la suma de X+Y) y la desviacin tpica desea-
da, que es 0,1 o bien 0,2, segn se ve en el Cuadro respectivo.
De manera anloga, el Cuadro 4.16 refleja los mismos datos, salvo que = 6.
Llama la atencin la gran diferencia que presentan las cifras de la columna (6)
cuando nos situamos en el caso de 3 grados de libertad y desviacin tpica 0,2,
con respecto a los dems supuestos contemplados. En la Figura 4.9117 se ve
con mayor claridad esta diferencia, as como la evolucin de la proporcin
VaR segn cpula/VaR segn aproximacin log-normal.
La escala de las ordenadas es, en esta Figura, muy distinta de la escala de las Figuras 4.8 y 4.10.
117
193
(1) (2) (3) (4) (5) (6)
VaR 99,5% s Desv VaR 99,5% 100 x (2)
distrib. (desv. estndar aprox. /(5)
exacta estndar cartera lognormal
de X y de Y)
0 2,429 0,1 0,1414 2,392 101,5
0,25 2,480 0,1 0,1581 2,444 101,5
0,50 2,525 0,1 0,1732 2,490 101,4
0,75 2,567 0,1 0,1871 2,533 101,3
0 2,893 0,2 0,2828 2,845 101,7
0,25 3,000 0,2 0,3162 2,961 101,3
0,50 3,090 0,2 0,3464 3,069 100,7
0,75 3,173 0,2 0,3742 3,170 100,
Cuadro 4.16: Comparacin del VaR al 99,5% segn el modelo t bivariante modificado con =6
y segn la aproximacin lognormal (Ejemplo 2).
En todos los casos contemplados en este ejemplo resulta superior el VaR de-
terminado mediante la cpula, al calculado bajo los criterios de Solvencia II.
Figura 4.9: Porcentaje que representa el VaR al 99,5% segn la distribucin t, sobre el calculado
utilizando la aproximacin lognormal
194
(1) (2) (3) (4) (5)
VaR 99,5% Desv estndar VaR 99,5% 100 x (2)
distrib. exacta cartera aprox. /(4)
lognormal
Como era de esperar el VaR segn la distribucin exacta es, en todos los ca-
sos, inferior a la aproximacin logartmica. La asimetra positiva propia de
esta distribucin tiende a exagerar el importe del VaR si nos situamos en el
caso de variables simtricas, cuya cpula carece de dependencia en colas. La
Figura 4.10 refleja estos resultados.
Figura 4.10: Porcentaje que representa el VaR al 99,5% segn la distribucin normal, sobre el
calculado utilizando la aproximacin lognormal.
195
Ejemplo 4: Cpula de Gumbel con mrgenes lognormales
196
(1) (2) (3) (4) (5) (6)
VaR 99,5% Desv VaR 99,5% 100 x (2)
distrib. estndar aprox. /(5)
exacta Cartera lognormal
1 2,392 0,00 0,1414 2,392 100,0
1,1 2,459 0,16 0,152 2,424 101,4
1,5 2,538 0,55 0,176 2,498 101,6
2 2,557 0,70 0,184 2,523 101,3
4 2,570 0,92 0,196 2,560 100,4
6 2,571 0,97 0,198 2,567 100,2
Cuadro 4.18: Comparacin del VaR al 99,5% segn el modelo Gumbel-lognormales y segn la
aproximacin lognormal (Ejemplo 4). Desviacin tpica marginal = 0,1.
197
Figura 4.11: Porcentaje que representa el VaR al 99,5% segn la distribucin gumbel-lognormal,
sobre el calculado utilizando la aproximacin lognormal.
198
CAPTULO V. APLICACIN DE LAS CPULAS AL CLCULO
DEL CAPITAL DE SOLVENCIA EXIGIDO A LAS EMPRESAS
ASEGURADORAS POR SOLVENCIA II. ESTUDIO DE UN
EJEMPLO NUMRICO
199
tericos que presenta118 y se adopt en su lugar el expected shortfall o tail conditio-
nal expectation119, que es:
ES(X) = E {XX>VaR(X)}
Por tanto, adoptar como criterio de solvencia el VaR al 99,5% resulta en la prcti-
ca aproximadamente similar al expected shortfall al 99%. Por otra parte, bajo la
hiptesis que adopta Solvencia II al configurar el modelo estndar para el capital
de solvencia bsico, el VaR es una medida del riesgo que tiene carcter subaditivo,
porque entre tales hiptesis est la de normalidad de los riesgos agregados121.
Estos trminos podran traducirse por dficit esperado o valor medio de la insuficiencia, el
119
A. Ferri, (2012).
121
200
2) Otra caracterstica de Solvencia II es que adopta el ao como amplitud de su
horizonte temporal.
3) Una tercera nota es el carcter finalista de Solvencia II; quiere ello decir que en
el sistema se marcan unos objetivos que pueden lograrse, bien con el modelo es-
tndar, que se describe con detalle en la propia normativa, o bien con los modelos
internos de carcter total o parcial que decida instaurar cada entidad asegurado-
ra, en su caso.
Esta breve descripcin del sistema Solvencia II se centrar slo en el Pilar 1, por
tener la mxima relevancia de los tres a nuestros efectos. Dentro de l, nos ceire-
mos a los aspectos que ms nos interesan sobre el capital de solvencia obligatorio
(en lo sucesivo, SCR, por su denominacin en ingls Solvency Capital Requirement) y
prescindiremos de las cuestiones que no sean necesarias para nuestro propsito,
que es comparar el importe del SCR a que conducen dichas normas con el que re-
sultara de la aplicacin de las cpulas. Diremos, no obstante, que otro concepto
importante es el de Minimum Capital Requirement o Capital Mnimo Obligatorio: se
trata del lmite por debajo del cual los asegurados y beneficiarios de las plizas de
seguro estaran expuestos a un nivel inaceptable de riesgo en caso de que el asegu-
rador o reasegurador continuara con su actividad.
201
b. El Capital Obligatorio para el Riesgo Operacional (este riesgo se refiere a la
posible prdida por procedimientos, personal o sistemas inadecuados o falli-
dos, o por acontecimientos externos; se incluyen los riesgos legales, pero no
los derivados de decisiones estratgicas ni los de reputacin).
d. Riesgo de mercado.
A cada uno de estos riesgos se asigna un capital obligatorio en virtud de las nor-
mas de clculo que se dictan. Y estos capitales se agregan, por va de los coefi-
cientes de correlacin que tambin se establecen, para llegar a dicho capital
122
El artculo 87 del Reglamento Delegado 2015/35 aade un capital de solvencia obligatorio fren-
te a los riesgos que afecten a los activos intangibles.
202
bsico, al que ha de sumarse el capital obligatorio para el riesgo operacional y el
ajuste (de signo negativo) por la capacidad de absorcin de prdidas de las provi-
siones tcnicas y de los impuestos diferidos.
b. Riesgo de catstrofes.
Es decir, que se suman todas las combinaciones posibles, (i,j) de los tres sub-
mdulos anteriores (primareserva, catstrofe y cada de cartera). En esta frmu-
la subyace la hiptesis de que las variables correspondientes a los riesgos con-
templados en ella se distribuyen normalmente, con media cero y desviacin tpica
igual al cociente entre el respectivo SCR y el cuantil al 99,5 por 100 segn la dis-
tribucin normal (que es, aproximadamente, 2,58).
Corr NL(i,j) = coeficiente de correlacin entre los riesgos i,j de los seguros no-
vida (NL= non-life).
203
las obligaciones de estos seguros, atendiendo a los eventos cubiertos y a los procesos
seguidos en el ejercicio de la actividad; este mdulo tendr en cuenta la incertidum-
bre de los resultados de las empresas de seguros y de reaseguros en relacin con las
obligaciones de seguros y reaseguros ya existentes y las nuevas actividades que se
espere realizar en los prximos doce meses. Es decir, que los riesgos propios de este
mdulo incluyen todas las desviaciones desfavorables de las obligaciones financiadas
con las primas y las provisiones tcnicas. Por lo que se refiere al riesgo de prima ello
incluye todas las obligaciones que sta ha de cubrir: No slo la siniestralidad, sino
tambin los gastos de administracin.
204
El importe mnimo de c, ser:
y por tanto:
El valor de N es, aproximadamente, 2,58 y resulta as una funcin c=f () que difie-
re muy poco123 de la recta c=3. De ah que se acepte como norma la expresin
(5.2), en la que subyacen, por tanto, dos hiptesis:
205
Admitir que la tasa de siniestralidad o factor de riesgo de una cartera se compor-
ta estadsticamente como una variable lognormal, sea cual fuere el nmero de
segmentos que integren esa cartera, supone una incongruencia porque la suma
de dos o ms lognormales no es una variable logartmiconormal124.
Para intentar justificar que esta incongruencia no tiene en la prctica efectos de-
masiado perturbadores, partimos de la frmula del coeficiente de asimetra de
una variable logartmico-normal de parmetros (m,s):
Hemos de insistir, de todas formas, en que este intento de justificacin slo tiene
alguna validez cuando las desviaciones tpicas son bajas.
124
Sobre la suma de variables log-normales puede consultarse Dufresne (2008).
Este hecho viene confirmado por la Figura 4.10 y la ltima columna del Cuadro 4.17. Ambos
125
muestran la cercana del VaR al 99,5% segn la distribucin normal bivariante y segn la aproxima-
cin lognormal, cuando la desviacin tpica marginal es 0,1.
206
A. Seguro novida, que comprende las lneas 1 a 12.
1. Ante todo, han de determinarse sus dos componentes, que son el volumen del
riesgo de prima y el de reserva, Vprem,s y Vres,s respectivamente, de cada segmen-
to. El primero, Vprem,s es la suma de las primas de doce meses (de los prximos
o de los ltimos, segn sea la cifra ms alta) y el valor actual esperado de las
207
que se devengarn despus de los doce meses siguientes, ya sea por los con-
tratos existentes, ya por aqullos cuya fecha de reconocimiento inicial se site
dentro de los doce meses siguientes. Las primas que entran en el clculo son
las devengadas o por devengar, con todos sus componentes (y no las primas de
riesgo). El volumen del riesgo de reserva, Vres,s es la mejor estimacin de las
provisiones para prestaciones pendientes.
2. Al 25 por 100 de la suma Vprem,s + Vres,s se aplica una reduccin, en su caso, por
diversificacin geogrfica. Pasamos por alto esta cuestin por no ser indispen-
sable a nuestros efectos. El resultado de aplicar tal reduccin es Vs, es decir, el
volumen de riesgo de prima y de reserva del segmento s.
Las desviaciones tpicas del sub-mdulo primas y reservas entre los distintos
segmentos son las especificadas en el Anexo II del Reglamento.
El ltimo paso es sumar las desviaciones tpicas de todos los segmentos, para
llegar a la desviacin tpica del riesgo de primas y reservas de todo el negocio no
vida. Esto se logra por medio de la expresin:
(s) y (t) son las desviaciones tpicas del riesgo total primas y reservas, de los
segmentos s y t respectivamente.
208
Vs y Vt representan los volmenes de riesgo de prima ms reserva de dichos seg-
mentos, s y t.
Por ltimo, nl y Vnl son las mismas medidas, referidas a toda la cartera novida.
1. Una vez determinados los volmenes de prima y de reserva para cada segmen-
to, se obtiene la desviacin tpica del riesgo conjunto de prima y de reserva de
cada segmento s, desviacin tpica a la que llamaremos s (para s =1,2,...).
209
En los ejemplos del Captulo IV, por razones de simplicidad, no hemos distinguido
entre riesgo de prima y de reserva.
En este Captulo nos proponemos ilustrar cmo las cpulas pueden ser un instru-
mento til en los modelos de riesgo para aseguradores sujetos a la normativa de
Solvencia II. Tal propsito vera muy limitado su alcance de no disponer de datos rea-
les que prestaran indiscutible conexin con la realidad al anlisis subsiguiente. Ahora
bien, una de las dificultades planteadas por las aplicaciones de la Teora de Cpulas
estriba en la obtencin de cifras empricas suficientemente numerosas y fiables.
Hemos obtenido los nmeros bsicos que vamos a utilizar en nuestro estudio del
documento del CEIOPS FS-14/07 titulado QIS3, Calibration of the Underwriting Risk,
Market Risk and MCR de fecha Abril de 2007. Su finalidad es describir las hiptesis y
lgica subyacentes en el calibrado de los antedichos riesgos (de suscripcin y de
mercado).
210
refiriendo reconoce que la escasez de datos disponibles impone un cierto grado de
juicio o valoracin subjetiva y aade algo importante para nosotros:
Sentados estos principios, llega al punto que nos ha proporcionado los datos em-
pricos que vamos a utilizar (pgina 21 del Documento CEIOPS): dice que, para
determinar la matriz de correlaciones para el riesgo de prima se llev a cabo un
anlisis de las correlaciones entre ramos a nivel de aseguradores individua-
lesen el mercado alemn. Se utilizaron ratios de siniestralidad del perodo com-
prendido entre 1988 y 2002 de forma que, para cada asegurador individual (en que
hubiera disponibles al menos 10 ratios en cada par de ramos), se estim el coefi-
ciente de correlacin y el coeficiente global como media de las correlaciones es-
pecficas de cada asegurador. Como ejemplo de lo dicho, muestra un diagrama de
dispersin (pgina 22 del Documento) en que se representan los residuos estn-
dar con respecto a los ratios de los seguros de RC automviles, en abscisas y RC
general, en ordenadas. Se trata de 1269 puntos representativos de otros tantos
pares de ratios de siniestralidad y gastos sobre primas, procedentes de 89 empre-
sas, para dichos segmentos de seguro. Adems de esos 1269 puntos, el grfico
contiene la recta de regresin de las ordenadas sobre las abscisas y cuya pen-
diente es 0,28, igual al coeficiente de correlacin lineal. Los datos a que nos refe-
rimos son residuos, es decir, diferencias entre los ratios individuales observados
y los generales, y con desviacin tpica igual a la unidad. El documento no es muy
explcito sobre esto; entendemos que los residuos positivos representan excesos
sobre la media y los negativos, defecto. Y suponemos que la media general de es-
tos ratios es la unidad, porque esta hiptesis subyace en la frmula que, segn
vimos antes126, conduce a la regla de que SCRnl es el triple del producto de la des-
viacin tpica por el volumen de riesgo.
Expresin (5.3)
126
211
Pues bien, es precisamente este grfico el origen de los datos que utilizaremos
a continuacin. Al no haber resultado posible obtener las coordenadas de
cada uno de estos puntos, no hemos tenido ms remedio que calcularlas por
medicin de la abscisa y la ordenada respectiva; esta labor, adems de tediosa,
implica el riesgo de errores, que hemos de aceptar al carecer de alternativas y
en la confianza de que no sean sustanciales en su cuanta ni sistemticos en su
signo.
212
Figura 5.1: Diagrama de dispersin de los ratios de siniestralidad sobre primas una vez tipificados
para que su media sea nula y su desviacin tpica la unidad, de los seguros de RC autos (abscisas)
y RC general (ordenadas).
Fuente: Elaboracin propia a partir de datos de: Documento CEIOPS FS-14/07 QIS3, Calibration of the
Underwriting Risk, Market Risk and MCR, abril 2007 y elaboracin propia segn lo explicado en el texto.
En definitiva, pues, admitimos que los datos representados en la Figura 5.1 tienen me-
dia nula y desviacin tpica igual a uno y proceden de una distribucin bivariante (X,Y), en
donde X tiene una desviacin tpica de 0,1 e Y de 0,14, mientras que la media es la uni-
dad en ambas variables. La Figura 5.2 contiene el diagrama de dispersin de (X,Y).
Figura 5.2: Diagrama de dispersin de los ratios de siniestralidad sobre primas, de los seguros de
RC autos, en abscisas y RC general, en ordenadas.
Elaboracin propia a partir de datos del Documento CEIOPS FS-14/07 QIS3, Calibration of the
Underwriting Risk, Market Risk and MCR, abril 2007 y elaboracin propia segn lo explicado en el texto.
213
Estos pares de valores (X,Y) son el dato bsico para nuestro trabajo. Sus caracte-
rsticas estadsticas principales se resumen en el Cuadro 5.1.
Las diferencias entre las medias y las desviaciones tpicas, por una parte, y los
valores que hemos decidido atribuir a las mismas segn lo ya indicado (1 para
ambas medias y 0,1 0,14 para las desviaciones) se deben a errores de redondeo.
El coeficiente de correlacin entre X e Y es 0,29.
214
cualesquiera de la distribucin emprica, (xi,xi+1) o bien (yj,yj+1) la distribucin de la
densidad es uniforme.
Una vez fijado el umbral, hay varias alternativas para determinar los parmetros
del modelo terico. Una solucin posible es la siguiente: Admitimos que el modelo
215
lognormal que vamos a utilizar tiene media igual a 1, con lo cual m=-s2/2, y deter-
minamos el parmetro s de forma que:
(a) (b)
Funciones de distribucin de X en (a) e Y en (b).
En rojo la distribucin emprica y en azul la curva logartmico-normal ajustada a cada cola derecha a
partir de x>1,2726 y de y>1,4095.
Figura 5.3.
5. LA CPULA EMPRICA
216
El Cuadro 5.2 ofrece en forma resumida los datos de esta cpula emprica para
nuestros datos bsicos, constituidos por los 1269 puntos representados en la
Figura 5.2. Cada fila del Cuadro representa un tramo de frecuencias acumuladas
para la variable U, de forma que la primera fila se refiere a valores u=FX(x) con
frecuencias no superiores a 0,1, la segunda a 0,2 y as sucesivamente hasta la d-
cima, referente a frecuencias mximas. Y lo mismo para las columnas, referidas
a la variable v=FY(y) (recordemos que las funciones FX y FY son las funciones em-
pricas y continuas, ajustadas segn lo descrito en el Epgrafe anterior).
Es decir, que el elemento (i,j) de la matriz significa la proporcin de pares (u,v) tales
que u<0,1.i, y, adems, v<0,1.j. Por ejemplo, para i=3, j=5 tenemos el nmero 0,2033,
que representa lo siguiente, recapitulando todo lo dicho para mayor claridad:
1. Partimos de nuestra base de datos constituida por los 1269 puntos (Figura 5.2)
representativos de nuestra distribucin bivariante de ratios RC autos (variable
X) y RC general (variable Y).
217
4. Computamos el nmero de dichos pares tales que u<0,3, v<0,5. La proporcin
de este nmero en relacin con el total de 1269 pares, es 0,2033. Dicha propor-
cin es, por tanto, el valor de la cpula emprica C(u,v)= C(0,3;0,5) o bien la fre-
cuencia de aquellos casos en que u<0,3 y v<0,5. Si las funciones de distribucin
marginales son continuas, se verificar que el valor de la cpula emprica en un
punto, C(uo,v0), significa tambin la proporcin de pares (x,y) (en donde x=FX-1(u),
y=FY-1(v)), tales que, simultneamente, x FX-1(uo), y FY-1(vo).
Hemos obtenido la idea de estos grficos como medio de evaluar a simple vista el signo y grado
127
218
se ve cmo las curvas de nivel de la cpula emprica estn comprendidas entre las
curvas de las cpulas mxima y mnima. El coeficiente de Kendall para la cpula
emprica es =0,20289.
Figura 5.5: Curvas de nivel de la cpula emprica y su comparacin con las de la cpula de
independencia.
219
Figura 5.6: Curvas de nivel de la cpula emprica y su comparacin con las de la mxima.
Figura 5.7: Curvas de nivel de la cpula emprica y su comparacin con las de la mnima.
Nos planteamos ahora la bsqueda del modelo de cpula que se adapte mejor a
nuestros datos.
220
En primer lugar, utilizaremos el mtodo grfico de Genest y Rivest, aplicable a las
cpulas arquimedianas y que nosotros proponemos extender, en las condiciones
que ms adelante veremos, a otras cpulas bivariantes. Segn qued expuesto en
III,2, este mtodo consiste en comparar la grfica de la funcin K(w)=w(w)/(w)
con las frecuencias del suceso C(u,v)w, en donde C(u,v) es la cpula emprica.
Ms concretamente: Recordemos que
w 0,05 0,10 0,15 0,20 0,25 0,30 0,35 0,40 0,45 0,50
K(w) 0,183 0,293 0,378 0,449 0,519 0,57 0,615 0,669 0,772 0,774
w 0,55 0,60 0,65 0,70 0,75 0,80 0,85 0,90 0,95 1,00
K(w) 0,812 0,852 0,885 0,911 0,931 0,953 0,967 0,984 0,998 1,00
Cpula de Gumbel
(w) = (LN(w))
221
de donde:
La representacin grfica de K(w), junto con los puntos (w, K(w)) de la cpula em-
prica, est en la Figura 5.8.
Figura 5.8: Funcin K(w) para la cpula de Gumbel, con parmetro =1,204 y su comparacin con
los valores de la cpula emprica.
Cpula de Clayton
(w)=(w1)/
222
y, por tanto, K(w)128:
Figura 5.9: Funcin K(w) para la cpula de Clayton, con parmetro =0,27 y su comparacin con
los valores de la cpula emprica.
Cpula de Frank
utilizada para realizar clculos y generar grficos pueden presentar algunas diferencias en cuanto
a tipo de letra y smbolos sobre el utilizado generalmente en este trabajo.
223
Por otra parte, recordemos que en el modelo de Frank,
224
Figura 5.10: Funcin K(w) para la cpula de Frank, con parmetro =1,946 y su comparacin con
los valores de la cpula emprica.
Cpula gaussiana
en donde x(u) y v(y) son las funciones de distribucin inversas de X e Y, es decir que
equivalen, respectivamente, a x= 1(u) e y= 1 (v), como ya vimos en su momento;
(.) es la funcin de distribucin normal estndar y es el coeficiente de correlacin.
Para obtener una aproximacin de la funcin K(w) para la abscisa w=w0 recurri-
mos a la simulacin estocstica mediante un programa que realiza los pasos si-
guientes:
225
1. Simula dos nmeros aleatorios (pseudoaleatorios) normales (0,1) con el coefi-
ciente de correlacin elegido.
4. Tras repetir los pasos anteriores muchas veces, da el cociente entre el nmero
de xitos y el de pares aleatorios generados.
Para =0,202894 resulta que =0,313. Elegimos, no obstante, el valor =0,29, que
es la estimacin mximo-verosmil, segn veremos ms adelante.
Si queremos obtener una curva que se adapte a los puntos obtenidos y que se re-
flejan en el Cuadro 5.5, obtenemos la siguiente por el mtodo de los mnimos
cuadrados:
w 0,05 0,10 0,15 0,20 0,25 0,30 0,35 0,40 0,45 0,50
K(w) 0,16 0,27 0,37 0,44 0,52 0,59 0,65 0,70 0,74 0,79
w 0,55 0,60 0,65 0,70 0,75 0,80 0,85 0,90 0,95 1,00
K(w) 0,83 0,86 0,88 0,91 0,94 0,96 0,97 0,99 0,996 1,00
Cuadro 5.5: Valores de la funcin K(w) para la cpula gaussiana.
226
En la Figura 5.11 quedan representados los puntos (w,K) segn la cpula empri-
ca, as como la funcin ajustada K(w).
Figura 5.11: Aproximacin de la funcin K(w) para la cpula gaussiana, con parmetro =0,29 y su
comparacin con los valores de la cpula emprica.
Cpula t
Su ecuacin es:
227
Figura 5.12a: Aproximacin de la funcin K(w) para la cpula t, con n=4 grados de libertad y
parmetro =0,29 y su comparacin con los valores de la cpula emprica.
Figura 5.12b: Ampliacin de las grficas de la Funcin K(w) ajustada por mnimos cuadrados de
las cpulas de Gauss y t, en el tramo de abscisas (0,040,16), aproximadamente.
228
En efecto, los dos nicos puntos que aparecen en esta ampliacin, correspondientes a
la funcin K(w) de la cpula emprica, se hallan distanciados de las curvas y la corres-
pondiente a la t de Student es la ms alejada de ambas. En el Epgrafe 8, en que deter-
minaremos la respectiva medida de Frees y Valdez, confirmaremos esta conclusin.
Hacemos ahora un breve apartado para referirnos a este tipo de grficos, que ya uti-
lizamos en IV,1. En lugar de las anteriores comparaciones entre la funcin K(w) de
cada modelo, representada como una curva, y varios valores de la misma funcin,
pero procedentes de la cpula emprica y que representamos por varios puntos, po-
dramos habernos servido de los grficos qq, que dan imagen del grado de cercana
entre dos funciones de distribucin, F1 y F2 por medio de los puntos (F1(xi), F2(xi)) (i=1,2,
... n). Cuanto ms prximos se encuentren dichos puntos a la recta y=x, ms cercanas
estn ambas funciones de distribucin. Aplicaremos este mtodo slo a dos casos de
los examinados: Las cpulas de Frank y de Clayton. El grfico obtenido representar
los puntos (KE(wi), K(wi)), siendo la abscisa KE(wi) de cada uno el valor de la funcin K
para el punto wi segn la cpula emprica, mientras que la ordenada significar el
mismo concepto, pero segn el respectivo modelo terico; distinguiremos a dicha or-
denada con el subndice F o C segn se trate del modelo de Frank o del de Clayton).
w 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0
KE(w) 0,293 0,449 0,57 0,669 0,774 0,852 0,911 0,953 0,984 1,0
KF(w) 0,274 0,439 0,57 0,679 0,771 0,849 0,911 0,959 0,989 1,0
KC(w) 0,271 0,461 0,608 0,725 0,816 0,886 0,938 0,973 0,993 1,0
Cuadro 5.6: Valores de la funcin K(w) para las cpulas emprica, de Frank (=1,9460) y de Clayton
(=0,27).
229
En la Figura 5.13 se representan los grficos qq para las cpulas de Frank y de
Clayton. Se observa que es mucho ms adecuada para nuestros datos la primera
de ellas, al estar ms prximos sus puntos a la lnea continua, expresin grfica
de la recta y=x.
(a) (b)
Figura 5.13: Grfico qq para comparar la cpula emprica con las de Frank (a) y de Clayton (b).
KE(w,n)P(t, ) = 7,5158w521,4353w4+22,5168w311,4751w2+3,8358w
El Cuadro 5.7 contiene la medida resultante para cada una de las cpulas consi-
deradas:
230
Cpula Parmetro Medida de Frees y Valdez x 1000
Frank 1,946 0,114
Gauss 0,29 0,277
Gumbel 1,204 0,294
t =0,29; n=4 0,440
Clayton 0,27 0,886
Cuadro 5.7: Medida de Frees y Valdez para los distintos modelos de cpula, ordenados de mejor a
peor.
9. LA FUNCIN DE VEROSIMILITUD
Otro mtodo para seleccionar el modelo de cpula ms adecuado a los datos con-
siste en calcular la verosimilitud y determinar el parmetro que la hace mxima.
Dicho mtodo se llama de la mxima pseudo-verosimilitud o mxima verosimi-
litud cannica, porque los valores de la cpula que entran en el cmputo no re-
sultan de una cpula emprica obtenida a partir de las distribuciones marginales
de las variables observadas, sino de una estimacin no paramtrica de dichas
distribuciones marginales.
Pues bien, el Cuadro siguiente refleja los valores del logaritmo de la mxima ve-
rosimilitud de cada modelo de cpula de los cinco que venimos examinando y
el valor correspondiente del parmetro130. El importe reflejado es, por tanto,
u=F1(x), v=F2(y).
231
Mximo
Cpula Parmetro Logaritmo
verosimilitud
Frank 1,946 61,596
Gauss 0,29 56,01
Gumbel 1,2040 48,627
t (4 g de l) 0,29 36,92
Clayton 0,27 25,14
Cuadro 5.8: Logaritmo del mximo de la funcin de verosimilitud de cada modelo de cpula, y
estimacin mximo verosmil del parmetro de dependencia. Los modelos estn ordenados de
mejor a peor.
El orden de los modelos coincide con el del Cuadro 5.7, basado en la medida de
Frees y Valdez. Y tambin los mtodos grficos utilizados en los Epgrafes 6 y 7
nos llevan a igual desenlace. Por otra parte, es evidente que el criterio de Akaike
conduce a idntica jerarqua; dicho criterio consiste en asignar a cada modelo el
valor:
232
(a) (b)
Figura 5.14: Funciones de distribucin emprica de X e Y (a) y (b) respectivamente (en rojo) y
lognormal mximo-verosmil (en azul).
m s
X 0,00695 0,0944
Y 0,01687 0,1465
Cuadro 5.9: Estimaciones mximo-verosmiles de los parmetros del modelo lognormal para X e Y.
a) Distribucin de X:
233
Las frecuencias tericas de cada intervalo (a,b) se obtienen con la integral:
El valor del estadstico 2 es ahora 5,68 y el nmero de grados de libertad es, igual-
mente, 5. El valor de las tablas para el nivel de significacin 5 por 100 es 11,07, por
lo que no rechazamos la hiptesis nula, consistente en admitir que la muestra de
valores de Y proviene de una poblacin logartmico-normal con parmetros m=
0,01687, s=0,1465. El valor p es 0,338.
234
En resumen: al nivel de significacin 5 por 100,
No hay evidencia suficiente para rechazar la hiptesis de que Y sigue dicho mo-
delo (m= 0,01687, s= 0,1465).
en donde K(.) es una funcin de densidad, que en nuestro caso ser la normal; xi,
el valor i-simo observado, n el nmero de elementos de la muestra y h es un pa-
rmetro que mide el grado de suavizado que logra la frmula: A mayor valor de h,
ms alta es la dispersin de la densidad representada por la funcin K(.); por tan-
to, ms amplio es cada sumando de la frmula y ms fuerte el suavizado que sta
opera. Si h es pequeo ocurre lo contrario: La dispersin de la densidad que apor-
ta cada sumando y que se centra en la respectiva abscisa xi queda ms concentra-
da en ella y el suavizado es menor. El valor de h no debe elegirse arbitrariamente,
sino que ha de guardar relacin con la desviacin tpica de la distribucin de X a la que
pretende aproximar. Una frmula prctica para elegir el valor de h es: h1,06/n0.2;
235
en nuestro caso, =0,10 y por tanto h0,0254. El resultado es una funcin de den-
sidad emprica y suavizada, que se representa en la Figura 5.15.
236
11. MEDIDAS DE RIESGO PARA UNA CARTERA DE SEGUROS COMPUESTA
POR RIESGOS X O POR RIESGOS Y
Empezaremos con el clculo del VaR al 99,5 por 100 (nivel de confianza exigido por
la normativa de Solvencia II) de cada uno de los riesgos X e Y. Se trata de resolver
la ecuacin
Fi(VaR)=0,995
ES(Z)=E(ZZ>VaR)
237
Para determinar ES99.5(Y) calculamos el valor de
Supondremos que la cartera est integrada por ambos tipos de riesgos, X e Y (ra-
tios de siniestralidad sobre primas en el seguro de RC automviles y RC general,
respectivamente) en proporciones respectivas p y q (p+q=1). Dichas proporciones
se refieren al volumen de primas.
238
Por tanto,
Si sustituimos las desviaciones tpicas por sus valores respectivos segn la nor-
mativa de Solvencia II (x=0,10; y =0,14) y tenemos en cuenta que q = 1p, resulta
que la desviacin tpica nl es una funcin nicamente de p, al igual que el capital
mnimo obligatorio por unidad de volumen de primas, cuya expresin ser:
239
En estas condiciones, la funcin de distribucin conjunta de las tasas o ratios de
siniestralidad y gastos sobre primas para los seguros de RC automviles (X) y RC
general (Y), ser:
240
Este mtodo tiene la ventaja de que es fcil de plantear y programar y que es flexi-
ble, al prestarse sin problemas a complejidades que exija el caso particular que
interese. Su inconveniente, en casos como el que ahora nos ocupa, es que requie-
re un nmero muy elevado de repeticiones; por ejemplo, cuando la probabilidad
de xito est en torno al 99,5 por 100 y deseemos una probabilidad muy alta (por
ejemplo, del 95 por 100) de que el error de la estimacin resultante del proceso
sea muy pequea (por ejemplo, una diezmilsima, de modo que si la verdadera
probabilidad de insolvencia es el 0,5 por 100, la estimacin quede comprendida
entre 0,49 y 0,51 por 100), el nmero de simulaciones no deber ser inferior a
1.911.196, que resulta de despejar n en la inecuacin 1.96 <
0,0001.
Se simula un nmero elevado (por ejemplo, 2000) de sumas X+Y mediante el pro-
cedimiento descrito en los prrafos anteriores. Obtenidos los n valores, se busca
una funcin de distribucin que se adapte a la muestra simulada y se toma esta
funcin como modelo de la suma T = X+Y.
241
Figura 5.17: Simulacin de 2000 puntos (u,v) de la cpula de Frank con =1,946.
m = 0,6848
s = 0,09724.
242
Figura 5.18: Los 2000 puntos (x,y) correspondientes a los representados en la Figura 5.17.
Coeficiente de
Valor mximo Valor mnimo Media Desviacin tpica
asimetra
2,5960 1,4010 1,993 0,194 +0,244
Cuadro 5.14: Algunas caractersticas de la muestra de T=X+Y obtenida por simulacin.
El valor del estadstico 2 es 9,75; para un nivel de significacin del 95 por 100, el
valor para 1112 = 8 grados de libertad, segn las tablas, es 15,507; el valor p
es 0,28; por tanto, no rechazamos la hiptesis de que el modelo lognormal con
parmetros m=0,6848; s=0,09724 sea adecuado para la variable T=X+Y.
Sobre esta base, calculamos el VaR al 99,5 por 100, de la variable T; dicho valor es
2,548 y el expected shortfall 2,63; pero si referimos los resultados a una cartera
que, en lugar de tener media 2 tuviera media unitaria quedara:
VaR0.995 = 1,274
ES0.995= 1,315
243
Figura 5.19: Funcin de distribucin emprica de las 2000 simulaciones de T.
244
c) Mtodo del clculo directo
F(x,y)=C[u(x), v(y)]
245
por tanto,
y, definitivamente,
ui* = FX(tyi)
vi = FY(yi).
Prob(X+Y2,532) = 0,995
Para una cartera con un volumen de primas igual a la unidad, el VaR es, por tanto,
2,532/2 = 1,266. Si la media de la distribucin de dicha cartera X+Y es la unidad,
resultar un capital de solvencia obligatorio de 0,266. El capital determinado
246
anteriormente por el mtodo de las correlaciones para la misma cartera es 0,313,
superior, por tanto, en un 18 por 100, aproximadamente.
Multiplicamos la derivada parcial C(u, v)/v en el punto (ui*, vi) por la densidad
de Y en el punto yi; multiplicamos el resultado por t/(n.q). La suma, cuando
introducimos un alto nmero de sub-intervalos, se acerca a la probabilidad
buscada.
As, calculamos el VaR al 99,5 por 100 para distintos valores de p y el resultado lo
tenemos en el Cuadro 5.16.
131
Bowers, Gerber, Hickman, Jones y Nesbit, (1986) pgina 62. Hay ediciones posteriores.
247
en donde:
***
En efecto,
Ahora bien, por hiptesis E(Tp)<. Con mayor razn, la integral desde t a infini-
to ser tambin finita y tender a cero al aumentar el lmite inferior de integra-
cin. Como quiera que t. (t;p)0 y, adems, es no superior a una cantidad que
se hace arbitrariamente pequea al tender t a infinito, queda demostrada la
expresin (5.7) y con ello la frmula (5.6), cuyo clculo es fcil de programar.
***
248
En el Cuadro 5.16 tenemos las dos medidas de riesgo para distintas composicio-
nes de la cartera (X,Y) cuyo ratio de siniestralidad sobre las primas sea la unidad:
En nuestro caso hemos llegado a determinar que, con la cpula de Frank y las
distribuciones marginales que admitimos como modelos para los riesgos X e Y
249
(coberturas del riesgo de responsabilidad de automviles y general respectiva-
mente), el capital de solvencia obligatorio (SRC) era de 0,266 por unidad de volu-
men de una cartera 50%-50% de ambos riesgos. Si queremos conocer el importe
del coeficiente que, segn el mtodo de las correlaciones auspiciado por
Solvencia II, conducira a un capital de solvencia equivalente a nuestro resultado
de 0,266, no tenemos ms que despejarlo en la ecuacin:
es decir, = 0,066, resultado muy alejado del coeficiente 0,5 segn la normativa de
Solvencia II para los ramos o segmentos que nos ocupan. Y si introdujramos en
el anterior clculo el coeficiente de correlacin de los datos originales, que, segn
vimos en el Epgrafe 2, es 0,28, el resultado sera:
Estos coeficientes de correlacin se han obtenido utilizando tan slo los puntos
(X,Y) cuya suma supera el cuantil correspondiente. El que sean negativos en todos
los casos sugiere que los datos que manejamos no justifican elevar la correlacin
para intentar modelizar dependencias en la cola que, en nuestro caso, no existen.
250
REFERENCIAS BIBLIOGRFICAS
CINTAS DEL RO, (2006): Teora de cpulas y control del riesgo financiero. Tesis
Doctoral. Universidad Complutense de Madrid, Facultad de Ciencias Matemticas.
Madrid.
DE MEY, Insurance and the Capital Markets, (2007): The Geneva Papers, 2007, 32
(35-41).
251
DONNELLY y EMBRECHTS, (2010): The Devil is in the tails: actuarial mathematics and
the subprime mortgage crisis.- http://www.macs.hw.ac.uk/~cd134/2010/donemb.pdf
252
INTERNATIONAL ACTUARIAL ASSOCIATION, A Global Framework for Insurer
Solvency Assessment. Research Report of the Insurer Solvency Assessment Working
Party, 2004.
JONES, The formula that felled Wall Street, Financial Times, 24 de abril de 2009.
253
STAUDT, (2010): Tail Risk, Systemic Risk and Copulas.- Casualty Actuarial Society
E- Forum, Volume 2, 2010.
254
BIBLIOGRAFA ADICIONAL
BARNETT, KREPS, MAJOR y VENTER, (2007) Multivariate copulas for financial mo-
deling, Variance, 2007, Vol. 1, No. 1, 103-119.
DEHEUVELS, (1981): A non parametric test for independence. Publ. Inst. Statist.
Univ. Paris 26: 29-50.
255
DESMEDT y WAHLIN, (2008) On the subadditivity of tail value at risk: An investiga-
tion with copulas, Variance, 2008, Vol. 2, No. 2, 231-252.
FRCHET, (1957): Les tableaux de corrlation dont les marges sont donnes.
Annals Univ. Lyon, Sci., 20:13-31.
GENEST y BOIES (2003): Detecting dependence with Kendall plots. The American
Statistician, 57, 275-284.
GENEST y FAVRE, (2007): Everything you always wanted to know about copula mo-
delling but were afraid to ask. Journal of Hydrologic Engineering, http://citeseerx.
256
ist.psu.edu/viewdoc/download;jsessionid=CC1B2A1F9164365AAF6BE126F579012
3?doi=10.1.1.160.8266&rep=rep1&type=pdf
GENEST y MCKEY, (1986): The joy of copulas: Bivariate distributions with uniform
marginals. American Statistics, 40: 280-285.
JOE, (1997): Multivariate models and dependence concepts. Chapman & Hall.
NELSEN, (2002): Concordance and copulas: a survey. Kluwer. Dordrecht, pgs 169-177.
PICKANDS, (1975): Statistical inference using extreme order statistics. The Annals
of Statistics 3: 199-131.
SAMARA y RANDLES, (1988). Kendalls tau. Comm. Statist. Theory Methods, 17:
3191-3205.
SHIH, y LOUIS, 1995): Inferences on the association parameter in copula models for
bivariate survival data. Biometrics, 51: 1384-1399.
257
SKLAR, (1959): Fonctions de rpartition n-dimensions et leur marges. Publ. Inst.
Statist. Univ. Paris. 8: 229-231.
SMITH, (2002): Risk management: Value at risk and beyond. Cambridge U.P.
TAJVIDI, (2004): Confidence intervals and accuracy estimation for heavy-tailed gene-
ralized Pareto distributions. Extremes 6: 111-123. Kluwer Academic Publishers.
Netherlands.
258
NDICE ALFABTICO
259
densidad: 78 Estimacin no paramtrica: 151
idea intuitiva: 65 y ss. Expected shortfall: 200, 237 y ss
implcitas: 38, 88, 124 Ferri: 34
mixtas: 163 y ss. Frchet y Hoeffding, cotas: 79, 82, 134
multivariantes: 172 y ss. Frees y Valdez: 32, 134
probabilidad condicionada: 78, 79 Frees y Valdez, medida: 147, 168, 230,
simulacin: 152 y ss, 158 231
Correlacin por rangos: 34, 62, 65 Funcin arco seno: 78
Densidad kernel: 235 Funcin completamente montona:
Dependencia: 43 y ss., 82 y ss. 173
Dependencia en colas: 87 (ver tambin Funcin de Debye: 110, 111
la de cada tipo de cpula) Funcin generador: 94 y ss.
Directiva 2009/138EC: 199 Generador estricto: 94, 95
Distribucin Genest y Rivest, mtodo grfico de
de Cauchy: 53 ajuste: 146, 221
logartmico-normal: 178, 182, 189, H-medida: 102
205, 233 y ss. H-volumen: 102
normal bivariante: 43 IFM (Inference function for margins):
Pareto: 135, 164, 176, 182 151, 214
t bivariante: 44 Iman y Conover, mtodo: 33
Distribuciones elpticas: 43, 47 y ss., Independencia de variables aleatorias:
53 y ss. 60
Distribuciones esfricas: 45 y ss., 51 y Industry loss warranties: 36
ss. Jacobiano: 57
Distribuciones marginales empricas: Kimberling, teorema: 174
150 y ss. Klugman, Panjer y Willmott: 167
Durn y Otero: 33 Li, David X.: 23
Embrechts: 24, 25 Matriz
Embrechts, Lindskog y McNeil: 173, de correlaciones: 47
175 de covarianzas: 47, 59
Embrechts, McNeil y Straumann: 43, ortogonal: 52
45, 60, 64 Mxima verosimilitud, mtodo: 149,
Estadstico de orden: 143, 151 231 y ss.
Estimacin mximo-verosmil McNeil, Frey y Embrechts: 77, 82
cannica o pseudo-verosmil: 151 y 152
260
Meta-distribuciones: 26, 113 de reserva: 203 y ss.
densidad: 114 Rivas, Prez y Cuesta: 36
Mtodo del tau acumulativo o tau Schmidt: 124, 128
parcial: 147, 148 Schweizer y Wolf, medida: 85
Mikosch: 149 Simulacin de cpulas: 152 y ss, 158
Modelo estndar (de Solvencia II): 200 Simulacin de nmeros aleatorios
Moscadelli: 51 normales: 52
National Association of Insurance Sklar, Teorema: 25, 75, 113
Commissioners (NAIC): 28 Solvencia II: 30, 31, 178 y ss., 188 y ss.,
Nelsen: 66, 75, 86, 93, 94, 96, 99, 100, 199 y ss.
102, 128, 131, 133, 137 hiptesis implcitas: 204, 205
Propiedades invariantes al cambio de Suma de variables aleatorias: 176 y
escala: 60, 75, 82 ss., 245 y ss.
Pseudo-inversa: 94 Staudt: 32, 151
Pseudo-observaciones: 151 Swiss Solvency Test: 28
Ratio de siniestralidad: 176 y ss., 188 Tail conditional expectation: 200, 237 y
Reglamento Delegado 2015/35 de la ss
Comisin: 199 The Economist: 25
Regresin mediana y cuantlica: 169 y Transformacin
ss afn: 53, 55
Riesgo montona: 76, 77
de catstrofes: 207, 208, 213, 214 ortogonal: 51, 52
de cada de cartera: 207, 208, 213, VaR (Value at risk): 185 y ss., 188 y ss.,
214 199, 237
de prima: 203 y ss. Venter: 32, 115, 122, 145
261
COLECCIN CUADERNOS DE LA FUNDACIN
205. Current Topics on Risk Analysis: ICRA6 and Risk 2015 Conference. 2015
202. Valoracin de los inmuebles del patrimonio histrico y los riesgos ssmicos en el
contrato de seguro: el caso de Lorca. 2014
183. Claves del Seguro Espaol: una aproximacin a la Historia del Seguro en
Espaa. 2012
180. Origen, situacin actual y futuro del seguro de Proteccin Jurdica. 2012
155. Los sistemas de salud en Latinoamrica y el papel del seguro privado. 2010
119. Anlisis del riesgo en seguros en el marco de Solvencia II: Tcnicas estadsticas
avanzadas Monte Carlo y Bootstrapping. 2008
110. El Fair Value de las provisiones tcnicas de los seguros de Vida. 2007
103. Las NIC y su relacin con el Plan Contable de Entidades Asegura-doras. 2006
72. Anlisis del proceso de exteriorizacin de los compromisos por pensiones. 2002
63. Anlisis de la repercusin fiscal del seguro de vida y los planes de pensiones.
Instrumentos de previsin social individual y empresarial. 2001
62. Seguridad Social: temas generales y rgimen de clases pasivas del Estado. 2001
27. La calidad total como factor para elevar la cuota de mercado en empresas de
seguros. 1995
26. El coste de los riesgos en la empresa espaola: 1993
LIBROS
The risk of longevity and its practical application to Solvency II. 2015
Desde 1994 se publican anualmente estudios que presentan una panormica concreta
de los mercados aseguradores europeos, de Espaa e Iberoamrica y que pueden con-
sultarse en formato electrnico desde la pgina Web:
www.fundacionmapfre.org
Seguros
Gerencia de riesgos
Prevencin
Centro de Documentacin
www.fundacionmapfre.org/documentacion
P.V.P.: 30
ISBN 978-84-9844-647-0