超幾何分布
確率質量関数 Hypergeometric PDF plot | |
累積分布関数 Hypergeometric CDF plot | |
母数 | [math]\begin{align}N&\in \left\{0,1,2,\dots\right\} \\ K&\in \left\{0,1,2,\dots,N\right\} \\ n&\in \left\{0,1,2,\dots,N\right\}\end{align}\,[/math] |
---|---|
台 | [math]k\, \in\, \left\{\max{(0,\, n+K-N)},\, \dots,\, \min{(n,\, K )}\right\}\,[/math] |
確率質量関数 | [math]{{{K \choose k} {{N-K} \choose {n-k}}}\over {N \choose n}}[/math] |
累積分布関数 | [math]1-{{{n \choose {k+1}}{{N-n} \choose {K-k-1}}}\over {N \choose K}} \,_3F_2\!\!\left[\begin{array}{c}1,\ k+1-K,\ k+1-n \\ k+2,\ N+k+2-K-n\end{array};1\right],[/math] [math]\,_pF_q[/math] は一般超幾何関数 |
期待値 | [math]n {K\over N}[/math] |
最頻値 | [math]\left \lfloor \frac{(n+1)(K+1)}{N+2} \right \rfloor[/math] |
分散 | [math]n{K\over N}{(N-K)\over N}{N-n\over N-1}[/math] |
歪度 | [math]\frac{(N-2K)(N-1)^\frac{1}{2}(N-2n)}{[nK(N-K)(N-n)]^\frac{1}{2}(N-2)}[/math] |
尖度 |
[math] \left.\frac{1}{n K(N-K)(N-n)(N-2)(N-3)}\cdot\right.[/math] [math]\Big[(N-1)N^{2}\Big(N(N+1)-6K(N-K)-6n(N-n)\Big)+[/math] [math]6 n K (N-K)(N-n)(5N-6)\Big][/math] |
モーメント母関数 | [math]\frac{{N-K \choose n} \scriptstyle{\,_2F_1(-n, -K; N - K - n + 1; e^{t}) } } {{N \choose n}} \,\![/math] |
特性関数 | [math]\frac{{N-K \choose n} \scriptstyle{\,_2F_1(-n, -K; N - K - n + 1; e^{it}) }} {{N \choose n}} [/math] |
超幾何分布(ちょうきかぶんぷ、英: hypergeometric distribution)とは、成功状態をもつ母集団から非復元抽出したときに成功状態がいくつあるかという確率を与える離散確率分布の一種である。 男女・合否などのように2種の排他的属性に分割できる有限母集団からの非復元抽出に適用される。 超幾何分布と対照的[1]な確率分布には二項分布がある。
定義
超幾何分布とは K 個の成功状態をもつ N 個の要素よりなる母集団から n 個の要素を非復元抽出したときに k 個の成功状態が含まれている確率を与える離散確率分布の一種である。 超幾何分布に従う確率変数 X の確率質量関数 fX は次で与えられる。
- [math] \operatorname{P}(X=k) = f_X(k;N,K,n) = \frac{\binom{K}{k}\binom{N-K}{n-k}}{\binom{N}{n}} = \frac{\binom{n}{k}\binom{N-n}{K-k}}{\binom{N}{K}} [/math]
確率質量関数は max(0, n + K - N) ≤ k ≤ min(K, n)のとき正となる。
超幾何分布は N が大きくなると、二項分布に近づく。また N と K の比が小さく、抽出数 n が大きいとき、ポアソン分布に近づく。
性質
- 期待値 [math]E(X) = \frac{nK}{N}[/math]
- 分散 [math]\operatorname{Var}(X) = \frac{(N-n)n(N-K)K}{(N-1)N^2}[/math]
- 最頻値 [math]\left \lfloor \frac{(n+1)(K+1)}{N+2} \right \rfloor[/math]
- 対称性 [math] \begin{align} f_X(k;N,K,n) &= f_X(k;N,n,K) \\ &= f_X(n-k;N,N-K,n) \\ &= f_X(K-k;N,K,N-n) \end{align} [/math]
例
例えば、赤い玉10個と白い玉20個を混ぜた、計30個の玉を入れた壺の中から5個の球を取りだすとき、赤い玉がちょうど1つである確率は
- [math]\frac{\binom{10}{1}\binom{30-10}{5-1} }{\binom{30}{5}}=\frac{8075}{23751}\approx 0.34[/math]
赤い玉の個数の期待値は
- [math]\frac{5\times 10}{30}\approx 1.67[/math]
多変量超幾何分布
定義
属性が 1 ≤ i ≤ c である要素を Ki 個含む N = K1 + … + Kc 個の要素よりなる母集団から n 個の要素を非復元抽出したとき、属性が i である要素を ki 個含んでいる確率を与える分布を多変量超幾何分布という。 超幾何分布と多変量超幾何分布の関係は、二項分布と多項分布の関係に相当する。
性質
多変量超幾何分布に従う確率変数を (X1, …, Xc) とする。
- 確率質量関数 [math]\operatorname{P}(X_1=k_1, \dots, X_c=k_c) = \frac{1}{\binom{N}{n}}\prod_{i=1}^{c} \binom{K_i}{k_i}[/math]
- 期待値 [math]E(X_i) = \frac{n K_i}{N}[/math]
- 分散 [math]\operatorname{Var}(X_i) = \frac{(N-n)n(N-K_i)K_i}{(N-1)N^2}[/math]
- 共分散 [math]\operatorname{Cov}(X_i,X_j) = -\frac{(N-n)nK_iK_j}{(N-1)N^2}[/math]
例
壺の中に黒い玉が5個、白い玉が10個、赤い玉が15個あるとする。その中から6個の玉を取り出すとき、各色2個ずつ取り出す確率は次の式で計算できる。
- [math]\frac{\binom{5}{2}\binom{10}{2}\binom{15}{2}}{\binom{30}{6}} \approx 0.0796[/math]
脚注
- ↑ 二項分布は超幾何分布の定義における「非復元抽出」を「復元抽出」に置き換えたものに相当する。
参考文献
- 蓑谷千凰彦, 統計分布ハンドブック, 朝倉書店 (2003).
- B. S. Everitt (清水良一訳), 統計科学辞典, 朝倉書店 (2002).
- M. Galassi et al. (富永大介訳), GNU Scientific Library リファレンスマニュアル ver. 1.8, p. 199 (2006).
関連項目
外部リンク
- Hypergeometric Probability Distribution Calculator (ALPHA) (超幾何分布の計算ができるウェブ・アプリケーション、英語)
- ちっぷす:超幾何分布を perl で計算 (日本語)
- Hypergeometric Probability Calculator: When Good Statistics Go Bad (ウェブ・アプリケーション、C++ および Ruby のソースコード、英語)
- Present Value Calculator Calculate the present value of future value sums.
- GSL reference manual Japanese version (GNU Scientific Library のマニュアルの超幾何分布を計算する関数のページ