公平である必要のないコイントスをモデル化した確率分布
確率論と統計学において、スイスの数学者ヤコブ・ベルヌーイにちなんで名付けられたベルヌーイ分布[1]は、確率 で値 1 をとり、確率 で値 0 をとるランダム変数の離散確率分布です。より形式的には、はい/いいえの質問をする単一の実験の可能な結果の集合のモデルと考えることができます。このような質問は、ブール値の結果、つまり、確率pで成功/はい/真/ 1 、確率qで失敗/いいえ/偽/ 0となる単一のビットの結果をもたらします。これは、1 と 0 がそれぞれ「表」と「裏」を表し、p がコインが表になる確率 (またはその逆で、1 が裏を表し、p が裏の確率) となる (おそらく偏りのある) コイントスを表すために使用できます。特に、不公平なコインは
ベルヌーイ分布は、1回の試行が行われる二項分布の特殊なケースです(したがって、このような二項分布ではnは1になります)。また、 2点分布の特殊なケースでもあり、可能な結果は0と1である必要はありません。
[2]
プロパティ
がベルヌーイ分布に従うランダム変数である
場合、次のようになります。

この分布の確率質量関数は、可能な結果k に対して、

[3]
これは次のようにも表現できる。

または

ベルヌーイ分布は二項分布の特殊なケースであり、 [4]
の値が高かったり低かったりすると尖度は無限大になりますが、ベルヌーイ分布を含む 2 点分布では、他のどの確率分布よりも過剰尖度が低く、つまり -2 になります。


のベルヌーイ分布は指数族を形成します。

ランダムサンプルに基づくの最大尤度推定値はサンプル平均です。

ベルヌーイ実験の確率質量分布関数とそれに対応する累積分布関数。
平均
ベルヌーイ確率変数の期待値は

![{\displaystyle \operatorname {E} [X]=p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0eb41a45634ab84b13b83cb1488b626aa2129285)
これは、ベルヌーイ分布の確率
変数に対して、


[3]
分散
ベルヌーイ分布の
分散は
![{\displaystyle \operatorname {Var} [X]=pq=p(1-p)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2d4e26d8a1fdfb90e91a2fafd5fb3841de88f1fb)
まず最初に
![{\displaystyle \operatorname {E} [X^{2}]=\Pr(X=1)\cdot 1^{2}+\Pr(X=0)\cdot 0^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/789d8c642c8c2a4da18ebe6833b6fdc7b97ce8ab)
![{\displaystyle =p\cdot 1^{2}+q\cdot 0^{2}=p=\operatorname {E} [X]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/22df95375a51f3068718cc64b926e4c8a517dabd)
これから
![{\displaystyle \operatorname {Var} [X]=\operatorname {E} [X^{2}]-\operatorname {E} [X]^{2}=\operatorname {E} [X]-\operatorname { E} [X]^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/78fc67baff72478c2f83e3ce1630b2b9769fc0cc)
[3]
この結果から、任意のベルヌーイ分布について、その分散が の範囲内の値を持つことが簡単に証明できます。
![{\displaystyle [0,1/4]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/604604d2122dc1c25141a841483b889d6832f261)
歪度
歪度はです。標準化されたベルヌーイ分布のランダム変数を取ると、このランダム変数は確率 で到達し、確率 で到達することが分かります。したがって、

![{\displaystyle {\frac {X-\operatorname {E} [X]}{\sqrt {\operatorname {Var} [X]}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3bdd3a134128b0517590174937df60485f0828d0)




![{\displaystyle {\begin{aligned}\gamma _{1}&=\operatorname {E} \left[\left({\frac {X-\operatorname {E} [X]}{\sqrt {\operatorname {Var} [X]}}}\right)^{3}\right]\\&=p\cdot \left({\frac {q}{\sqrt {pq}}}\right)^{3}+q\cdot \left(-{\frac {p}{\sqrt {pq}}}\right)^{3}\\&={\frac {1}{{\sqrt {pq}}^{3}}}\left(pq^{3}-qp^{3}\right)\\&={\frac {pq}{{\sqrt {pq}}^{3}}}(q^{2}-p^{2})\\&={\frac {(1-p)^{2}-p^{2}}{\sqrt {pq}}}\\&={\frac {1-2p}{\sqrt {pq}}}={\frac {qp}{\sqrt {pq}}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3ad7bf226e1980a4b3ed35773624b3e968d61933)
高次のモーメントとキュムラント
およびという事実により、生の瞬間はすべて等しくなります。


![{\displaystyle \operatorname {E} [X^{k}]=\Pr(X=1)\cdot 1^{k}+\Pr(X=0)\cdot 0^{k}=p\cdot 1+q\cdot 0=p=\operatorname {E} [X].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5cbafd464e73d482dc6c32d1c4f3eaedd5539952)
秩序の中心モーメントは次のように表される。


最初の6つの中心的な瞬間は

高次の中心モーメントは次のように簡潔に表現できる。

最初の6つのキュムラントは

エントロピ
エントロピーは、確率分布における不確実性またはランダム性の尺度です。成功確率と失敗確率を持つベルヌーイランダム変数の場合、エントロピーは次のように定義されます。




![{\displaystyle {\begin{aligned}H(X)&=\mathbb {E} _{p}\ln({\frac {1}{P(X)}})=-[P(X=0)\ln P(X=0)+P(X=1)\ln P(X=1)]\\H(X)&=-(q\ln q+p\ln p),\quad q=P(X=0),p=P(X=1)\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5639ed087a5db66ec3f0d967cb414d1a7caaa32d)
エントロピーは のときに最大になり、両方の結果が同等の可能性があるときに不確実性のレベルが最も高くなることを示します。またはのとき、エントロピーは 0 になり、一方の結果が確実になります。



フィッシャー情報量は、観測可能なランダム変数が、その確率が依存する未知のパラメータについて持つ情報量を測定するものです。ベルヌーイ分布の場合、パラメータに関するフィッシャー情報量は次のように表されます。





証拠:
- ベルヌーイ確率変数の尤度関数は次のとおりです。


これは、パラメータが与えられた場合に観測される確率を表します。



- スコア関数(対数尤度の一次導関数は次のようになります。



- フィッシャー情報は、対数尤度の 2 次導関数の負の期待値として計算されます。
![{\displaystyle {\begin{aligned}I(p)=-E\left[{\frac {\partial ^{2}}{\partial p^{2}}}\ln L(p;X)\right]=-\left(-{\frac {p}{p^{2}}}-{\frac {1-p}{(1-p)^{2}}}\right)={\frac {1}{p(1-p)}}={\frac {1}{pq}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/930c79a1f00b168593deb05fe32af5d9dcaed143)
これは、最大の不確実性と、したがってパラメータに関する最大の情報を反映しているときに最大になります。


- が独立した同一分布(iid)のランダム変数であり、すべてのベルヌーイ試行が成功確率 pである場合、それらの合計はパラメータnとp を持つ二項分布に従って分布します。

(二項分布)[3]
- ベルヌーイ分布は単純に とも表記される。


- カテゴリ分布は、任意の定数個の離散値を持つ変数に対するベルヌーイ分布の一般化です。
- ベータ分布はベルヌーイ分布の共役事前分布である。 [5]
- 幾何分布は、1 つの成功を得るために必要な独立した同一のベルヌーイ試行の数をモデル化します。
- の場合、 はラーデマッハ分布に従います。


参照
参考文献
- ^ ウスペンスキー、ジェームズ・ビクター(1937年)。数学的確率入門。ニューヨーク:マグロウヒル。p.45。OCLC 996937。
- ^ デッキング、フレデリック;クライカンプ、コーネリス。ロプハー、ヘンドリック。ルドルフ・メースター(2010年10月9日)。確率と統計の現代入門(第 1 版)。スプリンガーロンドン。 43~48ページ。ISBN 9781849969529。
- ^ abcd ベルツェカス、ディミトリ P. (2002)。確率の入門。チチクリス、ジョン N.、Τσιτσικλής、Γιάννης Ν。マサチューセッツ州ベルモント:Athena Scientific。ISBN 188652940X. OCLC 51441829.
- ^ McCullagh, Peter ; Nelder, John (1989).一般化線型モデル、第 2 版。ボカラトン: Chapman and Hall/CRC。セクション 4.2.2。ISBN 0-412-31760-5。
- ^ Orloff, Jeremy; Bloom, Jonathan. 「共役事前分布: ベータと正規分布」(PDF) . math.mit.edu . 2023年10月20日閲覧。
さらに読む
- Johnson, NL; Kotz, S.; Kemp, A. (1993).単変量離散分布(第 2 版). Wiley. ISBN 0-471-54897-9。
- ピートマン、ジョン・G. (1963)。応用統計学入門。ニューヨーク:ハーパー&ロウ。pp.162-171。
外部リンク
ウィキメディア・コモンズには、ベルヌーイ分布に関連するメディアがあります。