Probability distribution that has the most entropy of a class
統計学 および 情報理論 において 、 最大エントロピー確率分布は、特定の 確率分布 クラスの他のすべてのメンバーのエントロピーと少なくとも同じ大きさの エントロピー を持つ 。 最大エントロピーの原理 によれば、ある分布が特定のクラス(通常は特定の特性または尺度によって定義される)に属していること以外何も知られていない場合、最大エントロピーを持つ分布が最も情報量の少ないデフォルトとして選択されるべきである。その動機は2つある。1つは、エントロピーを最大化することで、分布に組み込まれる 事前情報 の量が最小化されること、2つ目として、多くの物理システムは時間の経過とともに最大エントロピー構成へと向かう傾向があることである。
エントロピーと微分エントロピーの定義 が確率密度 を持つ 連続確率変数 である 場合 、 の 微分エントロピーは 次のように定義される [1] [2] [3] X {\displaystyle X} p ( x ) {\displaystyle p(x)} X {\displaystyle X}
H ( X ) = − ∫ − ∞ ∞ p ( x ) log p ( x ) d x . {\displaystyle H(X)=-\int _{-\infty }^{\infty }p(x)\log p(x)\,dx~.}
が分布によって与えられる 離散確率変数 である 場合、 エントロピーは 次のように定義される。 X {\displaystyle X} Pr ( X = x k ) = p k for k = 1 , 2 , … {\displaystyle \Pr(X{=}x_{k})=p_{k}\qquad {\text{ for }}\quad k=1,2,\ldots } X {\displaystyle X} H ( X ) = − ∑ k ≥ 1 p k log p k . {\displaystyle H(X)=-\sum _{k\geq 1}p_{k}\log p_{k}\,.}
一見発散する項 は、常にゼロに置き換えられる。 p ( x ) log p ( x ) {\displaystyle p(x)\log p(x)} p ( x ) = 0 . {\displaystyle p(x)=0\,.}
これは、エントロピー(情報理論) 、 最大エントロピーの原理 、および微分エントロピーの各記事で説明されているより一般的な形式の特殊なケースです 。最大エントロピー分布に関しては、最大化によってより一般的な形式も最大化されるため、これが唯一必要な形式です 。 H ( X ) {\displaystyle H(X)}
対数 の底は 、一貫して同じ底が使われている限り重要ではありません。底の変更は、エントロピーのスケール変更に過ぎません。情報理論家は、エントロピーを ビット単位 で表すために底を2とすることを好む場合があります。一方、数学者や物理学者は、エントロピー の単位が 「nat」となる 自然対数を 好むことが多いです。
しかし、 ルベーグ測度 の典型的な使用は「自然な」選択として擁護されることが多いものの、選択された 測度 は非常に重要です。どの測度が選択されるかによって、エントロピーと、その結果として最大エントロピー分布が決まります。 d x {\displaystyle dx}
測定定数を持つ分布 適用可能な興味深い統計分布の多くは、 モーメントやその他の測定可能な量が定数となるように制約されている分布です。 ルートヴィヒ・ボルツマン による以下の定理は、 これらの制約の下での確率密度の形を示しています。
連続ケース が実数 の連続し た閉部分集合 である とし、 測定可能な関数 と 数 を指定することにします。 上でサポートされ (つまり、 の外側では密度関数がゼロになる )、モーメント条件を満たす すべての実数値ランダム変数の クラスを考えます 。 S {\displaystyle S} R {\displaystyle \mathbb {R} } n {\displaystyle n} f 1 , … , f n {\displaystyle f_{1},\ldots ,f_{n}} n {\displaystyle n} a 1 , … , a n . {\displaystyle a_{1},\ldots ,a_{n}.} C {\displaystyle C} S {\displaystyle S} S {\displaystyle S} n {\displaystyle n}
E [ f j ( X ) ] ≥ a j for j = 1 , … , n {\displaystyle \operatorname {E} [f_{j}(X)]\geq a_{j}\qquad {\text{for }}\quad j=1,\ldots ,n}
の どこにいても 密度関数 が正であるような要素があり 、 に対する最大エントロピー分布が存在する場合、 その確率密度は 次の形式になります。 C {\displaystyle C} S , {\displaystyle S,} C , {\displaystyle C,} p ( x ) {\displaystyle p(x)}
p ( x ) = exp ( ∑ j = 0 n λ j f j ( x ) ) for all x ∈ S {\displaystyle p(x)=\exp \left(\sum _{j=0}^{n}\lambda _{j}f_{j}(x)\right)\qquad {\text{ for all }}~x\in S}
ここで、定数 と ラグランジュ乗数は 制約付き最適化問題を解くと 仮定する (積分すると1になることを保証する ): [4] f 0 ( x ) = 1 . {\displaystyle f_{0}(x)=1\,.} λ 0 {\displaystyle \lambda _{0}} n {\displaystyle n} λ = ( λ 1 , … , λ n ) {\displaystyle {\boldsymbol {\lambda }}=(\lambda _{1},\ldots ,\lambda _{n})} a 0 = 1 {\displaystyle a_{0}=1} p {\displaystyle p}
max λ 0 ; λ { ∑ j = 0 n λ j a j − ∫ exp ( ∑ j = 0 n λ j f j ( x ) ) d x } subject to λ ≥ 0 {\displaystyle \max _{\lambda _{0};\,{\boldsymbol {\lambda }}}\left\{\sum _{j=0}^{n}\lambda _{j}a_{j}-\int \exp \left(\sum _{j=0}^{n}\lambda _{j}f_{j}(x)\right)dx\right\}\qquad ~{\text{ subject to }}~{\boldsymbol {\lambda }}\geq \mathbf {0} }
カルーシュ・キューン・タッカー条件 を用いると 、最適化問題には唯一の解が存在することが示せる。なぜなら最適化の目的関数は、 λ . {\displaystyle {\boldsymbol {\lambda }}\,.}
モーメント制約が不等式ではなく等式である場合、つまり、
E [ f j ( X ) ] = a j for j = 1 , … , n , {\displaystyle \operatorname {E} [f_{j}(X)]=a_{j}\qquad {\text{ for }}~j=1,\ldots ,n\,,}
制約条件 を削除すると、ラグランジュ乗数の最適化が制約されなくなります。 λ ≥ 0 {\displaystyle {\boldsymbol {\lambda }}\geq \mathbf {0} }
離散的なケース は 実数の(有限または無限の)離散部分集合であるとし、 関数 と 数を指定することにする。我々は 、上でサポートされ 、 モーメント条件を満たす すべての離散確率変数の クラスを考える。 S = { x 1 , x 2 , … } {\displaystyle S=\{x_{1},x_{2},\ldots \}} n {\displaystyle n} f 1 , … , f n {\displaystyle f_{1},\ldots ,f_{n}} n {\displaystyle n} a 1 , … , a n . {\displaystyle a_{1},\ldots ,a_{n}\,.} C {\displaystyle C} X {\displaystyle X} S {\displaystyle S} n {\displaystyle n}
E [ f j ( X ) ] ≥ a j for j = 1 , … , n {\displaystyle \operatorname {E} [f_{j}(X)]\geq a_{j}\qquad ~{\text{ for }}~j=1,\ldots ,n}
のすべてのメンバーに正の確率を割り当てる クラスのメンバーが存在し 、の最大エントロピー分布が存在する場合、 この分布は次の形状になります。 C {\displaystyle C} S {\displaystyle S} C , {\displaystyle C,}
Pr ( X = x k ) = exp ( ∑ j = 0 n λ j f j ( x k ) ) for k = 1 , 2 , … {\displaystyle \Pr(X{=}x_{k})=\exp \left(\sum _{j=0}^{n}\lambda _{j}f_{j}(x_{k})\right)\qquad {\text{ for }}~k=1,2,\ldots }
ここで、定数は 制約付き最適化問題を解くと 仮定する : [5] f 0 = 1 {\displaystyle f_{0}=1} λ 0 , λ ≡ ( λ 1 , … , λ n ) {\displaystyle \lambda _{0},\,{\boldsymbol {\lambda }}\equiv (\lambda _{1},\ldots ,\lambda _{n})} a 0 = 1 {\displaystyle a_{0}=1}
max λ 0 ; λ { ∑ j = 0 n λ j a j − ∑ k ≥ 1 exp ( ∑ j = 0 n λ j f j ( x k ) ) } for which λ ≥ 0 {\displaystyle \max _{\lambda _{0};\,{\boldsymbol {\lambda }}}\left\{\sum _{j=0}^{n}\lambda _{j}a_{j}-\sum _{k\geq 1}\exp \left(\sum _{j=0}^{n}\lambda _{j}f_{j}(x_{k})\right)\right\}\qquad {\text{ for which }}~{\boldsymbol {\lambda }}\geq \mathbf {0} }
上記と同様に、モーメント条件が不等式ではなく等式である場合、制約条件は 最適化に存在しません。 λ ≥ 0 {\displaystyle {\boldsymbol {\lambda }}\geq \mathbf {0} }
等式制約の場合の証明 等式制約の場合、この定理は 変分法 と ラグランジュ乗数 を用いて証明される。制約は次のように書ける。
∫ − ∞ ∞ f j ( x ) p ( x ) d x = a j {\displaystyle \int _{-\infty }^{\infty }f_{j}(x)p(x)\,dx=a_{j}}
機能的 な
J ( p ) = ∫ − ∞ ∞ p ( x ) ln p ( x ) d x − η 0 ( ∫ − ∞ ∞ p ( x ) d x − 1 ) − ∑ j = 1 n λ j ( ∫ − ∞ ∞ f j ( x ) p ( x ) d x − a j ) {\displaystyle J(p)=\int _{-\infty }^{\infty }p(x)\ln {p(x)}\,dx-\eta _{0}\left(\int _{-\infty }^{\infty }p(x)\,dx-1\right)-\sum _{j=1}^{n}\lambda _{j}\left(\int _{-\infty }^{\infty }f_{j}(x)p(x)\,dx-a_{j}\right)}
ここで 、 と はラグランジュ乗数である。ゼロ番目の制約は、 確率公理の第二公理 を保証する。他の制約は、関数の測定値が 次までの定数で与えられることである。 関数の微分 がゼロの とき、エントロピーは極値に達する。 η 0 {\displaystyle \eta _{0}} λ j , j ≥ 1 {\displaystyle \lambda _{j},j\geq 1} n {\displaystyle n}
δ J ( p ) δ p = ln p ( x ) + 1 − η 0 − ∑ j = 1 n λ j f j ( x ) = 0 {\displaystyle {\frac {\delta J(p)}{\delta p}}=\ln {p(x)}+1-\eta _{0}-\sum _{j=1}^{n}\lambda _{j}f_{j}(x)=0}
したがって、この場合の極値エントロピー確率分布は( )の形式でなければならない。 λ 0 := η 0 − 1 {\displaystyle \lambda _{0}:=\eta _{0}-1}
p ( x ) = e − 1 + η 0 e ∑ j = 1 n λ j f j ( x ) = exp ( ∑ j = 0 n λ j f j ( x ) ) , {\displaystyle p(x)=e^{-1+\eta _{0}}\,e^{\sum _{j=1}^{n}\lambda _{j}f_{j}(x)}=\exp \left(\sum _{j=0}^{n}\lambda _{j}f_{j}(x)\right),}
ということを覚えておいてください 。この解の周りの変化が常に負であることを確認することで、これが最大解であることが検証できます。 f 0 ( x ) = 1 {\displaystyle f_{0}(x)=1}
最大値の一意性 とが期待制約を満たす分布である と 仮定する。 分布を とすると、 この分布が期待制約を満たし、さらに が支持として存在することが明らかである。 エントロピーに関する基本的な事実から、 が成り立つ。それぞれ と の 極限をとる と 、 が得られる。 p {\displaystyle p} p ′ {\displaystyle p'} α ∈ ( 0 , 1 ) {\displaystyle \alpha \in (0,1)} q = α p + ( 1 − α ) p ′ {\displaystyle q=\alpha \,p+(1-\alpha )\,p'} supp ( q ) = supp ( p ) ∪ supp ( p ′ ) . {\displaystyle \operatorname {supp} (q)=\operatorname {supp} (p)\cup \operatorname {supp} (p')\,.} H ( q ) ≥ α H ( p ) + ( 1 − α ) H ( p ′ ) . {\displaystyle {\mathcal {H}}(q)\geq \alpha \,{\mathcal {H}}(p)+(1-\alpha )\,{\mathcal {H}}(p').} α → 1 {\displaystyle \alpha \to 1} α → 0 , {\displaystyle \alpha \to 0\,,} H ( q ) ≥ H ( p ) , H ( p ′ ) . {\displaystyle {\mathcal {H}}(q)\geq {\mathcal {H}}(p),{\mathcal {H}}(p')\,.}
したがって、期待値制約を満たし、エントロピーを最大化する分布は、必ず完全な支持を持つ、 つまり ほぼすべての点で正の分布となる必要がある。したがって、最大化分布は、期待値制約を満たす分布空間内の内部点、すなわち局所的な極値でなければならない。したがって、エントロピーを最大化する分布が唯一であることを示すには、局所的な極値が唯一であることを示すだけで十分である(そして、これは局所的な極値が全体的最大値であることも示す)。
とが局所的極値であると 仮定する 。上記の計算を再定式化すると、これらは と を介してパラメータによって特徴付けられる。ここでも 同様である。 ここで 一連の恒等式に注意する。期待値制約の充足と勾配/方向微分を用いることで、次式が得られる。 p {\displaystyle p} p ′ {\displaystyle p'} λ , λ ′ ∈ R n {\displaystyle {\boldsymbol {\lambda }},\,{\boldsymbol {\lambda }}'\in \mathbb {R} ^{n}} p ( x ) = exp ⟨ λ , f ( x ) ⟩ / C ( λ ) {\displaystyle p(x)={\exp \left\langle {\boldsymbol {\lambda }},\mathbf {f} (x)\right\rangle }/{C({\boldsymbol {\lambda }})}} p ′ , {\displaystyle p',} C ( λ ) = ∫ R exp ⟨ λ , f ( x ) ⟩ d x . {\textstyle C({\boldsymbol {\lambda }})=\int _{\mathbb {R} }\exp \left\langle {\boldsymbol {\lambda }},\mathbf {f} (x)\right\rangle \,dx\,.}
D log C ( ⋅ ) | λ = D C ( ⋅ ) C ( ⋅ ) | λ = E p [ f ( X ) ] = a {\displaystyle {\left.D\log C(\cdot )\right|}_{\boldsymbol {\lambda }}={\left.{\tfrac {DC(\cdot )}{C(\cdot )}}\right|}_{\boldsymbol {\lambda }}=\operatorname {E} _{p}\left[\mathbf {f} (X)\right]=\mathbf {a} } 同様にLetting についても 次のようになります。 λ ′ . {\displaystyle {\boldsymbol {\lambda }}'~.} u = λ ′ − λ ∈ R n {\displaystyle u={\boldsymbol {\lambda }}'-{\boldsymbol {\lambda }}\in \mathbb {R} ^{n}}
0 = ⟨ u , a − a ⟩ = D u log C ( ⋅ ) | λ ′ − D u log C ( ⋅ ) | λ = D u 2 log C ( ⋅ ) | γ {\displaystyle 0=\left\langle u,\mathbf {a} -\mathbf {a} \right\rangle ={\left.D_{u}\log C(\cdot )\right|}_{{\boldsymbol {\lambda }}'}-{\left.D_{u}\log C(\cdot )\right|}_{\boldsymbol {\lambda }}={\left.D_{u}^{2}\log C(\cdot )\right|}_{\boldsymbol {\gamma }}}
さらに計算を行う と 、 γ = θ λ + ( 1 − θ ) λ ′ {\displaystyle {\boldsymbol {\gamma }}=\theta {\boldsymbol {\lambda }}+(1-\theta ){\boldsymbol {\lambda }}'} θ ∈ ( 0 , 1 ) . {\displaystyle \theta \in (0,1).}
0 = D u 2 log C ( ⋅ ) | γ = D u ( D u C ( ⋅ ) C ( ⋅ ) ) | γ = D u 2 C ( ⋅ ) C ( ⋅ ) | γ − ( D u C ( ⋅ ) ) 2 C ( ⋅ ) 2 | γ = E q [ ⟨ u , f ( X ) ⟩ 2 ] − ( E q [ ⟨ u , f ( X ) ⟩ ] ) 2 = Var q [ ⟨ u , f ( X ) ⟩ ] {\displaystyle {\begin{aligned}0&={\left.D_{u}^{2}\log C(\cdot )\right|}_{\boldsymbol {\gamma }}\\[1ex]&={\left.D_{u}\left({\frac {D_{u}C(\cdot )}{C(\cdot )}}\right)\right|}_{\boldsymbol {\gamma }}={\left.{\frac {D_{u}^{2}C(\cdot )}{C(\cdot )}}\right|}_{\boldsymbol {\gamma }}-{\left.{\frac {{\left(D_{u}C(\cdot )\right)}^{2}}{C(\cdot )^{2}}}\right|}_{\boldsymbol {\gamma }}\\[1ex]&=\operatorname {E} _{q}\left[{\left\langle u,\mathbf {f} (X)\right\rangle }^{2}\right]-{\left(\operatorname {E} _{q}\left[\left\langle u,\mathbf {f} (X)\right\rangle \right]\right)}^{2}\\[2ex]&=\operatorname {Var} _{q}\left[\left\langle u,\mathbf {f} (X)\right\rangle \right]\end{aligned}}}
ここで 、は上記の分布に似ていますが、パラメータ化は によってのみ行われます。 観測量の非自明な線形結合は、 ほとんどどこでも (ae) 一定ではないと仮定します( 例えば 、 観測量が独立でaeが一定でない場合などに当てはまります)。 の 分散はゼロではありません。ただし、 の場合は、 上記の式から、後者が当てはまらなければならないことは明らかです。したがって、 局所的極値を特徴付けるパラメータは 同一であり、分布自体も同一であることを意味します。したがって、局所的極値は一意であり、上記の議論により、最大値も一意です(ただし、局所的極値が実際に存在する場合)。 q {\displaystyle q} γ , {\displaystyle {\boldsymbol {\gamma }}~,} ⟨ u , f ( X ) ⟩ {\displaystyle \langle u,\mathbf {f} (X)\rangle } u = 0 . {\displaystyle u=0~.} λ ′ − λ = u = 0 , {\displaystyle {\boldsymbol {\lambda }}'-{\boldsymbol {\lambda }}=u=0\,,} p , p ′ {\displaystyle p,\,p'}
注意点 すべての分布クラスが最大エントロピー分布を含むわけではないことに注意してください。クラスによっては、任意の大きさのエントロピーを持つ分布(例えば、平均が0で標準偏差が任意の R 上の連続分布のクラス)を含む場合や、エントロピーが上限で定義されていても最大エントロピーを達成する分布が存在しない場合があります。 [a] また、クラス Cの期待値制約により、 S の特定の部分集合において確率分布が0になる場合もあります。その場合、定理は適用されませんが、集合 S を縮小することでこの問題を回避できます 。
例 あらゆる確率分布は、その分布が自身のエントロピーを持つという制約の下では、自明に最大エントロピー確率分布となる。これを確認するには、密度を と書き直し 、上記の定理の式と比較する。 を測定可能な関数と とすることで、 p ( x ) = exp ( ln p ( x ) ) {\displaystyle p(x)=\exp {(\ln {p(x)})}} ln p ( x ) → f ( x ) {\displaystyle \ln {p(x)}\rightarrow f(x)}
∫ exp ( f ( x ) ) f ( x ) d x = − H {\displaystyle \int \exp {(f(x))}f(x)dx=-H}
定数となるのは、 制約条件の下での最大エントロピー確率分布である。 p ( x ) {\displaystyle p(x)}
∫ p ( x ) f ( x ) d x = − H . {\displaystyle \int p(x)f(x)\,dx=-H.}
非自明な例としては、エントロピーの割り当てとは異なる複数の制約を受ける分布が挙げられます。これらは、同じ手順から始めて 、 複数の部分に分割できることを発見することで見つかることがよくあります。 ln p ( x ) → f ( x ) {\displaystyle \ln {p(x)}\to f(x)} f ( x ) {\displaystyle f(x)}
最大エントロピー分布の例の表は、Lisman (1972) [6] とPark & Bera (2009) [7]に示されている。
区間 [ a , b ] 上の均一 分布は 、区間 [ a , b ] でサポートされるすべての連続分布の中で最大エントロピー分布であるため、区間外の確率密度は 0 になります。この均一密度は、ラプラスの 無差別原理 ( 不十分な理由の原理と呼ばれることもある) に関連付けることができます。より一般的には、 区間 [ a , b ] の細分 a = a 0 < a 1 < ... < a k = b と、合計が 1 になる 確率 p 1 、...、 p k が与えられている場合、すべての連続分布のクラスを次のように考えることができます
。このクラスの最大エントロピー分布の密度は、各区間 [ a j −1 、 a j ] で一定です。有限集合 { x 1 、...、 x n } (これらの値のそれぞれに 1/ n の確率を割り当てる ) 上の均一分布は、この集合でサポートされるすべての離散分布の中で最大エントロピー分布です。 Pr ( a j − 1 ≤ X < a j ) = p j for j = 1 , … , k {\displaystyle \Pr(a_{j-1}\leq X<a_{j})=p_{j}\quad {\text{ for }}j=1,\ldots ,k}
正の指定平均:指数分布 指数 分布 の密度関数は
p ( x | λ ) = { λ e − λ x x ≥ 0 , 0 x < 0 , {\displaystyle p(x|\lambda )={\begin{cases}\lambda e^{-\lambda x}&x\geq 0,\\0&x<0,\end{cases}}}
は、[0,∞)でサポートされている、指定された平均が1/λであるすべての連続分布の中で最大のエントロピー分布です。
[0,∞) でサポートされる分布の場合、最大エントロピー分布は第一モーメントと第二モーメントの関係に依存します。特定のケースでは、指数分布となる場合もあれば、他の分布となる場合もあり、あるいは定義不可能な場合もあります。 [8]
指定された平均と分散:正規分布 正規 分布 N ( μ , σ 2 )の密度関数は
p ( x | μ , σ ) = 1 σ 2 π e − ( x − μ ) 2 2 σ 2 , {\displaystyle p(x|\mu ,\sigma )={\frac {1}{\sigma {\sqrt {2\pi }}}}e^{-{\frac {(x-\mu )^{2}}{2\sigma ^{2}}}},}
は、指定された分散 σ 2 (特定の モーメント )を持つ、(−∞,∞) 上で支持されるすべての 実 数値分布の中で最大のエントロピーを持つ。 平均 μ と 分散 σ 2 が指定された場合(最初の2つのモーメント)も同様である 。これは、エントロピーが(−∞,∞) 上で並進不変であるためである。したがって、正規性の仮定は、これらのモーメントを超えると、最小の事前構造制約を課す。( 導出については 微分エントロピーの記事を参照のこと。)
指定された平均を持つ離散分布 指定された平均 μ を持つ 集合 { x 1 ,..., x n } でサポートされているすべての離散分布の中で、最大エントロピー分布は次の形状になります。
ここで、正の定数 C と r は 、すべての確率の合計が 1 で、期待値が μ であるという要件によって決定できます。 Pr ( X = x k ) = C r x k for k = 1 , … , n {\displaystyle \Pr(X{=}x_{k})=Cr^{x_{k}}\quad {\text{ for }}k=1,\ldots ,n}
例えば、 N 個のサイコロを投げ、出た目の合計が S であると告げられたとします。この情報だけに基づいて、1、2、…、6が出るサイコロの個数について、どのような仮定が妥当でしょうか?これは、{ x 1 ,…, x 6 } = {1,…,6}、 μ = S / N の場合の、上記で検討した状況の一例です。
最後に、平均 μ の無限集合上で支持されるすべての離散分布の中で 、最大エントロピー分布は次の形をとる。 ここでも定数 C と rは、すべての確率の合計が1で、期待値がμであるという要件によって決定される。例えば、 x k = k の場合 、次の式が成り立つ
。 { x 1 , x 2 , . . . } {\displaystyle \{x_{1},x_{2},...\}} Pr ( X = x k ) = C r x k for k = 1 , 2 , … , {\displaystyle \Pr(X{=}x_{k})=Cr^{x_{k}}\quad {\text{ for }}k=1,2,\ldots ,} C = 1 μ − 1 , r = μ − 1 μ , {\displaystyle C={\frac {1}{\mu -1}},\quad \quad r={\frac {\mu -1}{\mu }},}
それぞれの最大エントロピー分布は 幾何分布 となる。
循環確率変数 単位円の周りに分布する 連続ランダム変数の場合、 フォン・ミーゼス分布は、第1 円モーメント の実部と虚部が指定されている ときにエントロピーを最大化します [9] 。または、同等に、 円平均 と 円分散 が指定されています。 θ i {\displaystyle \theta _{i}}
角度の平均と分散を法として指定すると 、 ラップ された正規分布は エントロピーを最大化します。 [9] θ i {\displaystyle \theta _{i}} 2 π {\displaystyle 2\pi }
指定された平均、分散、歪度の最大化 指定された平均、分散、歪度を持つ 連続確率変数のエントロピーには上限が存在する。しかし、 この上限を満たす分布は存在しない 。なぜなら、のとき は非有界だ からである (Cover & Thomas (2006: 第12章)を参照)。 R {\displaystyle \mathbb {R} } p ( x ) = c exp ( λ 1 x + λ 2 x 2 + λ 3 x 3 ) {\displaystyle p(x)=c\exp {(\lambda _{1}x+\lambda _{2}x^{2}+\lambda _{3}x^{3})}} λ 3 ≠ 0 {\displaystyle \lambda _{3}\neq 0}
しかし、最大エントロピーは ε で達成可能です。つまり、分布のエントロピーは上限に任意に近い値を取ることができます。まず、指定された平均と分散を持つ正規分布から始めます。正の歪度を導入するには、正規分布を平均より σ 倍も大きい値で上方にわずかに摂動させます。歪度は3次モーメントに比例するため、低次のモーメントよりも影響を受けます。
これは、 x の任意の奇数次多項式の指数関数が 上で有界にならない一般的なケースの特殊なケースです 。例えば、 は 上で同様に有界ではありません が、台が有界または半有界の区間に限定されている場合、エントロピーの上限が達成される可能性があります(例えば、 x が 区間 [0,∞] にあり、 λ< 0 の場合、 指数分布 が成立します)。 R {\displaystyle \mathbb {R} } c e λ x {\displaystyle ce^{\lambda x}} R {\displaystyle \mathbb {R} }
指定された平均と偏差のリスク測定の最大化 対数凹密 度を持つすべての分布は 、指定された平均 μ と 偏差リスク尺度 D を持つ最大エントロピー分布である 。 [10]
特に、指定された平均 と偏差を持つ最大エントロピー分布は 次のようになります。 E ( X ) ≡ μ {\displaystyle E(X)\equiv \mu } D ( X ) ≡ d {\displaystyle D(X)\equiv d}
標準偏差 が の 場合の 正規 分布 ; N ( m , d 2 ) , {\displaystyle {\mathcal {N}}(m,d^{2}),} D ( X ) = E [ ( X − μ ) 2 ] {\textstyle D(X)={\sqrt {\operatorname {E} \left[{\left(X-\mu \right)}^{2}\right]}}} ラプラス 分布は 、 平均絶対偏差 が の場合 、 [6] D ( X ) = E { | X − μ | } {\displaystyle D(X)=\operatorname {E} \left\{\left|X-\mu \right|\right\}} の形の分布の密度は 標準下側半偏差であり、 は 定数であり、関数は 引数の負の値のみを返し、それ以外の場合はゼロを返す。 [10] f ( x ) = c exp ( a x + b ⇂ x − μ ⇃ − 2 ) {\displaystyle f(x)=c\exp \left(ax+b{\downharpoonright x-\mu \downharpoonleft }_{-}^{2}\right)} D ( X ) = E { ⇂ X − μ ⇃ − 2 } {\textstyle D(X)={\sqrt {\operatorname {E} \left\{{\downharpoonright X-\mu \downharpoonleft }_{-}^{2}\right\}}}} a , b , c {\displaystyle a,b,c} ⇂ y ⇃ − ≡ min { 0 , y } for any y ∈ R , {\displaystyle \downharpoonright y\downharpoonleft _{-}\equiv \min \left\{0,y\right\}~{\text{ for any }}y\in \mathbb {R} \,,}
その他の例 以下の表では、記載されている各分布は、3列目にリストされている特定の機能制約セットと、4列目にリストされている確率密度のサポートに含まれる制約のエントロピーを 最大化します。 [6] [7] x {\displaystyle x}
列挙されているいくつかの例( ベルヌーイ 、 幾何 、 指数 、 ラプラス 、 パレート )は、関連する制約がエントロピーの割り当てと等価であるため、自明に成り立ちます。これらの制約は、一般的な量、あるいは容易に測定できる量に関連しているため、いずれにせよここに含めています。
参考までに、 は ガンマ関数 、 は ディガンマ関数 、 は ベータ関数 、は オイラー・マスケロニ定数 です 。 Γ ( x ) = ∫ 0 ∞ e − t t x − 1 d t {\displaystyle \Gamma (x)=\int _{0}^{\infty }e^{-t}t^{x-1}\,dt} ψ ( x ) = d d x ln Γ ( x ) = Γ ′ ( x ) Γ ( x ) {\displaystyle \psi (x)={\frac {d}{dx}}\ln \Gamma (x)={\frac {\Gamma '(x)}{\Gamma (x)}}} B ( p , q ) = Γ ( p ) Γ ( q ) Γ ( p + q ) {\displaystyle B(p,q)={\frac {\Gamma (p)\,\Gamma (q)}{\Gamma (p+q)}}} γ E {\displaystyle \gamma _{\mathsf {E}}}
確率分布とそれに対応する最大エントロピー制約の表 配布名 確率密度/質量関数 最大エントロピー制約 サポート 均一(離散的) f ( k ) = 1 b − a + 1 {\displaystyle f(k)={\frac {1}{b-a+1}}} なし { a , a + 1 , . . . , b − 1 , b } {\displaystyle \{a,a+1,...,b-1,b\}} 均一(連続) f ( x ) = 1 b − a {\displaystyle f(x)={\frac {1}{b-a}}} なし [ a , b ] {\displaystyle [a,b]} ベルヌーイ f ( k ) = p k ( 1 − p ) 1 − k {\displaystyle f(k)=p^{k}{\left(1-p\right)}^{1-k}} E [ K ] = p {\displaystyle \operatorname {E} [K]=p} { 0 , 1 } {\displaystyle \{0,1\}} 幾何学的 f ( k ) = ( 1 − p ) k − 1 p {\displaystyle f(k)={\left(1-p\right)}^{k-1}p} E [ K ] = 1 p {\displaystyle \operatorname {E} [K]={\frac {1}{p}}} N ∖ { 0 } = { 1 , 2 , 3 , . . . } {\displaystyle \mathbb {N} \setminus \left\{0\right\}=\{1,2,3,...\}} 指数関数 f ( x ) = λ exp ( − λ x ) {\displaystyle f(x)=\lambda \exp \left(-\lambda x\right)} E [ X ] = 1 λ {\displaystyle \operatorname {E} [X]={\frac {1}{\lambda }}} [ 0 , ∞ ) {\displaystyle [0,\infty )} ラプラス f ( x ) = 1 2 b exp ( − | x − μ | b ) {\displaystyle f(x)={\frac {1}{2b}}\exp \left(-{\frac {|x-\mu |}{b}}\right)} E [ | X − μ | ] = b {\displaystyle \operatorname {E} [|X-\mu |]=b} R {\displaystyle \mathbb {R} } 非対称ラプラス f ( x ) = λ exp [ − ( x − m ) λ s κ s ] κ + 1 κ {\displaystyle f(x)={\frac {\lambda \exp \left[-(x-m)\,\lambda \,s\,\kappa ^{s}\right]}{\kappa +{\frac {1}{\kappa }}}}} どこ s ≡ sgn ( x − m ) {\displaystyle ~s\equiv \operatorname {sgn}(x-m)} E [ ( X − m ) s κ s ] = 1 λ {\displaystyle \operatorname {E} \left[(X-m)s\kappa ^{s}\right]={\frac {1}{\lambda }}} R {\displaystyle \mathbb {R} } パレート f ( x ) = α x m α x α + 1 {\displaystyle f(x)={\frac {\alpha x_{m}^{\alpha }}{x^{\alpha +1}}}} E [ ln X ] = 1 α + ln ( x m ) {\displaystyle \operatorname {E} [\ln X]={\frac {1}{\alpha }}+\ln(x_{m})} [ x m , ∞ ) {\displaystyle [x_{m},\infty )} 普通 f ( x ) = 1 2 π σ 2 exp ( − ( x − μ ) 2 2 σ 2 ) {\displaystyle f(x)={\frac {1}{\sqrt {2\pi \sigma ^{2}}}}\exp \left(-{\frac {(x-\mu )^{2}}{2\sigma ^{2}}}\right)} E [ X ] = μ , {\displaystyle \operatorname {E} [X]=\mu \,,} E [ X 2 ] = σ 2 + μ 2 {\displaystyle \operatorname {E} \left[X^{2}\right]=\sigma ^{2}+\mu ^{2}} R {\displaystyle \mathbb {R} } 切り捨て正規分布 (記事参照) E [ X ] = μ T , {\displaystyle \operatorname {E} \left[X\right]=\mu _{\mathsf {T}}\,,} E [ X 2 ] = σ T 2 + μ T 2 {\displaystyle \operatorname {E} \left[X^{2}\right]=\sigma _{\mathsf {T}}^{2}+\mu _{\mathsf {T}}^{2}} [ a , b ] {\displaystyle [a,b]} フォン・ミーゼス f ( θ ) = 1 2 π I 0 ( κ ) exp ( κ cos ( θ − μ ) ) {\displaystyle f(\theta )={\frac {1}{2\pi I_{0}(\kappa )}}\exp \left(\kappa \cos {(\theta -\mu )}\right)} E [ cos Θ ] = I 1 ( κ ) I 0 ( κ ) cos μ , {\displaystyle \operatorname {E} [\cos \Theta ]={\frac {I_{1}(\kappa )}{I_{0}(\kappa )}}\cos \mu \,,} E [ sin Θ ] = I 1 ( κ ) I 0 ( κ ) sin μ {\displaystyle \operatorname {E} [\sin \Theta ]={\frac {I_{1}(\kappa )}{I_{0}(\kappa )}}\sin \mu } [ 0 , 2 π ) {\displaystyle [0,2\pi )} レイリー f ( x ) = x σ 2 exp ( − x 2 2 σ 2 ) {\displaystyle f(x)={\frac {x}{\sigma ^{2}}}\exp \left(-{\frac {x^{2}}{2\sigma ^{2}}}\right)} E [ X 2 ] = 2 σ 2 , {\displaystyle \operatorname {E} \left[X^{2}\right]=2\sigma ^{2}\,,} E [ ln X ] = ln ( 2 σ 2 ) − γ E 2 {\displaystyle \operatorname {E} [\ln X]={\frac {\ln(2\sigma ^{2})-\gamma _{\mathrm {E} }}{2}}} [ 0 , ∞ ) {\displaystyle [0,\infty )} ベータ f ( x ) = x α − 1 ( 1 − x ) β − 1 B ( α , β ) {\displaystyle f(x)={\frac {x^{\alpha -1}(1-x)^{\beta -1}}{B(\alpha ,\beta )}}} のために 0 ≤ x ≤ 1 {\displaystyle 0\leq x\leq 1} E [ ln X ] = ψ ( α ) − ψ ( α + β ) , {\displaystyle \operatorname {E} [\ln X]=\psi (\alpha )-\psi (\alpha +\beta )\,,} E [ ln ( 1 − X ) ] = ψ ( β ) − ψ ( α + β ) {\displaystyle \operatorname {E} [\ln(1-X)]=\psi (\beta )-\psi (\alpha +\beta )} [ 0 , 1 ] {\displaystyle [0,1]} コーシー f ( x ) = 1 π ( 1 + x 2 ) {\displaystyle f(x)={\frac {1}{\pi \left(1+x^{2}\right)}}} E [ ln ( 1 + X 2 ) ] = 2 ln 2 {\displaystyle \operatorname {E} \left[\ln \left(1+X^{2}\right)\right]=2\ln 2} R {\displaystyle \mathbb {R} } チ f ( x ) = 2 2 k / 2 Γ ( k / 2 ) x k − 1 exp ( − x 2 2 ) {\displaystyle f(x)={\frac {2}{2^{k/2}\Gamma (k/2)}}x^{k-1}\exp \left(-{\frac {x^{2}}{2}}\right)} E [ X 2 ] = k , {\displaystyle \operatorname {E} \left[X^{2}\right]=k\,,} E [ ln X ] = 1 2 [ ψ ( k 2 ) + ln ( 2 ) ] {\displaystyle \operatorname {E} [\ln X]={\frac {1}{2}}\left[\psi {\left({\frac {k}{2}}\right)}+\ln(2)\right]} [ 0 , ∞ ) {\displaystyle [0,\infty )} カイ二乗 f ( x ) = 1 2 k / 2 Γ ( k / 2 ) x k 2 − 1 exp ( − x 2 ) {\displaystyle f(x)={\frac {1}{2^{k/2}\Gamma (k/2)}}x^{{\frac {k}{2}}\!-\!1}\exp \left(-{\frac {x}{2}}\right)} E [ X ] = k , {\displaystyle \operatorname {E} [X]=k\,,} E [ ln X ] = ψ ( k 2 ) + ln ( 2 ) {\displaystyle \operatorname {E} [\ln X]=\psi {\left({\frac {k}{2}}\right)}+\ln(2)} [ 0 , ∞ ) {\displaystyle [0,\infty )} アーラン f ( x ) = λ k ( k − 1 ) ! x k − 1 exp ( − λ x ) {\displaystyle f(x)={\frac {\lambda ^{k}}{(k-1)!}}x^{k-1}\exp(-\lambda x)} E [ X ] = k / λ , {\displaystyle \operatorname {E} [X]=k/\lambda \,,} E [ ln X ] = ψ ( k ) − ln ( λ ) {\displaystyle \operatorname {E} [\ln X]=\psi (k)-\ln(\lambda )}
[ 0 , ∞ ) {\displaystyle [0,\infty )} ガンマ f ( x ) = x k − 1 e − x / θ θ k Γ ( k ) {\displaystyle f(x)={\frac {x^{k-1}e^{-x/\theta }}{\theta ^{k}\Gamma (k)}}} E [ X ] = k θ , {\displaystyle \operatorname {E} [X]=k\theta \,,} E [ ln X ] = ψ ( k ) + ln θ {\displaystyle \operatorname {E} [\ln X]=\psi (k)+\ln \theta } [ 0 , ∞ ) {\displaystyle [0,\infty )} 対数正規分布 f ( x ) = 1 σ x 2 π exp ( − ( ln x − μ ) 2 2 σ 2 ) {\displaystyle f(x)={\frac {1}{\sigma x{\sqrt {2\pi }}}}\exp \left(-{\frac {{\left(\ln x-\mu \right)}^{2}}{2\sigma ^{2}}}\right)} E [ ln X ] = μ , {\displaystyle \operatorname {E} \left[\ln X\right]=\mu \,,} E [ ln ( X ) 2 ] = σ 2 + μ 2 {\displaystyle \operatorname {E} \left[\ln(X)^{2}\right]=\sigma ^{2}+\mu ^{2}} ( 0 , ∞ ) {\displaystyle (0,\infty )} マクスウェル・ボルツマン f ( x ) = 1 a 3 2 π x 2 exp ( − x 2 2 a 2 ) {\displaystyle f(x)={\frac {1}{a^{3}}}{\sqrt {\frac {2}{\pi }}}x^{2}\exp \left(-{\frac {x^{2}}{2a^{2}}}\right)} E [ X 2 ] = 3 a 2 , {\displaystyle \operatorname {E} \left[X^{2}\right]=3a^{2}\,,} E [ ln X ] = 1 + ln ( a 2 ) − γ E 2 {\displaystyle \operatorname {E} \left[\ln X\right]=1+\ln \left({\frac {a}{\sqrt {2}}}\right)-{\frac {\gamma _{\mathrm {E} }}{2}}} [ 0 , ∞ ) {\displaystyle [0,\infty )} ワイブル f ( x ) = k λ k x k − 1 exp ( − x k λ k ) {\displaystyle f(x)={\frac {k}{\lambda ^{k}}}x^{k-1}\exp \left(-{\frac {x^{k}}{\lambda ^{k}}}\right)} E [ X k ] = λ k , {\displaystyle \operatorname {E} \left[X^{k}\right]=\lambda ^{k},} E [ ln X ] = ln ( λ ) − γ E k {\displaystyle \operatorname {E} \left[\ln X\right]=\ln(\lambda )-{\frac {\gamma _{\mathrm {E} }}{k}}} [ 0 , ∞ ) {\displaystyle [0,\infty )} 多変量正規分布 f X ( x ) = exp [ − 1 2 ( x − μ ) T Σ − 1 ( x − μ ) ] ( 2 π ) N | Σ | {\displaystyle f_{X}(\mathbf {x} )={\frac {\exp \left[-{\frac {1}{2}}\left(\mathbf {x} -{\boldsymbol {\mu }}\right)^{\mathsf {T}}\Sigma ^{-1}\left(\mathbf {x} -{\boldsymbol {\mu }}\right)\right]}{\sqrt {{\left(2\pi \right)}^{N}\left|\Sigma \right|}}}} E [ x ] = μ , {\displaystyle \operatorname {E} \left[\mathbf {x} \right]={\boldsymbol {\mu }},} E [ ( x − μ ) ( x − μ ) T ] = Σ {\displaystyle \operatorname {E} \left[(\mathbf {x} -{\boldsymbol {\mu }})(\mathbf {x} -{\boldsymbol {\mu }})^{\mathsf {T}}\right]=\Sigma } R n {\displaystyle \mathbb {R} ^{n}} 二項式 f ( k ) = ( n k ) p k ( 1 − p ) n − k {\displaystyle f(k)={\binom {n}{k}}p^{k}{\left(1-p\right)}^{n-k}} E [ X ] = μ , {\displaystyle \operatorname {E} [X]=\mu \,,} f ∈ {\displaystyle f\in } n一般化二項分布 [11] { 0 , … , n } {\displaystyle \left\{0,{\ldots },n\right\}} ポアソン f ( k ) = λ k e − λ k ! {\displaystyle f(k)={\frac {\lambda ^{k}e^{-\lambda }}{k!}}} E [ X ] = λ , {\displaystyle \operatorname {E} [X]=\lambda ,} f ∈ ∞ {\displaystyle f\in \infty } 一般化二項分布 [11] N = { 0 , 1 , … } {\displaystyle \mathbb {N} =\left\{0,1,{\ldots }\right\}} ロジスティクス f ( x ) = e − x ( 1 + e − x ) 2 = e + x ( e + x + 1 ) 2 {\displaystyle f(x)={\frac {e^{-x}}{\left(1+e^{-x}\right)^{2}}}={\frac {e^{+x}}{\left(e^{+x}+1\right)^{2}}}} E [ X ] = 0 , {\displaystyle \operatorname {E} [X]=0,} E [ ln ( 1 + e − X ) ] = 1 {\displaystyle \operatorname {E} \left[\ln \left(1+e^{-X}\right)\right]=1} { − ∞ , ∞ } {\displaystyle \left\{-\infty ,\infty \right\}}
最大エントロピー原理は、統計的混合のエントロピーの上限を決定するために使用することができる。 [12]
参照
注記 ^ 例えば、 E( X )=0 かつ E( X2 )=E( X3 )=1を 満たす R 上 の すべての連続分布X のクラス(表紙、第12章を参照) 。
引用 ^ ウィリアムズ、D. (2001). 『Weighing the Odds 』 ケンブリッジ大学出版局 . pp. 197– 199. ISBN 0-521-00618-X 。 ^ ベルナルド, JM; スミス, AFM (2000). ベイズ理論 . ワイリー. pp. 209, 366. ISBN 0-471-49464-X 。 ^ O'Hagan, A. (1994), ベイズ推論 . Kendall's Advanced Theory of Statistics. 第2B巻. Edward Arnold . 1994年. セクション5.40. ISBN 0-340-52922-9 。 ^ Botev, ZI; Kroese, DP (2011). 「一般化クロスエントロピー法とその確率密度推定への応用」 (PDF) . 応用確率論における方法論と計算 . 13 (1): 1– 27. doi :10.1007/s11009-009-9133-7. S2CID 18155189. ^ Botev, ZI; Kroese, DP (2008). 「離散データの密度推定のための非漸近的帯域幅選択」. Methodology and Computing in Applied Probability . 10 (3): 435. doi :10.1007/s11009-007-9057-zv (2025年7月1日現在非アクティブ). S2CID 122047337. {{cite journal }}: CS1 maint: DOI inactive as of July 2025 (link )^ abc Lisman, JHC; van Zuylen, MCA (1972). 「最確頻度分布の生成に関する注記」. Statistica Neerlandica . 26 (1): 19– 23. doi :10.1111/j.1467-9574.1972.tb00152.x. ^ ab Park, Sung Y.; Bera, Anil K. (2009). 「最大エントロピー自己回帰条件付き異分散モデル」 (PDF) . Journal of Econometrics . 150 (2): 219– 230. CiteSeerX 10.1.1.511.9750 . doi :10.1016/j.jeconom.2008.12.014. オリジナル (PDF) から2016年3月7日にアーカイブ。 2011年6月2日 閲覧 。 ^ Dowson, D.; Wragg, A. (1973年9月). 「規定された第一および第二モーメントを持つ最大エントロピー分布」. IEEE Transactions on Information Theory (通信). 19 (5): 689– 693. doi :10.1109/tit.1973.1055060. ISSN 0018-9448. ^ ab Jammalamadaka, S. Rao; SenGupta, A. (2001). 循環統計のトピックス. ニュージャージー州: World Scientific. ISBN 978-981-02-3778-3 . 2011年5月15日 閲覧 。 ^ ab Grechuk, Bogdan; Molyboha, Anton; Zabarankin, Michael (2009). 「一般偏差尺度による最大エントロピー原理」. オペレーションズ・リサーチ数学 . 34 (2): 445– 467. doi :10.1287/moor.1090.0377 – researchgate.net経由. ^ ab Harremös, Peter (2001). 「最大エントロピー分布としての二項分布とポアソン分布」. IEEE Transactions on Information Theory . 47 (5): 2039– 2041. doi :10.1109/18.930936. S2CID 16171405. ^ Nielsen, Frank; Nock, Richard (2017). 「単変量連続分布の微分エントロピーのMaxEnt上限値」. IEEE Signal Processing Letters . 24 (4). IEEE : 402– 406. Bibcode :2017ISPL...24..402N. doi :10.1109/LSP.2017.2666792. S2CID 14092514.
参考文献 Cover, TM ; Thomas, JA (2006). 「第12章 最大エントロピー」 (PDF) . 情報理論の要素 (第2版). Wiley. ISBN 978-0471241959 。 F. Nielsen, R. Nock (2017)、 単変量連続分布の微分エントロピーのMaxEnt上限 、IEEE信号処理レター、24(4)、402–406 IJ Taneja (2001) 「一般化情報尺度とその応用」 第1章 Nader Ebrahimi, Ehsan S. Soofi, Refik Soyer (2008)、「多変量最大エントロピーの同定、変換、および依存性」、 Journal of Multivariate Analysis 99: 1217–1231、 doi :10.1016/j.jmva.2007.08.004
離散 一変数
連続 一変量
制限された間隔 でサポートされている 半無限 間隔 でサポートされている 実数直線 全体で サポートされている さまざまなタイプの サポート付き
混合 単変量
多変量 (ジョイント) 方向性 退化 と 特異性 家族