生存関数から期待値を求める方法を解説!!!

今回は、生存関数を用いて確率変数の期待値を求める方法について詳しく解説します。
数学的な導出だけでなく、直感的な理解や具体的な例も交えて説明していきます。

定義

確率変数\(X\)に対して、
\begin{align*} S(t) = P( t \leq X)\end{align*}
を\(X\)の生存関数といいます。

期待値は、生存関数の積分によって表現することができます。
どういうことかというと、

定理

\(X\)を確率変数とする。\(S\)を\(X\)の生存関数とすると、

\begin{align*} E(X) = \int_0^\infty S(t) dt \end{align*}

今回は、\(X\)が連続型で、確率密度関数\(f\)をもつ場合に限って話を進めていきます。
任意の\(n\)、部分集合\(A \subset \mathbb R^n\)に対して、\(\chi_A\)で, \(A\)上の定義関数を表すことにします。

任意の\((x, t) \in \mathbb R^2\)に対して、
\begin{align*} \chi_{[t, \infty)}(x) = \chi_{[0, x]}(t)\end{align*}
であることをリマインドしておきます(念の為ですが、上記式は関数として一致というわけではなく実数値としていつでも一致するということです)。


\begin{align*} \int_0^\infty S(t) dt &= \int_0^\infty \left(\int_t^\infty f(x) dx \right) dt
\\&= \int_0^\infty \left(\int_0^\infty \chi_{[t, \infty)}(x) f(x) dx \right) dt
\\&= \int_0^\infty \left(\int_0^\infty \chi_{[0, x]}(t) f(x) dx \right) dt
\\&= \int_0^\infty \left(\int_0^\infty \chi_{[0, x]}(t) f(x) dt \right) dx \quad (積分の順序を交換)
\\&= \int_0^\infty \left(\int_0^\infty \chi_{[0, x]}(t) dt \right) f(x)dx
\\&= \int_0^\infty \left(\int_0^x 1 dt \right) f(x)dx
\\&= \int_0^\infty x f(x)dx
\\&= E(X) \end{align*}


途中で行った積分の順序交換が正当化できるかどうかを確認しておきます。
\(\mathbb R\)にはボレル集合体\(\mathcal B (\mathbb R)\)を備えて(ボレル)可測空間\((\mathbb R, \mathcal B(\mathbb R))\)とします。
とりあえず気にせず積分の順序交換をしてやるぜーって人は気にしなくて大丈夫です。
\begin{align*} \chi_{[0, x]}(t) = \chi_{\{ (x, t) \in \mathbb R^2 \mid 0 \leq t \leq x \}} (x, t)\end{align*}
ですが、
\begin{align*} \{ (x, t) \in \mathbb R^2 \mid 0 \leq t \leq x \} \end{align*}
がボレル可測集合なので、
\begin{align*}\chi_{\{ (x, t) \in \mathbb R^2 \mid 0 \leq t \leq x \}} \end{align*}
について、任意の\(\mathbb R\)のボレル可測集合の逆像が\(\{ (x, t) \in \mathbb R^2 \mid 0 \leq t \leq x \}\)か\(\varnothing\)か\(\mathbb R\)なので、
\begin{align*}\chi_{\{ (x, t) \in \mathbb R^2 \mid 0 \leq t \leq x \}} \end{align*}
は\(\mathcal B(\mathbb R^2)\)可測です。なので、
\begin{align*} \mathcal B(\mathbb R) \otimes \mathcal B(\mathbb R) \end{align*}
可測です(ここは補足記事を用意してあります)。

ということで、フビニの定理が適用できて積分の順序交換ができます。
(最近数学をやっていないので細かいところで議論に致命傷があるかもしれません、コメント等で補足してもらえると助かります)

いずれにせよ、期待値は生存関数の積分と一致します。

記事をシェアして話のネタにする

コメント

コメントする

目次