20/09/02 05:49:23.58 62Rt+r6G.net
>>378
今、手元にPCないのでスマホで100回シミュレーションして平均値を出してみた。
sim <- function(n=900,p=0.02,m=100){
i = 0 # trial number
j = 0 # failure
k = 0# success
while(i<=n){
a=ifelse(j==(m-1),1,rbinom(1,1,p)) # ceiling or not
j=ifelse(a==1,0,j+1) # success or failure
k=k+a # total success
i=i+1
}
return(k)
}
mean(replicate(1e2,sim()))
実行結果
> mean(replicate(1e2,sim()))
[1] 21.26
>