院試hub

点推定と最尤推定不偏性・一致性・フィッシャー情報量・クラメール・ラオの不等式

最尤推定量は、観測値 x1,…,xnx_1, \dots, x_n を固定したときの尤度 L(θ)=∏if(xi;θ)L(\theta) = \prod_i f(x_i;\theta) を最大にする θ\theta である。ふつうは対数尤度 ℓ(θ)=log⁡L(θ)\ell(\theta) = \log L(\theta) を θ\theta で微分して 00 とおき、解が最大であることを確かめる。推定量の良さは偏り(バイアス)と分散、その和である平均2乗誤差で比べる。正則条件のもとでは、不偏推定量の分散はクラメール・ラオの下界 1/(nI(θ))1/(nI(\theta)) 以上で、最尤推定量は nn が大きいとき分散がこの下界に近い正規分布に従う。

確率・統計 の第5単元(全8単元)。先に読む単元:大数の法則と中心極限定理

院試数学 基礎完成シリーズ 確率統計編

基礎ノートと類題を答案の形で解いた問題集(2,480円)。収録章と見本ページを商品ページで確認できる。過去問の本文は収録していない。

確率統計編の内容を見る

推定量の良さを測る基準

定義

推定量・偏り・平均2乗誤差

母数 θ\theta を、標本 X1,…,XnX_1, \dots, X_n の関数 θ^=T(X1,…,Xn)\hat\theta = T(X_1, \dots, X_n) で推定する。θ^\hat\theta を推定量、観測値を代入した値を推定値という。

偏りを Bias⁡(θ^)=E⁡θ[θ^]−θ\operatorname{Bias}(\hat\theta) = \operatorname{E}_\theta[\hat\theta] - \theta と定め、すべての θ\theta で 00 なら不偏推定量という。n→∞n \to \infty で θ^\hat\theta が θ\theta に確率収束するとき一致推定量という。

平均2乗誤差の分解
MSE⁡(θ^)=E⁡θ[(θ^−θ)2]=Var⁡θ(θ^)+Bias⁡(θ^)2\operatorname{MSE}(\hat\theta) = \operatorname{E}_\theta[(\hat\theta - \theta)^2] = \operatorname{Var}_\theta(\hat\theta) + \operatorname{Bias}(\hat\theta)^2
標本平均
Xˉ=1n∑i=1nXi,E⁡[Xˉ]=μ,Var⁡(Xˉ)=σ2n\bar X = \frac{1}{n}\sum_{i=1}^{n} X_i, \quad \operatorname{E}[\bar X] = \mu, \quad \operatorname{Var}(\bar X) = \frac{\sigma^2}{n}
不偏分散
S2=1n−1∑i=1n(Xi−Xˉ)2,E⁡[S2]=σ2S^2 = \frac{1}{n-1}\sum_{i=1}^{n}(X_i - \bar X)^2, \quad \operatorname{E}[S^2] = \sigma^2

∑(Xi−Xˉ)2\sum (X_i - \bar X)^2 の期待値が (n−1)σ2(n-1)\sigma^2 になるため。nn で割ると偏りが −σ2/n-\sigma^2/n になる。

一致性は、偏りと分散がともに 00 に収束すれば成り立つ(チェビシェフの不等式から MSE⁡→0\operatorname{MSE} \to 0 が確率収束を導く)。不偏でも分散の大きい推定量より、少し偏っていても平均2乗誤差の小さい推定量のほうがよいことがある。

最尤推定量の求め方

定義

尤度と最尤推定量

独立同分布の標本で、密度(離散型なら確率関数)を f(x;θ)f(x;\theta) とする。観測値を固定し θ\theta の関数とみた L(θ)=∏i=1nf(xi;θ)L(\theta) = \prod_{i=1}^{n} f(x_i;\theta) を尤度関数、ℓ(θ)=log⁡L(θ)\ell(\theta) = \log L(\theta) を対数尤度という。L(θ)L(\theta) を最大にする θ^\hat\theta を最尤推定量という。

  1. 対数尤度 ℓ(θ)=∑ilog⁡f(xi;θ)\ell(\theta) = \sum_i \log f(x_i;\theta) を書く。
  2. 尤度方程式 ∂ℓ/∂θ=0\partial\ell/\partial\theta = 0 を解く。母数が複数なら連立方程式にする。
  3. 2階微分の符号や増減表、境界での値から、解が最大であることを確かめる。
  4. 台(f>0f > 0 となる範囲)が θ\theta に依存するときは微分を使わず、L(θ)L(\theta) の形を直接調べる。

計算例

正規分布の平均と分散

X1,…,XnX_1, \dots, X_n を N(μ,σ2)N(\mu, \sigma^2) からの標本とする。μ,σ2\mu, \sigma^2 の最尤推定量を求めよ。観測値が 2,4,4,5,52, 4, 4, 5, 5 のときの値も求めよ。

解答

ℓ(μ,σ2)=−n2log⁡(2πσ2)−12σ2∑i(xi−μ)2\ell(\mu, \sigma^2) = -\frac{n}{2}\log(2\pi\sigma^2) - \frac{1}{2\sigma^2}\sum_{i}(x_i - \mu)^2

∂ℓ/∂μ=∑(xi−μ)/σ2=0\partial\ell/\partial\mu = \sum(x_i - \mu)/\sigma^2 = 0 から μ^=xˉ\hat\mu = \bar x。∂ℓ/∂σ2=−n/(2σ2)+∑(xi−μ)2/(2σ4)=0\partial\ell/\partial\sigma^2 = -n/(2\sigma^2) + \sum(x_i-\mu)^2/(2\sigma^4) = 0 に μ=xˉ\mu = \bar x を入れて σ^2=1n∑(xi−xˉ)2\hat\sigma^2 = \dfrac{1}{n}\sum(x_i - \bar x)^2。

観測値では xˉ=20/5=4\bar x = 20/5 = 4、偏差の2乗和は 4+0+0+1+1=64 + 0 + 0 + 1 + 1 = 6 なので μ^=4\hat\mu = 4、σ^2=6/5=1.2\hat\sigma^2 = 6/5 = 1.2。

検算

σ^2\hat\sigma^2 は nn で割るので不偏でなく、E⁡[σ^2]=(n−1)σ2/n\operatorname{E}[\hat\sigma^2] = (n-1)\sigma^2/n。同じデータの不偏分散は 6/4=1.56/4 = 1.5 である。n→∞n \to \infty で両者の差は 00 に近づき、どちらも一致推定量である。

計算例

指数分布の率

X1,…,Xn∼Exp(λ)X_1, \dots, X_n \sim \mathrm{Exp}(\lambda)(密度 λe−λx\lambda e^{-\lambda x})の最尤推定量を求めよ。観測値が 1,2,3,61, 2, 3, 6 のときの値も求めよ。

解答

ℓ(λ)=nlog⁡λ−λ∑xi\ell(\lambda) = n\log\lambda - \lambda\sum x_i。ℓ′(λ)=n/λ−∑xi=0\ell'(\lambda) = n/\lambda - \sum x_i = 0 から λ^=n/∑xi=1/xˉ\hat\lambda = n/\sum x_i = 1/\bar x。ℓ′′(λ)=−n/λ2<0\ell''(\lambda) = -n/\lambda^2 < 0 なので最大である。観測値では λ^=4/12=1/3\hat\lambda = 4/12 = 1/3。

検算

平均 θ=1/λ\theta = 1/\lambda で表した最尤推定量は θ^=xˉ=3\hat\theta = \bar x = 3 で、1/λ^1/\hat\lambda と一致する(最尤推定量の不変性)。

フィッシャー情報量とクラメール・ラオの不等式

定義

フィッシャー情報量

1個の観測あたりのフィッシャー情報量を

I(θ)=E⁡θ[(∂∂θlog⁡f(X;θ))2]=−E⁡θ[∂2∂θ2log⁡f(X;θ)]I(\theta) = \operatorname{E}_\theta\left[\left(\frac{\partial}{\partial\theta}\log f(X;\theta)\right)^2\right] = -\operatorname{E}_\theta\left[\frac{\partial^2}{\partial\theta^2}\log f(X;\theta)\right]

と定める。2つ目の等号は下の正則条件のもとで成り立つ。独立同分布の nn 個の標本の情報量は nI(θ)nI(\theta) である。

定理・法則

クラメール・ラオの不等式

条件

X1,…,XnX_1, \dots, X_n は独立同分布。台 {x:f(x;θ)>0}\{x : f(x;\theta) > 0\} が θ\theta によらない。log⁡f\log f が θ\theta で微分でき、∫f dx\int f\,dx と ∫T∏f dx\int T\prod f\,dx の θ\theta での微分を積分記号の中に入れられる。0<I(θ)<∞0 < I(\theta) < \infty。TT は θ\theta の不偏推定量で分散が有限。

主張

Var⁡θ(T)≥1nI(θ)\operatorname{Var}_\theta(T) \ge \frac{1}{nI(\theta)}

要点

スコア U=∂ℓ/∂θU = \partial\ell/\partial\theta が E⁡[U]=0\operatorname{E}[U] = 0、Cov⁡(T,U)=1\operatorname{Cov}(T, U) = 1 を満たすことと、コーシー・シュワルツの不等式から出る。等号が成り立つ不偏推定量を有効推定量という。一様分布 U(0,θ)U(0,\theta) のように台が θ\theta に依存すると、この不等式は成り立たない。

計算例

ポアソン分布の標本平均は有効推定量

X1,…,Xn∼Po(λ)X_1, \dots, X_n \sim \mathrm{Po}(\lambda) のとき、I(λ)I(\lambda) とクラメール・ラオの下界を求め、Xˉ\bar X と比べよ。

解答

log⁡f(x;λ)=−λ+xlog⁡λ−log⁡x!\log f(x;\lambda) = -\lambda + x\log\lambda - \log x! で、∂2/∂λ2=−x/λ2\partial^2/\partial\lambda^2 = -x/\lambda^2。よって I(λ)=E⁡[X]/λ2=1/λI(\lambda) = \operatorname{E}[X]/\lambda^2 = 1/\lambda、下界は λ/n\lambda/n。Xˉ\bar X は不偏で Var⁡(Xˉ)=λ/n\operatorname{Var}(\bar X) = \lambda/n なので、下界に一致し有効推定量である。Xˉ\bar X は最尤推定量でもある。

検算

1階微分の2乗で計算しても E⁡[(X/λ−1)2]=Var⁡(X)/λ2=1/λ\operatorname{E}[(X/\lambda - 1)^2] = \operatorname{Var}(X)/\lambda^2 = 1/\lambda で同じ値になる。

最尤推定量の漸近的な性質

定理・法則

最尤推定量の一致性と漸近正規性

条件

独立同分布。母数が識別可能(θ≠θ′\theta \ne \theta' なら分布が異なる)。台が θ\theta によらない。log⁡f\log f が θ\theta について3回微分でき、その3階微分が可積分な関数で抑えられる。真の値 θ\theta が母数空間の内点で、0<I(θ)<∞0 < I(\theta) < \infty。θ^n\hat\theta_n は尤度方程式の一致解。

主張

θ^n\hat\theta_n は θ\theta の一致推定量で、

n(θ^n−θ)→dN(0,1I(θ))\sqrt{n}(\hat\theta_n - \theta) \xrightarrow{d} N\left(0, \frac{1}{I(\theta)}\right)

要点

尤度方程式を真の値のまわりでテイラー展開し、スコアの和に中心極限定理、2階微分の平均に大数の法則を使う。nn が大きいとき、分散がクラメール・ラオの下界 1/(nI(θ))1/(nI(\theta)) に近づく。

不変性
θ^ が θ の最尤推定量⇒g(θ^) は g(θ) の最尤推定量\hat\theta \text{ が } \theta \text{ の最尤推定量} \Rightarrow g(\hat\theta) \text{ は } g(\theta) \text{ の最尤推定量}

gg が1対1でないときは、g(θ)=ηg(\theta) = \eta となる θ\theta での尤度の上限を η\eta の尤度とみなして成り立つ。

近似の標準誤差
SE⁡(θ^)≈1nI(θ^)\operatorname{SE}(\hat\theta) \approx \frac{1}{\sqrt{nI(\hat\theta)}}

計算例

台が母数に依存する例:一様分布 U(0,θ)U(0,\theta)

X1,…,Xn∼U(0,θ)X_1, \dots, X_n \sim U(0,\theta) の最尤推定量を求め、その偏りと分散をモーメント推定量 2Xˉ2\bar X と比べよ。

解答

L(θ)=θ−nL(\theta) = \theta^{-n}(θ≥max⁡ixi\theta \ge \max_i x_i)、それ以外で 00。θ−n\theta^{-n} は減少関数なので、θ^=M=max⁡iXi\hat\theta = M = \max_i X_i で最大になる。尤度方程式 L′(θ)=0L'(\theta) = 0 には解がない。

P(M≤x)=(x/θ)nP(M \le x) = (x/\theta)^n から密度 nxn−1/θnnx^{n-1}/\theta^n で、E⁡[M]=nn+1θ\operatorname{E}[M] = \dfrac{n}{n+1}\theta、Var⁡(M)=nθ2(n+1)2(n+2)\operatorname{Var}(M) = \dfrac{n\theta^2}{(n+1)^2(n+2)}。不偏にした n+1nM\dfrac{n+1}{n}M の分散は θ2n(n+2)\dfrac{\theta^2}{n(n+2)}。2Xˉ2\bar X は不偏で分散 4⋅θ212n=θ23n4 \cdot \dfrac{\theta^2}{12n} = \dfrac{\theta^2}{3n}。

検算

n≥2n \ge 2 で θ2n(n+2)<θ23n\dfrac{\theta^2}{n(n+2)} < \dfrac{\theta^2}{3n} となり、MM にもとづく推定量の分散は 1/n21/n^2 の速さで小さくなる。正則な場合の 1/n1/n より速く、クラメール・ラオの不等式が使えない例になっている。n=1n = 1 では両者は一致する。

答案でよくある誤り

  1. 1

    尤度方程式の解を、最大であることを確かめずに最尤推定量とする。2階微分の符号や増減を調べ、解が複数あれば尤度を比べる。

  2. 2

    台が θ\theta に依存する分布で尤度を微分する。U(0,θ)U(0,\theta) では L(θ)L(\theta) の定義域 θ≥max⁡xi\theta \ge \max x_i を書き、端点で最大になることを示す。

  3. 3

    最尤推定量はいつも不偏だと考える。正規分布の σ^2\hat\sigma^2 や一様分布の max⁡Xi\max X_i のように偏りがある。不偏性が問われたら E⁡[θ^]\operatorname{E}[\hat\theta] を計算する。

  4. 4

    フィッシャー情報量を1個分の I(θ)I(\theta) のまま下界に使う。nn 個の標本では nI(θ)nI(\theta) で、下界は 1/(nI(θ))1/(nI(\theta)) になる。

  5. 5

    σ\sigma の最尤推定量を S2\sqrt{S^2}(不偏分散の平方根)とする。不変性から σ^=σ^2\hat\sigma = \sqrt{\hat\sigma^2} で、分母は nn である。

点推定と最尤推定が出た院試の過去問

院試hubが年度別解説の大問を分類した集計から、この単元の内容を扱う出題を挙げる。問題文は各大学の公式過去問で確かめ、解き方は各年度の解説ページで読める。

確率・統計の前後の単元

このページの本文を数式のLaTeXのまま読むには Markdown版。最終更新 2026-11-16。誤りに気づいたら 運営者 へ連絡してほしい。

院試数学 基礎完成シリーズ 確率統計編

基礎ノートと類題を答案の形で解いた問題集(2,480円)。収録章と見本ページを商品ページで確認できる。過去問の本文は収録していない。

確率統計編の内容を見る