最大似然vs最大后验

本文发布于 https://tswc.github.io/tswc.github.io/2018/01/10/statistic/
由于公式使用latex语法,在简书中好像不支持,所以这里作为草稿箱了,请移步上面的地址。

概率和统计

概率和统计通常被我们放在一起,虽然有联系,但它们的实际含义其实截然相反。
概率:我们不知道具体的数据,但是我有数据的模型。
比如说,我告诉你,a是一个整数,在$ (0,10] $平均分布。那么现在问你$a=3$的概率是多少?
$$P(a=3)=1/10$$
这里就是一个概率。 已知整体的分布,问你个体的概率。
统计:我们知道数据的出现概率,但是没有数据模型的知识。
现在,我告诉你,$a\in(0,10] & a is N$,且a是0到9的概率均为$1/10$,问你a是怎么分布的。
a在$(0,10]$上平均分布。
这里就是一个统计。我们知道个体出现的概率,问你整体的分布或者说模型是怎么样的。
总的来说就是,
概率是已知模型,用模型(方程)来计算概率。
统计是已知概率,用概率来推算模型(方程)。
对于机器学习领域来说,我们一般的流程是收集数据,希望通过已知的数据来学习规律,也就是推算出模型,用模型来处理未知的数据。
本文所侧重的极大似然和最大后验都属于统计的范畴。(用最大概率的方式来推算出模型的参数)

贝叶斯

先看看贝叶斯公式
$$P(A|B) = \frac{P(A,B)}{P(B)} \tag{1.1}$$
$$P(A|B) = \frac{P(B|A)P(A)}{P(B)} \tag{1.2}$$
<font color=red><center>贝叶斯公式就是在描述,你有多大把握能相信一件证据?(how much you can trust the evidence)</center></font>

这句话说的有点直击灵魂,棒!

例1【汽车被砸与警报器】

引用从http://blog.csdn.net/u011508640/article/details/72815981里看到的解释这份公式例子,十分通俗易懂。

$P(A|B)$表示B事件发生的情况下,A事件发生的概率,也就是说如果B发生了,我们有多大把握A事件也发生。
想想这个情况。一辆汽车(或者电瓶车)的警报响了,你通常是什么反应?有小偷?撞车了? 不。。 你通常什么反应都没有。因为汽车警报响一响实在是太正常了!每天都要发生好多次。本来,汽车警报设置的功能是,出现了异常情况,需要人关注。然而,由于虚警实在是太多,人们渐渐不相信警报的功能了。

$$P(A|B) = \frac{P(B|A)P(A)}{P(B|A)P(A) + P(B|\sim A)P(\sim A)}\tag{1.3}$$
我们假设响警报的目的就是想说汽车被砸了。把A计作“汽车被砸了”,B计作“警报响了”,带进贝叶斯公式里看。我们想求等式左边发生A|B的概率,这是在说警报响了,汽车也确实被砸了。汽车被砸引起(trigger)警报响,即B|A。但是,也有可能是汽车被小孩子皮球踢了一下、被行人碰了一下等其他原因(统统计作∼A),其他原因引起汽车警报响了,即B|∼A。那么,现在突然听见警报响了,这时汽车已经被砸了的概率是多少呢(这即是说,警报响这个证据有了,多大把握能相信它确实是在报警说汽车被砸了)?想一想,应当这样来计算。用警报响起、汽车也被砸了这事件的数量,除以响警报事件的数量(这是1.2)。进一步展开,即警报响起、汽车也被砸了的事件的数量,除以警报响起、汽车被砸了的事件数量加上警报响起、汽车没被砸的事件数量(这是1.3)。

再想想1.3,如果想让$P(A|B)=1$,也就是要将$P(B|\sim A)P(\sim A)$清0,也就是说要么杜绝小孩踢皮球,行人碰到车的可能性,要么就要让警报器在这些例外情况下不工作。这样的话,警报器的可行度大大提升。

做事之前,要综合考虑所以的可能性

当考虑$P(A|B)$时,即使$P(B|A)$很大,若是$P(A)$很小,$P(A|B)$还是很难变得较大,因为虽然A与B联系很紧,但是A的先验概率过小。

从这个角度思考贝叶斯公式:一个本来就难以发生的事情,就算出现某个证据和他强烈相关,也要谨慎。
证据很可能来自别的虽然不是很相关,但发生概率较高的事情。 发现刚才写的代码编译报错,可是我今天状态特别好,这语言我也很熟悉,犯错的概率很低。因此觉得是编译器出错了。 ————别,还是先再检查下自己的代码吧。

似然&后验

似然: $P(X;\theta)$
对于似然函数求解的问题,一般来说模型确定,但是参数未知,通过实验得到一组数据要我们来估计模型参数。其中X就是实验得到的数据,$\theta$就是我们要估计的参数。
$P(X|\theta)$表示的是在已知模型中,当参数为$\theta$时,X发生的概率。很自然的,我们要使这件事发生的概率最大,这就是我们求$\theta$的标准。
$$\theta_{MLE} = \mathop{\rm arg,max}\limits_{\theta} P(X;\theta)$$
对于这个式子中的分号,到底和表示条件概率的竖线有什么区别,这里有解释
http://blog.csdn.net/pipisorry/article/details/42715245
这里想想$P(X|\theta)$和$P(X;\theta)$的区别在哪?
(关于这个问题想了很久,下面是纯粹个人想法,不保证正确)
出现"|"时,其实表示的条件概率,等价于$P(X,\theta)/P(\theta)$;出现";"时,其实表示的是X的概率,只不过告诉我们有一个参数是$\theta$,等价于$P(X)$。
分歧点就在$\theta$上,$\theta$到底是不是一个随机变量?
如果我们将$P(X|\theta)$理解为条件概率,不考虑上面链接中说的等价情况。$P(X;\theta)$认为,$\theta$未知但确定,不存在概率的问题,只是具体的值需要听过上面的式子求得。或者说,我们对$\theta$不存在任何先验知识,所以不考虑$\theta$的分布情况,任何的$\theta$值对我们来说都是等价的。

后验:$P(X|\theta)P(\theta)$
几乎所以关于后验概率的解释都是
$$
P(\theta \vert X) = \frac{P(X \vert \theta) P(\theta)}{P(X)} \
\propto P(X \vert \theta) P(\theta)
$$
而我在这里的对于后验概率得到想法与这个过程恰恰相反。
所谓后验的意义,其实是我们有了先验的知识,所以才有了后验的概念。这个问题中,先验的知识其实就是$P(\theta)$。那么我们在极大似然中所用的$P(X;\theta)$应该写成$P(X|\theta)$,向我们前文中提到的,“一个本来就难以发生的事情,就算出现某个证据和他强烈相关,也要谨慎。”
我们希望最大化X的发生概率,不能仅仅最大化$P(X|\theta)$,因为即使X与$\theta$强烈相关,若是$P(\theta)$自身发生的概率太小,X还是不容易发生。
所以要将$P(\theta)$也考虑进去,才能保证X的发生概率。
$$
P(X \vert \theta) P(\theta) = \frac{P(X \vert \theta) P(\theta)}{P(X)}
\propto P(\theta \vert X)
$$
考虑到编程的例子,我这里很有可能还是错的,欢迎指正!

想法不一定正确,但是结论还是可以保证的
对于最大似然和最大后验的理解:两者的分歧点在于对于参数$\theta$的看法,将$\theta$视为未知定值时,是最大似然;将$\theta$视为随机变量时,考虑$P(\theta)$时是最大后验。

例2【掷硬币】

http://blog.csdn.net/u011508640/article/details/72815981
还是nebulaf91 博客里,有个掷硬币的例子来解释最大似然和最大后验。

参考:
关于具体例子:http://blog.csdn.net/pipisorry/article/details/42715245
关于$P(X|\theta)$和$P(X;\theta)$的理解:http://blog.csdn.net/u011508640/article/details/72815981

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 157,012评论 4 359
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 66,589评论 1 290
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 106,819评论 0 237
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 43,652评论 0 202
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 51,954评论 3 285
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,381评论 1 210
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 31,687评论 2 310
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,404评论 0 194
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,082评论 1 238
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,355评论 2 241
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 31,880评论 1 255
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,249评论 2 250
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 32,864评论 3 232
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,007评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,760评论 0 192
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 35,394评论 2 269
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,281评论 2 259