tecktom.com

授業で習った形との対応

2026年10月15日6分で読めますAI生成教科書『虎の巻』第1部 第7章 Fisher 情報と Cramér–Rao
目次
  1. 統計学の教科書の書き方
  2. 実験データの最小二乗法
  3. 木下の層

AI 生成・木下が未勉強

相対エントロピーから情報と物理を組み立てる 第1部 第7章 §1.7.6 ★★ 統計学の I(θ)I(\theta)・有効推定量・漸近正規性と、実験の最小二乗法の誤差を、この章の JJ と DD で並べる

統計学の教科書の書き方

数理統計学の授業や教科書では、確率密度を f(x;θ)f(x;\theta)、Fisher 情報を I(θ)I(\theta)(nn 個の標本では In(θ)=nI(θ)I_n(\theta)=nI(\theta))と書き、

I(θ)=E[(∂∂θlog⁡f(X;θ))2]=−E[∂2∂θ2log⁡f(X;θ)],V[θ^]≥1nI(θ)I(\theta)=E\Bigl[\Bigl(\frac{\partial}{\partial\theta}\log f(X;\theta)\Bigr)^2\Bigr]=-E\Bigl[\frac{\partial^2}{\partial\theta^2}\log f(X;\theta)\Bigr],\qquad V[\hat\theta]\ge\frac{1}{nI(\theta)}

の形で習う(log⁡\log は自然対数)。この章の式と次のように対応する。

統計学の教科書この章式
Fisher 情報 I(θ)I(\theta)J(θ)J(\theta)(1.7.5)
E[ ⋅ ]E[\,\cdot\,]、V[ ⋅ ]V[\,\cdot\,]⟨⋅⟩θ\langle \cdot \rangle _\theta、Varθ(⋅)\mathrm{Var}_\theta(\cdot)(1.7.1)
Cramér–Rao の下限 1/(nI)1/(nI)1/(nJ)1/(nJ)(1.7.11)
有効推定量(下限に達する不偏推定量)指数型分布族の平均パラメータでだけ存在(1.7.12)
最尤推定量の漸近正規性 n(θ^−θ)→N(0,1/I)\sqrt n(\hat\theta-\theta)\to N(0,1/I)同じ(1.7.15)
(教科書には出ない)D(pθ ∥ pθ+δ)≈12Jδ2D(p_\theta\,\|\,p_{\theta+\delta})\approx\frac12J\delta^2(1.7.6)

教科書の書き方では Fisher 情報は「スコアの分散」として天下りに定義され、なぜその量が推定の精度に関係するかは Cramér–Rao の証明を読むまで分からない。この章の書き方では、Fisher 情報は最初から「θ\theta を δ\delta ずらした2つの分布の見分けやすさ」で、推定の精度の限界はその帰結として出てくる。加法性もデータ処理による減少も、DD の性質をそのまま移したものになる。

実験データの最小二乗法

物理実験の授業では、測定値 yiy_i(i=1,…,ni=1,\dots,n)が理論値 μi(θ)\mu_i(\theta) のまわりに標準偏差 σi\sigma_i の正規分布で散らばるとして、χ2(θ)=∑i(yi−μi(θ))2/σi2\chi^2(\theta)=\sum_i(y_i-\mu_i(\theta))^2/\sigma_i^2 を最小にする θ\theta を選ぶ(最小二乗法)。このとき ln⁡p=−12χ2+const\ln p=-\frac12\chi^2+\text{const} なので、最小二乗法は正規分布の雑音のもとでの最尤推定である。Fisher 情報は

J(θ)=∑i1σi2(dμidθ)2(1.7.21)J(\theta)=\sum_i\frac{1}{\sigma_i^2}\Bigl(\frac{\mathrm{d} \mu_i}{\mathrm{d} \theta}\Bigr)^2\tag{1.7.21}

で、実験の授業で習う「パラメータの誤差 σθ2=[∑i(∂μi/∂θ)2/σi2]−1\sigma_\theta^2=\bigl[\sum_i(\partial\mu_i/\partial\theta)^2/\sigma_i^2\bigr]^{-1}」は Cramér–Rao の下限 1/J1/J そのものである。理論値 μi\mu_i が θ\theta について線形なら最小二乗の推定量は不偏で、この下限にちょうど達する。同じ量を測った重み付き平均 yˉ=∑i(yi/σi2)/∑i(1/σi2)\bar y=\sum_i(y_i/\sigma_i^2)/\sum_i(1/\sigma_i^2) の誤差 [∑i1/σi2]−1/2\bigl[\sum_i1/\sigma_i^2\bigr]^{-1/2} は、μi=θ\mu_i=\theta とした特別な場合である。

例題 1.7.6 原点を通る直線の傾きの誤差 ★★

オリジナル時刻 ti=1,2,3,4,5 st_i=1,2,3,4,5\,\mathrm{s} で物体の位置 yiy_i を測る。理論は y=θty=\theta t(θ\theta は速さ)で、各測定には標準偏差 σ=0.1 m\sigma=0.1\,\mathrm m の独立な正規分布の誤差が乗る。

  1. J(θ)J(\theta) と Cramér–Rao の下限を求めよ。
  2. 最小二乗の推定量 θ^=∑itiyi/∑iti2\hat\theta=\sum_it_iy_i/\sum_it_i^2 が不偏で、分散が下限に一致することを示せ。
  3. 速さの標準偏差を数値で求めよ。
解答

(1) 式 (1.7.21) で dμi/dθ=ti\mathrm{d}\mu_i/\mathrm{d}\theta=t_i なので J=∑iti2/σ2=55 s2/(0.01 m2)=5.5×103 s2 m−2J=\sum_it_i^2/\sigma^2=55\,\mathrm{s^2}/(0.01\,\mathrm{m^2})=5.5\times10^{3}\,\mathrm{s^2\,m^{-2}}。下限は 1/J=σ2/∑iti2=1.82×10−4 m2 s−21/J=\sigma^2/\sum_it_i^2=1.82\times10^{-4}\,\mathrm{m^2\,s^{-2}} である。

(2) ⟨yi⟩=θti\langle y_i \rangle =\theta t_i を代入すると ⟨θ^⟩=θ∑ti2/∑ti2=θ\langle \hat\theta \rangle =\theta\sum t_i^2/\sum t_i^2=\theta で不偏である。yiy_i が独立なので Var(θ^)=∑iti2σ2/(∑iti2)2=σ2/∑iti2\mathrm{Var}(\hat\theta)=\sum_it_i^2\sigma^2/(\sum_it_i^2)^2=\sigma^2/\sum_it_i^2 で、(1) の下限に等しい。

(3) 1.82×10−4=1.35×10−2 m s−1\sqrt{1.82\times10^{-4}}=1.35\times10^{-2}\,\mathrm{m\,s^{-1}} である。大きい tit_i の測定ほど JJ への寄与 ti2/σ2t_i^2/\sigma^2 が大きく、t=5 st=5\,\mathrm s の1点だけで JJ の 25/55=45%25/55=45\% を占める。したがって傾きを測る実験では、遠い時刻の測定を増やすほうが JJ が大きくなる。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ