二項分布

数学において、二項分布（にこうぶんぷ、英: binomial distribution）は、成功確率 p で成功か失敗のいずれかの結果となる試行（ベルヌーイ試行と呼ばれる）を独立に n 回行ったときの成功回数を確率変数Xとする離散確率分布である。 二項分布に基づく統計的有意性の検定は、二項検定と呼ばれている。 二項分布の典型例を次に示す。全住民の5%がある感染症に罹患しており、その全住民の中から無作為に500人を抽出する。ただし住民は500人よりずっと多いとする。このとき、抽出された集団の中に罹患者が30人以上いる確率はどれくらいだろうか。 500人のうちの感染症患者の分布は、大抵の場合は全住民のうちの患者の分布（真の分布）とおおよそ似通っていると考えられる。しかし、低確率ではあるが、選んだ500人の中に1人も患者が含まれないような真の分布とかけ離れた分布が得られる場合もある。直観的には、真の分布に近い分布が得られる確率は、真の分布から遠い分布が得られる確率より大きい。たとえば、500人中の患者の数が500×0.05=25人である確率は、24人や26人である確率より大きいだろうと思われる。しかし、その確率は定量的にどれほどだろうか。 これを定量的に表すことのできる分布が二項分布である。 抽出された集団の中に含まれる罹患者数を確率変数 X で表すとき、X は n = 500, p = 0.05 の二項分布に近似的に従う。ここで、罹患者が30人以上いる確率は Pr[X ≥ 30] である。 単純な定義としては、成功確率pの試行を独立にn回行い、成功回数を横軸にとってヒストグラムを作成した時のグラフの形である。このグラフの関数は下記の性質を持つ。 2つの母数p(0 ≤ p ≤ 1となる実数), n(自然数)に対して、0 以上の整数を値としてとる確率変数Xを定める。このとき、Xは試行の成功回数なので、(0 ≤ X ≤ n)である。 そして、X = xとなるような確率についての関数(確率質量関数)fX(x)について、 f X ( x ) = ( n x ) p x ( 1 − p ) n − x {\displaystyle f_{X}(x)={\binom {n}{x}}p^{x}(1-p)^{n-x}} となることが性質よりわかる。 また、fX(x)は確率であるため、 ∑ x = 0 n f X ( x ) = 1 {\displaystyle \sum _{x=0}^{n}f_{X}(x)=1} も明らかである。 上記を定義として、xをkに書き換え、 f X ( k ) = P ( X = k ) = ( n k ) p k ( 1 − p ) n − k {\displaystyle f_{X}(k)=P(X=k)={\binom {n}{k}}p^{k}(1-p)^{n-k}} で与えられるとき、Xは二項分布B(n, p)に従う、という。これはX ∼ B(n, p)と表記される。 ここで、 ( n k ) = n C k = n ! k ! ( n − k ) ! {\displaystyle {\binom {n}{k}}={}_{n}\!\mathrm {C} _{k}={\frac {n!}{k!\,(n-k)!}}} は n 個から k 個を選ぶ組合せの数、すなわち二項係数を表す。二項分布という名前は、この二項係数に由来している。 n = 1 の場合を特に、ベルヌーイ分布と呼ぶ。 上の定義式は次のように解釈することができる。1回の試行において成功する確率が p であるとき、pk は k 回成功する確率を表し、(1 − p)n−k は n − k 回失敗する確率を表している。ただし、k 回の成功は n 回の試行の中のどこかで発生したものであるから、nCk 通りの発生順序がある。これら全てを掛けると、n 回の独立な試行を行ったときの成功回数が k となる確率を求めることができる。 二項分布 B(n, p) に従う確率変数 X に対し、X の期待値 E[X] を求めると E [ X ] = ∑ k = 0 n k f X ( k ) = ∑ k = 0 n k ( n k ) p k ( 1 − p ) n − k = ∑ k = 1 n k ( n k ) p k ( 1 − p ) n − k = ∑ k = 1 n k n ! ( n − k ) ! k ! p k ( 1 − p ) n − k = ∑ k = 1 n n ! ( n − k ) ! ( k − 1 ) ! p k ( 1 − p ) n − k = ∑ k = 1 n n ( n − 1 ) ! ( n − k ) ! ( k − 1 ) ! p p k − 1 ( 1 − p ) n − k = n p ∑ k = 1 n ( n − 1 ) ! ( n − k ) ! ( k − 1 ) ! p k − 1 ( 1 − p ) n − k {\displaystyle {\begin{aligned}E[X]&=\sum _{k=0}^{n}kf_{X}(k)\\&=\sum _{k=0}^{n}k{\binom {n}{k}}p^{k}(1-p)^{n-k}\\&=\sum _{k=1}^{n}k{\binom {n}{k}}p^{k}(1-p)^{n-k}\\&=\sum _{k=1}^{n}k{\frac {n!}{(n-k)!\,k!}}p^{k}(1-p)^{n-k}\\&=\sum _{k=1}^{n}{\frac {n!}{(n-k)!\,(k-1)!}}p^{k}(1-p)^{n-k}\\&=\sum _{k=1}^{n}n{\frac {(n-1)!}{(n-k)!\,(k-1)!}}pp^{k-1}(1-p)^{n-k}\\&=np\sum _{k=1}^{n}{\frac {(n-1)!}{(n-k)!\,(k-1)!}}p^{k-1}(1-p)^{n-k}\\\end{aligned}}} ここで、k′ = k − 1とおくと、 = n p ∑ k = 1 n ( n − 1 ) ! ( n − k ) ! ( k − 1 ) ! p k − 1 ( 1 − p ) n − k = n p ∑ k ′ = 0 n − 1 ( n − 1 ) ! ( ( n − 1 ) − k ′ ) ! k ′ ! p k ′ ( 1 − p ) ( n − 1 ) − k ′ = n p ∑ k ′ = 0 n − 1 ( n − 1 k ′ ) p k ′ ( 1 − p ) ( n − 1 ) − k ′ = n p {\displaystyle {\begin{aligned}&{\hphantom {=}}np\sum _{k=1}^{n}{\frac {(n-1)!}{(n-k)!\,(k-1)!}}p^{k-1}(1-p)^{n-k}\\&=np\sum _{k'=0}^{n-1}{\frac {(n-1)!}{{\bigl (}(n-1)-k'{\bigr )}!\,k'!}}p^{k'}(1-p)^{(n-1)-k'}\\&=np\sum _{k'=0}^{n-1}{\binom {n-1}{k'}}p^{k'}(1-p)^{(n-1)-k'}\\&=np\\\end{aligned}}} これは全ての順序付けられた試行パターンについての平均値でもある。 続いて、分散 V[X] は V [ X ] = E [ X 2 ] − ( E [ X ] ) 2 {\displaystyle V[X]=E{\bigl [}X^{2}{\bigr ]}-{\bigl (}E[X]{\bigr )}^{2}} 上と同様に E [ X 2 ] = ∑ k = 0 n k 2 P ( X = k ) = ∑ k = 1 n k ( k − 1 ) ( n k ) p k ( 1 − p ) n − k + ∑ k = 1 n k ( n k ) p k ( 1 − p ) n − k = ∑ k = 1 n n ! ( n − k ) ! ( k − 2 ) ! p k ( 1 − p ) n − k + E [ X ] = ∑ k = 2 n n ! ( n − k ) ! ( k − 2 ) ! p k ( 1 − p ) n − k + n p = n ( n − 1 ) p 2 ∑ k = 2 n ( n − 2 ) ! ( n − k ) ! ( k − 2 ) ! p k − 2 ( 1 − p ) n − k + n p = n ( n − 1 ) p 2 ∑ k ′ = 0 n − 2 ( n − 2 ) ! ( n − 2 − k ′ ) ! k ′ ! p k ′ ( 1 − p ) n − 2 − k ′ + n p = n ( n − 1 ) p 2 + n p {\displaystyle {\begin{aligned}E{\bigl [}X^{2}{\bigr ]}&=\sum _{k=0}^{n}k^{2}P(X=k)\\&=\sum _{k=1}^{n}k(k-1){\binom {n}{k}}p^{k}(1-p)^{n-k}+\sum _{k=1}^{n}k{\binom {n}{k}}p^{k}(1-p)^{n-k}\\&=\sum _{k=1}^{n}{\frac {n!}{(n-k)!\,(k-2)!}}p^{k}(1-p)^{n-k}+E[X]\\&=\sum _{k=2}^{n}{\frac {n!}{(n-k)!\,(k-2)!}}p^{k}(1-p)^{n-k}+np\\&=n(n-1)p^{2}\sum _{k=2}^{n}{\frac {(n-2)!}{(n-k)!\,(k-2)!}}p^{k-2}(1-p)^{n-k}+np\\&=n(n-1)p^{2}\sum _{k'=0}^{n-2}{\frac {(n-2)!}{(n-2-k')!\,k'!}}p^{k'}(1-p)^{n-2-k'}+np\\&=n(n-1)p^{2}+np\\\end{aligned}}} 合計して、 V [ X ] = n ( n − 1 ) p 2 + n p − n 2 p 2 = n p ( 1 − p ) {\displaystyle {\begin{aligned}V[X]&=n(n-1)p^{2}+np-n^{2}p^{2}\\&=np(1-p)\\\end{aligned}}} となる。 二項分布 B(n, p) に従う確率変数 X の r 次モーメント E[Xr] は E [ X r ] = ∑ j = 0 r S ( r , j ) n ! ( n − j ) ! p j {\displaystyle E{\bigl [}X^{r}{\bigr ]}=\sum _{j=0}^{r}S(r,j){\frac {n!}{(n-j)!}}p^{j}} というやや複雑な表示をもつ。ここで S(r, j) は第二種スターリング数。低次から E [ X 1 ] = n p , E [ X 2 ] = n p + n ( n − 1 ) p 2 , … {\displaystyle E{\bigl [}X^{1}{\bigr ]}=np,\quad E[X^{2}]=np+n(n-1)p^{2},\dotsc } となる。一方 X の r 次階乗モーメント（英語版） E[(X)r] は E [ ( X ) r ] = ( n ) r p r = n ! ( n − r ) ! p r {\displaystyle E{\bigl [}(X)_{r}{\bigr ]}=(n)_{r}p^{r}={\frac {n!}{(n-r)!}}p^{r}} という単純な表示をもつ。ここで (n)r = n!/(n − r)! はポッホハマー記号。低次から E [ ( X ) 1 ] = n p , E [ ( X ) 2 ] = n ( n − 1 ) p 2 , … {\displaystyle E{\bigl [}(X)_{1}{\bigr ]}=np,\quad E{\bigl [}(X)_{2}{\bigr ]}=n(n-1)p^{2},\dotsc } となる。 二項分布は再生性を有する。すなわち B(n, p) に従う確率変数 X と B(m, p) に従う確率変数 Y が互いに独立であるとき、確率変数の和 X + Y は二項分布 B(n + m, p) に従う。 二項分布の近似として、以下の小節に挙げる分布などが知られている。 近似を用いることで計算の労力を削減できるという利点がある一方、各近似にはそれを適用可能とするための条件が存在する。 そのため、それらの条件や近似を用いることで生じる誤差が許容可能な範囲内に収まっていることの確認が必要となる。 特に、二項分布の母比率の信頼区間を求める際には、用いる近似と変数の値の組み合わせにより、厳密に求められた信頼区間との間に近似誤差が生じることになるため注意が必要である。 期待値 np および分散 np(1 − p) が 5 よりも大きい場合、二項分布 B(n, p) に対する良好な近似として正規分布がある。ただし、この近似を適用するにあたっては、変数のスケールに注意し、連続な分布への適切な処理がなされる必要がある。より厳密に述べれば、n が十分大きくかつ、期待値 np および 分散 np(1 − p) も十分大きい場合、期待値 np, 分散 np(1 − p) の正規分布 N(np, np(1 − p)) で近似することができ、期待値からの差 |k − np| が標準偏差 n p ( 1 − p ) {\textstyle {\sqrt {np(1-p)}}} と同程度となる k に対して P [ X = k ] ≃ 1 2 π n p ( 1 − p ) exp ⁡ ( − ( k − n p ) 2 2 n p ( 1 − p ) ) {\displaystyle P[X=k]\simeq {\frac {1}{\sqrt {2\pi np(1-p)}}}\exp {\biggl (}-{\frac {(k-np)^{2}}{2np(1-p)}}{\biggr )}} が漸近的に成り立つ。二項分布が一定の条件下で正規分布に近づく、この近似式は数学者アブラーム・ド・モアブルが1733年に著書 The Doctrine of Chances の中で紹介したのが最初であり、ド・モアブル=ラプラスの極限定理またはラプラスの定理と呼ぶことがある。これは、今日でいうところの中心極限定理の特別な場合に相当する。この正規分布への近似と標準正規分布表により、計算の労力を大きく削減することができる。 例えば、多数の住民の中から n 人を無作為に抽出し、ある質問について同意するかどうかを尋ねる場合を考える。同意する人数の割合は、もちろんサンプルに依存する。n 人を無作為に抽出する作業を何度も繰り返し行うとき、同意する人々の割合の分布は、実際の全住民の合意割合 p とほぼ等しい平均を持ち、標準偏差 σ = p ( 1 − p ) / n {\textstyle \sigma ={\sqrt {p(1-p)/n}}} である正規分布に近似される。未知の変数 p は、標準偏差が小さいほど正確な推定が可能である。そのため、抽出する人数 n は多い方が好ましい。 95%信頼区間ならば、正規分布で近似すると、その範囲は p − 1.959964 p ( 1 − p ) n ∼ p + 1.959964 p ( 1 − p ) n {\displaystyle p-1.959964{\sqrt {\frac {p(1-p)}{n}}}\sim p+1.959964{\sqrt {\frac {p(1-p)}{n}}}} となる。たとえば、p = 50% の場合、n = 100 なら40%–60%、n = 1000 ならば47%–53%、n = 10000 ならば49%–51%となる。n = 10 の場合、正規分布近似ではなく、本来の定義に従って計算すると、89%信頼区間で、30%–70%となる。 n が大きく p が十分小さい場合、np は適度な大きさとなるため、λ = np を母数とするポアソン分布が二項分布 B(n, p) の良好な近似を与える。すなわち、n が十分大きいとき、期待値 λ = np とおくと、 P [ X = k ] ≃ λ k e − λ k ! {\displaystyle P[X=k]\simeq {\frac {\lambda ^{k}e^{-\lambda }}{k!}}} が成り立つ（詳細はポアソン分布の項を参照）。この結果は数学者シメオン・ドニ・ポアソンが1837年に著書 Recherches sur la probabilite des jugements (Researches on the Probabilities) の中で与えており、ポアソンの極限定理と呼ばれる。 [脚注の使い方] ^ 藪 2012, p. 144. ^ 藪 2012, pp. 144–145. ^ Johnson, Kotz & Kemp 2005, p. 110, (3.6). ^ Johnson, Kotz & Kemp 2005, p. 109. ^ "EBCIC: Exact Binomial Confidence Interval Calculator" https://kazkobara.github.io/ebcic/README-jp.html ^ 伏見康治「確率論及統計論」第IV章　独立偶然量の和　27節 Bernoulliの定理, Laplaceの定理 p.452 ISBN 9784874720127 http://ebsa.ism.ac.jp/ebooks/ebook/204 ^ prob 3 <= x <= 7 for x binomial with n=10 and p=0.5 — Wolfram Alpha 藪友良『入門実践する統計学』東洋経済新報社、2012年。ISBN 978-4-492-47085-5。https://books.google.co.jp/books?id=EJWJ86L2bK4C。  Johnson, N. L.; Kotz, A. W.; Kemp, S. (2005). Univariate Discrete Distributions (Third ed.). Wiley. ISBN 0-471-27246-9. MR2163227. Zbl 1092.62010  負の二項分布 多項分布 ベルヌーイ過程 壺問題 ゴルトンボード