R言語でモンティ・ホール問題 - TadaoYamaokaの開発日記

モンティ・ホール問題をR言語でモンテカルロ法で解いてみた。

R言語ははじめて使ったので、もっとエレガントな書き方があるかもしれない。

choice <- function(a) {
  a[runif(1,1,length(a)+1)]
}
game <- function(change) {
  doors <- c("A", "B", "C")
  car <- choice(doors)
  you <- choice(doors)
  goat <- choice(doors[doors != car & doors != you])
  if (change) you <- choice(doors[doors != you & doors != goat])
  return(car == you)
}

change = TRUE # or FALSE
n = 10000
correct_count = 0
for(i in 1:n) {
  result <- game(change)
  correct_count <- correct_count + result
}

correct_count / n