Calculating probabilities

Avaliação de Produtos de Seguro de Vida em R

Katrien Antonio, Ph.D.

Professor, KU Leuven and University of Amsterdam

De probabilidade anual para multianual de sobrevivência

 

 

Avaliação de Produtos de Seguro de Vida em R

De probabilidade anual para multianual de sobrevivência

 

 

Avaliação de Produtos de Seguro de Vida em R

De probabilidade anual para multianual de sobrevivência

 

 

Avaliação de Produtos de Seguro de Vida em R

De probabilidade anual para multianual de sobrevivência

 

 

Avaliação de Produtos de Seguro de Vida em R

A regra da multiplicação

  • Reescrevendo as probabilidades de sobrevivência:

$$ _{t+u}p_x =\ _{u}p_x \cdot \ _{t}p_{x+u}. $$

Avaliação de Produtos de Seguro de Vida em R

A regra da multiplicação

  • Reescrevendo as probabilidades de sobrevivência:

    $$ _{t+u}p_x =\ _{u}p_x \cdot \ _{t}p_{x+u}. $$

  • Com $k$ um inteiro, obtemos:

$$ \begin{aligned} _kp_x &= p_x \cdot \ _{k-1}p_{x+1} \\ &= p_x \cdot p_{x+1} \cdots p_{x+k-1} \\ &= \prod_{l=0}^{k-1} p_{x+l} \end{aligned} $$

$\quad \,$ que é um produto de probabilidades anuais de sobrevivência.

Avaliação de Produtos de Seguro de Vida em R

Calculando probabilidades de sobrevivência no R

Calcule $_5p_{65} = p_{65} \cdot p_{66} \cdot p_{67} \cdot p_{68} \cdot p_{69}$.

# One-year survival probabilities
px <- 1 - life_table$qx 
px[(65 + 1):(69 + 1)]
0.98491 0.98320 0.98295 0.98091 0.97935
# Probability that (65) survives 5 more years
prod(px[(65 + 1):(69 + 1)])
0.9144015

Calcule $_5p_{65} = \frac{\ell_{70}}{\ell_{65}}$.

# Alternatively (difference due to 
                 rounding)
lx[70 + 1] / lx[65 + 1]
0.9143957
Avaliação de Produtos de Seguro de Vida em R

Produto cumulativo das probabilidades de sobrevivência no R

Calcule $_kp_{65}$ para $k = 1, 2, 3, 4, 5$.

# One-year survival probabilities
px[(65 + 1):(69 + 1)]
0.98491 0.98320 0.98295 0.98091 0.97935
# Multi-year survival probabilities of (65)
cumprod(px[(65 + 1):(69 + 1)])
0.9849100 0.9683635 0.9518529 0.9336820
0.9144015

Calcule $_kp_{65}$ para $k = 0, 1, 2, 3, 4, 5$.

# Multi-year survival probabilities of (65)
c(1, cumprod(px[(65 + 1):(69 + 1)]))
1.0000000 0.9849100 0.9683635 0.9518529 
0.9336820 0.9144015
Avaliação de Produtos de Seguro de Vida em R

Uma probabilidade de mortalidade diferida

  • Foque em uma probabilidade de mortalidade diferida específica:

  • $(x)$ sobrevive $k$ anos completos, mas morre antes de atingir a idade $x + k + 1$:

$$ \begin{aligned} _{k|}q_x &= {}_kp_x \cdot q_{x+k}. \end{aligned} $$

Avaliação de Produtos de Seguro de Vida em R

Uma probabilidade de mortalidade diferida no R

Calcule $_{5|}q_{65} = {}_5p_{65} \cdot q_{70}$.

# 5-year deferred mortality probability of (65) 
prod(px[(65 + 1):(69 + 1)]) * qx[70 + 1]
0.02086664

Calcule $_{5|}q_{65} = \frac{d_{70}}{\ell_{65}}$.

# Alternatively (difference due to rounding)
dx[70 + 1] / lx[65 + 1]
0.02086817
Avaliação de Produtos de Seguro de Vida em R

Probabilidades de mortalidade diferida no R

Calcule $_{k|}q_{65} = {}_kp_{65} \cdot q_{65 + k}$ para $k = 0, 1, 2, \ldots$

# Survival probabilities of (65)
kpx <- c(1, cumprod(px[(65 + 1):(length(px) - 1)]))
head(kpx)
1.0000000 0.9849100 0.9683635 0.9518529 0.9336820 0.9144015
# Deferred mortality probabilities of (65)
kqx <- kpx * qx[(65 + 1):length(qx)]
head(kpx)
0.01509000 0.01654649 0.01651060 0.01817087 0.01928053 0.02086664
Avaliação de Produtos de Seguro de Vida em R

Vamos praticar!

Avaliação de Produtos de Seguro de Vida em R

Preparing Video For Download...