tecktom.com

この章のまとめ

2026年10月15日9分で読めますAI生成教科書『虎の巻』第1部 第7章 Fisher 情報と Cramér–Rao

AI 生成・木下が未勉強

相対エントロピーから情報と物理を組み立てる 第1部 第7章 §1.7.7 ★★★ Fisher 情報は DD の2次の係数、Cramér–Rao はその見分けやすさが推定の精度を縛るという主張

この章では、パラメータの値だけが違う2つの分布 pθp_\theta と pθ+δp_{\theta+\delta} を比べた。D(pθ ∥ pθ+δ)D(p_\theta\,\|\,p_{\theta+\delta}) は δ\delta の2次から始まり、その係数の2倍を Fisher 情報 J(θ)J(\theta) と定めた(式 (1.7.6))。JJ はスコアの2乗の平均、または対数確率の2階微分の平均の符号を変えたものとしても書ける。DD の非負性・独立な観測での加法性・データ処理不等式を2次の係数に移すと、J≥0J\ge0、nn 回の観測で nJnJ、加工すると減り十分統計量なら保たれる、という3つの性質が出た。JJ はパラメータの付け替えで計量として変換し、分布 p\sqrt{p} を単位球面上の点とみたときの長さの要素になる。

Cramér–Rao の不等式は、不偏推定量の分散が 1/J1/J より小さくならないことを、推定量とスコアの共分散に Cauchy–Schwarz の不等式を当てて示した。等号に達するのは指数型分布族の平均パラメータに限られる。最尤推定は経験分布からの DD を最小にする推定で、観測が多いとき下限 1/(nJ)1/(nJ) に達する(B 型の極限の主張)。カノニカル分布では逆温度の Fisher 情報がエネルギーの分散になり、温度の推定精度は熱容量で決まる。

第1部第8章(第1部第8章)では、2つの分布が近くない場合に、観測を nn 回重ねると見分け損なう確率が e−nD\mathrm{e}^{-nD} の速さで減ること(Stein の補題)を示す。この章の D≈12Jδ2D\approx\frac12J\delta^2 は、δ\delta が 1/nJ1/\sqrt{nJ} 程度のとき nDnD が1程度になり見分けられるかどうかの境目に来る、という形でその章とつながる。

総合問題 指数分布の速さを推定する ★★★

オリジナル放射性の原子核が崩壊するまでの時間 XX は、崩壊の速さ λ>0\lambda>0 の指数分布 p(x;λ)=λe−λxp(x;\lambda)=\lambda\mathrm{e}^{-\lambda x}(x≥0x\ge0)に従う。nn 個の原子核の崩壊時間 X1,…,XnX_1,\dots,X_n を独立に測り、λ\lambda を推定する。和を S=∑iXiS=\sum_iX_i と書く(SS はガンマ分布に従い、⟨1/S⟩=λ/(n−1)\langle 1/S \rangle =\lambda/(n-1)、⟨1/S2⟩=λ2/((n−1)(n−2))\langle 1/S^2 \rangle =\lambda^2/((n-1)(n-2)) である)。

  1. (1) D(pλ ∥ pλ′)D(p_\lambda\,\|\,p_{\lambda'}) を閉じた式で求め、λ′=λ+δ\lambda'=\lambda+\delta で2次まで展開して J(λ)J(\lambda) を読み取れ。(使う道具:第1部第3章・1.7.1 節)
  2. (2) SS が十分統計量であることを示し、SS の分布から計算した Fisher 情報が nJ(λ)nJ(\lambda) に等しいことを確かめよ。(使う道具:第1部第4章・1.7.2 節)
  3. (3) 不偏推定量 λ^u=(n−1)/S\hat\lambda_{\mathrm u}=(n-1)/S の分散を求め、Cramér–Rao の下限との比を求めよ。(使う道具:1.7.3 節)
  4. (4) 最尤推定量 λ^ML\hat\lambda_{\mathrm{ML}} を求め、その偏りと2乗誤差の平均を求めよ。n→∞n\to\infty で nn 倍した2乗誤差の平均が 1/J1/J に近づくことを確かめよ。(使う道具:1.7.4 節)
  5. (5) 指数分布は、エネルギー x≥0x\ge0 を持つ状態が一様に並んだ系のカノニカル分布で λ=β\lambda=\beta としたものである。式 (1.7.18) を使って温度の推定の相対誤差の下限を求め、(1) の J(λ)J(\lambda) を付け替えた結果と一致することを確かめよ。(使う道具:1.7.5 節・解析力学から様々な理論分野を組み立てる 第3部第3章)
  6. (6) 崩壊時間を測らず、時刻 aa までに崩壊したかどうかだけを記録する。1個あたりの Fisher 情報 Ja(λ)J_a(\lambda) を求め、aa を選んで最大にしたとき J(λ)J(\lambda) の何%が残るかを数値で求めよ。(使う道具:1.7.2 節・例題 1.7.2)
解答

(1) ln⁡pλpλ′=ln⁡λλ′−(λ−λ′)x\ln\frac{p_\lambda}{p_{\lambda'}}=\ln\frac{\lambda}{\lambda'}-(\lambda-\lambda')x で ⟨X⟩λ=1/λ\langle X \rangle _\lambda=1/\lambda なので

D(pλ ∥ pλ′)=ln⁡λλ′+λ′λ−1D(p_\lambda\,\|\,p_{\lambda'})=\ln\frac{\lambda}{\lambda'}+\frac{\lambda'}{\lambda}-1

である。λ′=λ+δ\lambda'=\lambda+\delta とおき r=δ/λr=\delta/\lambda と書くと −ln⁡(1+r)+r=12r2−13r3+⋯-\ln(1+r)+r=\frac12r^2-\frac13r^3+\cdots なので D=δ22λ2+O(δ3)D=\frac{\delta^2}{2\lambda^2}+\mathcal{O}(\delta^3)、したがって J(λ)=1/λ2J(\lambda)=1/\lambda^2 である。DD は ln⁡(λ/λ′)\ln(\lambda/\lambda') と λ′/λ\lambda'/\lambda だけで書けるので、速さを何倍ずらしたかだけで決まり、J∝1/λ2J\propto1/\lambda^2 はその反映である。

(2) p(x1,…,xn;λ)=λne−λSp(x_1,\dots,x_n;\lambda)=\lambda^n\mathrm{e}^{-\lambda S} は SS と λ\lambda だけの関数 g(S;λ)=λne−λSg(S;\lambda)=\lambda^n\mathrm{e}^{-\lambda S} に h=1h=1 を掛けた形なので、Fisher–Neyman の分解により SS は十分統計量である。SS の密度は λnsn−1e−λs/(n−1)!\lambda^ns^{n-1}\mathrm{e}^{-\lambda s}/(n-1)! で、そのスコアは n/λ−sn/\lambda-s、2乗の平均は Var(S)=n/λ2\mathrm{Var}(S)=n/\lambda^2。これは nJ(λ)=n/λ2nJ(\lambda)=n/\lambda^2 に等しく、命題 1.7.2 の等号の場合になっている。

(3) ⟨λ^u⟩=(n−1)λ/(n−1)=λ\langle \hat\lambda_{\mathrm u} \rangle =(n-1)\lambda/(n-1)=\lambda で不偏である。⟨λ^u2⟩=(n−1)2λ2/((n−1)(n−2))=(n−1)λ2/(n−2)\langle \hat\lambda_{\mathrm u}^2 \rangle =(n-1)^2\lambda^2/((n-1)(n-2))=(n-1)\lambda^2/(n-2) より

Var(λ^u)=λ2n−2(n≥3)\mathrm{Var}(\hat\lambda_{\mathrm u})=\frac{\lambda^2}{n-2}\qquad(n\ge3)

である。下限は 1/(nJ)=λ2/n1/(nJ)=\lambda^2/n なので比は n/(n−2)>1n/(n-2)>1。有限の nn では下限に達しない。スコア n/λ−Sn/\lambda-S が 1/S1/S の1次式でないので、等号の条件 (1.7.12) が成り立たないからである。n→∞n\to\infty で比は 1 に近づく。

(4) ℓn(λ)=nln⁡λ−λS\ell_n(\lambda)=n\ln\lambda-\lambda S を最大にして λ^ML=n/S\hat\lambda_{\mathrm{ML}}=n/S を得る。⟨λ^ML⟩=nλ/(n−1)\langle \hat\lambda_{\mathrm{ML}} \rangle =n\lambda/(n-1) なので偏りは b=λ/(n−1)b=\lambda/(n-1)、分散は Var=n2λ2/((n−1)2(n−2))\mathrm{Var}=n^2\lambda^2/((n-1)^2(n-2))。2乗誤差の平均は

Var+b2=λ2(n2+(n−2))(n−1)2(n−2)=(n+2) λ2(n−1)(n−2)\mathrm{Var}+b^2=\frac{\lambda^2\bigl(n^2+(n-2)\bigr)}{(n-1)^2(n-2)}=\frac{(n+2)\,\lambda^2}{(n-1)(n-2)}

で、nn 倍すると λ2n(n+2)/((n−1)(n−2))→λ2=1/J\lambda^2n(n+2)/((n-1)(n-2))\to\lambda^2=1/J となる。偏りは O(1/n)\mathcal{O}(1/n) で、2乗誤差への寄与 b2=O(1/n2)b^2=\mathcal{O}(1/n^2) は分散の O(1/n)\mathcal{O}(1/n) より先に消える。

(5) この系では ⟨E⟩=kBT\langle E \rangle = k_{\mathrm B} T なので CV=kBC_V= k_{\mathrm B} であり、式 (1.7.19) から Var(T^)/T2≥1/n\mathrm{Var}(\hat T)/T^2\ge1/n となる。一方 (1) の J(λ)=1/λ2=β2J(\lambda)=1/\lambda^2=\beta^2 を付け替えの式 (1.7.9) で TT に移すと、dβ/dT=−1/(kBT2)\mathrm{d}\beta/\mathrm{d} T=-1/( k_{\mathrm B} T^2) より J(T)=β2/(kBT2)2=1/T2J(T)=\beta^2/( k_{\mathrm B} T^2)^2=1/T^2、したがって 1/(nJ(T))=T2/n1/(nJ(T))=T^2/n で一致する。1自由度の温度計を nn 回読む相対誤差は 1/n1/\sqrt n 以上である。

(6) 時刻 aa までに崩壊する確率は 1−e−λa1-\mathrm{e}^{-\lambda a} で、記録は Bernoulli 変数である。例題 1.7.1 と付け替えの式から

Ja(λ)=(ae−λa)2(1−e−λa)e−λa=a2e−λa1−e−λa=1λ2 y2ey−1(y=λa)J_a(\lambda)=\frac{\bigl(a\mathrm{e}^{-\lambda a}\bigr)^2}{(1-\mathrm{e}^{-\lambda a})\mathrm{e}^{-\lambda a}}=\frac{a^2\mathrm{e}^{-\lambda a}}{1-\mathrm{e}^{-\lambda a}}=\frac{1}{\lambda^2}\,\frac{y^2}{\mathrm{e}^{y}-1}\qquad(y=\lambda a)

となる。y2/(ey−1)y^2/(\mathrm{e}^y-1) の最大の条件は 2(ey−1)=yey2(\mathrm{e}^y-1)=y\mathrm{e}^y で、数値で解くと y=1.594y=1.594、そのとき y2/(ey−1)=0.648y^2/(\mathrm{e}^y-1)=0.648 である。平均寿命 1/λ1/\lambda の約1.6倍の時刻で区切ると、もとの Fisher 情報の64.8%が残る。これは命題 1.7.2 の Ja≤JJ_a\le J の1つの例で、どう aa を選んでも、少なくとも35%が失われる。

木下の勉強ステータス

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

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

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