年齢別の平均余命のプロット

シンシアは、年齢の関数として平均余命を可視化することにしました。このグラフを作成する方法を一緒に考えましょう。

R で効率的に行うには、まず指定した年齢と生命表から簡略平均余命を計算する関数を作成します。その後、生命表のすべての年齢にこの関数を適用します。sapply() （ドキュメント）を使うと、出力を最もシンプルなデータ構造に変換できます。

あらかじめ読み込まれている life_table オブジェクトには、ベルギーの女性に関する 1999 年の期間生命表が含まれています。

この演習はコースの一部です

R による生命保険商品の評価

コースを見る

演習の手順

curtate_future_lifetime() 関数を定義するコードを完成させましょう。
life_table の age 列を抽出して、ベクトル ages を作成します。
sapply() に ages、curtate_future_lifetime、および life_table を引数として指定し、生命表のすべての年齢における簡略平均余命を計算します。
ages に対して future_lifetimes をプロットします。

実践的なインタラクティブ演習

このサンプルコードを完成させて、この演習に挑戦してみましょう。

# Function to compute the curtate expected future lifetime for a given age and life table
curtate_future_lifetime <- function(age, life_table) {
  px <- ___
  kpx <- ___(px[(___):length(px)])
  ___(___)
}

# Vector of ages
ages <- ___

# Curtate future lifetimes for all ages
future_lifetimes <- ___(___, ___, ___)

# Future lifetime by age
plot(___, ___, type = 'l', lwd = 2, col = "green", xlab = "Age x", ylab = "Future lifetime", main = "Future lifetime by age")

コードを編集して実行