1年分の死亡率データの可視化

1999年の生命表についての理解が深まったところで、年齢 \(x\) に対する女性の死亡率 \(q_x\) の変化をグラフで表してみましょう。Cynthia の先生は、講義資料に死亡率の対数をプロットすることが多いそうです。これが有効なアプローチであることを、実際に確かめてみましょう。

前の演習で定義したオブジェクト life_table と、ベクトル age、qx はあらかじめ読み込まれています。

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

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

コースを見る

演習の手順

plot()（ドキュメント）を使って、1999年の女性の死亡率をプロットしてください。第1引数に age、第2引数に死亡率 qx を指定します。
もう一度 plot() を呼び出して、1999年の女性の死亡率を対数スケールで表示してください。前のコマンドとの違いは、第2引数に log()（ドキュメント）を適用して変換する点です。

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

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

# Plot the female mortality rates in the year 1999
plot(___, ___, 
    main = "Mortality rates (Belgium, females, 1999)", 
    xlab = "Age x", 
    ylab = expression(paste("Mortality rate ", q[x])), 
    type = "l")

# Plot the logarithm of the female mortality rates in the year 1999
plot(___, ___, 
    main = "Log mortality rates (Belgium, females, 1999)", 
    xlab = "Age x", 
    ylab = expression(paste("Log mortality rate ", log(q[x]))), 
    type = "l")

コードを編集して実行