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)
)
