重复二项分布的样本

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做到了。