1、什么是三门问题

三门问题（Monty Hall problem） 是一个关于概率的谜题。起源于一个电视节目 Let's Make a Deal。这个问题最开始是由Steve Selvin在1975年的一封信里提出并解决。
假设在一个游戏节目中，你可以选择三扇门：一扇门后面是一辆车，另两扇门后面是山羊。在你选择一扇门后，知道门后内容的主持人打开另一扇门，门后是山羊。然后他对你说：“你想更换选择吗吗？”改变选择对你有利吗？
------节选自维基百科

2、用R语言求解

关于三门问题的分析讲解网络上多如牛毛，也有直接应用Python模拟过程求近似解的，但应用R语言重采样的还很少，正好尝试一下。

需要用到的函数

sample()：用于随机实验取样
replicate()：重复实验
mean()：求概率
as.character():转换到字符串格式

2.1 计算坚持选择的概率

B <- 50000      #实验重复次数，越大越精确
stick <- replicate(B, {
  doors <- as.character(1:3)   #问题中的三门
  prize <- sample(c("car","goat","goat"))  
  car_door <- doors[prize == "car"]   #门后是车
  my_pick  <- sample(doors, 1)   #选择三门中的一扇
  show <- sample(doors[!doors %in% c(my_pick, car_door)],1)      #主持人展示
  stick <- my_pick     #坚持选择
  stick == car_door   # 坚持选择是车
})
mean(stick)   #求平均概率
[1] 0.33772    #结果约为三分之一

2.2 计算改变选择的概率

和上面的计算是一套逻辑，只是更改后的选择是另外两门中的其中之一。

B <- 50000
switch <- replicate(B, {
  doors <- as.character(1:3)
  prize <- sample(c("car","goat","goat"))
  car_door <- doors[prize == "car"]
  my_pick  <- sample(doors, 1)
  show <- sample(doors[!doors %in% c(my_pick, car_door)], 1)
  stick <- my_pick
  switch <- doors[!doors%in%c(my_pick, show)]
  switch == car_door
})
mean(switch)
[1] 0.66546