反復二項分布のサンプル

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です。