Động lực học

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

Katrien Antonio, Ph.D.

Professor, KU Leuven and University of Amsterdam

Nhiều dataframe với tỷ lệ tử vong

Xem $\{ q_x \} _{x\geq 0}$ hoặc $\{p _ x \} _{x \geq 0}$ qua nhiều giai đoạn thời gian:

Diễn biến tỷ lệ tử vong theo thời gian

$q_{x,t}$ là tỷ lệ tử vong của người $x$ tuổi trong giai đoạn $t$.

Cách tiếp cận theo giai đoạn (period)

Tính xác suất sống $_kp_{x,t}$ theo cách dọc:

$$ _kp_{x,t} = p_{x,t} \cdot p_{x+1,t} \cdot \ldots \cdot p_{x+k-1,t} $$

Cách tiếp cận theo thế hệ (cohort)

Tính xác suất sống $_kp_{x,t}$ theo cách chéo:

$$ _kp_{x,t} = p_{x,t} \cdot p_{x+1,t+1} \cdot \ldots \cdot p_{x+k-1,t+k-1} $$

Bảng sống theo thời gian trong R

head(life_table)

  year age      qx     lx    dx    ex
1 1841   0 0.16580 100000 16580 40.28
2 1841   1 0.07148  83420  5963 47.22
3 1841   2 0.03905  77457  3025 49.82
4 1841   3 0.02306  74432  1716 50.82
5 1841   4 0.01693  72716  1231 51.01
6 1841   5 0.01233  71485   881 50.88

tail(life_table)

      year age      qx lx dx   ex
19420 2015 105 0.49719 40 20 1.46
19421 2015 106 0.51410 20 10 1.40
19422 2015 107 0.52979 10  5 1.35
19423 2015 108 0.54425  5  3 1.31
19424 2015 109 0.55749  2  1 1.28
19425 2015 110 1.00000  1  1 1.26

Xác suất sống theo period của một người Bỉ nổi tiếng trong R

Jacques Brel là ca sĩ người Bỉ, sinh năm 1929 và mất ở tuổi 49.
Xác suất $_{49|}q_{0, 1929}$ rằng một trẻ sơ sinh, sinh năm 1929, sẽ mất ở tuổi 49 là bao nhiêu?
Cách period (dọc):

period_life_table <- subset(life_table, year == 1929)
qx <- period_life_table$qx
px <- 1 - qx
prod(px[1:(48 + 1)]) * qx[49 + 1]

0.008456378

Xác suất sống theo cohort của một người Bỉ nổi tiếng trong R

Jacques Brel là ca sĩ người Bỉ, sinh năm 1929 và mất ở tuổi 49.
Xác suất $_{49|}q_{0, 1929}$ rằng một trẻ sơ sinh, sinh năm 1929, sẽ mất ở tuổi 49 là bao nhiêu?
Cách cohort (chéo):

cohort_life_table <- subset(life_table, year - age == 1929)
qx <- cohort_life_table$qx
px <- 1 - qx
prod(px[(0 + 1):(48 + 1)]) * qx[49 + 1]

0.006410343

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

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