不偏分散 | MatometAI語彙unbiased-varianceバージョン 3
不偏分散
標本から母分散を見積もるとき、ずれの二乗和を n ではなく n−1 で割ったもの。n は標本の大きさ、xi は 番目の値、 は標本平均。
i
s2=n−11i=1∑n(xi−xˉ)2 標本を取り直すたびに s2 の値は変わるが、取り直しを繰り返したときの平均（期待値）がちょうど母分散 σ2 に一致する。この「平均的に真の値に一致する」性質を不偏という。n で割る標本分散 σ^2 は、平均的に母分散より小さくなる。
証明
示すこと
不偏分散 s2 は取り直しの平均で母分散に一致し、標本分散 σ^2 は母分散より nn−1 倍だけ小さくなること。
E[s2]=σ2,E[σ^2]=nn−1σ2(1) 準備
標本を取り直すたびに値が変わることを表すため、標本を確率変数 X1,…,Xn（n≥2）で書き、標本平均を Xˉ=n1∑ とする。
前提：X1,…,Xn は i.i.d. で、E[Xi]=μ、E[(Xi−μ)2]=σ2。σ2 は有限の値とする。分布によっては分散が無限大になる（存在しない）ものがあるので、それを除くための条件。
s2 と σ^2 はどちらも、標本平均から測ったずれの二乗和 ∑i(Xi−Xˉ)2 を n−1 か n で割ったもの。そこで、このずれの二乗和を取り直しながら平均した値（期待値）を求めればよい。示すのは次の式。
E[i=1∑n(Xi−Xˉ)2]=(n−1)σ2(2) 左辺は真の分散ではなく、「標本平均から測ったずれの二乗和」の取り直しの平均。真の分散 σ2 の n 倍ではなく n−1 倍にしかならない、というのがこの式の意味で、両辺を n−1 で割れば E[s2]=σ2、n で割れば E[σ^2]=nn−1 になる。
補題 1（分解）
示すこと
i=1∑n(Xi−c)2=i=1∑n(Xi−Xˉ)2+n(Xˉ−c)2(3) 証明
Xi−c=(Xi−Xˉ)+(Xˉ−c) と分けて二乗すると、
(Xi−c)2=(Xi−Xˉ)2+2(Xi−Xˉ)(Xˉ−c)+(Xˉ−c)2(4) (4) を i=1,…,n について足す。Xˉ−c は i によらないので和の外に出せ、最後の項は同じものを n 回足すことになる。
i=1∑n(Xi−c)2=i=1∑n(Xi−Xˉ)2+2(Xˉ−c)i=1∑n(Xi−Xˉ)+n(Xˉ−c)2(5) (5) の真ん中の項の和は、Xˉ の定義から 0 になる。
i=1∑n(Xi−Xˉ)=i=1∑nXi−nXˉ=nXˉ−nXˉ=0(6) (6) を (5) に入れると真ん中の項が消え、(3) が成り立つ。(3) の右辺第 2 項は 0 以上なので、Xˉ はずれの二乗和を最小にする点でもある。
補題 2（標本平均のばらつき）
示すこと
標本平均 Xˉ の母平均 μ からのずれの二乗は、取り直しの平均で σ2/n になる。
E[(Xˉ−μ)2]=nσ2(7) 証明
Xˉ の定義から、Xˉ−μ は各データのずれ Xi−μ の平均で書ける。
Xˉ−μ=n1i=1∑nXi−μ=n1i=1∑n(Xi−μ)(8) (8) を二乗する。和の二乗は、すべての組 (i,j) の積の和になる。
(Xˉ−μ)2=n21(i=1∑n(Xn21i=1∑nj=1μ)(Xj−μ)(9) (9) の組を、i=j のものと i=j のものに分ける。
(Xˉ−μ)2=n21i=1(10) (10) の両辺の期待値をとる。期待値の線形性から、和の期待値は期待値の和になる。
E[(Xˉ−μ)2]=n21i=1(11) (11) の第 2 項を考える。i=j のとき Xi と Xj は独立なので、Xi−μ と Xj−μ も独立。独立な確率変数の積の期待値は期待値の積に分かれる（期待値を参照）ので、
E[(Xi−μ)(Xj−μ)]=E[Xi−μ]E[Xj−μ]=(μ−μ)(μ−μ)=0(12) (11) の第 1 項は、前提から各項が σ2 なので nσ2。これと (12) を (11) に入れると、
E[(Xˉ−μ)2]=n21(nσ2+0)=nσ2(13) 本題
(3) で c=μ とおいて移項する。
i=1∑n(Xi−Xˉ)2=i=1∑n(Xi−μ)2−n(Xˉ−μ)2(14) (14) の両辺の期待値をとり、右辺第 1 項に前提、第 2 項に (7) を使う。
E[i=1∑n(Xi−Xˉ)2]=i=1∑nE[(Xi−μ)2]−nE[(Xˉ−μ)2]=nσ2−n⋅nσ2=(n−1)σ2(15) (15) は (2) そのものなので、両辺を n−1 と n で割って (1) が成り立つ。
- 使った仮定は「期待値と分散が等しい」「i=j で E[(Xi−μ)(Xj−μ)]=0」の 2 つだけなので、独立でなくても無相関なら成り立つ
- 直感的には、n 個のデータのうち 1 個分の情報を xˉ を決めるのに使ったので、ばらつきを測るのに残るのは n−1 個分。これを自由度という
- 不偏は「取り直しの平均が真の値に一致する」という性質で、1 回ごとの見積もりが標本分散より真の値に近いことまでは意味しない
- s2 は不偏でも、その平方根 s は母標準偏差 σ の不偏な見積もりにはならない
例
母集団が「0 か 2 が確率 1/2 ずつ出る」（μ=1、σ2=1）とき、2 個取る標本は 4 通りで、どれも確率 1/4。
| 標本 | xˉ | 標本分散 σ^2 | 不偏分散 s2 |
|---|
| (0, 0) | 0 | 0 | 0 |
| (0, 2) | 1 | 1 | 2 |
| (2, 0) | 1 | 1 | 2 |
| (2, 2) | 2 | 0 | 0 |
| 4 通りの平均 | | 0.5 | 1 |
標本分散は平均すると 0.5 で真の値の半分、不偏分散は平均するとちょうど 1 になる。
i=1
n
Xi
σ2
i
−
μ
)
)
2
=
∑
n
(
Xi
−
∑
n
(
Xi
−
μ
)2
+
i=j∑
(
Xi
−
μ
)
(
Xj
−
μ
)
∑
n
E
[
(
Xi
−
μ
)2
]
+
i=j∑
E
[
(
Xi
−
μ
)
(
Xj
−
μ
)
]