1.3 相互情報量
1.2 節の終わりで,「𝑌を知って減る𝑋の不確かさ」𝐻(𝑋) −𝐻(𝑋 ∣𝑌)が𝑋,𝑌について対称な量であることを見た.対称であるという事実は,それ自体が驚きに値する.𝐻(𝑋) −𝐻(𝑋 ∣𝑌)は𝑋の側から書いた式なのに,𝑋と𝑌を入れ替えても値が変わらない.つまりこれは「𝑋が𝑌から受け取る量」でも「𝑌が𝑋から受け取る量」でもなく,二つのあいだで共有されている量である.共有量であるならば,𝑋を主語にした𝐻(𝑋) −𝐻(𝑋 ∣𝑌)という書き方をやめ,𝑋と𝑌を対等に扱う形で直接定義したい.本節ではその定義を与え,それが確かに𝐻(𝑋) −𝐻(𝑋 ∣𝑌)に一致することを確かめる.
定義
定義 1.3.1(相互情報量). 同時分布𝑝(𝑥,𝑦),周辺分布𝑝(𝑥),𝑝(𝑦)をもつ(𝑋,𝑌)について,
𝐼(𝑋;𝑌)=∑𝑥,𝑦𝑝(𝑥,𝑦)log𝑝(𝑥,𝑦)𝑝(𝑥)𝑝(𝑦).和は𝑝(𝑥,𝑦) >0の項についてとる(𝑝(𝑥,𝑦) =0の項は 0 と約束する.なお𝑝(𝑥,𝑦) >0なら𝑝(𝑥) ≥𝑝(𝑥,𝑦) >0,𝑝(𝑦) >0なので,対数の中身はつねに正である).
この式は次のように読める.もし𝑋と𝑌が独立なら同時分布は𝑝(𝑥)𝑝(𝑦)に一致し,対数の中身は全点で1,したがって𝐼(𝑋;𝑌) =0になる.一般には比𝑝(𝑥,𝑦)/(𝑝(𝑥)𝑝(𝑦))が1からどれだけずれるかが各点での「独立からの隔たり」で,その対数を同時分布そのもので平均したのが相互情報量である.つまり𝐼(𝑋;𝑌)は「𝑋と𝑌が独立からどれだけ隔たっているか」を一つの数にまとめた量である.
定義式をもう少し読み込んでおこう.対数の中身は𝑝(𝑥,𝑦)𝑝(𝑥)𝑝(𝑦) =𝑝(𝑥∣𝑦)𝑝(𝑥)とも書ける.これは「𝑦を見たことで,𝑥が起きるという見込みが何倍になったか」である.1 倍(見込みが変わらない)なら対数は 0 で寄与なし,見込みが上がれば正,下がれば負の寄与をする.相互情報量は,この見込みの更新率の対数を実際に起きる(𝑥,𝑦)で平均した量にほかならない.1.1 節の「驚き−log𝑝」の言葉でいえば,𝑦を知る前の𝑥への驚き−log𝑝(𝑥)から,知った後の驚き−log𝑝(𝑥 ∣𝑦)を引いた「減った驚き」の平均である.
個々の(𝑥,𝑦)での寄与は負にもなりうる(誤解を招く観測はありうる)ことに注意したい.それでも平均すれば必ず非負になる,というのが次の命題で,1.2 節で「平均的に」を強調したのと同じ構図である.
記法の先取り. この右辺は,1.6 節で導入する相対エントロピー(KL ダイバージェンス)𝐷( ⋅ ‖ ⋅)を使えば𝐼(𝑋;𝑌) =𝐷(𝑝(𝑥,𝑦) ∥ 𝑝(𝑥)𝑝(𝑦))と一行で書ける.本節はこの記法を使わずに進むので,いまは読み飛ばしてよい.
基本性質
命題 1.3.2. 𝐼(𝑋;𝑌) ≥0,かつ𝐼(𝑋;𝑌) =0は𝑋と𝑌が独立であることと同値.
証明. ここでも道具は 1.1 節の有限 Jensen(補題 1.1.9)ひとつである.今度は凹関数としてlogをとる(補題 1.1.6).𝑆 :={(𝑥,𝑦) :𝑝(𝑥,𝑦) >0}とおき,𝑆上で
𝑡(𝑥,𝑦):=𝑝(𝑥)𝑝(𝑦)𝑝(𝑥,𝑦)>0と定める.重み𝑝(𝑥,𝑦)は𝑆上で総和 1 だから,補題 1.1.9 を重み𝑝(𝑥,𝑦)・点𝑡(𝑥,𝑦)・凹関数logに適用して
−𝐼(𝑋;𝑌)=∑𝑆𝑝(𝑥,𝑦)log𝑡(𝑥,𝑦)≤log(∑𝑆𝑝(𝑥,𝑦)𝑡(𝑥,𝑦))=log(∑𝑆𝑝(𝑥)𝑝(𝑦))≤log1=0.最後の不等号は∑𝑆𝑝(𝑥)𝑝(𝑦) ≤∑𝑥,𝑦𝑝(𝑥)𝑝(𝑦) =1による.よって𝐼(𝑋;𝑌) ≥0.
等号を調べる.𝐼(𝑋;𝑌) =0なら上の二つの≤がともに等号である.logは狭義凹だから,第 1 の等号は補題 1.1.9 の等号条件により「𝑆上で𝑡(𝑥,𝑦)が定数」を意味する.その定数を𝑐とすると𝑆上で𝑝(𝑥)𝑝(𝑦) =𝑐 𝑝(𝑥,𝑦),𝑆上で足して∑𝑆𝑝(𝑥)𝑝(𝑦) =𝑐.第 2 の等号は∑𝑆𝑝(𝑥)𝑝(𝑦) =1を意味するので𝑐 =1,すなわち𝑆上で𝑝(𝑥,𝑦) =𝑝(𝑥)𝑝(𝑦).さらに∑𝑆𝑝(𝑥)𝑝(𝑦) =1 =∑𝑥,𝑦𝑝(𝑥)𝑝(𝑦)から,𝑆の外では𝑝(𝑥)𝑝(𝑦) =0であり,そこでは𝑝(𝑥,𝑦) =0でもあるから両辺は一致する.以上より全点で𝑝(𝑥,𝑦) =𝑝(𝑥)𝑝(𝑦),すなわち𝑋と𝑌は独立.逆に独立なら対数の中身が全点で 1 なので𝐼(𝑋;𝑌) =0.◼
非負性は「観測は平均的には損にならない」ということである.個々の観測は誤解を招きうる(さきほど見たとおり寄与が負の(𝑥,𝑦)はある)が,平均すれば必ず得になる.等号条件はそれを裏から言ったもので,「まったく得にならない」のは「本当に無関係だった」場合に限る.つまり𝐼(𝑋;𝑌)は独立性の定量的な指標である:0 なら独立,大きいほど強く結びついている.
命題 1.3.3(対称性). 𝐼(𝑋;𝑌) =𝐼(𝑌;𝑋).
証明. 定義式の𝑝(𝑥,𝑦)𝑝(𝑥)𝑝(𝑦)は𝑥,𝑦の入れ替えで不変だから.◼
エントロピーとの関係
定理 1.3.4.
𝐼(𝑋;𝑌)=𝐻(𝑋)−𝐻(𝑋∣𝑌)=𝐻(𝑌)−𝐻(𝑌∣𝑋)=𝐻(𝑋)+𝐻(𝑌)−𝐻(𝑋,𝑌).
これが相互情報量の「正体」を示す中心的な等式群である.最初の表現は「𝑌を知ることで減る𝑋の不確かさ」,対称形𝐻(𝑋) +𝐻(𝑌) −𝐻(𝑋,𝑌)は「別々に測った不確かさの和から,まとめて測った不確かさを引いた重複分」と読める.定義 1.3.1 という一見別物の式が,1.2 節で導いた𝐻(𝑋) −𝐻(𝑋 ∣𝑌)とぴったり一致する.これで,節の冒頭に掲げた「共有量を対等な形で定義する」という目標が達成されたことになる.
面積として読む. 二つの集合の面積として描くと見通しがよい.𝐻(𝑋)と𝐻(𝑌)を二つの円の面積,𝐻(𝑋,𝑌)を和集合の面積とみなすと,𝐼(𝑋;𝑌)はちょうど共通部分にあたる.𝐻(𝑋 ∣𝑌)は𝑋の円から共通部分を除いた部分,𝐻(𝑌 ∣𝑋)はその逆である.チェイン則𝐻(𝑋,𝑌) =𝐻(𝑌) +𝐻(𝑋 ∣𝑌)は「和集合 = 一方 + 他方の残り」,𝐼 =𝐻(𝑋) +𝐻(𝑌) −𝐻(𝑋,𝑌)は包除原理そのものになっている.この絵は 2 変数では正確で,実際に等式をすべて言い当てる(ただし 3 変数以上では「共通部分」が負になりうるので,絵として素直には延長できない).
両端で確かめる. 𝑋と𝑌が独立なら𝐻(𝑋 ∣𝑌) =𝐻(𝑋)なので𝐼(𝑋;𝑌) =0,共通部分なし.逆に𝑌 =𝑋なら𝐻(𝑋 ∣𝑋) =0なので𝐼(𝑋;𝑋) =𝐻(𝑋),二つの円が完全に重なる.相互情報量の値は必ずこの0とmin{𝐻(𝑋),𝐻(𝑌)}のあいだにおさまる.
証明. 定義 1.3.1 の対数を𝑝(𝑥,𝑦) =𝑝(𝑥 ∣𝑦)𝑝(𝑦)を使って分解する:
𝐼(𝑋;𝑌)=∑𝑥,𝑦𝑝(𝑥,𝑦)log𝑝(𝑥∣𝑦)𝑝(𝑥)=−∑𝑥,𝑦𝑝(𝑥,𝑦)log𝑝(𝑥)⏟____⏟____⏟=𝐻(𝑋)−(−∑𝑥,𝑦𝑝(𝑥,𝑦)log𝑝(𝑥∣𝑦))⏟______⏟______⏟=𝐻(𝑋∣𝑌).第 1 項は𝑦について先に和をとって𝐻(𝑋),第 2 項は定義 1.2.2 により𝐻(𝑋 ∣𝑌).よって𝐼(𝑋;𝑌) =𝐻(𝑋) −𝐻(𝑋 ∣𝑌).対称性(命題 1.3.3)から𝐼(𝑋;𝑌) =𝐻(𝑌) −𝐻(𝑌 ∣𝑋)も従う.最後にチェイン則(定理 1.2.3)で𝐻(𝑋 ∣𝑌) =𝐻(𝑋,𝑌) −𝐻(𝑌)を代入すれば𝐼(𝑋;𝑌) =𝐻(𝑋) +𝐻(𝑌) −𝐻(𝑋,𝑌).◼
いま見た𝐼(𝑋;𝑋) =𝐻(𝑋)という特別な場合から,相互情報量はしばしば「自己情報量」としてのエントロピーを含む一般化とみなされる.エントロピーは「自分自身と共有している情報の量」だ,というわけである.
有限性
命題 1.3.5. 定義 1.1.1 の設定(有限アルファベット)で𝐼(𝑋;𝑌) <∞.
証明. 定理 1.3.4 より𝐼(𝑋;𝑌) =𝐻(𝑋) −𝐻(𝑋 ∣𝑌)で,𝐻(𝑋)は有限個の項の有限和だから有限,𝐻(𝑋 ∣𝑌)は非負(定義 1.2.2 の直後で見た).よって𝐼(𝑋;𝑌) ≤𝐻(𝑋) <∞.◼
拡張実数値をとりうる量を扱うとき,有限性はいちいち確かめる必要がある.以降ではこの命題により,相互情報量を実数として足し引きしてよい.
InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.