20/09/06 09:31:01 tN16f4S+.net
>>924
１０万回のシミュレーション結果と照合
> sim <- function(p=0.0833,s=17){
+ Su=0 # 成功数
+ i=0 # 試行数
+ while(Su<17){
+ Su=Su+rbinom(1,1,p)
+ i=i+1
+ }
+ return(i)
+ }
> re=replicate(1e5,sim())
> quantile(re,c(0.025,0.975))
2.5% 97.5%
122 308
> HDInterval::hdi(re)
lower upper
115 297
attr(,"credMass")
[1] 0.95