def test_binomial(rnd, n, p): # bimonial random variate generator binomial_dist = RVGs.Binomial(n, p) # obtain samples samples = get_samples(binomial_dist, rnd) # report mean and variance print_test_results('Binomial', samples, expectation=n*p, variance=n*p*(1-p))
def test_multinomial(rnd, n, pvals): # multinomial random variate generator multinomial_dist = RVGs.Binomial(n, pvals) # obtain samples samples = get_samples(multinomial_dist, rnd) # report mean and variance print_test_results('Multinomial', samples, expectation=n*pvals, variance=n*pvals*(1-pvals) )