シンシアが100歳まで生きる確率はどのくらい?
R統計好きのシンシアは、1999年の生命表を使った最初の探索を楽しみました。次は、Rを使ってさらに深く分析してみましょう。シンシアが特に関心を持っているのは、100歳(センテナリアン)になる確率です。1999年のデータを使うと、0歳の人がこの確率はどのくらいでしょうか?18歳の場合はどうでしょうか?これらの確率の計算方法を確認したうえで、\((k,\ _kp_{18})\) を可視化しましょう。
以降の演習では、ベルギーの女性を対象とした1999年の期間生命表があらかじめ読み込まれ、life_table オブジェクトに格納されています。また、lx 列はすでに抽出済みです。
この演習はコースの一部です
R による生命保険商品の評価
演習の手順
lxを使って、0歳の人が100歳(センテナリアン)に達する確率を計算してください。インデックスで抽出するようにしてください。- 同じ確率を18歳の人について計算してください。
- 18歳の人が18歳、19歳、…100歳まで生存する確率を計算し、プロットしてください。
kを0から82までのベクトルとして定義し、kを使って確率を計算してください。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Compute the probabilty for (0) to reach the age 100
___
# Compute the probabilty for (18) to reach the age 100
___
# Plot the survival probabilties for (18) up to age 100
k <- ___
plot(k, ___,
pch = 20,
xlab = "k",
ylab = expression(paste(""[k], "p"[18])),
main = "Survival probabilities for (18)")