10.1 最大エントロピー問題

第1章はエントロピーを分布の不確かさの尺度として立て,定理 1.1.5 でその最大値を求めた.すなわち,アルファベットの大きさしか分かっていないとき,いちばん不確かなのは一様分布である.だが分布について分かっていることが,それだけとは限らない.サイコロの目の平均が測定で分かっている,粒子の平均エネルギーが分かっている,文書の中で母音が占める割合だけが分かっている.こうした部分的な知識をもって,では分布そのものを何と定めればよいか.

候補は無数にある.目の平均がである分布は 1 つではないからである.そこで選び方の原則を一つ置く.知っていることは反映し,知らないことは決めつけない.測定が言っているのは平均がだということだけで,それ以上は何も言っていない.にもかかわらず「の目は出ない」と決めた分布を選べば,測定が支持していない主張を勝手に足したことになる.逆に,条件を満たす分布のうちもっとも不確かなものを選べば,条件が強いること以外は何も決めつけていない.この選び方を定式化したのが,本章の最大エントロピー問題である.

本章ではの底を自然対数にとる.以下で主役になる分布はで書かれた形をしており,が互いに逆でなければ本章の恒等式が閉じないからである.単位はナットで,ビットで読みたいときはで割る.アルファベットは空でない有限集合とし,その上の分布,すなわち各点に非負の値を与えて総和が 1 になる関数を,一般には大文字で書き,あとで主役になる特定の分布だけを小文字に星印を付けて書く.

形式化上の注記(本章共通). 形式化もを自然対数にとるので,本章は本文と形式化で単位が一致する.

命題 10.1.2 の形式化は確率変数の像測度に対して述べられており,10.2 節以降が紐付ける宣言はアルファベット上の関数としての分布に対して述べられている.相対エントロピーも,前者は測度に対する形,後者は pmf に対する形で,単独の宣言としては別のものである.

制約と実行可能集合

定義 10.1.1(モーメント制約と実行可能集合). 空でない有限アルファベット上の実数値関数𝑘 0)を特徴関数(feature function),実数をその制約の値と呼ぶ.上の分布モーメント制約 を満たすとは,

𝔼𝑃[𝑓𝑖]:=𝑥X𝑃(𝑥)𝑓𝑖(𝑥)=𝑐𝑖(𝑖=1,,𝑘)

がすべてので成り立つことをいう.モーメント制約を満たす分布を 実行可能 と呼び,実行可能な分布全体の集合を 実行可能集合 という.実行可能集合の上でを最大にする分布を求める問題を,最大エントロピー問題 という.

制約はについて線形である.は値の並びの一次結合だから,番目の制約は分布全体のなす集合を超平面で切る(が恒等的に 0 でという退化した場合だけは,切り口が空になる).実行可能集合はその枚の切り口の共通部分であり,𝑘 =0(制約を何も課さない場合)なら分布全体そのものである.何を知っているかは特徴関数の選び方に入る.なら平均を知っていることになり,を足せば 2 次モーメントまで,ある部分集合の指示関数をとればその部分集合の確率を知っていることになる.

形式化上の注記. 定義 10.1.1 の実行可能集合を名指す宣言はない.モーメント制約は,形式化では定理 10.2.5定理 10.2.6 の仮定として直接現れる.

一様分布からの隔たりとして読む

命題 10.1.2. を空でない有限アルファベット,とし,上の一様分布とする.上の任意の分布に対して

𝐷(𝑃𝑈)=log𝑀𝐻(𝑃).

証明. はすべてので正だから,1.6 節の相対エントロピーの和はどの項も有限で,

𝐷(𝑃𝑈)=𝑥𝑃(𝑥)log𝑃(𝑥)1/𝑀=𝑥𝑃(𝑥)log𝑃(𝑥)+log𝑀𝑥𝑃(𝑥)

と 2 つの和に分けられる(の項はどちらの形でも 0 なので,分け方に注意は要らない).第 1 の和は定義 1.1.1 より,第 2 の和はよりである.

形式化: klDiv_uniformOn_univ_toReal_eq (ソース)

この恒等式が本章の全体を貫く読み替えを与える.に依らない定数だから,を大きくすることとを小さくすることは,同じ一つのことである.エントロピーが最大の分布とは,一様分布からの隔たりがもっとも小さい分布にほかならない.1.6 節が情報不等式から定理 1.1.5 を再現してみせたのも,この恒等式にを当てただけだった.等式そのものはそこで一度使っている.本章ではこれが出発点になる.

系 10.1.3. を空でない有限アルファベット,をその上の一様分布とする.特徴関数と制約の値の定める実行可能集合をと書く.に属する分布について,次の二つは同値である.

  1. 任意のに対してである.
  2. 任意のに対してである.

証明. とおく.の 2 元をとり,命題 10.1.2 をそれぞれに当てるとだから

𝐷(𝑄𝑈)𝐷(𝑃𝑈)=𝐻(𝑃)𝐻(𝑄)

が成り立つ.左辺が 0 以上であることと右辺が 0 以上であることは同じであり,これはが同値だということである.全体で動かせば,1 と 2 の同値を得る.

形式化上の注記. 系 10.1.3 に対応する単独の宣言はない.命題 10.1.2 の形式化を 2 つの分布で使い,差をとった形になる.

基準を取り替える. この読み替えの値打ちは,基準にとる分布を取り替えられるところにある.一様分布は,制約が何もないときの「何も決めつけていない分布」だった.だが制約があるとき一様分布が実行可能とは限らない.例 10.1.4 がまさにそうである.一様分布が実行可能でなければ,実行可能集合の外から隔たりを測っていることになる.もし基準を実行可能集合の中の分布に取り替えてもこの読み替えが生き延びるなら,隔たりが 0 になる点(基準そのもの)がそのまま最大化子になり,問題は「基準をどう選ぶか」に移る.どんな分布ならそれが本当に言えるのかが 10.2 節の主題であり,命題 10.1.2 はその特別な場合,制約が空で基準が一様分布である場合にあたる.

解きたい問題

例 10.1.4(平均の分かったサイコロ). X ={1,2,3,4,5,6}𝑘 =1𝑓1(𝑥) =𝑥とする.すなわち「目の平均がである」という制約を 1 つだけ課す.このとき実行可能集合は空ではなく,一様分布を含まない.さらに実行可能な分布はすべてを満たす.

証明.と目に確率ずつを置いた分布は平均がだから,実行可能集合は空でない.いっぽう一様分布の平均はと異なるので,一様分布は実行可能でない.最後にを実行可能な分布とすると,定理 1.1.5 よりであり,その等号が成り立つのはが一様分布のときに限る.は一様分布ではないので,不等号は狭義である.

形式化上の注記. 例 10.1.4 に対応する宣言は形式化されていない.本節の計算がその保証のすべてである.

制約が実際に効いていることが,この二つの事実に出ている.どの実行可能な分布もナット(ビット)には届かない.では実際にどこまで行けるのか.つまり,上限が真に下がっているのか,下がるとしてどこまでか.実行可能集合は無限個の分布を含むので,しらみつぶしにはできない.10.2 節がこの問いに閉じた答えを与える.

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