用贝叶斯定理解决三门问题并用Python进行模拟(Bayes' Rule Monty Hall Problem Simulation Python)

三门问题(Monty Hall problem)也称为蒙提霍尔问题或蒙提霍尔悖论,出自美国的电视游戏节目《Let’s Make a Deal》。问题名字来自该节目的主持人蒙提·霍尔(Monty Hall)。

这个游戏的玩法是:参赛者会看见三扇关闭的门,其中一扇门后面有一辆汽车,选中后面有车的那扇门就可以赢得该汽车,而另外两扇门后面则各藏有一只山羊。当参赛者选定了一扇门,但未去开启它的时候,节目主持人会开启剩下两扇门中的一扇,露出其中一只山羊。主持人其后会问参赛者要不要更换其初始的选择,选另一扇仍然关上的门。

那么问题来了,参赛者到底要不要更换其初始的选择呢?

解决这个问题需要用到贝叶斯定理:用贝叶斯定理解决三门问题并用Python进行模拟(Bayes' Rule Monty Hall Problem Simulation Python)

让我们选一个特定的例子来看看:假设三扇门分别为Door A,Door  B,Door C,并且参赛者初始选定了Door A,然后主持人展示了Door  B。那么参赛者是坚持选择Door A还是更换成Door C呢?这就要根据Door A和Door C哪个门后汽车出现的概率较大决定了。

也就是说,我们需要解决P(Door A=car|Door A is selected, Door  B is revealed)P(Door C=car|Door A is selected, Door  B is revealed)哪个大的问题。

用贝叶斯定理解决三门问题并用Python进行模拟(Bayes' Rule Monty Hall Problem Simulation Python)

用贝叶斯定理解决三门问题并用Python进行模拟(Bayes' Rule Monty Hall Problem Simulation Python)

首先,每个门后有车的概率都是1/3:用贝叶斯定理解决三门问题并用Python进行模拟(Bayes' Rule Monty Hall Problem Simulation Python)

其次,如果Door A门后有汽车,那么Door A被选择的几率是1/3,假设初始选择了Door A,那么Door B被主持人打开的几率是1/2:用贝叶斯定理解决三门问题并用Python进行模拟(Bayes' Rule Monty Hall Problem Simulation Python)

再次,普通情况下,Door A被选择的几率是1/3,Door B被主持人打开的几率是1/2(因为已经有一扇门被选择了,选择的门不能被打开):用贝叶斯定理解决三门问题并用Python进行模拟(Bayes' Rule Monty Hall Problem Simulation Python)

同理,如果Door C门后有汽车,那么Door A被选择的几率是1/3,假设初始选择了Door A,那么Door B被主持人打开的几率是1:用贝叶斯定理解决三门问题并用Python进行模拟(Bayes' Rule Monty Hall Problem Simulation Python)

因此,我们可以看到,P(Door C=car|Door A is selected, Door  B is revealed)是P(Door A=car|Door A is selected, Door  B is revealed)的两倍。也就是说,更换初始的选择将会使我们的获胜几率提高2倍!

可以用probability tree来帮助理解一下:

用贝叶斯定理解决三门问题并用Python进行模拟(Bayes' Rule Monty Hall Problem Simulation Python)

如果对上面的计算公式还有疑问,那么让我们用计算机来模拟一下:

from random import randint
from random import choice N = 1000 def simulate(N):
m=0 #设置不更换初始选择赢得汽车的次数
n=0 #设置更换初始选择赢得汽车的次数
for i in range(N): #模拟1000次游戏
win=randint(1,3) #设置藏有汽车的门,在1-3之间随机选出
bet1=randint(1,3) #设置初始选择的门,在1-3之间随机选出
remain=[i for i in range(1,4) if i!=win and i!=bet1] #剩余可选的门(除去初始选择的门和藏有汽车的门)
monty_reveal=choice(remain) #monty会在剩余可选的门中选择一扇门打开
bet2=6-bet1-monty_reveal #bet2表示更换初始选择(用6减是因为三扇门加起来等于6)
if bet1==win: #如果初始选择和藏有汽车的门吻合,那么初始选择的获胜次数+1
m+=1
if bet2==win: ##如果更换初始选择的bet2和藏有汽车的门吻合,那么bet2的获胜次数+1
n+=1
return n/m print(simulate(N))
2.0211480362537766

最后的结果: 更换初始选择获胜的次数差不多是不更换初始选择获胜次数的两倍。

三门问题是有些反直觉的,我们可以这样来理解:当参赛者选择Door A时,他的获胜概率是1/3,当主持人展示了Door B门后没有汽车以后,这个信息并没有给参赛者的初始选择带来任何有用的信息 ,选择Door A获胜的概率仍然是1/3,但是鉴于选择Door B获胜的概率降为了0,因此选择Door C获胜的概率变为1-1/3,也就是2/3。

参考:https://classroom.udacity.com/courses/st101/lessons/48744119/concepts/484806120923

上一篇:PHP生成器Generators


下一篇:SET操作符