相対エントロピーから情報と物理を組み立てる 第1部 第7章 §1.7.1 ★★★ D ( p θ ∥ p θ + δ ) ≈ 1 2 J ( θ ) δ 2 D(p_\theta\,\|\,p_{\theta+\delta})\approx\frac12J(\theta)\delta^2 D ( p θ ∥ p θ + δ ) ≈ 2 1 J ( θ ) δ 2 が Fisher 情報の定義であり、章の出発点
この章で比べる2つの分布は、同じ形の分布の族の中で、パラメータの値だけをわずかに変えた p θ p_\theta p θ と p θ + δ p_{\theta+\delta} p θ + δ である。第1部第3章(第1部第3章 )で定めた相対エントロピー D ( p ∥ q ) D(p\,\|\,q) D ( p ∥ q ) は、2つの分布がどれだけ見分けやすいかを測る量であった。δ \delta δ が小さいとき D ( p θ ∥ p θ + δ ) D(p_\theta\,\|\,p_{\theta+\delta}) D ( p θ ∥ p θ + δ ) は δ 2 \delta^2 δ 2 に比例し、その係数の2倍が Fisher 情報 J ( θ ) J(\theta) J ( θ ) になる。つまり Fisher 情報は、パラメータを少し動かしたときに分布がどれだけ見分けやすくなるかを表す数である。
D D D の性質を当てると、Fisher 情報の性質が順に出てくる。非負性からは J ( θ ) ≥ 0 J(\theta)\ge0 J ( θ ) ≥ 0 が、独立な観測で D D D が足し算になることからはn n n 回の観測で Fisher 情報が n n n 倍になることが、データ処理不等式(第1部第4章 )からは、データを加工するとFisher 情報が減り、十分統計量なら減らないことが従う。そのうえで、パラメータを推定する人がどれだけ精度よく推定できるかに下限があること(Cramér–Rao の不等式)、観測の回数が多いとき最尤推定がその下限に達すること、統計力学のカノニカル分布では Fisher 情報がエネルギーの揺らぎ、したがって熱容量そのものになることを示す。最後に、統計学の教科書や実験データの解析で習う形と、この章の形を並べる。
パラメータの付いた分布の族
コインの表が出る確率 θ \theta θ 、正規分布の平均 θ \theta θ 、指数分布の速さ θ \theta θ のように、1つの実数 θ \theta θ を決めると分布が1つ決まる族 { p θ } \{p_\theta\} { p θ } を考える。確率変数 X X X の値 x x x における確率(連続なら確率密度)を
p θ ( x ) ≡ p ( x ; θ ) (1.7.1) p_\theta(x)\equiv p(x;\theta)\tag{1.7.1} p θ ( x ) ≡ p ( x ; θ ) ( 1.7.1 )
と書き、p θ p_\theta p θ のもとでの平均を ⟨ f ⟩ θ = ∑ x p ( x ; θ ) f ( x ) \langle f \rangle _\theta=\sum_xp(x;\theta)f(x) ⟨ f ⟩ θ = ∑ x p ( x ; θ ) f ( x ) 、分散を V a r θ ( f ) \mathrm{Var}_\theta(f) Var θ ( f ) と書く。和は X X X が連続なら積分 ∫ d x \int\mathrm{d} x ∫ d x に読み替える。θ \theta θ を知らない人が X X X の値を見て θ \theta θ を当てる問題を推定 と呼ぶ。
この章を通して次の正則条件 を仮定する。第一に、p ( x ; θ ) > 0 p(x;\theta)>0 p ( x ; θ ) > 0 となる x x x の範囲(台)が θ \theta θ によらない。第二に、p ( x ; θ ) p(x;\theta) p ( x ; θ ) は θ \theta θ について3回微分でき、θ \theta θ での微分と x x x についての和(積分)の順序を入れ替えてよい。第三に、後で定める J ( θ ) J(\theta) J ( θ ) が正で有限である。第一の条件が成り立たない例(一様分布 [ 0 , θ ] [0,\theta] [ 0 , θ ] )は 1.7.3 節で扱う。
スコア
対数をとった確率の θ \theta θ 微分
s θ ( x ) ≔ ∂ ∂ θ ln p ( x ; θ ) = ∂ θ p ( x ; θ ) p ( x ; θ ) (1.7.2) s_\theta(x)\coloneqq\frac{\partial }{\partial \theta}\ln p(x;\theta)=\frac{\partial_\theta p(x;\theta)}{p(x;\theta)}\tag{1.7.2} s θ ( x ) : = ∂ θ ∂ ln p ( x ; θ ) = p ( x ; θ ) ∂ θ p ( x ; θ ) ( 1.7.2 )
をスコア と呼ぶ。x x x を観測した人から見て、θ \theta θ を増やすと x x x の確率が相対的にどれだけ増えるかを表す。スコアの平均は 0 になる。実際、規格化 ∑ x p ( x ; θ ) = 1 \sum_xp(x;\theta)=1 ∑ x p ( x ; θ ) = 1 を θ \theta θ で微分すると
0 = ∑ x ∂ θ p ( x ; θ ) = ∑ x p ( x ; θ ) s θ ( x ) = ⟨ s θ ⟩ θ (1.7.3) 0=\sum_x\partial_\theta p(x;\theta)=\sum_xp(x;\theta)\,s_\theta(x)=\langle s_\theta \rangle _\theta\tag{1.7.3} 0 = x ∑ ∂ θ p ( x ; θ ) = x ∑ p ( x ; θ ) s θ ( x ) = ⟨ s θ ⟩ θ ( 1.7.3 )
となる。ここで正則条件の第二(微分と和の交換)と第一(和の範囲が θ \theta θ によらない)を使った。
D D D を2次まで展開する
D ( p θ ∥ p θ + δ ) = ∑ x p ( x ; θ ) [ ln p ( x ; θ ) − ln p ( x ; θ + δ ) ] D(p_\theta\,\|\,p_{\theta+\delta})=\sum_xp(x;\theta)\bigl[\ln p(x;\theta)-\ln p(x;\theta+\delta)\bigr] D ( p θ ∥ p θ + δ ) = ∑ x p ( x ; θ ) [ ln p ( x ; θ ) − ln p ( x ; θ + δ ) ] の ln p ( x ; θ + δ ) \ln p(x;\theta+\delta) ln p ( x ; θ + δ ) をδ \delta δ について Taylor 展開すると
ln p ( x ; θ + δ ) = ln p ( x ; θ ) + δ s θ ( x ) + δ 2 2 ∂ θ s θ ( x ) + O ( δ 3 ) \ln p(x;\theta+\delta)=\ln p(x;\theta)+\delta\,s_\theta(x)+\frac{\delta^2}{2}\,\partial_\theta s_\theta(x)+\mathcal{O}(\delta^3) ln p ( x ; θ + δ ) = ln p ( x ; θ ) + δ s θ ( x ) + 2 δ 2 ∂ θ s θ ( x ) + O ( δ 3 )
である。これを代入し、式 (1.7.3) で1次の項を消すと
D ( p θ ∥ p θ + δ ) = − δ 2 2 ⟨ ∂ θ s θ ⟩ θ + O ( δ 3 ) (1.7.4) D(p_\theta\,\|\,p_{\theta+\delta})=-\frac{\delta^2}{2}\langle \partial_\theta s_\theta \rangle _\theta+\mathcal{O}(\delta^3)\tag{1.7.4} D ( p θ ∥ p θ + δ ) = − 2 δ 2 ⟨ ∂ θ s θ ⟩ θ + O ( δ 3 ) ( 1.7.4 )
を得る。0次の項は D ( p θ ∥ p θ ) = 0 D(p_\theta\,\|\,p_\theta)=0 D ( p θ ∥ p θ ) = 0 で、1次の項は D D D が δ = 0 \delta=0 δ = 0 で最小値 0 をとることと合っている。2次の係数に現れた − ⟨ ∂ θ s θ ⟩ θ -\langle \partial_\theta s_\theta \rangle _\theta − ⟨ ∂ θ s θ ⟩ θ はスコアの2乗の平均に等しい。式 (1.7.3) をもう一度 θ \theta θ で微分すると
0 = ∂ θ ∑ x p s θ = ∑ x ( ∂ θ p ) s θ + ∑ x p ∂ θ s θ = ⟨ s θ 2 ⟩ θ + ⟨ ∂ θ s θ ⟩ θ 0=\partial_\theta\sum_xp\,s_\theta=\sum_x\bigl(\partial_\theta p\bigr)s_\theta+\sum_xp\,\partial_\theta s_\theta=\langle s_\theta^2 \rangle _\theta+\langle \partial_\theta s_\theta \rangle _\theta 0 = ∂ θ x ∑ p s θ = x ∑ ( ∂ θ p ) s θ + x ∑ p ∂ θ s θ = ⟨ s θ 2 ⟩ θ + ⟨ ∂ θ s θ ⟩ θ
となるからである(∂ θ p = p s θ \partial_\theta p=p\,s_\theta ∂ θ p = p s θ を使った)。
定義 1.7.1(Fisher 情報)
正則な族 { p θ } \{p_\theta\} { p θ } のFisher 情報 を
J ( θ ) ≔ ⟨ s θ 2 ⟩ θ = ⟨ ( ∂ ∂ θ ln p ( X ; θ ) ) 2 ⟩ θ = − ⟨ ∂ 2 ∂ θ 2 ln p ( X ; θ ) ⟩ θ (1.7.5) J(\theta)\coloneqq\langle s_\theta^2 \rangle _\theta=\langle \Bigl(\frac{\partial }{\partial \theta}\ln p(X;\theta)\Bigr)^2 \rangle _\theta=-\langle \frac{\partial^2}{\partial\theta^2}\ln p(X;\theta) \rangle _\theta\tag{1.7.5} J ( θ ) : = ⟨ s θ 2 ⟩ θ = ⟨ ( ∂ θ ∂ ln p ( X ; θ ) ) 2 ⟩ θ = − ⟨ ∂ θ 2 ∂ 2 ln p ( X ; θ ) ⟩ θ ( 1.7.5 )
で定める。スコアの平均が 0 なので、J ( θ ) = V a r θ ( s θ ) J(\theta)=\mathrm{Var}_\theta(s_\theta) J ( θ ) = Var θ ( s θ ) でもある。
Fisher 情報は D D D の2次の係数
式 (1.7.4) と定義 1.7.1 を合わせると
D ( p θ ∥ p θ + δ ) = 1 2 J ( θ ) δ 2 + O ( δ 3 ) (1.7.6) D(p_\theta\,\|\,p_{\theta+\delta})=\frac12J(\theta)\,\delta^2+\mathcal{O}(\delta^3)\tag{1.7.6} D ( p θ ∥ p θ + δ ) = 2 1 J ( θ ) δ 2 + O ( δ 3 ) ( 1.7.6 )
である。逆向きの D ( p θ + δ ∥ p θ ) D(p_{\theta+\delta}\,\|\,p_\theta) D ( p θ + δ ∥ p θ ) も同じ 1 2 J ( θ ) δ 2 \frac12J(\theta)\delta^2 2 1 J ( θ ) δ 2 で始まり、2つの向きは δ 3 \delta^3 δ 3 の項で初めて違う(例題 1.7.1 )。D D D は対称でないが、近くの2点の間では対称な2次形式になる。
名前は R. A. Fisher に由来する。Fisher はスコアの2乗の平均を推定の精度を測る量として導入した(R. A. Fisher (1925), Proc. Cambridge Philos. Soc. 22 , 700–725)。D D D を2つの分布の見分けやすさとして定式化したのは S. Kullback と R. A. Leibler である(S. Kullback and R. A. Leibler (1951), Ann. Math. Statist. 22 , 79–86)。統計学の教科書では Fisher 情報を I ( θ ) I(\theta) I ( θ ) と書くが、この巻では I I I を相互情報量に使うので J ( θ ) J(\theta) J ( θ ) と書く。
図 1.7.1 Bernoulli 分布(表の確率 θ = 0.3 \theta=0.3 θ = 0.3 )で、表の確率を θ + δ \theta+\delta θ + δ にずらした分布との相対エントロピー。点線と点は厳密な値、実線は2次の近似 1 2 J δ 2 \frac12J\delta^2 2 1 J δ 2 (J = 1 / ( 0.3 × 0.7 ) J=1/(0.3\times0.7) J = 1/ ( 0.3 × 0.7 ) )。∣ δ ∣ ≲ 0.1 \lvert \delta\rvert \lesssim0.1 ∣ δ ∣ ≲ 0.1 では2つはほぼ重なり、離れると3次の項のために左右が非対称になる。
近似の裏付け
型 B(極限) | 小さいもの δ J ( θ ) → 0 \delta\sqrt{J(\theta)}\to0 δ J ( θ ) → 0 (パラメータのずれが小さい) | 証明 条件付き(p p p が θ \theta θ について3回微分でき、微分と和を交換できる。剰余は O ( δ 3 ) \mathcal{O}(\delta^3) O ( δ 3 ) ) | 行き先 近似の数理 第1章
例題 1.7.1 Bernoulli 分布の Fisher 情報と D D D の向き ★★
オリジナル 序章の例題 1.1.3 で J = 1 / ( θ ( 1 − θ ) ) J=1/(\theta(1-\theta)) J = 1/ ( θ ( 1 − θ )) と 1 2 J δ 2 \frac12J\delta^2 2 1 J δ 2 の一致を見た。ここでは3次の係数まで調べる。表が確率 θ \theta θ で出るコイン(x = 1 x=1 x = 1 が表、x = 0 x=0 x = 0 が裏)を考える。
スコア s θ ( x ) s_\theta(x) s θ ( x ) と Fisher 情報 J ( θ ) J(\theta) J ( θ ) を求めよ。
D ( p θ ∥ p θ + δ ) D(p_\theta\,\|\,p_{\theta+\delta}) D ( p θ ∥ p θ + δ ) と D ( p θ + δ ∥ p θ ) D(p_{\theta+\delta}\,\|\,p_\theta) D ( p θ + δ ∥ p θ ) を δ \delta δ の3次まで展開し、2次の係数が等しく3次の係数が違うことを示せ。
θ = 0.3 \theta=0.3 θ = 0.3 、δ = 0.05 \delta=0.05 δ = 0.05 で2つの D D D と 1 2 J δ 2 \frac12J\delta^2 2 1 J δ 2 の値を比べよ。
解答 (1) ln p ( x ; θ ) = x ln θ + ( 1 − x ) ln ( 1 − θ ) \ln p(x;\theta)=x\ln\theta+(1-x)\ln(1-\theta) ln p ( x ; θ ) = x ln θ + ( 1 − x ) ln ( 1 − θ ) なので s θ ( x ) = x θ − 1 − x 1 − θ = x − θ θ ( 1 − θ ) s_\theta(x)=\dfrac{x}{\theta}-\dfrac{1-x}{1-\theta}=\dfrac{x-\theta}{\theta(1-\theta)} s θ ( x ) = θ x − 1 − θ 1 − x = θ ( 1 − θ ) x − θ である。⟨ ( X − θ ) 2 ⟩ θ = θ ( 1 − θ ) \langle (X-\theta)^2 \rangle _\theta=\theta(1-\theta) ⟨( X − θ ) 2 ⟩ θ = θ ( 1 − θ ) を使うと
J ( θ ) = θ ( 1 − θ ) θ 2 ( 1 − θ ) 2 = 1 θ ( 1 − θ ) J(\theta)=\frac{\theta(1-\theta)}{\theta^2(1-\theta)^2}=\frac{1}{\theta(1-\theta)} J ( θ ) = θ 2 ( 1 − θ ) 2 θ ( 1 − θ ) = θ ( 1 − θ ) 1
となる。θ \theta θ が 0 か 1 に近いほど J J J は大きい。表がほとんど出ないコインでは、表の確率がわずかに変わっただけでも表の出る頻度の比が大きく変わり、見分けやすいからである。
(2) D ( p θ ∥ p θ + δ ) = θ ln θ θ + δ + ( 1 − θ ) ln 1 − θ 1 − θ − δ D(p_\theta\,\|\,p_{\theta+\delta})=\theta\ln\frac{\theta}{\theta+\delta}+(1-\theta)\ln\frac{1-\theta}{1-\theta-\delta} D ( p θ ∥ p θ + δ ) = θ ln θ + δ θ + ( 1 − θ ) ln 1 − θ − δ 1 − θ を展開すると
D ( p θ ∥ p θ + δ ) = δ 2 2 θ ( 1 − θ ) − ( 1 − 2 θ ) δ 3 3 θ 2 ( 1 − θ ) 2 + O ( δ 4 ) D(p_\theta\,\|\,p_{\theta+\delta})=\frac{\delta^2}{2\theta(1-\theta)}-\frac{(1-2\theta)\,\delta^3}{3\theta^2(1-\theta)^2}+\mathcal{O}(\delta^4) D ( p θ ∥ p θ + δ ) = 2 θ ( 1 − θ ) δ 2 − 3 θ 2 ( 1 − θ ) 2 ( 1 − 2 θ ) δ 3 + O ( δ 4 )
となり、逆向きは
D ( p θ + δ ∥ p θ ) = δ 2 2 θ ( 1 − θ ) − ( 1 − 2 θ ) δ 3 6 θ 2 ( 1 − θ ) 2 + O ( δ 4 ) D(p_{\theta+\delta}\,\|\,p_\theta)=\frac{\delta^2}{2\theta(1-\theta)}-\frac{(1-2\theta)\,\delta^3}{6\theta^2(1-\theta)^2}+\mathcal{O}(\delta^4) D ( p θ + δ ∥ p θ ) = 2 θ ( 1 − θ ) δ 2 − 6 θ 2 ( 1 − θ ) 2 ( 1 − 2 θ ) δ 3 + O ( δ 4 )
となる。2次の係数はどちらも 1 2 J ( θ ) \frac12J(\theta) 2 1 J ( θ ) で、3次の係数は2倍違う(θ = 1 2 \theta=\frac12 θ = 2 1 のときだけ3次の項が両方消える)。
(3) J = 1 / 0.21 = 4.762 J=1/0.21=4.762 J = 1/0.21 = 4.762 なので 1 2 J δ 2 = 5.95 × 10 − 3 \frac12J\delta^2=5.95\times10^{-3} 2 1 J δ 2 = 5.95 × 1 0 − 3 である。厳密な値は D ( p 0.3 ∥ p 0.35 ) = 5.63 × 10 − 3 D(p_{0.3}\,\|\,p_{0.35})=5.63\times10^{-3} D ( p 0.3 ∥ p 0.35 ) = 5.63 × 1 0 − 3 、D ( p 0.35 ∥ p 0.3 ) = 5.78 × 10 − 3 D(p_{0.35}\,\|\,p_{0.3})=5.78\times10^{-3} D ( p 0.35 ∥ p 0.3 ) = 5.78 × 1 0 − 3 で、どちらも2次の近似と5%程度の差に収まる。差は (2) の3次の項でほぼ説明でき、たとえば前者は 5.95 × 10 − 3 − 3.02 × 0.05 3 = 5.57 × 10 − 3 5.95\times10^{-3}-3.02\times0.05^3=5.57\times10^{-3} 5.95 × 1 0 − 3 − 3.02 × 0.0 5 3 = 5.57 × 1 0 − 3 となる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com