16.7 エントロピーべきと Fisher 情報量
本節からは第7章の微分エントロピーの上に乗る.微分エントロピーℎ(𝑋)(定義 7.1.2)は目盛の取り替えで動く量で,𝑋を𝑐倍するとlog|𝑐|だけ増えた(命題 7.3.3).散らばりを測るもう一つの量である分散は,𝑐倍すると𝑐2倍になる(これを,学部の確率で扱う事実として既知とする).二つを同じ物差しに載せるには,ℎを2倍して指数に持ち上げればよい.そうして作った量なら𝑐倍で𝑐2倍になるはずで,実際そうなる(命題 16.7.2).本節はその量を定めてガウス分布と一様分布で計算し,分散との大小を決める.続いて,密度の傾きの大きさを測る Fisher 情報量を定め,こちらもガウス分布で計算する.二つの量は,独立なガウス雑音を少しずつ足していく道筋の上で結びつく.
16.4 節 で既知とした,logが正の実数の全体から実数の全体への狭義単調な全単射であることから,どの実数についても,その値を対数にもつ正の実数はただ一つ定まる.
定義 16.7.1(エントロピーべき). 𝑋を密度をもつ実数値確率変数とし,ℎ(𝑋)(定義 7.1.2)が定まるとする.log𝑁(𝑋) =2 ℎ(𝑋)を満たすただ一つの正の実数𝑁(𝑋)を,𝑋の エントロピーべき と呼ぶ.
𝑁はエントロピーべきを表す.量そのものを名指すときのほかは引数をとり,第8章 が雑音の分散に使う裸の𝑁とは別の記号である.定める式は,底を𝑏にとれば𝑁(𝑋) =𝑏2ℎ(𝑋)と書き直せる.底を取り替えても𝑁(𝑋)の値は動かない.取り替えによってℎもlogも同じ倍率で変わり,定める式の両辺が同じ倍率で変わるからである.
命題 16.7.2(定数倍でのふるまい). 𝑋を密度をもつ実数値確率変数,𝑐 ≠0を実数とする.ℎ(𝑋)が定まるならばℎ(𝑐𝑋)も定まり,𝑁(𝑐𝑋) =𝑐2 𝑁(𝑋)である.
証明. 命題 7.3.3 よりℎ(𝑐𝑋)は定まってℎ(𝑐𝑋) =ℎ(𝑋) +log|𝑐|である.よって
log𝑁(𝑐𝑋)=2ℎ(𝑐𝑋)=2ℎ(𝑋)+2log|𝑐|=log𝑁(𝑋)+log𝑐2であり,右辺はlog(𝑐2𝑁(𝑋))に等しい(|𝑐|2 =𝑐2である).定義 16.7.1 は𝑁(𝑐𝑋)を,対数が2ℎ(𝑐𝑋)に等しいただ一つの正の実数として定めているから,𝑁(𝑐𝑋) =𝑐2𝑁(𝑋)である.◼
命題 16.7.2 が,ℎを2倍して指数に持ち上げた理由である.分散も𝑐倍で𝑐2倍になるので,𝑁と分散は同じ倍率で動く.二つの大小を比べる意味が出てくるのはそのためで,比べた結果は 命題 16.7.5 で出る.
例 16.7.3(一様分布のエントロピーべき). 𝑤 >0とし,𝑋が長さ𝑤の区間の上の一様分布に従うとすると,ℎ(𝑋)は定まり𝑁(𝑋) =𝑤2である.
証明. 例 7.1.3 よりℎ(𝑋)は定まってℎ(𝑋) =log𝑤である.よって2ℎ(𝑋) =log𝑤2であり,定義 16.7.1 より𝑁(𝑋) =𝑤2である.◼
命題 16.7.4(ガウス分布のエントロピーべき). 𝜇を実数,𝜎2 >0とし,𝑋がN(𝜇,𝜎2)(定義 7.2.1)に従うとすると,ℎ(𝑋)は定まり𝑁(𝑋) =2𝜋𝑒 𝜎2である.
証明. 定理 7.2.2 よりℎ(𝑋)は定まってℎ(𝑋) =12log(2𝜋𝑒𝜎2)である.よって2ℎ(𝑋) =log(2𝜋𝑒𝜎2)であり,定義 16.7.1 より主張を得る.◼
幅の2乗として読む. 7.1 節 は2ℎ(𝑋)をその分布の実効的な幅と読んだ.底を2にとれば𝑁(𝑋) =22ℎ(𝑋)はその幅の2乗である.例 16.7.3 の一様分布では幅が𝑤そのもので𝑁(𝑋) =𝑤2になり,7.2 節 が幅を√2𝜋𝑒 𝜎と読んだガウス分布では,命題 16.7.4 のとおり𝑁(𝑋) =2𝜋𝑒𝜎2になる.長さの2乗の次元をもつという点で,𝑁は分散と同じ物差しの上にある.ただしガウス分布では,𝑁(𝑋)は分散そのものではなく2𝜋𝑒倍になる.定義 16.7.1 がℎを指数に持ち上げるだけで,ガウス分布に値をそろえる定数を入れていないからである.Cover–Thomas は2𝜋𝑒で割った量を同じ名前で呼ぶ流儀で,そちらではガウス分布での値が分散そのものになる.どちらをとっても 命題 16.7.2 の𝑐2倍のふるまいは変わらないので,本書は定義に定数を持ち込まないほうをとる.
分散を固定したときにガウス分布が微分エントロピーを最大にする(定理 7.2.4)ことは,指数に持ち上げてもそのまま残る.
命題 16.7.5(分散による上界). 𝜎2 >0とし,𝑋を密度をもつ実数値確率変数とする.ℎ(𝑋)が定まり,𝑋が平均をもち,その分散が定まって𝜎2以下であるならば𝑁(𝑋) ≤2𝜋𝑒 𝜎2である.
証明. 定理 7.2.4 よりℎ(𝑋) ≤12log(2𝜋𝑒𝜎2)である.両辺を2倍すると,定義 16.7.1 より左辺はlog𝑁(𝑋)であり,右辺はlog(2𝜋𝑒𝜎2)である.logが狭義単調(16.4 節)だから主張を得る.◼
この上界が達成されることを見るために,第7章 7.2 節 が借りた ガウス分布の基本性質 をそのまま引く.すなわち,定義 7.2.1 の関数は実軸の上での積分が1であり,それを密度とする確率変数の平均は𝜇,分散は𝜎2であって,∫ℝ(𝑥 −𝜇)2 𝑔𝜇,𝜎2(𝑥) 𝑑𝑥 =𝜎2である.最後の形は,7.2 節 が 定理 7.2.2 の証明で使ったものと同じである.当てる対象は 定義 7.2.1 のガウス分布に従う実数値確率変数だけで,7.2 節 が宣言した射程からは広がらない.依存するのは 系 16.7.6 と 命題 16.7.8 の証明の二つである.
系 16.7.6. 𝜇を実数,𝜎2 >0とし,𝑋がN(𝜇,𝜎2)(定義 7.2.1)に従うとすると,𝑋は 命題 16.7.5 の仮定をこの𝜎2について満たし,その不等式は等号で成り立つ.
一様分布では等号にならない.例 7.2.6 の証明が求めたとおり,長さ𝑤の区間の上の一様分布の分散は𝑤2/12である.命題 16.7.5 をこの値でとると上界は2𝜋𝑒 𝑤2/12 =(𝜋𝑒/6) 𝑤2 =1.4232… 𝑤2で,例 16.7.3 の𝑁(𝑋) =𝑤2を確かに上から抑える.比𝜋𝑒/6は,例 7.2.6 が出した微分エントロピーの差12log𝜋𝑒6を指数に持ち上げたものである.
ここまでは,密度がどれだけ広がっているかを測ってきた.次に測るのは,密度がどれだけ急に変わるかである.
定義 16.7.7(Fisher 情報量). 𝑋を密度𝑓をもつ実数値確率変数とし,𝑓は実軸のすべての点で微分可能かつ正であるとする.∫ℝ(𝑓′(𝑥)𝑓(𝑥))2𝑓(𝑥) 𝑑𝑥が有限であるとき,その値を𝑋の Fisher 情報量 と呼んで𝐽(𝑋)と書き,このとき𝐽(𝑋)が 定まる という.
𝐽は Fisher 情報量を表す.この文字は既出の章がそれぞれ別の意味に使っているが,本節のものはそのどれとも別である.定義 16.7.7 は密度を一つ選んだうえでの量だが,値は選び方によらない.定義 7.1.1 より同じ𝑋の二つの密度はどの有界閉区間の上でも積分が等しく,定義 16.7.7 の条件を満たす密度は微分可能だから連続であって,どの有界閉区間の上でも積分が等しい二つの連続関数は全点で一致する(最後の一歩を,学部の解析で扱う事実として既知とする).被積分関数に現れる𝑓′/𝑓は,密度の増え方をその点の密度自身で割ったもので,密度が急に立ち上がっている点ほど絶対値が大きい.その2乗を,やはり密度自身で平均している.広がりを測る 定義 16.7.1 と違って,𝐽が見ているのは密度の形の細かさである.この量を情報量と呼ぶのは,密度を横に少しずらしたときの見分けやすさを測るからである.形は同じで位置だけが分かっていない分布から標本を得て,その位置を当てにいくとしよう.位置を少し動かしたときに密度が大きく変わる分布ほど,どこにいるかは当てやすい.統計でその当てやすさを測るのに現れるのが同じ量で,定義 16.7.7 はそれを,位置がずれるだけの場合について密度の形で書いたものである.
計算のために 微積分の計算規則 を既知とする.積分の線形性,合成関数の微分と指数関数の導関数,log𝑡の導関数が(log𝑒)/𝑡であること,1次式による置換積分,部分積分,𝑡2log𝑡が𝑡 ↓0で0に近づくこと,連続関数の定積分が原始関数の差に等しいこと(微積分の基本定理),および有限な値に収束する広義積分の扱いである.このうち1次式による置換積分は 第7章 7.3 節 が,積分の線形性・log𝑡の導関数・部分積分・𝑡2log𝑡の極限・広義積分の扱いは 7.1 節 が,同じ資格ですでに借りている.本節が新たに加えるのは,合成関数の微分と指数関数の導関数と微積分の基本定理である.当てる対象は,命題 16.7.8 と 例 16.8.2 の証明に現れる密度と,それに多項式や対数を掛けて作る関数,および 例 16.7.9 の証明が雑音の分散を変数として書く,1次式の対数とその逆数だけである.
命題 16.7.8(ガウス分布の Fisher 情報量). 𝜇を実数,𝜎2 >0とし,𝑋がN(𝜇,𝜎2)(定義 7.2.1)に従うとすると,𝐽(𝑋)は定まり𝐽(𝑋) =1/𝜎2である.
証明. 定義 7.2.1 の密度𝑔𝜇,𝜎2は実軸のすべての点で正であり,借りた微積分の計算規則より微分可能で
𝑔′𝜇,𝜎2(𝑥)=−𝑥−𝜇𝜎2𝑔𝜇,𝜎2(𝑥)である.よって 定義 16.7.7 の被積分関数は(𝑥−𝜇)2𝜎4 𝑔𝜇,𝜎2(𝑥)であり,借りたガウス分布の基本性質より
∫ℝ(𝑥−𝜇)2𝜎4𝑔𝜇,𝜎2(𝑥)𝑑𝑥=𝜎2𝜎4=1𝜎2である.これは有限だから𝐽(𝑋)は定まり,その値は1/𝜎2である.◼
二つの量を掛けてみる. ガウス分布では 命題 16.7.4 より𝑁(𝑋) =2𝜋𝑒𝜎2,命題 16.7.8 より𝐽(𝑋) =1/𝜎2だから,積は𝜎2に依らず2𝜋𝑒である.散らばりが大きいほど𝑁は大きく𝐽は小さい,という向きが,ガウス分布では正確に釣り合っている.一般の分布でも,二つが定まるかぎり,この積は2𝜋𝑒を下回らないことが知られている(エントロピーの等周不等式と呼ばれる).ガウス分布での値がちょうど2𝜋𝑒だから,ガウス分布は積を最小にする.
等周とあるのは,16.5 節 と同じ形の問いだからである.あちらは大きさにあたる要素数を決めて,境界にあたる辺の本数を最小にする問いで,要素数が2のべきなら,最小は部分立方体でとられた(系 16.5.9).こちらは大きさにあたるエントロピーべきを決めて,境界にあたる Fisher 情報量を最小にする問いで,最小はガウス分布でとられる.本書はこの不等式を示さず,形式化もしていない.本章のどの主張もこれを使わないので,借りる道具には数えず,名前を書くにとどめる.
二つの量が結びつくのは,分布を動かしたときである.動かし方として,独立なガウス雑音を少しずつ足していくもの,すなわち𝑠を0から増やしながら𝑋 +√𝑠 𝑍を見ていくものをとり,これを以下,熱流と呼ぶ.名前は,𝑠を時間とみると,この道筋の密度が熱の伝わり方と同じ形の方程式に従うことによる(本書はその方程式を扱わない).ここから次節でエントロピーべき不等式を借りるところまでは,その不等式がどこから来るかを見せるのが目的である.途中で借りる de Bruijn の恒等式と Stam の不等式は,本書のどの主張の証明も使わない.用意するものが二つある.
熱流を書くには,ガウス雑音を足したものがまたガウス分布に従うことが要る.第9章 9.4 節 が確率変数の言葉で借り直した ガウス分布の畳み込み を引く.独立でN(𝜇1,𝜎21)とN(𝜇2,𝜎22)(定義 7.2.1)に従う二つの実数値確率変数の和がN(𝜇1 +𝜇2, 𝜎21 +𝜎22)に従う,という事実である.9.4 節 は差についても借りたが,本章が使うのは和だけで,射程はそのぶん狭い.あわせて,独立な二つの実数値確率変数の一方に0でない定数を掛けても独立性が保たれることを既知とする.当てる対象は,本節と次節に現れるガウス分布に従う実数値確率変数だけである.依存するのは 例 16.7.9 と 命題 16.8.1 の証明の二つである.
de Bruijn の恒等式 を借りる.𝑋を密度をもつ実数値確率変数で2乗の平均をもつもの,𝑍をN(0,1)(定義 7.2.1)に従う確率変数とし,𝑋と𝑍は独立であるとする.このとき𝑠 >0について𝑋 +√𝑠 𝑍は密度をもち𝐽(𝑋 +√𝑠 𝑍)(定義 16.7.7)は定まって,次の二つが成り立つ,という主張である.第一に,𝑡 >0において𝑠 ↦ℎ(𝑋 +√𝑠 𝑍)は微分可能で,その𝑡での導関数はlog𝑒2 𝐽(𝑋 +√𝑡 𝑍)に等しい.第二に,さらにℎ(𝑋)が定まるならば,𝑡 ≥0について
ℎ(𝑋+√𝑡𝑍)−ℎ(𝑋)=log𝑒2∫𝑡0𝐽(𝑋+√𝑠𝑍)𝑑𝑠である.微分の形を0から𝑡まで積分すれば積分の形が出そうに見えるが,微分の形が言えているのは𝑡 >0だけなので,𝑠 ↓0の端をつなぐ段が別に要る.そこで二つを別々に借りる.当てる対象は,本節の𝑋 +√𝑠 𝑍と,16.8 節 がエントロピーべき不等式の出どころを述べる筋書きの中で当てる,密度をもつ独立な二つの実数値確率変数とその和である.本書のどの主張の証明もこれを使わない.例 16.7.9 がガウス分布の場合にこれを確かめ,16.8 節 がその筋書きの中で引く.証明を載せないのは,積分と微分の順序を入れ替える段が要り,本書がその道具を用意していないからである.
増え方が細かさで決まると読む. 恒等式が言っているのは,雑音を少し足したときに微分エントロピーが増える速さが,そのときの密度の細かさ𝐽で決まるということである.でこぼこの多い密度は,雑音を少し足すだけで大きく平らになるのでℎが大きく増え,すでに滑らかな密度は増えにくい.
例 16.7.9(熱流に沿った微分エントロピー). 𝜇を実数,𝜎2 >0とし,𝑋がN(𝜇,𝜎2)(定義 7.2.1)に従い,𝑍がN(0,1)に従い,𝑋と𝑍は独立であるとする.このとき𝑠 ≥0についてℎ(𝑋 +√𝑠 𝑍)は定まって12log(2𝜋𝑒(𝜎2 +𝑠))に等しく,次の二つが成り立つ.
- 𝑡 >0において,𝑠 ↦ℎ(𝑋 +√𝑠 𝑍)は微分可能で,その𝑡での導関数はlog𝑒2 𝐽(𝑋 +√𝑡 𝑍)に等しい.
- 𝑡 ≥0についてℎ(𝑋 +√𝑡 𝑍) −ℎ(𝑋) =log𝑒2∫𝑡0𝐽(𝑋 +√𝑠 𝑍) 𝑑𝑠である.
証明. 𝑠 >0とする.例 7.3.6 を,掛ける定数√𝑠,足す定数0ととって当てると√𝑠 𝑍はN(0,𝑠)に従い,ガウス分布の畳み込みとあわせて既知とした定数倍での独立性の保存より,𝑋と√𝑠 𝑍は独立である.借りたガウス分布の畳み込みより𝑋 +√𝑠 𝑍はN(𝜇, 𝜎2 +𝑠)に従い,定理 7.2.2 よりℎ(𝑋 +√𝑠 𝑍)は定まって12log(2𝜋𝑒(𝜎2 +𝑠))である.𝑠 =0のときは𝑋 +√𝑠 𝑍 =𝑋で,同じ式が 定理 7.2.2 そのものである.
1 を示す.借りた微積分の計算規則より,𝑠 ↦12log(2𝜋𝑒(𝜎2 +𝑠))は微分可能で,その𝑡での導関数はlog𝑒2(𝜎2+𝑡)である.一方 命題 16.7.8 をN(𝜇, 𝜎2 +𝑡)に当てると𝐽(𝑋 +√𝑡 𝑍) =1/(𝜎2 +𝑡)だから,二つは一致する.
2 を示す.借りた微積分の計算規則より𝑠 ↦1log𝑒log(𝜎2 +𝑠)の導関数は1𝜎2+𝑠であり,これは[0,𝑡]の上で連続だから,微積分の基本定理より
∫𝑡0𝑑𝑠𝜎2+𝑠=1log𝑒(log(𝜎2+𝑡)−log𝜎2)である.命題 16.7.8 より被積分関数は𝐽(𝑋 +√𝑠 𝑍)に等しいので,両辺にlog𝑒2を掛けると,右辺は12log(𝜎2 +𝑡) −12log𝜎2になる.いっぽう 2 の左辺は,上で求めた値より12log(2𝜋𝑒(𝜎2 +𝑡)) −12log(2𝜋𝑒𝜎2)であり,logの中の共通の因子2𝜋𝑒が引き算で消えて同じ値になる.◼
恒等式の両辺を数で突き合わせる. 𝜎2 =1,𝑡 =3とし,底を2にとる.例 16.7.9 の 2 の左辺は12log2(2𝜋𝑒 ⋅4) −12log2(2𝜋𝑒) =12log24 =1ビットである.右辺は,𝐽(𝑋 +√𝑠 𝑍) =1/(1 +𝑠)で∫30𝑑𝑠1+𝑠 =1log2𝑒log24 =2log2𝑒だから,log2𝑒2 ⋅2log2𝑒 =1ビットである.雑音の分散を0から3まで増やすあいだに,微分エントロピーはちょうど1ビット増える.
de Bruijn の恒等式は,微分エントロピーの増え方を Fisher 情報量で書いた.増え方が書けたということは,Fisher 情報量の側で成り立つ不等式を,積分して微分エントロピーの側へ運べるということである.次節はその運び先を見る.
InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.