18/06/30 09:59:55.93 j2CU1Lw0.net
>>316
シミュレーションスクリプトを書いてみた。
A=c(1,1,0)
B=c(rep(1,101),rep(0,100))
f<-function(){
idx=rbinom (1,1,0.5)
if(idx)AB=A else AB=B
AB
ball=sample (AB,2, replace=FALSE)
c(ball[1],ball[2],idx)
}
re=replicate (1e3,f())
sum(re[1,]&re[2,]&re[3,])/sum(re[1,]&re[2,]) #P(A|RR)
sum(re[1,]&!re[2,]&re[3,])/sum(re[1,]&!re[2,]) #P(A|RW)
sum(!re[1,]&re[2,]&re[3,])/sum(!re[1,]&re[2,]) #P(A|WR)
sum(!re[1,]&!re[2,]&re[3,])/sum(!re[1,]&!re[2,]) #P(A|WW)
sum(re[1,]&re[2,]&!re[3,])/sum(re[1,]&re[2,]) #P(B|RR)
sum(re[1,]&!re[2,]&!re[3,])/sum(re[1,]&!re[2,]) #P(B|RW)
sum(!re[1,]&re[2,]&!re[3,])/sum(!re[1,]&re[2,]) #P(B|WR)
sum(!re[1,]&!re[2,]&!re[3,])/sum(!re[1,]&!re[2,]) #P(B|WW)