正規分布の分散を計算する

2017/12/31

カテゴリ：数理統計

概要

正規分布 $N(\mu,\sigma^2)$ の分散は $\sigma^2$ です。

公式のように覚えてしまうのもよいのですが、実際に計算して $\sigma^2$ であることを確かめましょう。

ただし本記事では、意図的に式展開や計算の根拠などの行間を省いている部分があります。練習問題のように実際に計算して、行間を埋めてください。行間は、来月の記事で解説します。

この記事の目安

統計検定でのおおよその目安

統計検定2級の内容をある程度理解しており、統計検定準１級や統計検定１級に向けて勉強したい人

数学的な内容の目安

理系で高校を卒業した人、または次の内容を簡単に覚えている人

合成関数の微分
指数関数の微分
部分積分
積分の変数変換（置換積分）

計算

$N(\mu,\sigma^2)$ の確率密度関数は次式で表されます。

$\displaystyle f(x)= \frac{1}{\sqrt{2\pi \sigma^2}} \exp \left(-\frac{(x- \mu)^2} {2\sigma^2} \right)$

平均は $\mu$ であるため、分散 $V$ は次の式を計算すればよいことがわかります。

$\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$

一見、とても難しそうです。しかしこの計算は、変形によって簡単に計算することができます。

最初に、変数変換を行います。 $\frac{x-\mu}{\sigma}=y$ とおきましょう。

このとき、 $\frac{dy}{dx}=\frac{1}{\sigma}$ となることから次のように変形できます。 $y$ の積分区間についても忘れないようにしましょう。

$\displaystyle V = \sigma^2 \int_{-\infty }^{\infty} y^2 \frac{1}{\sqrt{2\pi}} \exp \left(-\frac{y^2} {2} \right) dy$

直接計算するにはまだ少し難しいので、部分積分を行います。どのように変形すればよいかを検討するため、 $\exp \left(-\frac{y^2}{2} \right)$ を微分しましょう。

$\displaystyle \frac{d}{dy} \exp \left(-\frac{y^2}{2} \right)= -y \exp \left(-\frac{y^2}{2} \right)$

これを利用すると、与式は次のように書き換えることができます。

$\displaystyle \sigma^2 \int_{-\infty }^{\infty} y \left(- \frac{1}{\sqrt{2\pi}} \exp \left(-\frac{y^2} {2} \right) \right)' dy$

部分積分を行うと、次のように計算できます。

$\begin{eqnarray*} \displaystyle \lefteqn{\sigma^2 \int_{-\infty }^{\infty} y \left(- \frac{1}{\sqrt{2\pi}} \exp \left(-\frac{y^2} {2} \right) \right)' dy} \\ &=& \sigma^2 \left[ -y \frac{1}{\sqrt{2\pi}} \exp \left(-\frac{y^2} {2} \right) \right]_{-\infty }^{\infty} + \sigma^2 \int_{-\infty }^{\infty} \frac{1}{\sqrt{2\pi}} \exp \left(-\frac{y^2} {2} \right) dy \\ &=& \sigma^2 \times 0 + \sigma^2 \times 1 \\ &=& \sigma^2 \end{eqnarray*}$

以上のことから、正規分布の分散が $\sigma^2$ と求められました。

ブログ