7.5 連続な相対エントロピーと相互情報量
7.3 節で見たとおり,微分エントロピーの数値そのものは目盛の取り方で動く.意味をもつのは同じ目盛で測った量の差である.第1章で相対エントロピー𝐷(𝑝 ‖ 𝑞)と相互情報量𝐼(𝑋;𝑌)を導入したとき,どちらも二つの分布の比の対数を平均する形をしていた(1.6 節と定義 1.3.1).比の形は,分母と分子に同じ係数がかかる取り替えを見なくなる.本節はこの二つの量を密度で書き直し,離散のときと同じ性質が成り立つことを見る.
書き直したあとの相対エントロピーは,補題 7.2.3 が示した「ℎ(𝑋)と−∫𝑓log𝑞の差」にほかならない.したがって非負性はすでに手元にある.新しく要るのは定義と,離散のときの読み方がそのまま通ることの確認である.
相対エントロピー
定義 7.5.1(相対エントロピー). 𝑛 ≥1とし,𝑓と𝑞をℝ𝑛上の密度とする.𝑆 :={𝑥 ∈ℝ𝑛 :𝑓(𝑥) >0 かつ 𝑞(𝑥) >0}とおく.∫𝑆𝑓(𝑥) 𝑑𝑥 =1であり∫𝑆𝑓(𝑥)∣log𝑓(𝑥)𝑞(𝑥)∣ 𝑑𝑥 <∞であるとき,𝑓の𝑞に対する 相対エントロピー を
𝐷(𝑓‖𝑞):=∫𝑆𝑓(𝑥)log𝑓(𝑥)𝑞(𝑥)𝑑𝑥で定め,このとき𝐷(𝑓 ‖ 𝑞)が 定まる という.
∫𝑆𝑓 =1という条件は,𝑞が0になる点に𝑓が確率を置いていない,ということである.離散のときはその場合に𝐷(𝑝 ‖ 𝑞) = +∞と読んだ(1.6 節)が,本書は密度についてその場合の𝐷を定めない.+∞を値として認めると,以下の主張がどれも「両辺が有限なら」という但し書きを抱えることになり,得るものより手間のほうが大きいからである.𝑞を「起こらないと決めつけた想定」と読むなら,決めつけが外れたときの罰が無限大になるという事情は離散のときと同じである.定めない場合は実際にある.例 7.2.6 の二つの分布を逆向きにとって,𝑓をガウス分布の密度,𝑞を同じ分散をもつ一様分布の密度にすると,区間の外では𝑞 =0なのに𝑓 >0だから∫𝑆𝑓 <1であり,𝐷(𝑓 ‖ 𝑞)は定まらない.
定理 7.5.2(情報不等式(連続版)). 𝑛 ≥1とし,𝑋をℝ𝑛に値をとり密度𝑓をもつ確率変数,𝑞をℝ𝑛上の密度とし,𝑆を定義 7.5.1 のとおりとする.∫𝑆𝑓(𝑥) 𝑑𝑥 =1であり,ℎ(𝑋)が定まり,∫𝑆𝑓(𝑥)|log𝑞(𝑥)| 𝑑𝑥 <∞であるならば,𝐷(𝑓 ‖ 𝑞)は定まり
𝐷(𝑓‖𝑞)=−ℎ(𝑋)−∫𝑆𝑓(𝑥)log𝑞(𝑥)𝑑𝑥≥0である.
証明. まず𝐷(𝑓 ‖ 𝑞)が定まることを見る.𝑆の上でlog𝑓𝑞 =log𝑓 −log𝑞であり,∫𝑆𝑓|log𝑓| ≤∫ℝ𝑛|𝑓log𝑓|はℎ(𝑋)が定まることから有限,∫𝑆𝑓|log𝑞|は仮定より有限だから,∫𝑆𝑓|log𝑓𝑞|は有限である.同時に,𝐷(𝑓 ‖ 𝑞) =∫𝑆𝑓log𝑓 −∫𝑆𝑓log𝑞と分けてよい.
第 1 項を−ℎ(𝑋)に直す.本定理の仮定は補題 7.2.3 の仮定そのものだから,補題 7.2.3 の証明と同じく,7.2 節で借りた零集合の上での積分により∫𝑆𝑓log𝑓 =∫ℝ𝑛𝑓log𝑓 = −ℎ(𝑋)である.よって等式が従う.
非負性も補題 7.2.3 から出る.同じ仮定のもとでℎ(𝑋) ≤ −∫𝑆𝑓log𝑞が成り立ち,移項すれば右辺が0以上であることがわかる.◼
定理 7.5.2 は,離散の情報不等式(定理 1.6.1)を密度で書き直したものである.離散のときは−∑𝑥𝑝(𝑥)log𝑞(𝑥)と𝐻(𝑝)をそれぞれ平均符号長として読めたが,連続では−∫𝑆𝑓log𝑞もℎ(𝑋)も負になりうる(例 7.1.3)ので,その読み方は個々の項には当てられない.当てられるのは差のほうで,𝐷(𝑓 ‖ 𝑞)が「密度が𝑞だと思い込んだために余計に払うぶん」にあたる.思い込みで得はしない,という主張になっている.離散のときに述べた等号条件(𝑝 =𝑞のとき,かつそのときに限る)にあたるものを,本書は密度について述べない.定理 7.2.4 の等号のときと同じ理由で,密度どうしの一致は確率0の集合の上での食い違いを許す形でしか言えないからである.
命題 7.5.3(目盛の取り替えでの不変性). 𝑛 ≥1とし,𝑓と𝑞をℝ𝑛上の密度,𝑐1,…,𝑐𝑛を0でない実数,𝑏1,…,𝑏𝑛を実数とし
𝑓∗(𝑥):=1|𝑐1⋯𝑐𝑛|𝑓(𝑥1−𝑏1𝑐1,…,𝑥𝑛−𝑏𝑛𝑐𝑛),𝑞∗(𝑥):=1|𝑐1⋯𝑐𝑛|𝑞(𝑥1−𝑏1𝑐1,…,𝑥𝑛−𝑏𝑛𝑐𝑛)とおく.𝐷(𝑓 ‖ 𝑞)が定まるならば𝐷(𝑓∗ ‖ 𝑞∗)も定まり,𝐷(𝑓∗ ‖ 𝑞∗) =𝐷(𝑓 ‖ 𝑞)である.
証明. 以下,𝑥 ∈ℝ𝑛に対して𝑦 :=(𝑥1−𝑏1𝑐1,…,𝑥𝑛−𝑏𝑛𝑐𝑛)と書き,7.3 節で借りた置換積分をこの置き換えで使う.
𝑓∗と𝑞∗がℝ𝑛上の密度であることを見る.どちらも非負であり,置換積分より∫ℝ𝑛𝑓∗(𝑥) 𝑑𝑥 =∫ℝ𝑛𝑓(𝑦) 𝑑𝑦 =1となり,𝑞∗についても同じ計算である.𝑓∗(𝑥) >0と𝑓(𝑦) >0は同値,𝑞∗についても同様だから,定義 7.5.1 の集合は𝑆∗ ={𝑥 ∈ℝ𝑛 :𝑦 ∈𝑆}である.また𝑆∗の上で
log𝑓∗(𝑥)𝑞∗(𝑥)=log𝑓(𝑦)𝑞(𝑦)であり,1/|𝑐1⋯𝑐𝑛|は分母と分子で打ち消し合う.置換積分により積分域𝑆∗は𝑆に移り,𝑑𝑥 =|𝑐1⋯𝑐𝑛| 𝑑𝑦だから
∫𝑆∗𝑓∗(𝑥)∣log𝑓∗(𝑥)𝑞∗(𝑥)∣𝑑𝑥=∫𝑆𝑓(𝑦)∣log𝑓(𝑦)𝑞(𝑦)∣𝑑𝑦であり,同じ置き換えで∫𝑆∗𝑓∗ =∫𝑆𝑓 =1である.よって𝐷(𝑓∗ ‖ 𝑞∗)は定まり,絶対値を外した同じ計算でその値が𝐷(𝑓 ‖ 𝑞)に等しいことがわかる.◼
命題 7.5.3 が,本節の量が微分エントロピーと違うところである.𝑛 =1で見ると,ℎは引き伸ばしでlog|𝑐|だけ動いた(命題 7.3.3)のに対し,𝐷は動かない.比べる二つの密度に同じ変換をかければ,比が変わらないからである.目盛を変えても意味が変わらないので,𝐷の値はそれ自体で読んでよい量になっている.
相互情報量
定義 7.5.4(相互情報量). 𝑋と𝑌を実数値確率変数とし,(𝑋,𝑌)が結合密度𝑓をもつとする.定義 7.4.1 を𝑓に当てて得られる𝑋と𝑌の周辺密度を𝑓1,𝑓2,その積を𝑞(𝑥,𝑦) :=𝑓1(𝑥)𝑓2(𝑦)とする(補題 7.4.2 より𝑞はℝ2上の密度である).𝐷(𝑓 ‖ 𝑞)が定まるとき,𝑋と𝑌の 相互情報量 を
𝐼(𝑋;𝑌):=𝐷(𝑓‖𝑞)で定める.
離散のときの定義 1.3.1 は,相互情報量を「同時分布と周辺分布の積との隔たり」として読めるものだった(1.6 節).定義 7.5.4 はその読み方をそのまま定義に据えている.なおℎのときと同じく,𝐷(𝑓 ‖ 𝑞)と𝐼(𝑋;𝑌)の値も密度の選び方によらない.これも測度論を使って示される事実で,本書では証明しない.
独立なときの値を先に見ておく.
命題 7.5.5(独立なときの相互情報量). 𝑋と𝑌を独立な実数値確率変数とし,どちらも密度をもつとする.このとき(𝑋,𝑌)は結合密度をもち,𝐼(𝑋;𝑌)は定まって𝐼(𝑋;𝑌) =0である.
証明. 𝑓1,𝑓2を𝑋,𝑌の密度とする.7.4 節で借りた独立と積の密度により𝑞(𝑥,𝑦) :=𝑓1(𝑥)𝑓2(𝑦)は(𝑋,𝑌)の結合密度であり,𝑓1と𝑓2はその周辺密度でもある.結合密度としてこの𝑞をとれば,定義 7.5.4 の𝑓も𝑞に等しい.どちらも同じ関数だから𝑆 ={𝑞 >0}で∫𝑆𝑞 =1であり,𝑆の上でlog𝑞𝑞 =0だから,𝐷(𝑞 ‖ 𝑞)は定まって0である.◼
定理 7.5.6(相互情報量とエントロピーの関係). 𝑋と𝑌を実数値確率変数とし,(𝑋,𝑌)が結合密度をもつとする.ℎ(𝑋,𝑌)が定まり,ℎ(𝑋)とℎ(𝑌)が定まるならば,𝐼(𝑋;𝑌)は定まり
𝐼(𝑋;𝑌)=ℎ(𝑋)+ℎ(𝑌)−ℎ(𝑋,𝑌)≥0である.
証明. 𝑓を結合密度とし,𝑓1,𝑓2,𝑞,𝑆を補題 7.4.2 の𝑛 =2の場合のとおりにとる.補題 7.4.2 より𝑞はℝ2上の密度で,∫𝑆𝑓 =1と∫𝑆𝑓|log𝑞| <∞が成り立つ.よって定理 7.5.2 の仮定が満たされ,𝐷(𝑓 ‖ 𝑞)すなわち𝐼(𝑋;𝑌)は定まって
𝐼(𝑋;𝑌)=−ℎ(𝑋,𝑌)−∫𝑆𝑓log𝑞≥0である.補題 7.4.2 より−∫𝑆𝑓log𝑞 =ℎ(𝑋) +ℎ(𝑌)だから,主張を得る.◼
定理 7.5.6 は,離散のときのエントロピー表現𝐼(𝑋;𝑌) =𝐻(𝑋) +𝐻(𝑌) −𝐻(𝑋,𝑌)(定理 1.3.4)と同じ形をしている.離散のエントロピーと微分エントロピーは別物だったのに,それを組み合わせた相互情報量は同じ式で書ける.右辺が三つの差でできていて,微分エントロピーの数値そのものが表に出ないからである.7.1 節の終わりに予告した「差だけが現れる」というのは,このことである.
定理 7.4.3 との関係も読める.𝑛 =2の場合の劣加法性ℎ(𝑋,𝑌) ≤ℎ(𝑋) +ℎ(𝑌)は,定理 7.5.6 の𝐼(𝑋;𝑌) ≥0を移項したものにほかならない.離散のときに 1.6 節で情報不等式が既出の不等式を束ねたのと同じことが,ここでも起きている.例 7.4.5 の三角形の上の一様分布では𝐼(𝑋;𝑌) =log𝑒2で,底を2にとれば約0.4427ビットである.𝑋と𝑌の依存の強さが,それだけあるということである.
7.3 節が残した問いにも,ここで答えが付く.微分エントロピーの数値は目盛の取り替えで動いたが,相互情報量は動かない.
系 7.5.7(目盛の取り替えでの相互情報量の不変性). 𝑋と𝑌を実数値確率変数とし,(𝑋,𝑌)が結合密度をもち𝐼(𝑋;𝑌)が定まるとする.𝑐と𝑐′を0でない実数,𝑏と𝑏′を実数とすると,𝐼(𝑐𝑋 +𝑏; 𝑐′𝑌 +𝑏′)は定まって𝐼(𝑋;𝑌)に等しい.
証明. 𝑓を(𝑋,𝑌)の結合密度,𝑓1,𝑓2を定義 7.4.1 の周辺密度,𝑞(𝑥,𝑦) :=𝑓1(𝑥)𝑓2(𝑦)とする.補題 7.4.2 より𝑞はℝ2上の密度であり,定義 7.5.4 より𝐷(𝑓 ‖ 𝑞) =𝐼(𝑋;𝑌)は定まっている.命題 7.5.3 を𝑛 =2,(𝑐1,𝑐2) :=(𝑐,𝑐′),(𝑏1,𝑏2) :=(𝑏,𝑏′)ととって読み,そこで作られる密度を𝑓∗,𝑞∗とする.
𝑓∗が(𝑐𝑋 +𝑏, 𝑐′𝑌 +𝑏′)の結合密度であることを見る.任意の直方体𝑅について,(𝑐𝑋 +𝑏, 𝑐′𝑌 +𝑏′)が𝑅に入ることと,(𝑋,𝑌)が𝑅の二つの辺を𝑡 ↦(𝑡 −𝑏)/𝑐と𝑡 ↦(𝑡 −𝑏′)/𝑐′で送った直方体に入ることとは同値である.定義 7.1.1 をその直方体に当て,7.3 節で借りた置換積分を𝑦1 =(𝑥1 −𝑏)/𝑐,𝑦2 =(𝑥2 −𝑏′)/𝑐′で使うと,Pr[(𝑐𝑋 +𝑏, 𝑐′𝑌 +𝑏′) ∈𝑅] =∫𝑅𝑓∗(𝑥) 𝑑𝑥を得る.𝑓∗は非負で,同じ置き換えにより∫ℝ2𝑓∗ =∫ℝ2𝑓 =1だからℝ2上の密度であり,したがって定義 7.1.1 の意味の結合密度である.
𝑓∗の周辺密度を求める.定義 7.4.1 の積分で,第1周辺密度については第2変数について,第2周辺密度については第1変数について,借りた置換積分を使うと,前者は𝑥 ↦1|𝑐|𝑓1(𝑥−𝑏𝑐),後者は𝑥 ↦1|𝑐′|𝑓2(𝑥−𝑏′𝑐′)である.この二つの積は𝑞∗にほかならない.
よって定義 7.5.4 より𝐼(𝑐𝑋 +𝑏; 𝑐′𝑌 +𝑏′) =𝐷(𝑓∗ ‖ 𝑞∗)であり,命題 7.5.3 よりこれは定まって𝐷(𝑓 ‖ 𝑞) =𝐼(𝑋;𝑌)に等しい.◼
𝑋と𝑌を別々の目盛で測り直しても,二つのあいだの依存の強さは変わらないということである.
例 7.5.8(二つのガウス分布の相対エントロピー). 𝜇1,𝜇2を実数,𝜎21,𝜎22 >0とし,𝑓をN(𝜇1,𝜎21)の密度,𝑞をN(𝜇2,𝜎22)の密度とする(定義 7.2.1).このとき𝐷(𝑓 ‖ 𝑞)は定まり
𝐷(𝑓‖𝑞)=12log𝜎22𝜎21+log𝑒2(𝜎21+(𝜇1−𝜇2)2𝜎22−1)である.
証明. 𝑋がN(𝜇1,𝜎21)に従うとする.どちらの密度も実軸のすべての点で正だから定義 7.5.1 の𝑆はℝ全体で,∫𝑆𝑓 =1である.定義 7.2.1 の二つの密度について対数の差をとると
log𝑓(𝑥)𝑞(𝑥)=12log𝜎22𝜎21+log𝑒((𝑥−𝜇2)22𝜎22−(𝑥−𝜇1)22𝜎21)である.右辺は𝑥の2次式で,借りたガウス分布の基本性質より𝔼[𝑋] =𝜇1と𝔼[(𝑋 −𝜇1)2] =𝜎21はどちらも有限だから,𝑓を掛けた積分は絶対収束し𝐷(𝑓 ‖ 𝑞)は定まる.
平均をとる.𝔼[(𝑋 −𝜇1)2] =𝜎21であり,(𝑋 −𝜇2)2 =(𝑋 −𝜇1)2 +2(𝜇1 −𝜇2)(𝑋 −𝜇1) +(𝜇1 −𝜇2)2の平均をとると𝔼[𝑋 −𝜇1] =0より𝔼[(𝑋 −𝜇2)2] =𝜎21 +(𝜇1 −𝜇2)2である.よって
𝐷(𝑓‖𝑞)=12log𝜎22𝜎21+log𝑒(𝜎21+(𝜇1−𝜇2)22𝜎22−𝜎212𝜎21)となり,右辺は主張の式である.◼
両端で確かめる. 𝜎21 =𝜎22 =𝜎2のとき,例 7.5.8 の式はlog𝑒2 ⋅(𝜇1−𝜇2)2𝜎2になる.平均が離れているほど大きく,二つの密度が同じもの(𝜇1 =𝜇2)になったところで0である.定理 7.5.2 の非負性と整合している.
距離ではない. 離散のときに 1.6 節で見たのと同じく,𝐷(𝑓 ‖ 𝑞)と𝐷(𝑞 ‖ 𝑓)は一致しない.N(0,1)とN(0,4)で例 7.5.8 の式を両向きに計算すると,底を2にとって前者から後者への隔たりが約0.4590ビット,後者から前者への隔たりが約1.1640ビットである.狭い分布を広い分布だと思い込む損と,広い分布を狭い分布だと思い込む損は,別物だということである.
ガウス通信路の章は,本章で組み立てた量を通信路に当てる.雑音がガウス分布に従い,送る信号の電力に制約のかかる通信路では,容量が本節の相互情報量の最大値として書け,その最大を決めるのが定理 7.2.4 の最大エントロピー性である.本章が条件付き微分エントロピーを用意しなかったぶんは,そこで補うことになる.
InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.