Function rgsl::randist::hypergeometric::hypergeometric

pub fn hypergeometric(r: &mut Rng, n1: u32, n2: u32, t: u32) -> u32

This function returns a random integer from the hypergeometric distribution. The probability distribution for hypergeometric random variates is,

p(k) = C(n_1, k) C(n_2, t - k) / C(n_1 + n_2, t)

where C(a,b) = a!/(b!(a-b)!) and t <= n_1 + n_2. The domain of k is max(0,t-n_2), ..., min(t,n_1).

If a population contains n_1 elements of “type 1” and n_2 elements of “type 2” then the hypergeometric distribution gives the probability of obtaining k elements of “type 1” in t samples from the population without replacement.