来自重复二项分布的样本

Question

我有一个重复的二项式检验（p = 成功概率，n = 试验次数）。我想随机生成成功次数。

我可以只生成 n 个随机数，但这太低效了。即不是这样：

def successes(p, n):
    count = 0
    for i in xrange(n):
        if random.random() < p:
            count += 1
    return count

我希望这是一个精确的样本（不是正态或其他近似值），因为这将被重复用于高度偏斜的分布。

我觉得这可以使用密度函数等 (http://en.wikipedia.org/wiki/Binomial_distribution#Probability_mass_function) 来完成，但我无法完全理解它。

Answer 1

numpy.random.binomial 做到了。