Thí nghiệm nhị thức

Định giá Sản phẩm Bảo hiểm Nhân thọ bằng R

Roel Verbelen, Ph.D.

Statistician, Finity Consulting

Bảng sống trong R

  • life_table chứa bảng sống giai đoạn cho nam giới ở Bỉ năm 2013.
head(life_table, 10)

   age      qx     lx  dx    ex
1    0 0.00381 100000 381 77.95
2    1 0.00047  99619  47 77.24
3    2 0.00019  99572  19 76.28
4    3 0.00015  99553  15 75.30
5    4 0.00013  99538  13 74.31
6    5 0.00010  99525  10 73.32
7    6 0.00011  99514  11 72.32
8    7 0.00008  99504   8 71.33
9    8 0.00011  99496  11 70.34
10   9 0.00008  99485   8 69.34

# Các biến dùng trong video này
qx <- life_table$qx
px <- 1 - qx
lx <- life_table$lx
dx <- life_table$dx
Định giá Sản phẩm Bảo hiểm Nhân thọ bằng R

Thí nghiệm nhị thức: sống sót thêm một năm

  • Tập trung vào $\ell_x$ trong life_table.
lx[0 + 1]

1e+05
Định giá Sản phẩm Bảo hiểm Nhân thọ bằng R

Thí nghiệm nhị thức: sống sót thêm một năm

  • Số người sống đến tuổi $x+1$ tuân theo BIN($\ell_x$,$\ p_{x}$). 
lx[72 + 1]

73977

px[72 + 1]

0.97369

rbinom(n = 1, size = lx[72 + 1], prob = px[72 + 1])

72022
Định giá Sản phẩm Bảo hiểm Nhân thọ bằng R

Thí nghiệm nhị thức: sống sót thêm một năm

  • Giờ làm theo kiểu vector hóa!
sims <- rbinom(n = length(lx), size = lx, prob = px)
head(sims)

99637 99567 99553 99546 99525 99515
Định giá Sản phẩm Bảo hiểm Nhân thọ bằng R

Thí nghiệm nhị thức: sống sót $k$ năm

  • Số người sống sót 1 năm tuân theo BIN($\ell_x$,$\ p_{x}$).
    $\quad \,$ Kỳ vọng:
    $$ \ell_{x+1} = \ell_x \cdot p_x. $$

  • Số người sống sót $k$ năm tuân theo BIN($\ell_x$,$\ _{k}p_{x}$).

$\quad \,$ Kỳ vọng:

$$ \ell_{x+k} = \ell_{x} \cdot \ _{k}p_{x}. $$

$\quad \,$ Do đó:

$$ _{k}p_{x} = \frac{\ell_{x+k}}{\ell_{x}}. $$

Định giá Sản phẩm Bảo hiểm Nhân thọ bằng R

Thí nghiệm nhị thức: số người tử vong

  • Số người tử vong tuân theo BIN($\ell_x$,$\ q_{x}$).

$\quad \,$ Kỳ vọng: $$ \large \begin{aligned} d_x &= \ell_x \cdot q_x \\ &= \ell_x \cdot (1-p_x) \\ &= \ell_x -\ell_{x+1}. \end{aligned} $$ 

dx[72 + 1]

1946

lx[72 + 1] - lx[73 + 1]

1946
Định giá Sản phẩm Bảo hiểm Nhân thọ bằng R

Xác suất sống sót trong R

Tính $_5p_{65} = \dfrac{\ell_{70}}{\ell_{65}}$.

# Xác suất (65) sống thêm 5 năm
lx[age == 70] / lx[age == 65]

0.9143957

# Cách khác
lx[70 + 1] / lx[65 + 1]

0.9143957
Định giá Sản phẩm Bảo hiểm Nhân thọ bằng R

Minh họa xác suất sống sót trong R

# xác suất (65) sống đến tuổi 65 + k
k <- 0:45
plot(k, 
     lx[65 + k + 1] / lx[65 + 1], 
     pch = 20, 
     xlab = "k", 
     ylab = expression(paste(""[k], "p"[65])))

Định giá Sản phẩm Bảo hiểm Nhân thọ bằng R

Cùng luyện tập nào!

Định giá Sản phẩm Bảo hiểm Nhân thọ bằng R

Preparing Video For Download...