2012-02-13 3 views
0

나는 반복 된 이항 테스트를 가진다 (p = 성공 확률, n = 시도 없음). 성공 횟수를 무작위로 생성하고 싶습니다.반복 된 이항 분포의 샘플

나는 임의의 숫자를 생성 할 수는 있지만 너무 비효율적입니다. 즉이 싫어 :

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

I는 이것이 정확한 샘플 (정상이 아닌 다른 근사)으로 원하는이 매우 비대칭 분포를 반복적으로 사용하는 것 때문이다.

밀도 함수 등을 사용하여 수행 할 수 있다고 생각하지만 (http://en.wikipedia.org/wiki/Binomial_distribution#Probability_mass_function) 머리 주위를 감쌀 수 없습니다.

+0

당신은 너무 비효율적이라고 말하지만 그것이 비효율적이라는 것을 어떻게 알 수 있습니까? 테스트 했습니까? 아니면 그냥 추측입니까? –

+0

이것은 "반복 된 이항 분포"가 아니며 단지 이항 분포에 불과합니다. – ShreevatsaR

답변