不偏推定量とは|定義・公式とアクチュアリー試験の関連問題 | acpass統計・用語・点推定
ひとことで言うと
推定量を何度も繰り返し計算したときの「平均的なズレ」がゼロであることです。1回の標本では真の値からズレることがあっても、長期的に見れば偏りなく真の値の周りに散らばる推定量です。
こんなデータが従う
標本平均は母平均の不偏推定量標本分散をn-1で割ると母分散の不偏推定量になる回帰係数の最小二乗推定量比率pの標本比率p̂最尤推定量は一般に不偏とは限らない例
推定量を選ぶときの基本的な良さの基準の一つです。不偏性だけでなく分散の小ささ(有効性)も合わせて評価するのが実務的です。
真の値θに分布の中心が一致する推定量(不偏、実線)と、中心がズレている推定量(偏りあり、破線)の比較。
数式で表すと
E[θ^]=θ 期待値が真値に等しい推定量。標本分散は n−1 で割ると不偏になる。
推定量 θ^ が不偏であるとは、E[θ^]=θ(真の値)が常に成り立つことです。偏り(バイアス)は Bias(θ^)=E[θ^]−θ
試験に出る性質
不偏性の定義
E[θ^]=θ
バイアス
Bias(θ^)=E[θ^]−θ。
例で見る
n=5の標本で ∑(xi−xˉ)2=80 のとき、母分散の不偏推定量は u2=80/(5−1)=20。
つまずきポイント
- 母分散の不偏推定量を求めるのに n で割ってしまう(n 割りは標本分散 s2 で有偏。不偏なのは n−1 割りの u2)
定着クイズ
Q1n=6の標本で偏差二乗和が50のとき、母分散の不偏推定量u2は?
Q2推定量 θ^ が不偏であるとはどういう意味か?
Q3標本分散をnで割った場合、母分散と比べてどうなるか?
Q4U(0,θ) の X(n) の期待値は?
Q5X(n) を不偏にする補正は?
Q6n=5,θ=1 でこの補正がMSEに与える効果は?
と定義され、不偏推定量はバイアスがゼロの推定量にあたります。
代表例が不偏分散です。標本平均からの偏差の二乗和を単純に n で割った標本分散 s2 は母分散を少し過小評価してしまいます(標本平均自体がデータから推定された値であるため、自由度が1つ減るのが原因)。これを補正するため、n−1 で割った u2=n−11∑i=1n(Xi−Xˉ)2 が母分散σ²の不偏推定量として使われます。
不偏性は「良い推定量」の必要条件ではありません。実務では、不偏性に加えて分散の小ささ(有効性)や一致性(標本数を増やすと真の値に収束する性質)も合わせて評価します。最尤推定量(MLE)は不偏とは限りませんが、一致性をもち漸近的に効率的という別の強みがあります。
バイアス補正は、推定量 θ^ の系統的なずれ Bias(θ^)=E[θ^]−θ(不偏推定量 で定義)を、定数倍や定数加算で消す操作です。バイアスが θ に比例する形で出るときは、定数倍するだけで不偏化できるのがうれしい点です。題材として、非正則モデルで登場した U(0,θ) の最尤推定量 X(n)=maxXi を補正します。
まず X(n) の期待値を求めます。U(0,θ) の最大値統計量は ベータ分布 Beta(n,1) をスケールした分布にしたがい、その期待値は
E[X(n)]=n+1nθ
です。よってバイアスは Bias=E[X(n)]−θ=n+1nθ−θ=−n+1θ<0 で、つねに負=下方バイアス(真値より小さめ)です。これは X(n)≤θ が構造的に成り立つことの当然の帰結です。補正は簡単で、E[X(n)]=n+1nθ を θ に戻すために逆数 nn+1 を掛ければよく、
θ^unbiased=nn+1X(n)
とすると E[θ^unbiased]=nn+1⋅n+1nθ=θ となり、ぴったり不偏になります。
ここで強調したいのは、この補正がバイアスを消すだけでなくMSE(平均二乗誤差、平均二乗誤差/concepts/%E5%B9%B3%E5%9D%87%E4%BA%8C%E4%B9%97%E8%AA%A4%E5%B7%AE)まで小さくすることです。一般にはバイアスを消すと分散が増えてMSEが悪化することもありますが、この非正則ケースでは『不偏になり、かつMSEも下がる』という一石二鳥が起きます。n=5、真値を比較のため θ=1 とすると、補正前 X(5) の MSE=(n+1)(n+2)2θ2=6×72=21θ2≈0.0476。補正後 nn+1X(5) の MSE=n(n+2)θ2=5×71=35θ2≈0.0286。補正後のほうが MSE が小さく、補正が純粋に推定を改善しています。なお補正係数 nn+1 は n が大きいと 1 に近づくので、大標本ではバイアスは自然に小さくなり、補正の効果も相対的に小さくなります(バイアスのオーダーは 1/n)。
不偏ならBias=0。
不偏分散はn-1で割る
u2=n−11∑(Xi−Xˉ)2 がσ²の不偏推定量。nで割ると過小評価になる(自由度の補正)。
不偏性は唯一の基準ではない
分散の小ささ(有効性)・一致性も重要で、不偏性だけでは推定量の良さを決められない。
不偏推定量の線形和
不偏推定量どうしの線形結合(重みの合計が1)も不偏になる。
順序統計量による不偏推定量の構成
X∼U(−a,a)からの標本で、S=C1Min(Xi)、T=C2Max(Xi)、U=C3Max(∣Xi∣)を不偏推定量にする定数は、順序統計量の期待値から逆算してC1=−n−1n+1、C2=n−1n+1、C3=nn+1となる(2019年度問題1(6))。
バイアス補正:バイアスの定義
Bias(θ^)=E[θ^]−θ
θ に比例する形で出るなら定数倍で不偏化できる。
バイアス補正:X(n) の期待値
U(0,θ) で E[X(n)]=n+1nθ。Beta(n,1) をスケールした分布から。
バイアス補正:下方バイアス
Bias=−θ/(n+1)<0
X(n)≤θ が構造的に成り立つため必ず真値より小さめ。
バイアス補正:補正で不偏化
θ^unbiased=nn+1X(n) とすると E[θ^]=θ。逆数を掛けるだけ。
バイアス補正:MSEも改善(一石二鳥)
この非正則ケースでは補正で MSE も縮む(21θ2→35θ2)。不偏化が悪化を伴わない好例。
n で割った標本分散 s2=80/5=16 は過小評価になる。
▶ バイアス補正:n=5、比較のため真値 θ=1 とする。E[X(5)]=65≈0.833、バイアス =−61。
補正 θ^=56X(5) は E[θ^]=56⋅65=1=θ で不偏。
MSE比較: MSE(X(5))=6×72=211≈0.0476、MSE(θ^)=5×71=351≈0.0286。補正後は不偏かつMSEも小さい。
不偏性=良い推定量と思い込み、分散の大小を考慮しない1回の標本での推定値と真の値が一致しないことを「不偏性が崩れている」と誤解する(不偏性は期待値の話であり、1回の結果が一致するとは限らない)最大値・最小値など順序統計量に基づく推定量を不偏化する係数を求める問題で、期待値の式をそのまま係数だと思い込む(正しくは期待値の逆数を掛けて真値に一致させる)。絶対値の最大値Max(∣Xi∣)を使う変形(台が[0,a]になる)にも対応できるようにする。 バイアス補正をすると必ずMSEも下がると一般化する(一般にはバイアスを消すと分散が増えMSEが悪化することもある。U(0,θ) は両方改善する特別な例) 補正係数を n+1n と逆にする(E[X(n)]=n+1nθ を θ に戻すには逆数 nn+1 を掛ける) X(n) を上方バイアスと誤る(最大値だから大きめに思えるが、X(n)≤θ なので必ず真値以下=下方バイアス) 「バイアス補正」という言葉が単独の論点として使われる過去問は10年分で確認できなかった。ただし非正則モデルのMLE(本ページのU(0,θ)や2022年度問題1(6)のU(0,2θ)のθ̂=Max(Xi)/2など)は構造的にバイアスを持つため、点推定・区間推定の複合問題の中で暗黙的に関係しうる論点である。 モンテカルロ法 13-1
モデリング