初心者データサイエンティストの備忘録

調べたことは全部ここに書いて自分の辞書を作る

標準正規分布の期待値

はじめに

 確率変数 X \sim N(0, 1)について、 E\left[ X ^ n \right]を求める方法を知りたかったので、導出方法を少し整理しました。

導出の方針と導出結果

 結論だけを知りたい方のために、まず導出結果を書いておきます。


E\left[ X ^ n \right] = 
\begin{cases}
\dfrac{n!}{2 ^ {\frac{n}{2}} \left( \frac{n}{2}\right)!} & (n {\rm が偶数}) \\
0 & (n {\rm が奇数})
\end{cases} \tag{1}

 本記事では、(1)式をモーメント母関数とマクローリン展開を用いて導出します。

導出

 まず、一般的に e ^ sマクローリン展開


e ^ s = \displaystyle{\sum _{k=0} ^ \infty} \frac{s ^ k}{k!} \tag{2}

です。また、標準正規分布のモーメント母関数は


M(t) = e ^ {\frac{t ^ 2}{2}}

です。したがって、(2)式に s = \dfrac{t ^ 2}{2}を代入することで、標準正規分布のモーメント母関数のマクローリン展開は、


\begin{eqnarray}
  e ^ {\frac{t ^ 2}{2}} &=& \displaystyle{ \sum _{k=0} ^ \infty} \frac{1}{k!} \left( \frac{t ^ 2}{2} \right) ^ k \\
  &=& \displaystyle{ \sum _{k=0} ^ \infty} \frac{1}{2 ^ k k!} t ^ {2k}
\end{eqnarray} \tag{3}

となります。

 また、一般的にモーメント母関数のマクローリン展開


\begin{eqnarray*}
  M(t) = \displaystyle{ \sum _{n=0} ^ \infty} \frac{M ^ {(n)}(0)}{n!} t ^ n \tag{4}
\end{eqnarray*}

です。したがって、(3)・(4)式の係数を比較することで nが偶数のとき、 n = 2kとすることで、


\begin{eqnarray*}
  \dfrac{M ^ {(n)}(0)}{n!} = \dfrac{1}{2 ^ {n/2} (n/2)!}
\end{eqnarray*}

です。以上より


M ^ {(n)}(0) = \dfrac{n!}{2 ^ {n/2} (n/2)!}

です。同様に(3)・(4)式の係数を比較することで nが奇数のとき


\begin{eqnarray*}
  \dfrac{M ^ {(n)}(0)}{n!} = 0
\end{eqnarray*}

です。これと E\left[ X ^ n \right] = M ^ {(n)} (0)を組み合わせると、(1)式が得られます。

導出のポイント

 (3)式を導出する際に、(2)式を経由することで、標準正規分布のモーメント母関数のマクローリン展開が簡単にできます。私は、標準正規分布のモーメント母関数をそのままマクローリン展開しようとしたので、てこずりました。