確率の問題をシミュレーションで解く

斎藤さんには二人の子供がいる。日曜日生まれの女の子はいるかと聞くと、いると言う。
では、もう一人も女の子である確率は？

分かりやすくするために，冗長に書く

n = 100000000
d = data.frame(mf1=sample(2, n, replace=TRUE), mf2=sample(2, n, replace=TRUE),
               dw1=sample(7, n, replace=TRUE), dw2=sample(7, n, replace=TRUE))
d2 = d[(d$mf1 == 2 & d$dw1 == 1) | (d$mf2 == 2 & d$dw2 == 1),]
sum(d2$mf1 == 2 & d2$mf2 == 2)/nrow(d2)
13/27 # 理論値

mf1，mf2 は性別（2 が女），dw1，dw2 は生まれ曜日（1 が日曜日）