20/09/02 06:10:57.22 62Rt+r6G.net
>>386
求めるのは16回以上当たる確率だった
１万回のシミュレーション
sim <- function(n=900,p=0.02,m=100){
i = 0 # 試行回数
j = 0 # 現在の外れの連続数
k = 0 # 当たった総数
while(i<=n){
a=ifelse(j==(m-1),1,rbinom(1,1,p)) # 天井なら１でそれ以外はpで当たる
j=ifelse(a==1,0,j+1) # 当たりか否かでjをリセット
k=k+a # 現在の当たりの総数
i=i+1 # 試行回数
}
return(k) # 当たりの数を返す
}
K=replicate(1e4,sim())
mean(K>=16) # 16回以上の割合
> mean(K>=16) # 16回以上の割合
[1] 0.9376