離散型確率変数のモーメント母関数
確率空間\(\left( \Omega ,\mathcal{F},P\right) \)に加えて離散型の確率変数\(X:\Omega \rightarrow \mathbb{R} \)が与えられているものとします。加えて、\(X\)の確率分布が確率質量関数\(f_{X}:\mathbb{R} \rightarrow \mathbb{R} \)によって記述されているものとします。つまり、確率変数\(X\)が値\(x\in \mathbb{R} \)をとる確率は、\begin{equation*}P\left( X=x\right) =f_{X}\left( x\right)
\end{equation*}であり、確率変数\(X\)の値が集合\(A\subset \mathbb{R} \)に属する確率は、\begin{equation*}P\left( X\in A\right) =\sum_{x\in A}f_{X}\left( x\right)
\end{equation*}であるということです。
実数\(t\in \mathbb{R} \)を任意に選んだ上で、変数\(x\)に関する関数\(e^{tx}:\mathbb{R} \rightarrow \mathbb{R} \)を定義します。この関数\(e^{tx}\)の定義域は\(\mathbb{R} \)であるため、これと先の確率変数\(X\)を用いることにより、それぞれの\(\omega \in \Omega \)に対して、\begin{equation*}e^{tX\left( \omega \right) }\in \mathbb{R} \end{equation*}を定める新たな確率関数\begin{equation*}
e^{tX}:\Omega \rightarrow \mathbb{R} \end{equation*}が定義可能です。つまり、\(e^{tX}\)は\(X\)と\(e^{tx}\)の合成関数として定義される確率変数です。この確率変数の期待値は、\begin{equation*}E\left( e^{tX}\right) =\sum_{x\in X\left( \Omega \right) }e^{tx}f_{X}\left(
x\right) \quad \because \text{LOTUS}
\end{equation*}であることに注意してください。\(t\)の値が変われば確率変数\(e^{tX}\)の形状も変化し、それに応じて期待値\(E\left(e^{tX}\right) \)もまた変化します。また、そもそも、確率変数\(e^{tX}\)の期待値は有限な実数として定まるとは限りません。ただ、\(t\)がとり得る値の範囲を点\(0\)を中心とする近傍\(N_{\varepsilon}\left( 0\right) =\left( -\varepsilon ,\varepsilon \right) \)に制限したとき、十分小さい半径\(\varepsilon >0\)のもとでは、この近傍に属する任意の値\(t\in N_{\varepsilon }\left(0\right) \)について期待値\(E\left(e^{tX}\right) \)が有限な値として定まることが保証されるのであれば、それぞれの\(t\in N_{\varepsilon }\left(0\right) \)に対して、\begin{equation*}M_{X}\left( t\right) =E\left( e^{tX}\right)
\end{equation*}を定める関数\begin{equation*}
M_{X}:\mathbb{R} \supset N_{\varepsilon }\left( 0\right) \rightarrow \mathbb{R} \end{equation*}が定義可能です。このような関数\(M_{X}\)が存在する場合、これを確率変数\(X\)のモーメント母関数(moment generating function of \(X\))や積率母関数などと呼びます。
\end{equation*}であるとともに、確率質量関数\(f_{X}:\mathbb{R} \rightarrow \mathbb{R} \)がそれぞれの\(x\in \mathbb{R} \)に対して定める値が、\begin{equation*}0\leq p\leq 1
\end{equation*}を満たす定数\(p\in \mathbb{R} \)を用いて、\begin{equation*}f_{X}\left( x\right) =\left\{
\begin{array}{cl}
p & \left( if\ x=1\right) \\
1-p & \left( if\ x=0\right) \\
0 & \left( otherwise\right)
\end{array}\right.
\end{equation*}と表されるものとします。\(t\in \mathbb{R} \)を任意に選んだとき、確率変数\(e^{tX}:\Omega \rightarrow \mathbb{R} \)の期待値は、\begin{eqnarray*}E\left( e^{tX}\right) &=&\sum_{x\in X\left( \Omega \right)
}e^{tx}f_{X}\left( x\right) \quad \because \text{LOTUS} \\
&=&e^{t\cdot 0}f_{X}\left( 0\right) +e^{t\cdot 1}f_{X}\left( 1\right) \\
&=&e^{0}\left( 1-p\right) +e^{t}p \\
&=&1-p+e^{t}p
\end{eqnarray*}です。任意の\(t\in \mathbb{R} \)について\(E\left( e^{tX}\right) \)は有限な実数として定まるためモーメント母関数\begin{equation*}M_{X}:\mathbb{R} \rightarrow \mathbb{R} \end{equation*}が存在して、これはそれぞれの\(t\in \mathbb{R} \)に対して、\begin{equation*}M_{X}\left( t\right) =1-p+e^{t}p
\end{equation*}を定めます。
モーメント母関数とモーメントの関係
モーメント母関数はその名の通りモーメントと深い関係があります。モーメントについて簡単に復習した上で、モーメント母関数との関係を以下で解説します。
離散型の確率変数\(X:\Omega\rightarrow \mathbb{R} \)の確率分布が確率質量関数\(f_{X}:\mathbb{R} \rightarrow \mathbb{R} \)として記述されているものとします。自然数\(m\in \mathbb{N} \)を任意に選んだとき、それぞれの\(\omega \in \Omega \)に対して、\begin{equation*}X^{m}\left( \omega \right) =\left[ X\left( \omega \right) \right] ^{m}
\end{equation*}を定める新たな確率変数\(X^{m}:\Omega \rightarrow \mathbb{R} \)が定義可能です。この確率変数\(X^{m}\)の期待値は、\begin{equation*}E\left( X^{m}\right) =\sum_{x\in X\left( \Omega \right) }x^{m}f_{X}\left(
x\right) \quad \because \text{LOTUS}
\end{equation*}となりますが、これを確率変数\(X\)の\(m\)次のモーメントと呼びます。
確率変数\(X\)のモーメント母関数\(M_{X}:\mathbb{R} \supset N_{\varepsilon }\left( 0\right) \rightarrow \mathbb{R} \)が存在する場合、任意の自然数\(m\)について、確率変数\(X\)の\(m\)次のモーメント\(E\left( X^{m}\right) \)が有限な実数として定まることが保証されるとともに、それはモーメント母関数\(M_{X}\)の点\(0\)における\(m\)次の微分係数と一致することが保証されます。つまり、以下の関係\begin{equation*}E\left( X^{m}\right) =\left. \frac{d^{m}M_{X}\left( t\right) }{dt^{m}}\right\vert _{t=0}
\end{equation*}が成り立つということです。
\end{equation*}が成り立つ。
\end{equation*}であるとともに、確率質量関数\(f_{X}:\mathbb{R} \rightarrow \mathbb{R} \)がそれぞれの\(x\in \mathbb{R} \)に対して定める値が、\begin{equation*}0\leq p\leq 1
\end{equation*}を満たす定数\(p\in \mathbb{R} \)を用いて、\begin{equation*}f_{X}\left( x\right) =\left\{
\begin{array}{cl}
p & \left( if\ x=1\right) \\
1-p & \left( if\ x=0\right) \\
0 & \left( otherwise\right)
\end{array}\right.
\end{equation*}と表されるものとします。先に示したように、\(X\)のモーメント母関数\(M_{X}:\mathbb{R} \rightarrow \mathbb{R} \)が存在して、これはそれぞれの\(t\in \mathbb{R} \)に対して、\begin{equation*}M_{X}\left( t\right) =1-p+e^{t}p
\end{equation*}を定めます。したがって、先の命題より、自然数\(m\in \mathbb{N} \)を任意に選んだとき、\(X\)の\(m\)次のモーメント\(E\left( X^{m}\right) \)が有限な実数として定まるとともに、以下の関係\begin{eqnarray*}E\left( X^{m}\right) &=&\left. \frac{d^{m}M_{X}\left( t\right) }{dt^{m}}\right\vert _{t=0} \\
&=&\left. \frac{d^{m}}{dt^{m}}\left( 1-p+e^{t}p\right) \right\vert _{t=0} \\
&=&\left. e^{t}p\right\vert _{t=0} \\
&=&e^{0}p \\
&=&p
\end{eqnarray*}が成立します。他方で、定義にもとづいてモーメントを求めると、\begin{eqnarray*}
E\left( X^{m}\right) &=&\sum_{x\in X\left( \Omega \right) }x^{m}f_{X}\left(
x\right) \quad \because \text{LOTUS} \\
&=&0^{m}f_{X}\left( 0\right) +1^{m}f_{X}\left( 1\right) \\
&=&0\left( 1-p\right) +1\left( p\right) \\
&=&p
\end{eqnarray*}となりますが、これは先の結果と整合的です。
\end{equation*}が成り立ちます。つまり、\(X\)の原点まわりの\(1\)次のモーメントは\(X\)の期待値と一致します。さらに、先の命題より、\begin{equation*}E\left( X^{1}\right) =\left. \frac{dM_{X}\left( t\right) }{dt}\right\vert
_{t=0}
\end{equation*}が成り立つため、\begin{equation*}
E\left( X\right) =\left. \frac{dM_{X}\left( t\right) }{dt}\right\vert _{t=0}
\end{equation*}を得ます。モーメント母関数から期待値を導出できるということです。
\end{equation*}が成り立ちます。また、確率変数\(X\)の期待値とモーメントの間には以下の関係\begin{equation*}E\left( X\right) =E\left( X^{1}\right)
\end{equation*}が成り立ちます。加えて、先の命題より、\begin{eqnarray*}
E\left( X^{1}\right) &=&\left. \frac{d^{2}M_{X}\left( t\right) }{dt}\right\vert _{t=0} \\
E\left( X^{2}\right) &=&\left. \frac{dM_{X}\left( t\right) }{dt^{2}}\right\vert _{t=0}
\end{eqnarray*}を得るため、\begin{equation*}
\mathrm{Var}\left( X\right) =\left. \frac{dM_{X}\left( t\right) }{dt^{2}}\right\vert _{t=0}-\left[ \left. \frac{d^{2}M_{X}\left( t\right) }{dt}\right\vert _{t=0}\right] ^{2}
\end{equation*}を得ます。モーメント母関数から分散を導出できるということです。
モーメント母関数による確率分布の特徴づけ
確率変数の確率分布はモーメント母関数によって特徴づけられます。つまり、2つの確率変数が同一のモーメント母関数を持つことと、それらの確率変数が同一の確率分布にしたがうことは必要十分です。順番に解説します。
2つの離散型の確率変数\(X,Y:\Omega \rightarrow \mathbb{R} \)の確率分布が確率質量関数\(f_{X},f_{Y}:\mathbb{R} \rightarrow \mathbb{R} \)としてそれぞれ表現されているものとします。加えて、これらの確率変数\(X,Y\)がともにモーメント母関数\(M_{X},M_{Y}:\mathbb{R} \supset N_{\varepsilon }\left( 0\right) \rightarrow \mathbb{R} \)を持つものとします。このとき、\(f_{X}\)と\(f_{Y}\)が一致することと\(M_{X}\)と\(M_{Y}\)が一致することは必要十分になります。
&&\left( b\right) \ \forall t\in N_{\varepsilon }\left( 0\right)
:M_{X}\left( t\right) =M_{Y}\left( t\right)
\end{eqnarray*}は必要十分である。
離散型確率変数の確率分布は分布関数によって表現することもできるため、上の命題より以下を得ます。
&&\left( b\right) \ \forall t\in N_{\varepsilon }\left( 0\right)
:M_{X}\left( t\right) =M_{Y}\left( t\right)
\end{eqnarray*}は必要十分である。
多くの場合、2つの確率変数の確率分布が一致することを示すことよりも、それらのモーメント母関数が一致することを示すことの方が容易です。上の命題は、2つの確率変数が同一の確率分布にしたがうことを示す代わりに、それらのモーメント母関数が一致することを示してもよいことを保証します。
演習問題
\end{equation*}であるとともに、確率質量関数\(f_{X}:\mathbb{R} \rightarrow \mathbb{R} \)がそれぞれの\(x\in \mathbb{R} \)に対して定める値が、\begin{equation*}f_{X}\left( x\right) =\left\{
\begin{array}{cl}
\frac{1}{3} & \left( if\ x=1\right) \\
\frac{2}{3} & \left( if\ x=2\right) \\
0 & \left( otherwise\right)
\end{array}\right.
\end{equation*}を定めるものとします。モーメント母関数\(M_{X}\)は存在するでしょうか。存在する場合には具体的に求めてください。
プレミアム会員専用コンテンツです
【ログイン】【会員登録】