11.1 型と型類

公平なコインを回投げて,出たの割合がちょうどになる確率は,が増えるとどれだけ速く落ちるだろうか.この落ちる速さを指数で書き下すことが本章の主題であり,その指数には相対エントロピーが現れる.手がかりは第2章にある.第2章は長さの系列を,確率という一つの数だけで分類した.定義 2.2.1 の典型集合がそれで,そこに入る確率がに近づくこと(定理 2.2.3)と,要素数がおよそであること(定理 2.2.5)の二つから,圧縮の限界が出た.ところがこの分類は粗く,まったく違う構成の系列どうしが同じ集合に入りうる(例 2.4.2).2.4 節はその粗さを埋めるために,各文字の出現回数を直接見る強典型性へ移り,そこで系列の型(経験分布)を 定義 2.4.1 で置いた.本章はこの型による分類を最後まで押し進める.

型で分類すると,二つの勘定が立つ.一つは,長さの系列を型で類別したときの類の個数であり,これはの多項式で抑えられる.もう一つは,一つの類に入る系列の個数であり,こちらは指数で増えて,その指数は類の型のエントロピーになる.本節はこの二つを示す.類の個数が多項式でしかないことは,あとで有限個の項の和を最大の項一つで置き換えるときに効く.というのも,多項式倍の差はをとると消えるからである.この道具立てを型の方法と呼ぶ.めったに起きない事象の確率がとともにどれだけ速く落ちるか(大偏差),二つの分布のどちらが真かを当てる問題(仮説検定)で誤りがどれだけ速く落ちるか,本章が問うのはこの二つである.

本章ではの底を自然対数にとり,エントロピーと相対エントロピーの単位をナットで測る.型類の要素数や確率がの形で書けるところが本章の眼目なので,指数と対数の底がそろっていると式が読みやすい.第2章が底をにとってと書いたのに対し,単位の名前が変わるだけである.

以下,アルファベット上の長さの系列をと書く.上の分布に対し,各成分が独立にに従う長さのブロックの分布をと書く.すなわちとし,部分集合についてはとする(第6章が通信路の回の使用をと書いたのと同じ,肩に長さを乗せる書き方である).また本章では,最適化して選んだものに星印を付け,分布にはを,値にはを使い分ける.対数が狭義単調増加であること(補題 8.2.5)は,以下ことわりなく使う.

定義

定義 11.1.1(型と型類). を空でない有限アルファベット,とする.系列に含まれる文字の個数を定義 2.4.1)と書き,

ˆ𝑃𝑥(𝑎):=𝑁(𝑎𝑥)𝑛(𝑎X)

で定める.上の分布長さの型 であるとは,となるが存在することをいう.長さの型に対し,その 型類

T𝑛(𝑃):={𝑥X𝑛:ˆ𝑃𝑥=𝑃}

で定める.

は,系列を一度ばらばらにして,どの文字が何回出たかだけを残したものである.順序の情報は捨てられているので,の並べ替えはすべて同じ型を与える.値は非負で,だから総和はであり,は確かに上の分布である.定義 2.4.1 が名前だけ与えていたこの分布に,本章は記号を与えたことになる.長さの型がとる値はに限られる.すなわちを固定すると,分布のうちごく一部だけが型として現れる.

型類は,型がちょうどである長さの系列の全体である.どのもただ一つの型をもつので,長さの型すべての型類はを重なりなく覆う.または「各文字について」と同じことだから,型類は文字ごとの個数を指定して切り出した集合でもある.とくにが長さの型ならはどの文字でも非負整数である.

形式化: 文字の個数 typeCount (ソース),型類 typeClassByCount (ソース)

例 11.1.2(長さの二値系列). X ={0,1}とし,定義 11.1.1 のとおりとする.分布をの形で書くと,長さの型はちょうど

(1,0),(23, 13),(13, 23),(0,1)

個である.系列と並べて書くと,型類はこの順に

T3(1,0)={000},T3(23, 13)={001, 010, 100},T3(13, 23)={011, 101, 110},T3(0,1)={111}

であり,要素数はである.

証明. 長さの系列についてのいずれかであり,だから,型だけで決まる.に対応する型は,の順に上の個である.上の並びは本をの値で分けたものであり,どの型類も空でないから,個はすべて長さの型であって,型はちょうど個である.要素数は並びを数えて順にである.

型が個なのに対し,系列は本ある.一般のアルファベットと一般ので型がいくつあるかを,次に数える.

型は多項式個しかない

命題 11.1.3. を空でない有限アルファベット,とする.長さの型(定義 11.1.1)の個数は以下である.

証明. を長さの型とすると,となるがあり,各文字についてである.したがっては,個数の組で割ったものとして書ける.個数以上以下の整数だから,文字ごとに通りの値しかとれず,組は全部で通り以下である.異なる型は異なる組を与えるので,型の個数もこれ以下である.

形式化: numTypes_le (ソース)

この数え方は粗い.個数の組の総和がでなければならないという条件を使っておらず,型にならない組まで数えている.それでも以降の評価に入るのはの多項式なので,をとる水準では消える.系列のほうはの要素数だけあり,とともに指数で増える.本節の残りは,一つの型類がどれだけ大きいかを数える.

型類の中では確率が一定である

型が分かると,その系列の確率まで決まってしまう.長さのブロックの確率が文字ごとの確率の積であり,同じ文字については何回出たかだけが効くからである.次の命題はその値を書き下したものである.

命題 11.1.4. を空でない有限アルファベット,とし,上の全点で正の分布,をその重の積分布とする.を長さの型(定義 11.1.1),とすると

𝑄𝑛({𝑥})=exp(𝑛(𝐻(𝑃)+𝐷(𝑃𝑄)))

である.ここで定義 1.1.1 のエントロピー,1.6 節の相対エントロピーである.

証明. の積を,同じ文字ごとにまとめる.文字にちょうど回現れ,よりだから

𝑄𝑛({𝑥})=𝑎X𝑄(𝑎)𝑁(𝑎𝑥)=𝑎X𝑄(𝑎)𝑛𝑃(𝑎)

である.は全点で正だからどの因子も正で,対数をとると

log𝑄𝑛({𝑥})=𝑛𝑎X𝑃(𝑎)log𝑄(𝑎)

になる(の文字では因子がであり,右辺の項もである).

右辺の和を自身の対数を経由して書き直す.の文字では

𝑃(𝑎)log𝑄(𝑎)=𝑃(𝑎)log𝑃(𝑎)𝑃(𝑎)log𝑃(𝑎)𝑄(𝑎)

であり,の文字では,の約束(定義 1.1.1)により三つの項がすべてである.について足すと,第項の和は,第項の和はだから

log𝑄𝑛({𝑥})=𝑛(𝐻(𝑃)+𝐷(𝑃𝑄))

となる.両辺の指数をとれば主張を得る.

右辺にが現れないことが,この命題の言っていることのすべてである.型類の中の点はどれも同じ確率をもち,しかもその確率は,型のエントロピーと,型がからどれだけ隔たっているかの二つだけで決まる.第2章の経験エントロピーはという一つの数だったが,命題 11.1.4 の指数はの二つに分かれている.この分かれ方を 定義 2.1.1 の言葉に戻したものが 系 11.1.9 である.

形式化上の注記. 命題 11.1.4 に対応する単独の宣言は無い.typeClassByCount_prod_eq (InformationTheory/Shannon/Sanov/LDP.lean) が同じ等式を与えているが,第因子がの形で残っている.これをに書き換えるのが prod_div_pow_eq_prod_pow_div_npow_of_sum (InformationTheory/Shannon/Sanov/MultinomialLowerBound.lean) と pow_div_prod_pow_eq_exp_n_entropyByCount (InformationTheory/Shannon/TypeClassLowerBound.lean) であり,命題 11.1.4 はこの三つの合成として得られる.

系 11.1.5. を空でない有限アルファベット,とし,を長さの型(定義 11.1.1),をその重の積分布とする.とするとである(定義 1.1.1 のエントロピー).とくにの点はどれも同じ確率をもつ.

証明. を同じ文字ごとにまとめると,よりだからである.の文字ではなので因子はであり,積はの文字だけをわたる.その対数はであり,の約束(定義 1.1.1)のもとでこれはに等しい.両辺の指数をとれば等式を得る.右辺はに依らないから,後半も従う.

形式化上の注記. 系 11.1.5 に対応する単独の宣言も無い.が全点で正であれば,命題 11.1.4 と同じ三つの宣言にとして自身をとった合成で得られる(相対エントロピーの項が消える).

型類の要素数は多項係数である

型類は文字ごとの個数を指定して切り出した集合だった(定義 11.1.1).したがってその要素数は,位置を文字ごとに割り振る割り振り方の数であり,組合せの数として書き下せる.

命題 11.1.6. を空でない有限アルファベット,とし,を長さの型(定義 11.1.1)とすると

T𝑛(𝑃)=𝑛!𝑎X(𝑛𝑃(𝑎))!

である.

証明. は長さの型だから,どの文字でもは非負整数であり,その総和はである.の系列は,位置の全体を,各文字に大きさの集合を割り当てて重なりなく分ける分け方とに対応する(割り当てた位置にその文字を書けばよい).の文字に順番を付け,先頭の文字から順に,まだ使っていない位置の中から必要な個数を選んでいくと,選び方の総数は二項係数の積になり,約分すると主張の右辺になる.

形式化上の注記. 命題 11.1.6 に対応する単独の宣言は無い.形式化は InformationTheory/Shannon/Sanov/MultinomialLowerBound.lean の中でこの等式を両向きの不等式として証明しているが,どちらもファイルの外に公開されていない補助補題なので名前で指すことができない.

右辺の分母は,同じ文字どうしの並べ替えを打ち消す因子である.例 11.1.2 で数えたは,この式のの場合であり,たとえば型ではになる.要素数がこの形で書けると,二つの型類の重さを比べるのは階乗の比を比べることになる.

いちばん重いのは自分の型類

命題 11.1.4 は一点の確率を与えた.型類全体の確率は,これに要素数を掛けたものである.要素数の指数評価に移る前に,で測ったときどの型類がいちばん重いかを見ておく.型類はを重なりなく覆い,その個数は多項式で抑えられる(命題 11.1.3)から,全体の確率のうち多項式の逆数ぶんは,いちばん重い型類が担っている.どれがいちばん重いかをここで決めておけば,系 11.1.5 で確率を要素数に読み替えて 定理 11.1.8 の下界が出る.

補題 11.1.7. を空でない有限アルファベット,とする.を長さの型(定義 11.1.1),重の積分布とすると

𝑃𝑛(T𝑛(𝑃))𝑃𝑛(T𝑛(𝑃))

である.

証明. まず型類全体の確率を書き下す.長さの型をとると,ならであり,これはに依らない.よって型類全体の確率は要素数にこの値を掛けたものであり,命題 11.1.6 より

𝑃𝑛(T𝑛(˜𝑃))=𝑛!𝑎(𝑛˜𝑃(𝑎))! 𝑎𝑃(𝑎)𝑛˜𝑃(𝑎)

である.この式をについて比べる.

ある文字かつである場合を先に片づける.このときの式は因子をもつので値はであり,は確率だから非負で,主張が成り立つ.

以下,どの文字でもならばである場合を見る.このときの文字ではでもあるから,二つの式の対応する因子はどちらもで一致する.また 系 11.1.5 よりで,は長さの型だからは空でなく,この値は正である.よって割ることができ,も約分されて

𝑃𝑛(T𝑛(𝑃))𝑃𝑛(T𝑛(𝑃))=𝑎:𝑃(𝑎)>0(𝑛𝑃(𝑎))!(𝑛𝑃(𝑎))! 𝑃(𝑎)𝑛𝑃(𝑎)𝑛𝑃(𝑎)

である.

各因子を階乗の比で抑える.非負整数を満たすときである.実際,ならからまでの個の整数の積で,どの因子も以下だから以下である.ならからまでの個の整数の積で,どの因子もより大きいから以上であり,逆数をとると同じ不等式になる.

の文字ごとにととる.は長さの型だからこれらは非負整数で,よりである.に注意して,正の数を両辺に掛けると

(𝑛𝑃(𝑎))!(𝑛𝑃(𝑎))! 𝑃(𝑎)𝑛𝑃(𝑎)𝑛𝑃(𝑎)𝑢𝑢𝑣(𝑢𝑛)𝑣𝑢=𝑛𝑢𝑣

を得る.これをの文字について掛け合わせると,の肩に乗るのは

𝑎:𝑃(𝑎)>0(𝑛𝑃(𝑎)𝑛𝑃(𝑎))=𝑛𝑛=0

である(第の和は全文字にわたる和に等しく,第の和も,いま見ている場合にはの文字でだからに等しい).よって比は以下であり,主張を得る.

この不等式は,大数の法則を型類の言葉で言い直したものである.に従う情報源から見ると,型がの系列は各文字の出現回数が期待どおりの系列であり,型がずれた系列はそのぶん確率を失う.いちばん重い型類が自分の型類になるのはそのためである.

形式化上の注記. 補題 11.1.7 の不等式は multinomial_mul_prod_ratio_pow_le (InformationTheory/Shannon/Sanov/MultinomialLowerBound.lean) が与えるが,型類の要素数を多項係数で書いた形になっている.要素数と多項係数が一致すること(命題 11.1.6)のほうはファイルの外に公開されていないので,本文の形そのままの単独の宣言は無い.

要素数はの前後にある

定理 11.1.8(型類の要素数). を空でない有限アルファベット,とし,を長さの型(定義 11.1.1)とすると

(𝑛+1)|X|𝑒𝑛𝐻(𝑃)T𝑛(𝑃)𝑒𝑛𝐻(𝑃)

である(定義 1.1.1 のエントロピー).

証明. 上界を示す.重の積分布とすると,系 11.1.5 よりの各点の確率はだから

𝑃𝑛(T𝑛(𝑃))=T𝑛(𝑃)𝑒𝑛𝐻(𝑃)

である.左辺は確率だから以下であり,両辺にを掛ければ上界を得る.

下界に移る.どのもただ一つの型をもつ(定義 11.1.1)から,長さの型の型類はを重なりなく覆う.よってが長さの型の全体をわたる和として

1=𝑃𝑃𝑛(T𝑛(𝑃))

が成り立つ.右辺の各項は 補題 11.1.7 より以下であり,項の個数は 命題 11.1.3 より以下だから

1(𝑛+1)|X|𝑃𝑛(T𝑛(𝑃))=(𝑛+1)|X|T𝑛(𝑃)𝑒𝑛𝐻(𝑃)

である.両辺にを掛ければ下界を得る.

型類の要素数はの前後にあり,上下の隔たりは倍である.をとるとこの隔たりはになる.第2章 定理 2.2.5 も典型集合の要素数をのまわりで評価したが,そこには典型集合を切り出すときに選んだ幅が指数に残っていた.典型集合が「経験エントロピーがに近い」という幅のある条件で切り出した集合であるのに対し,型類は一つの型で切り出した集合なので,指数に幅が残らない.

形式化上の注記. 定理 11.1.8 の上界にも下界にも,対応する単独の宣言は無い.typeClassByCount_card_letypeClassByCount_card_ge (InformationTheory/Shannon/Sanov/MultinomialLowerBound.lean) が同じ挟み込みをの形で与え,pow_div_prod_pow_eq_exp_n_entropyByCount (InformationTheory/Shannon/TypeClassLowerBound.lean) がこの量をに書き換える.二つの合成である.

経験エントロピーを型で読む

系 11.1.9. を空でない有限アルファベットとし,定義 2.1.1 の設定(i.i.d. 情報源,各に値をとり,その分布は)で,上の全点で正の分布とし,とする.ブロックの型を定義 11.1.1)と書くと,定義 2.1.1 の経験エントロピー

ˆ𝐻𝑛=𝐻(ˆ𝑃𝑋𝑛)+𝐷(ˆ𝑃𝑋𝑛𝑝)

を満たす(定義 1.1.1 のエントロピー,1.6 節の相対エントロピー).

証明. 補題 2.1.3 よりである.対数は積を和に変えるので,和を対数の中へ戻すと,重の積分布としてである.

と置く.は長さの型であり,である.は全点で正だから,命題 11.1.4に当てて

𝑝𝑛({𝑋𝑛})=exp(𝑛(𝐻(𝑃)+𝐷(𝑃𝑝)))

を得る.両辺の対数をとってを掛ければ主張になる.

系 11.1.9 は,第2章が一つの数として扱っていた経験エントロピーを二つに割る.第は,実際に出た系列の型そのもののエントロピーであり,第は,その型が真の分布からどれだけ隔たっているかである.からずれる原因が,型のエントロピーのずれと,型と真の分布の隔たりの二つに分かれる,ということでもある.

形式化上の注記. 第2章 2.1 節の形式化は経験エントロピーを文字ごとの対数尤度の相加平均として定義しており,型を経由した 系 11.1.9 の形に対応する単独の宣言は無い.

本節は型類を,中の点がどれも同じ確率をもち,が全点で正ならその値がであること(命題 11.1.4)と,要素数がの前後にあること(定理 11.1.8)の二つで押さえた.次節はこの二つを掛け合わせる.掛けるとエントロピーの項が打ち消し合い,相対エントロピーだけが指数に残る.

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