20/09/30 22:33:45.08 fmOlkQ9c.net
シミュレーションしてみた。
> sim <- function(){
+ x=rbinom(2,1,1/2) # コインを2枚投げる
+ y=sample(x,1) # そのうち1枚を選ぶ(チラ見コイン)
+ c(y==1,sum(x)==2) # チラ見コインが表か?2枚が表か?を返す
+ }
> k=1e7
> z=replicate(k,sim()) # k回コインを投げる
> n=ncol(z[,z[1,]==TRUE]) # チラ見コインが表の試行数
> m=ncol(z[,z[1,]==TRUE & z[2,]==TRUE]) # そのうち2枚が表の試行数
> m/n
[1] 0.5000268