logo

標本分散をn-1で割る理由 📂数理統計学

標本分散をn-1で割る理由

なぜ n-1で割るのか?

Xi(μ,σ2)X_{i} \sim \left( \mu , \sigma^{2} \right) とすると、標本分散 S2S^{2} は次のようになる。 S2:=1n1i=1n(XiX)2 S^{2} := {{1} \over {n-1}} \sum_{i=1}^{n} \left( X_{i} - \overline{X} \right)^{2}

ご存知の通り、標本平均と異なり、標本分散では偏差の平方をすべて足した後、標本サイズの nn ではなく、n1n-1 で割る。これがおかしいとは言わないが、式に対する一般的な感覚があるなら、nn 個を足して n1n-1 で割ることに強い嫌悪感を感じるのも普通だ。これについて、高校生、統計学部の新入生、統計学部の卒業生のレベルで簡単に説明してみよう。

ちなみに、統計学部の卒業生レベルの説明を除いて、残りは数学的な虚構が混じっているので、納得できなかったり、おかしなところがあっても問題ない。また、統計学の分野によっては、目的に応じてn1n-1 ではなくnn で標本分散を求める場合もあるため、正しいか間違っているかを議論するよりは、正確な理由を理解することが大切だ。

高校生レベルの説明

n=1n=1 のとき、分母が 00 になって、定義できなくするためだ。基本的に、nnn1n-111 しか違わないから、標本サイズが大きくなれば、nn で割ろうが n1n-1 で割ろうが、大した違いはない。問題は標本のサイズが小さいときで、極端な話、n=1n=1 だと、標本が一つしかない場合になる。しかし、分散はデータの分布を話しており、一つしかないデータがどれだけ散らばっているかは、そもそも定義できないというのが常識だ。

ただ一つのデータ x1x_{1} が与えられた場合、その標本平均も x=x1\overline{x} = x_{1} となり、標本分散は s2=0s^{2} = 0 となって、「全く散らばっていない」という意味になるかもしれない。しかし、これは標本サイズが大きくなっても同じだ。すべてのデータが正確に同じであれば、x=x1==xn\overline{x} = x_{1} = \cdots = x_{n} となり、結局標本分散は s2=0s^{2} = 0 となることができる。数値的に全く散らばっていないことと、実際に散らばることができないことは、この感覚で本質的に異なる。標本分散を求めるときに n1n-1 で割ると、このような長い話が短い式に要約される。

そもそも分散が「偏差の平方の平均」として定義されたことはなく、期待値という表現を使わなければ、実際はnn で割る理由も特になく、正式な定義に従うなら「偏差の平方の期待値」となれば、n+1n+1n7n-7 などのよく分からない数で割られても構わない。ただ、上述の説明のように「一つの標本に対して分散が定義されるのはおかしい」と主張したいなら、その数は正確に n1n-1 であるべきだ。

統計学部新入生のレベルの説明

nn は標本サイズに過ぎず、正確には n1n-1自由度だ。一見、XiX_{i}nn 回足すから、データも nn 種類を使っているように見えるが、x=i=1nxi/n\overline{x} = \sum_{i=1}^{n} x_{i} / n が与えられた時点で、x1,,xn1x_{1} , \cdots , x_{n-1} を知っていれば xn=n(xi=1n1xi/n) x_{n} = n \left( \overline{x} - \sum_{i=1}^{n-1} x_{i} / n \right) のように逆算して定数 xnx_{n} を特定できる。つまり、式の見た目とは異なり、実際の分散の計算で使うデータの数は nn 個ではなく、そこから一つが失われた n1n-1 個が正しい。これを自由度といい、この文脈では「実際に使う標本サイズ」と受け取ってもよい。だから、標本分散を求めるときは、意味のない標本サイズ nn ではなく、自由度 n1n-1 で割らなければならない。

統計学部の卒業生レベルの説明

簡単に言えば、標本分散の期待値が母分散になるようにするためだ。難しく言えば、標本分散は n1n-1 で割ることによって不偏推定量になる。不偏推定量は偏りのない推定量で、実際に式を展開してみると、n1n-1 で割ることが X\overline{X} を扱う過程で生じる問題を適切に解決してくれる。簡単にまとめられた式は以下の通り。 ES2=(n1)1Ei=1n(XiX)2=(n1)1[i=1n(σ2+μ2)n(μ2+σ2n)]=(n1)1(n1)σ2=σ2 \begin{align*} E S^{2} =& (n-1)^{-1} E \sum_{i=1}^{n} \left( X_{i} - \overline{X} \right)^{2} \\ =& (n-1)^{-1} \left[ \sum_{i=1}^{n} \left( \sigma^{2} + \mu^{2} \right) - n \left( \mu^{2} + {{\sigma^{2}} \over {n}} \right) \right] \\ =& (n-1)^{-1} (n-1) \sigma^{2} \\ =& \sigma^{2} \end{align*} この式の展開は非常に端折られているので、正しく理解したい場合は、詳細な証明を参照することをお勧めする。このような数理的な議論によれば、標本分散を求めるときに n1n-1 で割るのは、nn で割るときの改善策ではなく、必然としてそうすべきだという正当性を持っている。計算した標本分散が実際の母分散と期待されないなら、それを標本分散と呼ぶ理由がない。

もし、卒業生レベルまで納得した後に、高校生、新入生レベルの説明を再び見れば、説明が非常に不十分であることが分かる。高校生の立場から、標本が一つの時に分散が定義されないようにn1n-1で割るという説明は、一つの標本で定義されなければならないという前提に賛同する必要があり、新入生が自由度n1n-1であることを理解したとしても、なぜ自由度で割る必要があるのかについての明確な説明がない。これは単に、できるだけ式に頼らずに、聞く人が感覚的に理解できるように説明する方法に過ぎない。

逆に言えば、自分自身が統計学を勉強して理解したと主張するには、少なくともそのような式を扱うだけの数学的な素養を備える必要がある。自分のレベルに満足していて、十分だと感じるなら、自分の主要な分野に集中しても良いが、何か物足りなければ、時間と努力を投資して、式を自然に理解できるようにしよう。