有一枚硬币(不知道它是否公平),假如抛了三次,三次都是“花”:
能够说明它两面都是“花”吗?
按照传统的算法,抛了三次得到三次“花”,那么“花”的概率应该是:
但是抛三次实在太少了,完全有可能是运气问题。我们应该怎么办?
托马斯·贝叶斯(1702-1761),18世纪英国数学家,1742年成为英国皇家学会会员。
贝叶斯认为在实验之前,应根据不同的情况对硬币有所假设。不同的假设会得到不同的推断。
比如和滑不溜手的韦小宝玩。韦小宝可能拿出各种做过手脚的硬币,让我们猜不透,只能假设对硬币一无所知。这种假设之下,我们就只能根据实验结果来猜测。
因此,实验结果是“扔三次,三次花”,倾向于认为韦小宝有可能作弊:
大侠陈近南用的可能是公平硬币:
而憨坏的多隆,真的有可能用两面“花”来和你玩:
各种假设称为先验分布,结合刚才“扔三次,三次花”的实验数据,推断出硬币的后验分布,这就是贝叶斯推断:
这里补充一下,可能大家觉得再多抛几次硬币就可以了,何必弄什么贝叶斯推断。不过现实生活中有一些事件不是能够多“抛”几次的,比如地震、彗星撞击地球等等。这里只是借着硬币来讨论问题。
那么问题来了,“先验分布”,“后验分布”用数学怎么表示:
对于扔硬币,
根据
这个特性非常适合用来做先验分布。比如,在韦小宝面前,我们对硬币一无所知。
贝叶斯说一无所知也就是意味着任何概率都是一样的,都是有可能的,所以选用均匀分布(所谓的无信息先验,可以参看这篇文章):
正直的陈近南,可能用的是公平硬币,也就是说概率在0、1之间(0表示“字”,1表示“花”),
而憨坏的多隆,可能用了两面花,也就是说概率可能集中到1附近,
也就是说可以用
用
具体到这里:
再具体到韦小宝的情况就是:
其中,用
图像上的变化就是:
可以看到,作弊的可能性还是比较大的。
陈近南的情况:
结合实验数据之后,图像的中心从0.5往0.6方向移动了,作弊可能性有所增加,不过总体来看应该还是公平硬币的可能性大。
多隆的情况:
更向1集中,作弊的可能性非常高。
贝叶斯推断:
的应用到二项式分布的数学细节如下。假设实验数据
上面的式子根据贝叶斯定理(离散贝叶斯可以参看怎样用非数学语言讲解贝叶斯定理(Bayes theorem)?,连续贝叶斯可以参看这里)可以表示为:
其中
因此,写成下面这样更容易看清楚重点(其中
其中,
随着
对于二项式分布,用
这种特性称为共轭先验。
并且:
关于这点的证明请参看这里,需要科学上网。