死力のソースを表示

{{Dablink|icon=Wiktfavicon en.svg|icon_size=16px|icon_alt=ウィクショナリー|1=<!--この項目では、について記述しています。-->その他の「'''死力'''」の語義については、[[ウィクショナリー]]の「[[wikt:死力|死力]]」の項目をご覧ください。}}
'''死力'''（しりょく、英: force of mortality）とは、[[保険数理]]で用いられる用語で、''X'' 歳に達した人が次の瞬間に死亡する確率を統計的に表している。[[自然人]]だけでなく、企業が倒産する確率や夫婦が離婚する確率なども死力と呼ぶことがある。

== 概要 ==
[[保険金]]の掛け金を推測するための重要な概念であり、保険においては'''死力'''を正確に求めることに膨大な労力を費やしている。一般的には年齢とともに上昇していくが、条件によって死亡する可能性が高くなる年代なども存在するため、必ずしも滑らかな曲線にはならない。

== 計算法 ==
[[生命表]]では、<math>x</math>歳になった人間が<math>x+1</math>歳までに死亡する率を、<math>q_x</math>で表現する。これを、ある瞬間の時間の死亡率として捉え直すため、<math>x</math>歳になった人間が<math>x+\Delta x</math>歳までに死亡する率を、下記のように条件付き確率 <math>P_x(\Delta x)</math> として表現する。

:<math>P_{x}(\Delta x)=P(x<X<x+\Delta\;x\mid\;X>x)=\frac{F_X(x+\Delta\;x)-F_X(x)}{(1-F_X(x))}</math>

ここで、<math>F_X(x)</math> は、死亡年齢を[[確率変数]]<math>X</math>で表すとき、<math>x</math>歳までに死亡する確率を示す[[累積分布関数]]である。

上記式を、勾配を求めるために<math>\Delta x</math>で除算し、<math>\Delta x</math> を 0 に近づけることによって、'''死力''' <math>\mu\,(x)</math> を得ることができる。

:<math>\mu\,(x)= \lim_{\Delta x \rightarrow 0} \frac{F_X(x+\Delta\;x)-F_X(x)}{\Delta x (1-F_X(x))} = \frac{F'_X(x)}{1-F_X(x)}</math>

ここで、<math>f_X(x) = F'_X(x)</math> とし、[[生存関数]]を <math>S(x) = 1-F_X(x)</math> とすると、死力は生存関数で以下のように表現できる。

:<math>\mu\,(x)=\frac{f_X(x)}{1-F_X(x)}=-\frac{S'(x)}{S(x)}=-{\frac{d}{dx}}\ln[S(x)]</math>

死力 <math>\mu\,(x)</math> は、死亡年齢の確率変数 <math>X</math>の条件付き[[確率密度関数]]である一方で、 <math>f_X(x)</math> は条件のない確率密度関数である<ref name="MQR">R. Cunningham, T. Herzog, R. London (2008). ''Models for Quantifying Risk, 3rd Edition'', Actex.</ref>。そのため、<math>x</math> 歳までの生存関数 <math>S(x)</math>が与えられたとき、<math>f_X(x)</math>は、条件付き確率<math>\mu(x)</math>と<math>S(x)</math>の積となるため、死力 <math>\mu(x)</math> は、

:<math>\mu(x) = \frac{f_X(x)}{S(x)}</math>

と表現される。

死力 <math>\mu</math> を、<math>x</math>から<math>x+t</math>まで積分すると、

:<math> \int_x^{x+t}\mu(y)\, dy\, = \int_{x}^{x+t} -\frac{d}{dy} \ln[S(y)]\, dy = -\ln[S(x + t)] + \ln[S(x)]</math>

となる。

<math>x</math>歳に到達した人間が、そこから<math>t</math>年生存する確率を、<math>S_x(t) = \frac{S(x+t)}{S(x)}</math> と定義し、上式の両辺の負の指数をとると、

:<math> S_x(t) = e^{-\int_x^{x+t}\mu(y)\, dy\,}</math>

となる。

== 脚注 ==
{{脚注ヘルプ}}
{{reflist}}

{{Applied-math-stub}}
{{DEFAULTSORT:しりよく}}
[[Category:保険]]
[[Category:保険数理]]
[[Category:数学に関する記事]]