(a, b, c, d)
| 686 | # Hypergeometric distribution. |
| 687 | # (a+b)!(c+d)!(a+c)!(b+d)! / a!b!c!d!n! for n=a+b+c+d |
| 688 | def p(a, b, c, d): |
| 689 | return C(a + b, a) * C(c + d, c) / C(a + b + c + d, a + c) |
| 690 | # Binomial coefficient. |
| 691 | # n! / k!(n-k)! for 0 <= k <= n |
| 692 | def C(n, k): |
no test coverage detected