白玉4個、赤玉5個を準備して5個取り出す作業を1000万回やって、同じ色が2個でた回数を数えれば確率の概算がでる。

> (balls=rep(c('白','赤'),c(4,5)))
[1] "白" "白" "白" "白" "赤" "赤" "赤" "赤" "赤"

> sim <- function(){
+ picked=sample(balls,5)
+ sum(picked=="白")==2 | sum(picked=="赤")==2
+ }
> mean(replicate(1e7,sim()))
[1] 0.7936473

数を数えるに還元できた。