tecktom.com

当てる DD の性質:非負性・加法性・データ処理

2026年10月15日7分で読めますAI生成教科書『虎の巻』第1部 第7章 Fisher 情報と Cramér–Rao
目次
  1. 非負性と加法性
  2. データ処理不等式と十分統計量
  3. パラメータの付け替えと幾何
  4. 木下の層

AI 生成・木下が未勉強

相対エントロピーから情報と物理を組み立てる 第1部 第7章 §1.7.2 ★★★ DD の3つの性質を2次の係数に移すと、J≥0J\ge0・nn 倍・データ処理で減ることがそのまま出る

式 (1.7.6) は「DD について成り立つ不等式や等式を δ2/2\delta^2/2 で割って δ→0\delta\to0 とすれば、JJ について成り立つ」ことを意味する。この節では、第1部第3章と第4章で示した DD の性質を順にこの方法で JJ へ移す。

非負性と加法性

D≥0D\ge0(第1部第3章)を δ2/2\delta^2/2 で割れば J(θ)≥0J(\theta)\ge0 である。定義 J=⟨sθ2⟩J=\langle s_\theta^2 \rangle からも明らかだが、J(θ)=0J(\theta)=0 となるのは sθ(x)=0s_\theta(x)=0 がほとんどすべての xx で成り立つとき、つまり θ\theta を少し動かしても分布がまったく変わらない(θ\theta が識別できない)ときに限る。

独立な2つの確率変数の組 (X,Y)(X,Y) で、pθ(x,y)=pθ(x) rθ(y)p_\theta(x,y)=p_\theta(x)\,r_\theta(y) と分かれるとき、DD は足し算になる:D(pθrθ ∥ pθ+δrθ+δ)=D(pθ ∥ pθ+δ)+D(rθ ∥ rθ+δ)D(p_\theta r_\theta\,\|\,p_{\theta+\delta}r_{\theta+\delta})=D(p_\theta\,\|\,p_{\theta+\delta})+D(r_\theta\,\|\,r_{\theta+\delta})。両辺の δ2\delta^2 の係数を比べて、Fisher 情報も足し算になる。特に同じ分布から独立に nn 回観測した Xn=(X1,…,Xn)X^n=(X_1,\dots,X_n) では

JXn(θ)=n J(θ)(1.7.7)J_{X^n}(\theta)=n\,J(\theta)\tag{1.7.7}

である。以後、どの確率変数の Fisher 情報かを区別したいときは JX(θ)J_X(\theta) のように添字を付ける。

データ処理不等式と十分統計量

観測値 XX をそのまま使わず、何らかの加工 T=t(X)T=t(X)(ビンにまとめる、符号だけ記録する、平均だけ残す、雑音の入った通信路を通す)をしたとする。加工は確率行列 WW を掛ける操作で、θ\theta によらない。第1部第4章のデータ処理不等式 D(Wp ∥ Wq)≤D(p ∥ q)D(Wp\,\|\,Wq)\le D(p\,\|\,q) をp=pθp=p_\theta、q=pθ+δq=p_{\theta+\delta} に当てて δ2/2\delta^2/2 で割ると、次の命題になる。決定的な加工 T=t(X)T=t(X) については第1部第4章の式 (1.4.24) で先に見た。ここでは確率的な加工まで含めて命題にまとめる。

命題 1.7.2(Fisher 情報のデータ処理不等式 ★★★)

θ\theta によらない加工 T=t(X)T=t(X)(確率的な加工でもよい)について

JT(θ)≤JX(θ)(1.7.8)J_T(\theta)\le J_X(\theta)\tag{1.7.8}

が成り立つ。TT が θ\theta の十分統計量なら、すべての θ\theta で等号が成り立つ。

等号の条件は、十分統計量では DD が保たれる(第1部第4章)ことから従う。十分統計量に直接当てはめると次のように確かめられる。Fisher–Neyman の分解 p(x;θ)=g(t(x);θ) h(x)p(x;\theta)=g(t(x);\theta)\,h(x) が成り立てばsθ(x)=∂θln⁡g(t(x);θ)s_\theta(x)=\partial_\theta\ln g(t(x);\theta) となり、スコアは t(x)t(x) だけの関数である。このときスコアの分布は TT だけから計算でき、JT=JXJ_T=J_X となる。

加工で情報が増えないこと

命題 1.7.2 は、「データを後からどう処理しても、パラメータについての情報は増えない」という直観を、Fisher 情報という数で述べたものである。θ\theta が少し違う2つの分布 pθp_\theta と pθ+δp_{\theta+\delta} は、加工の後のほうが見分けにくい。推定量 θ^(X)\hat\theta(X) 自体も XX の加工の1つなので、推定量の分布が持つ Fisher 情報も JXJ_X を超えない。これが次の節の Cramér–Rao の不等式の背景になる。

パラメータの付け替えと幾何

同じ族を別のパラメータ η\eta で表す(θ=θ(η)\theta=\theta(\eta))と、分布 pθp_\theta と pθ+δp_{\theta+\delta} そのものは変わらないので DD も変わらない。η\eta を ϵ\epsilon だけずらすと θ\theta は δ=θ′(η)ϵ\delta=\theta'(\eta)\epsilon だけずれるから、2次の係数を比べて

J(η)(η)=J(θ)(θ)(dθdη)2(1.7.9)J^{(\eta)}(\eta)=J^{(\theta)}(\theta)\Bigl(\frac{\mathrm{d} \theta}{\mathrm{d} \eta}\Bigr)^2\tag{1.7.9}

となる。JJ はパラメータの変化の2乗に掛かる係数なので、長さの2乗を測る計量と同じ変換をする。

この計量の意味は、値が有限個の分布で見るとはっきりする。pip_i(i=1,…,mi=1,\dots,m)についてJ=∑i(∂θpi)2/pi=4∑i(∂θpi)2J=\sum_i(\partial_\theta p_i)^2/p_i=4\sum_i(\partial_\theta\sqrt{p_i})^2 と書き直せる。ベクトル (p1,…,pm)(\sqrt{p_1},\dots,\sqrt{p_m}) は ∑i(pi)2=1\sum_i(\sqrt{p_i})^2=1 を満たすので、半径1の球面上の点である。したがって 14J(θ) dθ2\frac14J(\theta)\,\mathrm{d}\theta^2 は、分布の族が球面上に描く曲線の長さの要素の2乗に等しい。この計量で測った2つの分布 pp、qq の間の最短距離は球面上の大円の弧で、2arccos⁡∑ipiqi2\arccos\sum_i\sqrt{p_iq_i} になる(係数 2 は ds2=J dθ2\mathrm{d} s^2=J\,\mathrm{d}\theta^2 と測るため)。この距離は、係数 2 を除いて W. K. Wootters が「統計的距離」と呼んだもので、量子状態の間の距離へそのまま拡張された(W. K. Wootters (1981), Phys. Rev. D 23, 357–362)。量子版は第3部第8章(第3部第8章)で扱う。

例題 1.7.2 符号だけ記録すると Fisher 情報はどれだけ減るか ★★★

オリジナルXX は平均 θ\theta、分散 1 の正規分布に従う。XX の値は記録せず、T=1T=1(X>0X>0 のとき)、T=0T=0(X≤0X\le0 のとき)だけを残す。

  1. JX(θ)J_X(\theta) を求めよ。
  2. JT(θ)J_T(\theta) を標準正規分布の密度 ϕ\phi と分布関数 Φ\Phi で表し、θ=0\theta=0 での値を求めよ。
  3. JT(θ)≤JX(θ)J_T(\theta)\le J_X(\theta) を θ=0, 1, 2\theta=0,\,1,\,2 で数値で確かめ、θ\theta が大きいとき符号の記録がほとんど役に立たない理由を述べよ。
解答

(1) ln⁡p(x;θ)=−12(x−θ)2−12ln⁡(2π)\ln p(x;\theta)=-\frac12(x-\theta)^2-\frac12\ln(2\pi) なので sθ(x)=x−θs_\theta(x)=x-\theta、JX=⟨(X−θ)2⟩θ=1J_X=\langle (X-\theta)^2 \rangle _\theta=1 である。

(2) TT は表の確率 Pr⁡(X>0)=Φ(θ)\Pr(X>0)=\Phi(\theta) の Bernoulli 変数である。例題 1.7.1 の結果と付け替えの式 (1.7.9)(Bernoulli の確率 Φ(θ)\Phi(\theta) を θ\theta で微分すると ϕ(θ)\phi(\theta))から

JT(θ)=ϕ(θ)2Φ(θ)(1−Φ(θ))J_T(\theta)=\frac{\phi(\theta)^2}{\Phi(\theta)\bigl(1-\Phi(\theta)\bigr)}

となる。θ=0\theta=0 では ϕ(0)2=1/(2π)\phi(0)^2=1/(2\pi)、Φ(0)=12\Phi(0)=\frac12 なので JT(0)=1/(2π)1/4=2π=0.637J_T(0)=\dfrac{1/(2\pi)}{1/4}=\dfrac{2}{\pi}=0.637 である。

(3) JT(0)=0.637J_T(0)=0.637、JT(1)=0.439J_T(1)=0.439、JT(2)=0.131J_T(2)=0.131 で、どれも JX=1J_X=1 より小さい。θ\theta が大きいと X>0X>0 がほとんど確実に起こり、θ\theta を少し動かしても TT の分布がほとんど変わらない。符号という1ビットの記録は、平均が 0 の近くにあるときだけ、もとの観測の Fisher 情報の6割強を残す。

木下の勉強ステータス

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

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

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