コホート生存確率
シンシアが以前に計算した、1999年における18歳の人が23歳まで生存する確率は、実際には期間生存確率でした。生存確率の時間的な変化を考慮するために、計算を修正してみましょう。乗算則で使用する1年生存確率は、生命表から対角方向に取り出す必要があります。
$$ _5p_{18, 1999} = p_{18, 1999} \cdot p_{19,2000} \cdot p_{20,2001} \cdot p_{21,2002} \cdot p_{22,2003} $$
前の演習で使用した life_table オブジェクトはそのまま読み込まれています。
この演習はコースの一部です
R による生命保険商品の評価
演習の手順
subset()を使って、\(1999 - 18 = 1981\) 年生まれの女性のコホート生命表を作成し、表示してください。life_table_1981のqx列から、コホートの1年生存確率pxを定義してください。- (18) に対する5年コホート生存確率を計算してください。
- 1世紀前の1881年生まれについて、
with()とsubset()を両方使った1行のコマンドで同じ計算を繰り返せますか？
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Construct and print the cohort life table for birth year 1981
life_table_1981 <- ___(___, ___)
life_table_1981
# 1981 cohort one-year survival probabilities
px <- ___
# 1981 cohort survival probability that (18) survives 5 more years
___(px[(___):(___)])
# 1881 cohort survival probability that (18) survives 5 more years
___(___(life_table, year - age == ___), prod(1 - qx[(18 + 1):(22 + 1)]))