11.3 Sanov の定理
前節は型を一つ指定したときの確率を挟み込んだ(定理 11.2.1).指数に残るのは相対エントロピーだけで,𝑄から見て型が𝑃の系列が出る確率は𝑒−𝑛𝐷(𝑃‖𝑄)の前後にあった.本節は指定するものを,型一つから分布の集合に替える.X上の分布の集合Eを先に決めておき,長さ𝑛の系列を𝑄から引いたとき,その経験分布がEに落ちる確率を問う.
この確率が小さいことは,前節までの二つから見当がつく.𝐷(𝑃 ‖ 𝑄)が0になるのは𝑃 =𝑄のときに限る(定理 1.6.1)ので,𝑄自身から離れたところにEを置けば,Eに属するどの型𝑃でも𝐷(𝑃 ‖ 𝑄)は正であり,その型類の確率は 定理 11.2.1 より𝑛とともに指数で落ちる.落ちる速さは型ごとに違うが,型の個数は𝑛の多項式で抑えられている(命題 11.1.3).有限個の正の数の和は,最大の項以上で,最大の項に項の個数を掛けたもの以下である.したがって1𝑛logをとって𝑛を大きくすると多項式倍の差は消え,残るのはいちばん大きい項の指数,すなわち𝐷( ⋅ ‖ 𝑄)がいちばん小さい型の指数だけになる.これが本節の主張で,Sanov の定理と呼ばれる.
示し方は上下に分ける.上界は,型の個数の多項式評価と 定理 11.2.1 の上界を合わせた数え上げで出る.下界は,Eの中の分布を一つ選び,その型類ぶんだけで確率を下から押さえて出す.二つを合わせると,指数が一点で決まるという主張になる.
集合に属する型
定義 11.3.1(集合に属する型). Xを空でない有限アルファベットとし,EをX上の分布の集合とする.𝑛 ≥1に対し,Eに属する長さ𝑛の型(定義 11.1.1)の全体を
E𝑛:={𝑃∈E:𝑃 は長さ 𝑛 の型}と書く.
Eは分布の集合であって,系列の集合ではない.系列の側でこれに対応するのは,型がEに入る系列の全体{𝑥 ∈X𝑛 :ˆ𝑃𝑥 ∈E}である.どの系列もただ一つの型をもつ(定義 11.1.1)から,この集合はE𝑛に属する型の型類を重なりなく合わせたものにほかならない.Eが無限に多くの分布を含んでいてもよいが,長さ𝑛の型は有限個しかない(命題 11.1.3)のでE𝑛は有限集合である.Eの中に長さ𝑛の型が一つも無いこともあり,そのときこの系列の集合は空になる.本節が測るのは,X上の分布𝑄とその𝑛重の積分布𝑄𝑛に対する,この集合の確率である.
上界
定理 11.3.2 の𝛾はinf𝑃∈E𝐷(𝑃 ‖ 𝑄)と思って読めばよい.条件をE𝑛に属する型についてだけ課しているのは,証明が型の上の和しか使わないからである.
定理 11.3.2(Sanov の上界). Xを空でない有限アルファベット,𝑄をX上の全点で正の分布,𝑄𝑛をその𝑛重の積分布とし,EをX上の分布の集合,E𝑛を 定義 11.3.1 のとおりとする.𝑛 ≥1とし,実数𝛾が,E𝑛に属するどの型𝑃についても𝛾 ≤𝐷(𝑃 ‖ 𝑄)を満たすとすると
𝑄𝑛({𝑥∈X𝑛:ˆ𝑃𝑥∈E})≤(𝑛+1)|X|𝑒−𝑛𝛾であり,左辺が正ならば
1𝑛log𝑄𝑛({𝑥∈X𝑛:ˆ𝑃𝑥∈E})≤−𝛾+|X|log(𝑛+1)𝑛である(ˆ𝑃𝑥は 定義 11.1.1 の型,𝐷は 1.6 節の相対エントロピー).
証明. どの系列もただ一つの型をもつ(定義 11.1.1)から,{𝑥 ∈X𝑛 :ˆ𝑃𝑥 ∈E}はE𝑛に属する型𝑃の型類T𝑛(𝑃)を重なりなく合わせたものであり
𝑄𝑛({𝑥:ˆ𝑃𝑥∈E})=∑𝑃∈E𝑛𝑄𝑛(T𝑛(𝑃))である.右辺の各項に 定理 11.2.1 の上界を当てると𝑄𝑛(T𝑛(𝑃)) ≤𝑒−𝑛𝐷(𝑃‖𝑄)であり,仮定より𝛾 ≤𝐷(𝑃 ‖ 𝑄)であり,指数関数は単調増加(logの単調性とlog𝑒𝑠 =𝑠から出る)だから𝑄𝑛(T𝑛(𝑃)) ≤𝑒−𝑛𝛾である.項の個数はEに属する長さ𝑛の型の個数だから,長さ𝑛の型の総数以下であり,命題 11.1.3 より(𝑛 +1)|X|以下である.よって和は(𝑛 +1)|X|𝑒−𝑛𝛾以下である.
第2の不等式は第1の両辺の対数をとって𝑛で割ったものである.左辺が正なら対数がとれ,logは単調だから
log𝑄𝑛({𝑥:ˆ𝑃𝑥∈E})≤|X|log(𝑛+1)−𝑛𝛾であり,両辺を𝑛で割れば主張を得る.◼
型の個数の多項式評価は,ここで初めて本来の役目を果たす.型ごとの上界𝑒−𝑛𝛾を項の個数だけ足し合わせても,掛かるのは𝑛の多項式だけなので,1𝑛logをとると0に向かう項しか足されない.そのことを極限の形にしたのが次の系である.
系 11.3.3. Xを空でない有限アルファベット,𝑄をX上の全点で正の分布,𝑄𝑛をその𝑛重の積分布とし,EをX上の分布の集合,E𝑛を 定義 11.3.1 のとおりとする.実数𝛾が,すべての𝑛 ≥1とE𝑛に属するどの型𝑃についても𝛾 ≤𝐷(𝑃 ‖ 𝑄)を満たすとする.𝑄𝑛({𝑥 ∈X𝑛 :ˆ𝑃𝑥 ∈E})が正であるような𝑛に限って考えることにすると,どの𝜀 >0についても,ある𝑛0があって,𝑛 ≥𝑛0を満たすそのような𝑛のすべてについて
1𝑛log𝑄𝑛({𝑥∈X𝑛:ˆ𝑃𝑥∈E})≤−𝛾+𝜀が成り立つ.
証明. 定理 11.3.2 の第2の不等式より,左辺は−𝛾 +|X|log(𝑛+1)𝑛以下である.したがって|X|log(𝑛+1)𝑛が0に収束することを見れば,与えられた𝜀に対してその収束から𝑛0がとれる.
補題 11.2.2 よりlog(𝑛+1)𝑛 →0であり,|X|は𝑛に依らない有限の数だから,|X|log(𝑛+1)𝑛 →0である.◼
𝛾としてどれだけ大きい値がとれるかで,系 11.3.3 の強さが決まる.とれるのは,すべての𝑛のE𝑛を通して𝐷( ⋅ ‖ 𝑄)を下から抑える値までである.Eの上で𝐷( ⋅ ‖ 𝑄)を最小にする分布があれば,その値は条件を満たす𝛾の一つになる.同じ値が下からの評価でも現れることを,次に見る.
下界
下界は,Eの中の分布を一つ選び,その型類ぶんだけを数えて出す.ただし選んだ分布がそのまま長さ𝑛の型であるとは限らない.長さ𝑛の型がとる値は0,1/𝑛,…,1に限られる(定義 11.1.1)からである.そこで,選んだ分布に近い型を作って代用する.各文字について𝑛倍した値を整数に切り下げ,切り下げで足りなくなったぶんを一つの文字にまとめて押し付ければ,個数の総和が𝑛になって型になる.以下,⌊𝑡⌋で実数𝑡以下の最大の整数を表す(⌊𝑡⌋ ≤𝑡 <⌊𝑡⌋ +1である).選ぶ分布は˜𝑃と書く.ここで選ぶ分布に最適性は要らないので星印を付けない.
定義 11.3.4(丸め型). Xを空でない有限アルファベット,˜𝑃をX上の分布,𝑎0をXの文字とし,𝑛 ≥1とする.
˜𝑃𝑛(𝑎):=⌊𝑛˜𝑃(𝑎)⌋𝑛(𝑎≠𝑎0),˜𝑃𝑛(𝑎0):=1−∑𝑎≠𝑎0⌊𝑛˜𝑃(𝑎)⌋𝑛で定まるX上の関数˜𝑃𝑛を,𝑎0を端数の引き受け手とする˜𝑃の 丸め型 と呼ぶ(⌊𝑡⌋は𝑡以下の最大の整数).
補題 11.3.5(丸め型は長さ𝑛の型であり,各文字で元の分布に収束する). Xを空でない有限アルファベット,˜𝑃をX上の分布,𝑎0をXの文字とし,各𝑛 ≥1について˜𝑃𝑛を𝑎0を端数の引き受け手とする˜𝑃の丸め型(定義 11.3.4)とする.このとき,どの𝑛 ≥1でも˜𝑃𝑛は長さ𝑛の型(定義 11.1.1)であり,どの文字𝑎についても˜𝑃𝑛(𝑎) →˜𝑃(𝑎)(𝑛 →∞)である.
証明. まず˜𝑃𝑛が長さ𝑛の型であることを見る.𝑎 ≠𝑎0については𝑛˜𝑃𝑛(𝑎) =⌊𝑛˜𝑃(𝑎)⌋で,これは非負整数である.𝑎0については𝑛˜𝑃𝑛(𝑎0) =𝑛 −∑𝑎≠𝑎0⌊𝑛˜𝑃(𝑎)⌋で,これも整数であり,⌊𝑡⌋ ≤𝑡と∑𝑎˜𝑃(𝑎) =1から
∑𝑎≠𝑎0⌊𝑛˜𝑃(𝑎)⌋≤∑𝑎≠𝑎0𝑛˜𝑃(𝑎)≤𝑛なので非負である.よって(𝑛˜𝑃𝑛(𝑎))𝑎∈Xは総和が𝑛の非負整数の組であり,各文字をその個数だけ並べた系列の型は˜𝑃𝑛だから,˜𝑃𝑛は長さ𝑛の型である.
次に˜𝑃𝑛が各文字で˜𝑃に収束することを見る.𝑎 ≠𝑎0については⌊𝑡⌋ ≤𝑡 <⌊𝑡⌋ +1より0 ≤˜𝑃(𝑎) −˜𝑃𝑛(𝑎) <1/𝑛である.𝑎0については,˜𝑃𝑛と˜𝑃の総和がどちらも1だから˜𝑃𝑛(𝑎0) −˜𝑃(𝑎0) =∑𝑎≠𝑎0(˜𝑃(𝑎) −˜𝑃𝑛(𝑎))であり,右辺の各項は0以上1/𝑛未満だから,差の絶対値は|X|/𝑛以下である.よってどの文字でも˜𝑃𝑛(𝑎) →˜𝑃(𝑎)である.◻
定理 11.3.6(Sanov の下界). Xを空でない有限アルファベット,𝑄をX上の全点で正の分布,𝑄𝑛をその𝑛重の積分布とし,EをX上の分布の集合,E𝑛を 定義 11.3.1 のとおりとする.˜𝑃をX上の分布,𝑎0をXの文字とし,各𝑛 ≥1について˜𝑃𝑛を𝑎0を端数の引き受け手とする˜𝑃の丸め型(定義 11.3.4)とする.十分大きいすべての𝑛について˜𝑃𝑛 ∈E𝑛であるならば,十分大きい𝑛で𝑄𝑛({𝑥 ∈X𝑛 :ˆ𝑃𝑥 ∈E})は正であり
lim inf𝑛→∞1𝑛log𝑄𝑛({𝑥∈X𝑛:ˆ𝑃𝑥∈E})≥−𝐷(˜𝑃‖𝑄)である(ˆ𝑃𝑥は 定義 11.1.1 の型,𝐷は 1.6 節の相対エントロピー).
証明. 補題 11.3.5 より,どの𝑛 ≥1でも˜𝑃𝑛は長さ𝑛の型であり,どの文字でも˜𝑃𝑛(𝑎) →˜𝑃(𝑎)である.
仮定より,ある𝑛0があって𝑛 ≥𝑛0のとき˜𝑃𝑛 ∈E𝑛,とくに˜𝑃𝑛 ∈Eである.型が˜𝑃𝑛である系列は型がEに属するから
T𝑛(˜𝑃𝑛)⊆{𝑥∈X𝑛:ˆ𝑃𝑥∈E}であり,包含している側の集合の確率のほうが小さくないから𝑄𝑛(T𝑛(˜𝑃𝑛)) ≤𝑄𝑛({𝑥 :ˆ𝑃𝑥 ∈E})である.左辺は 定理 11.2.1 の下界より(𝑛 +1)−|X|𝑒−𝑛𝐷(˜𝑃𝑛‖𝑄)以上で,これは正だから,𝑛 ≥𝑛0では右辺も正である.
logは単調だから,𝑛 ≥𝑛0について
1𝑛log𝑄𝑛({𝑥:ˆ𝑃𝑥∈E})≥1𝑛log𝑄𝑛(T𝑛(˜𝑃𝑛))である.˜𝑃𝑛は長さ𝑛の型で各文字で˜𝑃に収束するから,系 11.2.3 より右辺は−𝐷(˜𝑃 ‖ 𝑄)に収束する.よって左辺の下極限は−𝐷(˜𝑃 ‖ 𝑄)以上である.◼
二つを合わせる
定理 11.3.7(Sanov の定理). Xを空でない有限アルファベット,𝑄をX上の全点で正の分布,𝑄𝑛をその𝑛重の積分布とし,EをX上の分布の集合,E𝑛を 定義 11.3.1 のとおりとする.𝑃⋆をX上の分布,𝑎0をXの文字とし,𝑃⋆𝑛を𝑎0を端数の引き受け手とする𝑃⋆の丸め型(定義 11.3.4)とする.次の二つが成り立つとする.
- すべての𝑛 ≥1とE𝑛に属するどの型𝑃についても𝐷(𝑃⋆ ‖ 𝑄) ≤𝐷(𝑃 ‖ 𝑄)である.
- 十分大きいすべての𝑛について𝑃⋆𝑛 ∈E𝑛である.
このとき,十分大きい𝑛で𝑄𝑛({𝑥 ∈X𝑛 :ˆ𝑃𝑥 ∈E})は正であり
1𝑛log𝑄𝑛({𝑥∈X𝑛:ˆ𝑃𝑥∈E})⟶−𝐷(𝑃⋆‖𝑄)(𝑛→∞)である(ˆ𝑃𝑥は 定義 11.1.1 の型,𝐷は 1.6 節の相対エントロピー).
証明. 第2の仮定から 定理 11.3.6 が使えて,十分大きい𝑛でこの確率は正であり,その下極限は−𝐷(𝑃⋆ ‖ 𝑄)以上である.
上からの評価には 系 11.3.3 を𝛾 :=𝐷(𝑃⋆ ‖ 𝑄)ととって当てる.第1の仮定がその𝛾についての条件そのものだから,どの𝜀 >0についても,ある𝑛0があって,𝑛 ≥𝑛0でこの確率が正であるかぎり
1𝑛log𝑄𝑛({𝑥:ˆ𝑃𝑥∈E})≤−𝐷(𝑃⋆‖𝑄)+𝜀である.十分大きい𝑛では確率は正なのだから,上極限は−𝐷(𝑃⋆ ‖ 𝑄) +𝜀以下であり,𝜀 >0は任意だから,上極限は−𝐷(𝑃⋆ ‖ 𝑄)以下である.
下極限が−𝐷(𝑃⋆ ‖ 𝑄)以上で,上極限が−𝐷(𝑃⋆ ‖ 𝑄)以下だから,この数列は収束して極限は−𝐷(𝑃⋆ ‖ 𝑄)である.◼
読み方は素直である.経験分布がEに落ちる確率の指数は,Eに属する型のうち𝑄にいちばん近いもの一つで決まり,残りの型は指数の水準では何も寄与しない.第1の仮定は𝑃⋆がその一つであること,すなわちすべての𝑛のE𝑛を通した最小化子であることを求めている.第2の仮定は,Eが型で近づけられる形をしていることを求めている.たとえばEが一点だけからなり,その一点がどの𝑛でも長さ𝑛の型でなければ,どの𝑛でもE𝑛は空になり,第2の仮定は成り立たない.このとき系列の集合も空で,確率は0である.
第2の仮定を確かめる手筋は決まっている.Eが凸で,境界から離れた内側の分布が一つとれるなら,最小化子をその分布のほうへ少しだけ寄せればよい.寄せた先も内側にあるから,その丸め型は十分大きい𝑛でEに入り,第2の仮定が成り立つ.寄せたぶん𝐷( ⋅ ‖ 𝑄)の値は上がるが,寄せ幅を0に近づければその損は消える.この手筋は 定理 11.6.9 の証明で実際に使う.
数値で見る
例 11.3.8(コインの表が7割以上出る確率). X ={0,1}とし,1を表と読む.𝑄を𝑄(0) =𝑄(1) =1/2で定まる分布,𝑄𝑛をその𝑛重の積分布とし,E :={ 𝑃 :𝑃 は X 上の分布で 𝑃(1) ≥0.7 }とする.ˆ𝑃𝑥を 定義 11.1.1 の型,𝐻𝑏を 例 1.1.2 の二値エントロピー関数とし,𝑃⋆を𝑃⋆(0) =0.3,𝑃⋆(1) =0.7で定まる分布とすると,次の四つが成り立つ.
- Eに属するどの分布𝑃についても𝐷(𝑃⋆ ‖ 𝑄) ≤𝐷(𝑃 ‖ 𝑄)である(𝐷は 1.6 節の相対エントロピー).
- 𝐷(𝑃⋆ ‖ 𝑄) =log2 −𝐻𝑏(0.7)であり,その値は約0.0823ナットである.
- 1𝑛log𝑄𝑛({𝑥 ∈X𝑛 :ˆ𝑃𝑥 ∈E}) ⟶ −(log2 −𝐻𝑏(0.7))(𝑛 →∞)である.
- 端数の引き受け手を𝑎0 :=0ととった𝑃⋆の丸め型(定義 11.3.4)は,𝑛が10の倍数でないかぎりEに属さない.とくに,その選び方では 定理 11.3.7 の第2の仮定は成り立たない.
証明.
-
補題 11.2.5 より,X上のどの分布𝑃についても𝐷(𝑃 ‖ 𝑄) =log2 −𝐻𝑏(𝑃(0))である.𝑃 ∈Eとすると𝑃(1) ≥0.7だから𝑃(0) ≤0.3であり,0 ≤𝑃(0) ≤0.3 ≤1/2だから 補題 9.3.1 の第3の主張より𝐻𝑏(𝑃(0)) ≤𝐻𝑏(0.3)である.𝑃⋆(0) =0.3だから,log2から引く向きに直すと𝐷(𝑃⋆ ‖ 𝑄) ≤𝐷(𝑃 ‖ 𝑄)を得る.
-
補題 11.2.5 を𝑃⋆に当てると𝐷(𝑃⋆ ‖ 𝑄) =log2 −𝐻𝑏(𝑃⋆(1)) =log2 −𝐻𝑏(0.7)である.𝐻𝑏(0.7) = −0.7log0.7 −0.3log0.3 =0.61086…,log2 =0.69314…だから,差は0.08228…である.
-
定理 11.3.7 の仮定を,端数の引き受け手を𝑎0 :=1ととって確かめる.第1の仮定は,E𝑛 ⊆Eだから第1の主張から従う.第2の仮定を見る.丸め型は𝑃⋆𝑛(0) =⌊0.3𝑛⌋/𝑛,𝑃⋆𝑛(1) =1 −⌊0.3𝑛⌋/𝑛である.𝑛𝑃⋆𝑛(0) =⌊0.3𝑛⌋と𝑛𝑃⋆𝑛(1) =𝑛 −⌊0.3𝑛⌋はどちらも非負整数で和は𝑛だから,0を⌊0.3𝑛⌋個並べたあと1を並べた系列の型は𝑃⋆𝑛であり,𝑃⋆𝑛は長さ𝑛の型である.また⌊0.3𝑛⌋ ≤0.3𝑛より𝑃⋆𝑛(0) ≤0.3だから𝑃⋆𝑛(1) ≥0.7であり,𝑃⋆𝑛 ∈Eである.よってどの𝑛 ≥1でも𝑃⋆𝑛 ∈E𝑛で,第2の仮定も成り立つ.𝑄は全点で正だから 定理 11.3.7 が使えて,第2の主張と合わせて結論を得る.
-
端数の引き受け手を𝑎0 :=0にとると丸め型は𝑃⋆𝑛(1) =⌊0.7𝑛⌋/𝑛である.0.7𝑛 =7𝑛/10は,7と10が互いに素だから𝑛が10の倍数のときに限り整数であり,整数でなければ⌊0.7𝑛⌋ <0.7𝑛だから𝑃⋆𝑛(1) <0.7で,𝑃⋆𝑛 ∉Eである.10の倍数でない𝑛はいくらでも大きくとれるから,「十分大きいすべての𝑛で𝑃⋆𝑛 ∈E𝑛」は成り立たない.
◼
第3と第4の主張は,端数の引き受け手の選び方が仮定の成否を左右することを示している.制約𝑃(1) ≥0.7が緩む側の文字に端数を押し付ければ丸め型はEに入り,きつくなる側に押し付ければ入らない.𝑛 =11が後者の例で,⌊7.7⌋/11 =7/11は0.7より小さい.
指数の値0.0823が言っているのは,1𝑛logをとった量が−0.0823に近づくということであって,有限の𝑛での確率そのものを与えるものではない.公平なコインを投げて表が7割以上出るのはめったに起きないが,その「めったに」の速さは,表の割合を0.7に固定した分布が公平なコインからどれだけ隔たっているかだけで決まる.
本節は,経験分布が指定した集合に落ちる確率を測った.測ったのは𝑄から引いた系列についてであり,集合Eは𝑄とは関わりなく先に決めておいた.次節はEにあたるものを,二つの分布のどちらが真かを当てるという目的から決める.そこでも指数に現れるのは相対エントロピーで,型の方法がそのまま効く.
InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.