正規分布の分散を計算する の解説
2018/01/22
カテゴリ:数理統計
タグ:
概要
ブログ記事「正規分布の分散を計算する」の解説記事です。
次の式で表される正規分布の分散(Vとします)がとなることについて、その計算の過程を補ったものです。
![Rendered by QuickLaTeX.com \displaystyle V= \int_{-\infty }^{\infty} (x-\mu)^2 \frac{1}{\sqrt{2\pi \sigma^2}} \exp \left(-\frac{(x-\mu)^2}{2\sigma^2} \right) dx = \sigma^2](https://bellcurve.jp/statistics/wp-body/wp-content/ql-cache/quicklatex.com-00e7992710c325f7ef056346c6e6d5dc_l3.png)
なお、数学的な記述についてはあまり込み入った解説はせず、高校卒業程度の知識で理解できる範囲に留めます。
ポイント
正規分布の分散を計算する上で、式変形のポイントが二つあります。
- 簡単な式に変形する
- 部分積分で計算する
置換積分により、を
に変形します。
1.で変形した式に部分積分を適用し、が
となることを計算します。
簡単な式に変形する
正規分布の分散を求めるには、次の式の計算が必要です。
(1)
まず変数変換を行うのですが、正規分布の指数部分が複雑な形をしています。
このまま計算するのは煩雑です。そこで、次のような変換を行いについての積分を考えます。
(2)
次の3つの点に注意しながら、元の積分を書き換えます。
- 積分区間
- 被積分関数
- dxとdyの関係
元々の積分区間はがが
から
までを考えていたのですが、
という変換によって、
がどのように変化するか考えます。
![Rendered by QuickLaTeX.com \begin{array}{c|ccc} x & -\infty & \cdots & \infty \\ \hline \downarrow & \downarrow & & \downarrow \\ \hline y & -\infty & \cdots & \infty \\ \end{array}](https://bellcurve.jp/statistics/wp-body/wp-content/ql-cache/quicklatex.com-dbdeb3ce9a4eda9ab4c07a27491063e4_l3.png)
も
から
まで値を取ることが分かります。つまり、積分区間は
です。
被積分関数について、を含む部分を全て
に置き換えます。
最後にdxです。(2)式をについての式に変形します。
![Rendered by QuickLaTeX.com \displaystyle y=\frac{x-\mu}{\sigma} \longleftrightarrow x=\sigma y + \mu](https://bellcurve.jp/statistics/wp-body/wp-content/ql-cache/quicklatex.com-dad4f8aa50f9d3791f3aad95dc6a9b51_l3.png)
これを、yについて微分します。
![Rendered by QuickLaTeX.com \displaystyle \frac{dx}{dy}=\frac{d}{dy} \left( \sigma y + \mu \right) = \sigma](https://bellcurve.jp/statistics/wp-body/wp-content/ql-cache/quicklatex.com-2369dff94b2dbfcd87b7f7fbe4e97ddb_l3.png)
これを変形するとと書き換えられます。
以上のことから、(1)式の積分は次のように書き換えられます。
(3)
このような変数変換による積分の書き換えを置換積分と呼びます。 現在の学習指導要領では高校の数学IIIの範囲です。
部分積分で計算する
次に、部分積分による計算です。
部分積分は次の式で表される積分公式の一つです。や
はそれぞれの関数を微分したものを表します。
![Rendered by QuickLaTeX.com \displaystyle \int_a^b f(x)g(x)' dx = \left[ f(x)g(x) \right]_a^b - \int_a^b f(x)'g(x)dx](https://bellcurve.jp/statistics/wp-body/wp-content/ql-cache/quicklatex.com-f11dc152d2c8861973082bc6c1b23ab5_l3.png)
まず、(3)式の被積分関数を二つの関数の積
で表せるか考えます。
やや唐突ですが、について次の関係が成り立ちます。
![Rendered by QuickLaTeX.com \displaystyle \frac{d}{dx}\exp \left(-\frac{x^2}{2} \right)= -x \exp \left(-\frac{x^2}{2} \right) \\ \longleftrightarrow \left( - \exp \left(-\frac{x^2}{2} \right) \right)' = x \exp \left(-\frac{x^2}{2} \right)](https://bellcurve.jp/statistics/wp-body/wp-content/ql-cache/quicklatex.com-958936714b07d704fac4577d79e91418_l3.png)
これを用いると、の形をした関数を積分することができます。
の積分
とおき、部分積分によって変形しましょう。
最後の式の第1項が分かればと
は何らかの数列で表すことができそうです。
ここで、第1項は次のような極限となります。
![Rendered by QuickLaTeX.com \begin{eqnarray*} \displaystyle \lefteqn{ \left[ -x^{n-1} \exp \left(-\frac{x^2}{2} \right) \right]_{-\infty }^{\infty } } \\ &=&-\lim_{x \to \infty} x^{n-1} \exp \left( -\frac{x^2}{2} \right) + \lim_{ x \to -\infty} x^{n-1} \exp \left( -\frac{x^2}{2} \right) \\ &=& 0](https://bellcurve.jp/statistics/wp-body/wp-content/ql-cache/quicklatex.com-9b6b6dc6de67b24d058340bfad399fb9_l3.png)
最後の変形には、であることを用いました。
つまり、
![Rendered by QuickLaTeX.com I_n = (n-1) I_{n-2}](https://bellcurve.jp/statistics/wp-body/wp-content/ql-cache/quicklatex.com-561f93176a2fd79eeafa535c5d6060de_l3.png)
となります。
以上のことから、式(3)は次のように変形できます。
最後にの値です。
ガウス積分でも計算できますが、確率分布についての知識があれば計算は不要です。
について、具体的にどのような式になっているのでしょうか。
よく見ると、この式の被積分関数は標準正規分布の確率密度関数と全く同じであることが分かります。
ここで、確率分布の性質を思い出してみましょう。
統計学の時間でも解説していますが、確率密度関数はその定義から、全域で積分すると1となることが分かっています。
つまりです。これを(4)に代入します。
これで、正規分布の分散がと求められました。