21/01/06 15:24:30.42 CkBfrSdp.net
>>213-214
当たる確率を増やしてシミュレーションして検算
n=10 # 引けるクジの枚数
p1=20/100 ; p2=30/100 ; p3=40/100 # 当たる確率
q1=2 ; q2=4 ; q3=1 # 当たって枚数
p0=1-p1-p2-p3 # ハズレの確率
q0=n-q1-q2-q3 # ハズレの枚数
sim <- function(){
lot=sample(c(0,1,2,3),n,rep=T,prob=c(p0,p1,p2,p3))
sum(lot==1)==q1 & sum(lot==2)==q2 & sum(lot==3)==q3
}
mean(replicate(1e8,sim()))
calc <- function(){
factorial(n)/(factorial(q0)*factorial(q1)*factorial(q2)*factorial(q3))*
p0^q0*p1^q1*p2^q2*p3^q3
}
calc()
> mean(replicate(1e8,sim()))
[1] 0.00162789
>
> calc()
[1] 0.00163296