Cálculo de Seguros de Vida Continuos
Autor: Santiago Yael Morales Torres
Materia: Análisis Numérico
Resumen: Los seguros de vida corresponden al cálculo de todas las obligaciones futuras que
tiene una compañía al momento de pactar un contrato de seguro. Estas obligaciones son
estimaciones de la probabilidad de pagar la Suma Asegurada pactada lo largo de la duración del
contrato de seguro. Para calcular estos flujos, se tienen 2 enfoques: el enfoque discreto, en el cual
se estima la probabilidad de que el asegurado sobreviva cierta cantidad de años y muera al siguiente
año (e.g. muera en el primer año, o que sobreviva 1 año y muera en el 2do año, o que sobreviva 2
años y muera en el 3er año, etc.), donde el valor del seguro es la suma (discreta) de todos estos flujos
posibles, y es el más usado actualmente por la facilidad de los cálculos, ya que a lo mucho, una
persona sobrevive entre 100 a 110 años. Sin embargo, uno de los supuestos de este enfoque discreto
es que el pago de la Suma Asegurada es al final del año de la muerte, lo cual en la realidad eso no
pasa casi nunca, y es por esto que existe otro enfoque más preciso pero mucho menos trabajado: el
enfoque continuo, en el cual en lugar de considerar periodos anuales se considera cada instante del
tiempo, por lo cual tendríamos infinitos flujos, y en lugar de una suma discreta tenemos una integral
para el cálculo, y es por esto que este enfoque no es muy utilizado. Por lo tanto, el objetivo de este
proyecto es utilizar las herramientas que tenemos para el cálculo de seguros de vida discretos, y
usarlas para construir las funciones necesarias para hacer la aproximación de un seguro de vida
continuo y realizar una comparación entre ambos enfoques para observar cuál de ellos es mejor
utilizar tanto numéricamente como computacionalmente.
Marco Teórico
Para la construcción de los seguros de vida en el caso discreto, se definen las siguientes variables:
x -> Edad de la persona que contrata el seguro.
n -> Temporalidad del seguro.
-> Edad máxima que un asegurado puede sobrevivir (en nuestro caso consideraremos = 110).
S.A. -> Suma Asegurada que paga por muerte el seguro.
Sx (t) = 𝑡𝑝𝑥 <- Probabilidad de que una persona de edad x sobreviva t años. (función de
supervivencia)
Fx (t) = 𝑡q 𝑥 <- Probabilidad de que una persona de edad x muera en los siguientes t años, es decir,
que la persona muera entre edades x y x+t. (Función de Distribución)
i <- Tasa de interés compuesto con la que descontarán los flujos.
1
v t = (1 + 𝑖)-t = (1+𝑖)t
<- Factor de descuento al tiempo t con interés compuesto.
Las funciones 𝑡𝑝𝑥 y 𝑡q 𝑥 son complementos, pues si una persona de edad x no sobrevive t años, es
porque murió entre edad x y x+t, y viceversa, por lo que existe la relación 𝑡𝑝𝑥 = 1 - 𝑡q 𝑥 . Por otro
lado, en caso de que t = 1, la notación suele omitir a la letra t, por lo que Sx (1) = 1𝑝𝑥 = 𝑝𝑥 es la
probabilidad que una persona de edad x sobreviva 1 año y análogamente Fx (1) = 1q 𝑥 = q 𝑥 es la
probabilidad de que una persona de edad x muera entre edades x y x+1. También debemos aclarar
que 0𝑝𝑥 se refiere a que la persona de edad x sobreviva 0 años, es decir, que sobreviva a este preciso
instante, lo cual obviamente siempre es 1, por lo que su complemento 0q 𝑥 es 0.
Por otro lado, para que una persona de edad x sobreviva a edad x+t, es necesario que primero la
persona de edad x sobreviva a edad x+1, despúes la persona de edad x+1 sobreviva a edad x+2 y así
sucesivamente hasta que una persona de edad x+t-1 sobreviva a edad x+t, por lo que también
tenemos la relación:
t-1
𝑡 p𝑥 = px px+1 px+2 . . . . px+t-2 px+t-1 = ∏ px+k
k=0
Y que por la formula complemento 𝑡𝑝𝑥 = 1 - 𝑡q 𝑥 , tenemos entonces que:
t-1
𝑡 p𝑥 = (1-q x )(1- q x+1 )(1- q x+2 ). . . . (1-q x+t-2 )(1- q x+t-1 ) = ∏(1-q x+k )
k=0
Una vez construidas las probabilidades, ahora si tenemos los elementos necesarios para construir
un seguro de vida. El seguro de vida para una persona de edad x y temporalidad n se define como
todos aquellos flujos futuros que la aseguradora espera pagar. En general la aseguradora pagará la
suma asegurada si la persona sobrevive t años, para 0 < t < n, y muere al siguiente año. Este flujo de
pago se expresa como:
S. A. 𝑡p𝑥 q x+t v t
Donde primeramente 𝑡p𝑥 q x+t expresa la probabilidad de que una persona de edad x sobreviva t
años y luego, teniendo edad x+t muera al siguiente año, es decir, en edades x+t y x+t+1, después
multiplicamos ese factor por la suma asegurada representando que en ese año esa sería la
proporción esperada de suma asegurada a pagar, y v t nos permite traernos ese flujo a valor presente
considerando que el pago de suma asegurada es al final del año de la muerte.
Por ejemplo, suponiendo un seguro de vida para una persona de edad 20 temporal 5 años, con una
suma asegurada de 100000 pesos y tasa de interés 5%, los flujos de la aseguradora serían:
100000 0p20 q 20 (1.05)-0 <- Flujo correspondiente a que el asegurado muera entre edades 20 y 21
100000 1p20 q 21 (1.05)-1 <- Flujo correspondiente a que el asegurado muera entre edades 21 y 22
100000 2p20 q 22 (1.05)-2 <- Flujo correspondiente a que el asegurado muera entre edades 22 y 23
100000 3p20 q 23 (1.05)-3 <- Flujo correspondiente a que el asegurado muera entre edades 23 y 24
100000 4p20 q 24 (1.05)-4 <- Flujo correspondiente a que el asegurado muera entre edades 24 y 25
Y para calcular por lo tanto el valor del seguro, es necesario sumar todos estos flujos.
Por lo tanto, definimos un seguro de vida discreto con suma asegurada S.A., para una persona de
edad x temporal n años con tasa de interés i como:
n-1 n-1
t -t
́̅ = ∑ S. A. 𝑡p𝑥 q x+t v = ∑ S. A. 𝑡p𝑥 q x+t (1 + i)
Ax:n|
t=0 t=0
Solo falta la cuestión de saber cómo calcular 𝑡p𝑥 y q x+t para cada tiempo t. La Circular Única de
Seguros y Fianzas nos proporciona una tabla para encontrar el valor de q x (probabilidad de que una
persona de edad x muera entre x y x+1) para todas las edades, desde x = 0 hasta x = 110.
CNSF M CNSF M CNSF M
Edad Edad Edad
2013 qx 2013 qx 2013 qx
0 0.000433 37 0.000896 74 0.015235
1 0.000433 38 0.000938 75 0.017009
2 0.000434 39 0.000983 76 0.019024
3 0.000434 40 0.001033 77 0.021312
4 0.000435 41 0.001087 78 0.023915
5 0.000436 42 0.001145 79 0.026879
6 0.000438 43 0.001208 80 0.030257
7 0.00044 44 0.001278 81 0.03411
8 0.000443 45 0.001353 82 0.038509
9 0.000446 46 0.001435 83 0.043533
10 0.000449 47 0.001525 84 0.049274
11 0.000453 48 0.001623 85 0.055833
12 0.000457 49 0.00173 86 0.063329
13 0.000463 50 0.001848 87 0.071889
14 0.000468 51 0.001977 88 0.08166
15 0.000475 52 0.002119 89 0.092798
16 0.000482 53 0.002274 90 0.105476
17 0.000489 54 0.002446 91 0.119875
18 0.000498 55 0.002635 92 0.136184
19 0.000507 56 0.002844 93 0.154594
20 0.000517 57 0.003074 94 0.175291
21 0.000528 58 0.003329 95 0.198441
22 0.00054 59 0.003612 96 0.224184
23 0.000553 60 0.003926 97 0.252613
24 0.000567 61 0.004275 98 0.28376
25 0.000582 62 0.004664 99 0.317576
26 0.000598 63 0.005096 100 0.353919
27 0.000616 64 0.005579 101 0.39254
28 0.000635 65 0.006119 102 0.433078
29 0.000656 66 0.006723 103 0.475068
30 0.000678 67 0.0074 104 0.517949
31 0.000703 68 0.00816 105 0.561099
32 0.000729 69 0.009015 106 0.603861
33 0.000757 70 0.009977 107 0.645589
34 0.000788 71 0.011061 108 0.685682
35 0.000821 72 0.012285 109 0.72362
36 0.000857 73 0.013668 110 0.758991
Con estos valores es suficiente para calcular las probabilidades necesarias para el cálculo del
seguro: q x+t (basta con buscarlo en la tabla) y 𝑡p𝑥 (mediante la formula ya descrita):
t-1
𝑡 p𝑥 = (1-q x )(1- q x+1 )(1- q x+2 ). . . . (1-q x+t-2 )(1- q x+t-1 ) = ∏(1-q x+k )
k=0
para toda edad x y todo tiempo t siempre que 0 <= x+t <=110. Nótese que conforme la edad aumenta,
q x es mayor, pues a mayor edad mayor es la probabilidad de fallecer en el siguiente año.
¿Cuál es el caso continuo?
En el caso continuo, buscamos ahora cubrir cada instante de tiempo. Para esto vamos a definir la
siguiente función:
𝜇x (t) = 𝜇x+t -> Fuerza de mortalidad para una persona de edad x+t. Es una tasa que mide que tan
propensa es una persona de edad x+t en un instante inmediato.
Con probabilidad condicional podemos definir la fuerza de mortalidad, la cual es definir la
probabilidad de que el tiempo de vida de una persona sea x+t dado que esta persona tiene al menos
edad x+t. Sea X la edad en la que la persona de edad x morirá:
P[X = x+t, X ≥ x+t ] P[X = x+t ]
𝜇x (t) = 𝜇x+t = P[X = x + t| X ≥ x + t] = P[X ≥ x]
= P[X ≥ x+t]
=
𝑑 𝑑 𝑑
P[X ≤ x+t] F (t) q
𝑑𝑥 𝑑𝑥 x 𝑑t t x
= =
P[X ≥ x+t] Sx (t) tpx
Esta fuerza de mortalidad 𝜇x+t reemplazará a q x+t en el caso continuo, pues ahora en cada tiempo t
en lugar de considerar la muerte entre edades x+t y x+t+1, considerará la muerte instantánea que le
puede ocurrir a la persona de edad x en cada punto t del tiempo, por lo que t ya no es necesariamente
un valor discreto, sino un valor continuo y para barrer todo tiempo t, en lugar de una suma, es
necesario usar una integral.
Por lo tanto, bajo estos supuestos, definimos un seguro de vida continuo con suma asegurada S.A.,
para una persona de edad x temporal n años con tasa de interés i como:
n n
̅ ́̅ = ∫ S. A. 𝑡p𝑥 𝜇x+t v t dt = ∫ S. A. 𝑡p𝑥 𝜇x+t (1 + i)-t dt
A x:n|
0 0
Justificación
El poder calcular seguros de vidas continuos brinda una alternativa adicional a la manera tradicional
de calcular los seguros de vida (el método discreto), además de que evita los sesgos que pueden
producirse por suponer que la suma asegurada será pagada al final del año, cuando la realidad es
que ningún asegurado esperará tanto tiempo para recibir el beneficio. La razón por la que son
necesarios métodos numéricos para el cálculo de seguros continuos, una vez conocida su
expresión, son principalmente 3:
- En las funciones Sx (t) = 𝑡𝑝𝑥 y Fx (t) = 𝑡q 𝑥 , si bien ya tenemos la tabla dada por la CUSF
para cada valor 1q 𝑥 para cualquier edad x y a partir de estas probabilidades se puede
construir 𝑡𝑝𝑥 y 𝑡q 𝑥 , los valores de t que podemos tener solamente son valores discretos (por
ejemplo, no tenemos forma de calcular 3.54𝑝𝑥 , la probabilidad de que una persona de edad
x sobreviva 3.54 años). Para el seguro continuo es necesario tener las funciones Sx (t) = 𝑡𝑝𝑥
y 𝜇x+t definida para toda t entre 0 y n.
𝑑
q
𝑑t t x
- Para la función 𝜇x+t = tenemos el cálculo de una derivada para la función Fx (t) =
p
t x
𝑡 q 𝑥 , la cual requerimos que sea continua también para poder hacer una aproximación de su
derivada en cada punto t y poder evaluar 𝜇x+t .
- La integral que calcula el seguro involucra el producto de 3 funciones 𝑡𝑝𝑥 , 𝜇x+t y (1 + i) -t,
con i constante, la cual es difícil de calcular exactamente considerando que las expresiones
𝑡 𝑝𝑥 , 𝜇x+t no tienen una forma exacta de definirse aún.
Desarrollo
Se busca resolver los 3 puntos planteados en la justificación de la siguiente manera:
- Dada la tabla de la CUSF con los valores de q x , x=0,1,2,…109,110, construimos una función
tal que reciba como parámetro una edad x y nos devuelva la tabla con todos los valores de
tpx , para t = 0,1,2,…110-x, calculados mediante la formula:
t-1
𝑡 p𝑥 = (1-q x )(1- q x+1 )(1- q x+2 ). . . . (1-q x+t-2 )(1- q x+t-1 ) = ∏(1-q x+k )
k=0
Y mediante esta tabla, fácilmente calculamos otra tabla para tq x como 1 - tpx .
- Una vez teniendo ambas tablas, las de Sx (t) = 𝑡𝑝𝑥 y Fx (t) = 𝑡q 𝑥 , obtendremos un total de
111-x puntos por tabla, a saber (0, 0𝑝𝑥 ), (1, 1𝑝𝑥 ), (2, 2𝑝𝑥 ), . . . , (109-x, 109-x𝑝𝑥 ), (110-x, 110-x𝑝𝑥 ).
Para crear las funciones continuas, realizamos una interpolación de estos puntos. Todavía
no defino que tipo de interpolación usar, pero tengo solo 2 opciones: o splines lineales o
Polinomio de Lagrange (revisaré cual es la más sencilla o incluso veré si puedo hacerlo por
ambos métodos).
- Teniendo una vez las funciones interpoladas Sx (t) = 𝑡𝑝𝑥 y Fx (t) = 𝑡q 𝑥 , mediante
𝑑
diferencias finitas calculamos la función q para cada punto t y una vez calculada esta
𝑑t t x
𝑑
q
𝑑t t x
función, obtenemos 𝜇x+t = para cada punto t.
p
t x
- Finalmente, una vez teniendo las funciones 𝑡𝑝𝑥 y 𝜇x+t , se realizará el cálculo de la integral
del seguro continuo como:
n n
̅ ́̅ = ∫ S. A. 𝑡p𝑥 𝜇x+t v t dt = ∫ S. A. 𝑡p𝑥 𝜇x+t (1 + i)-t dt
A x:n|
0 0
Por lo que se creará una función que reciba como parámetros la edad x, la suma asegurada
S.A., la temporalidad n y la tasa de interés i y calcule la integral por alguno de los métodos
de cuadratura numérica.
Resultados
Además del seguro continuo, se pretende realizar el cálculo por el método discreto (el cual
es más fácil de calcular pues solo requiere de las tablas). Por lo que esperamos obtener en
ambos cálculos, para una suma asegurada S.A., edad x, tasa de interés i y una temporalidad
n fijas, resultados similares tanto para la metodología continua como para la metodología
discreta.
