1.3 相互情報量

1.2 節の終わりで,「を知って減るの不確かさ」について対称な量であることを見た.対称であるという事実は,それ自体が驚きに値する.の側から書いた式なのに,を入れ替えても値が変わらない.つまりこれは「から受け取る量」でも「から受け取る量」でもなく,二つのあいだで共有されている量である.共有量であるならば,を主語にしたという書き方をやめ,を対等に扱う形で直接定義したい.本節ではその定義を与え,それが確かにに一致することを確かめる.

定義

定義 1.3.1(相互情報量). 同時分布,周辺分布をもつについて,

𝐼(𝑋;𝑌)=𝑥,𝑦𝑝(𝑥,𝑦)log𝑝(𝑥,𝑦)𝑝(𝑥)𝑝(𝑦).

和はの項についてとる(の項は 0 と約束する.なおならなので,対数の中身はつねに正である).

この式は次のように読める.もしが独立なら同時分布はに一致し,対数の中身は全点で,したがってになる.一般には比からどれだけずれるかが各点での「独立からの隔たり」で,その対数を同時分布そのもので平均したのが相互情報量である.つまりは「が独立からどれだけ隔たっているか」を一つの数にまとめた量である.

定義式をもう少し読み込んでおこう.対数の中身はとも書ける.これは「を見たことで,が起きるという見込みが何倍になったか」である.1 倍(見込みが変わらない)なら対数は 0 で寄与なし,見込みが上がれば正,下がれば負の寄与をする.相互情報量は,この見込みの更新率の対数を実際に起きるで平均した量にほかならない.1.1 節の「驚き」の言葉でいえば,を知る前のへの驚きから,知った後の驚きを引いた「減った驚き」の平均である.

個々のでの寄与は負にもなりうる(誤解を招く観測はありうる)ことに注意したい.それでも平均すれば必ず非負になる,というのが次の命題で,1.2 節で「平均的に」を強調したのと同じ構図である.

記法の先取り. この右辺は,1.6 節で導入する相対エントロピー(KL ダイバージェンス)を使えばと一行で書ける.本節はこの記法を使わずに進むので,いまは読み飛ばしてよい.

形式化上の注記. 本ライブラリは相互情報量を,いま先取りした「同時分布と周辺積との相対エントロピー」の形でそのまま定義する(Mathlib の klDiv を用いる).この量は(拡張非負実数)値なので,相互情報量も同じ型をとる.エントロピー表現(定理 1.3.4)との橋渡しでは .toReal で実数に落とす一手間が入るが,これは型をまたぐ技術処理で,数学的内容は変わらない.

形式化: mutualInfo (ソース)

基本性質

命題 1.3.2. 𝐼(𝑋;𝑌) 0,かつが独立であることと同値.

証明. ここでも道具は 1.1 節の有限 Jensen(補題 1.1.9)ひとつである.今度は凹関数としてをとる(補題 1.1.6).とおき,上で

𝑡(𝑥,𝑦):=𝑝(𝑥)𝑝(𝑦)𝑝(𝑥,𝑦)>0

と定める.重み上で総和 1 だから,補題 1.1.9 を重み・点・凹関数に適用して

𝐼(𝑋;𝑌)=𝑆𝑝(𝑥,𝑦)log𝑡(𝑥,𝑦)log(𝑆𝑝(𝑥,𝑦)𝑡(𝑥,𝑦))=log(𝑆𝑝(𝑥)𝑝(𝑦))log1=0.

最後の不等号はによる.よって

等号を調べる.なら上の二つのがともに等号である.狭義凹だから,第 1 の等号は補題 1.1.9 の等号条件により「上でが定数」を意味する.その定数をとすると上で上で足して.第 2 の等号はを意味するので,すなわち上で.さらにから,の外ではであり,そこではでもあるから両辺は一致する.以上より全点で,すなわちは独立.逆に独立なら対数の中身が全点で 1 なので

非負性は「観測は平均的には損にならない」ということである.個々の観測は誤解を招きうる(さきほど見たとおり寄与が負のはある)が,平均すれば必ず得になる.等号条件はそれを裏から言ったもので,「まったく得にならない」のは「本当に無関係だった」場合に限る.つまりは独立性の定量的な指標である:0 なら独立,大きいほど強く結びついている.

形式化上の注記. 本文と形式化で難所が入れ替わる例である.形式化では相互情報量が拡張非負実数値なので非負性は型から従い,本文が手をかけた凹性の議論は現れない.代わりに実質的な内容が等号条件の側独立)に集まる.mutualInfo_eq_zero_iff_indep は,相対エントロピーの等号条件(klDiv_eq_zero_iff)と独立の特徴づけ(indepFun_iff_map_prod_eq_prod_map_map)を貼り合わせるだけになる.非自明さは Mathlib 側のこの二つが担っているからである.

形式化: 非負性 mutualInfo_nonneg,独立との同値 mutualInfo_eq_zero_iff_indep (ソース)

命題 1.3.3(対称性). 𝐼(𝑋;𝑌) =𝐼(𝑌;𝑋)

証明. 定義式のの入れ替えで不変だから.

形式化上の注記. 見れば自明な対称性が,形式化では測度同型MeasurableEquiv.prodComm に沿った像測度の移送になり,「相対エントロピーは測度同型で値を保つ」という自作補題 klDiv_map_measurableEquiv を要する.

形式化: mutualInfo_comm (ソース)

エントロピーとの関係

定理 1.3.4.

𝐼(𝑋;𝑌)=𝐻(𝑋)𝐻(𝑋𝑌)=𝐻(𝑌)𝐻(𝑌𝑋)=𝐻(𝑋)+𝐻(𝑌)𝐻(𝑋,𝑌).

これが相互情報量の「正体」を示す中心的な等式群である.最初の表現は「を知ることで減るの不確かさ」,対称形は「別々に測った不確かさの和から,まとめて測った不確かさを引いた重複分」と読める.定義 1.3.1 という一見別物の式が,1.2 節で導いたとぴったり一致する.これで,節の冒頭に掲げた「共有量を対等な形で定義する」という目標が達成されたことになる.

面積として読む. 二つの集合の面積として描くと見通しがよい.を二つの円の面積,を和集合の面積とみなすと,はちょうど共通部分にあたる.の円から共通部分を除いた部分,はその逆である.チェイン則は「和集合 = 一方 + 他方の残り」,は包除原理そのものになっている.この絵は 2 変数では正確で,実際に等式をすべて言い当てる(ただし 3 変数以上では「共通部分」が負になりうるので,絵として素直には延長できない).

両端で確かめる. が独立ならなので,共通部分なし.逆にならなので,二つの円が完全に重なる.相互情報量の値は必ずこののあいだにおさまる.

証明. 定義 1.3.1 の対数をを使って分解する:

𝐼(𝑋;𝑌)=𝑥,𝑦𝑝(𝑥,𝑦)log𝑝(𝑥𝑦)𝑝(𝑥)=𝑥,𝑦𝑝(𝑥,𝑦)log𝑝(𝑥)________=𝐻(𝑋)(𝑥,𝑦𝑝(𝑥,𝑦)log𝑝(𝑥𝑦))____________=𝐻(𝑋𝑌).

第 1 項はについて先に和をとって,第 2 項は定義 1.2.2 により.よって.対称性(命題 1.3.3)からも従う.最後にチェイン則(定理 1.2.3)でを代入すれば

いま見たという特別な場合から,相互情報量はしばしば「自己情報量」としてのエントロピーを含む一般化とみなされる.エントロピーは「自分自身と共有している情報の量」だ,というわけである.

形式化上の注記. 形式化ではが拡張非負実数値,が実数値なので,橋渡し定理は左辺に .toReal を付けた等式として述べる.数学的内容は定理 1.3.4 と同一.本文の「について和をとって」は,共通注記のとおりの像測度の上の積分補題に対応する.

形式化: mutualInfo_eq_entropy_sub_condEntropy (ソース),対称形mutualInfo_eq_entropy_add_entropy_sub_jointEntropy (ソース)

有限性

命題 1.3.5. 定義 1.1.1 の設定(有限アルファベット)で

証明. 定理 1.3.4 よりで,は有限個の項の有限和だから有限,は非負(定義 1.2.2 の直後で見た).よって

拡張実数値をとりうる量を扱うとき,有限性はいちいち確かめる必要がある.以降ではこの命題により,相互情報量を実数として足し引きしてよい.

形式化上の注記. 本文一行の有限性が,形式化では結合分布の周辺積への絶対連続性と対数尤度比の可積分性を経て Mathlib の klDiv_ne_top に帰着する測度論の補題になる.

形式化: mutualInfo_ne_top (ソース)

InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.