확률 목록으로 이항 확률을 효율적으로 계산하는 방법은 무엇입니까? R

Question

에서 나는 독립적 인 이벤트에 대한 성공의 확률의 목록을 가지고 :

probs <- c(.2,.3,.4) 

나는이 질문에 대한 대답 함수를 작성하려면 :이 독립적 인 사건에서 적어도 X의 성공을 얻기의 가능성은 무엇인가?

현재 나의 접근 방식은 for 루프를 사용하여 다음 코드를 수천 번 실행하는 것입니다. 더 효율적인 방법이 있습니까?

sim_results <- probs - runif(length(probs), 0, 1) 
length(sim_results[sim_results > 0]) 

Answer 1

한 샷으로 모든 난수를 생성하고 루프를 피하십시오.

# number of random probabilities 
p = 100 
# number of simulations 
N = 1000 
# random probabilities 
probs = runif(p) 
M = matrix(runif(p*N), N, p) 
# outcome 
y = apply(M,1,function(x) sum(probs>x)) 
hist(y)