贝叶斯定理为什么这么强大?这篇关于它的原理与应用的书终于讲清楚了!

0cd2cceaf9b52f0ba0b3b497cabd211d.png

相信很多人都或多或少的听过贝叶斯定理,但它到底怎么用?又怎么学?许多的书都无法讲清楚。但《趣学贝叶斯统计》却独辟蹊径,用一个个趣味十足、脑洞大开的例子,将贝叶斯统计的原理和用途娓娓道来。

在书中你可以评估UFO出现在自家后院中的可能性、《星球大战》中汉•索罗穿越小行星带幸存下来的可能性、抓鸭子中大奖游戏的公平性,并学会用乐高积木理解贝叶斯定理。

简直是最好的贝叶斯入门书!今天我们就研究如何将贝叶斯定理当作一种概率工具,对不确定性进行逻辑推理。将利用贝叶斯定理来计算和量化在给定数据的情况下,信念有多大的可能性为真。为此,需要使用该定理的三要素——后验概率、似然和先验概率。这3个要素将在这场贝叶斯统计和概率探险之旅中频繁出现。

来源 | 《趣学贝叶斯统计:橡皮鸭、乐高和星球大战中的统计学》

作者 | [美] 威尔·库尔特(Will Kurt)

译者 | 王凌云

1 贝叶斯定理三要素

贝叶斯定理可以准确地量化所观察到的数据改变我们信念的概率。这也就是。简单来说,我们想量化的是:在所观察到的数据下,自己对信念的坚信程度。在贝叶斯公式中,这个要素的术语是后验概率(posterior probability,简称为“后验”),也就是将通过贝叶斯定理所求出的解。

为了得到后验概率,还需要用到下一个要素:似然(likelihood)。它表示在给定信念的情况下,观察到某一数据的概率,也就是(数据|信念)。

最后,需要量化初始信念的概率,即 (信念)。这一要素在贝叶斯定理中被称为先验概率(prior probability,简称为“先验”),它表示我们在看到数据之前的信念强度。似然和先验结合在一起就会形成后验。通常情况下,我们需要使用数据的概率 (数据)对后验归一化,从而使其值介于0和1之间。然而在实践中并不总是需要 (数据),所以这个值没有特殊的名字。

你已经知道,我们将信念称为假设 ,并用变量  来表示数据。图8-1展示了贝叶斯定理的各个要素。

cec63023147ff10b0b99c5e42c5ebc7a.jpeg

图8-1 贝叶斯定理的要素

在下文中,我们将调查一起犯罪案件,并结合这些要素进行推理。

2 调查犯罪现场

假设,一天你下班回家后,发现家里的窗户玻璃碎了,前门开着,你的笔记本计算机也不见了。你的第一反应可能是:“家里被盗了!”但你是如何得出这个结论的?更重要的是,你如何量化这个信念呢?

你的第一反应是家里被盗了,所以这里 e71af6dc3f1e3b6e38bbc25fa650e2a2.png 被盗。我们需要一个概率来描述家里被盗的可能性有多大,所以根据现有的数据,想要求解的后验是:

ddcd5faf80923f8e9ac88e762dcab5ed.png(被盗 | 窗户玻璃碎了, 前门开着, 笔记本计算机不见了)

为了解决这个问题,我们需要补充贝叶斯定理中缺失的部分。

1 求解似然

首先,需要求解似然,具体到这个例子也就是,如果家里真的被盗了,同样的数据会被观察到的概率。换句话说,也就是数据与假设的吻合程度:

2809b5d14c53a243962be2dc49ad931d.png(窗户玻璃碎了, 前门开着, 笔记本计算机不见了 | 被盗)

这里,我们问的是:“如果发生了盗窃,那么你看到当前这些数据的概率有多大?”你可以想象一下,被盗时以上所有数据并非都存在的任何场景。例如,聪明的小偷可能撬开了你家的锁,偷走笔记本计算机之后再把门锁上,这不需要打破窗户玻璃。或者他可能只打破了窗户玻璃,拿走笔记本计算机之后再通过窗户爬出来。从直觉上来说,我们看到的场景在盗窃现场很常见,所以我们假定,如果家里被盗,你回家后有 7df251fe9c6fba6009a47a29d04b8913.png 的概率会发现这些数据。

值得注意的是,尽管在这个例子中,我们只是猜测有哪些数据,但其实也可以通过一些调查来获得更好的判断。比如,我们可以去当地警察局询问盗窃案件犯罪现场的统计数据,或者浏览最近关于盗窃案的新闻报道。这样就可以得到更准确的似然估计:如果被盗,你看到这些数据的概率。

贝叶斯定理的不可思议之处在于,我们既可以用它来衡量信念,也可以用它来处理具有精确概率的大数据集。即使认为 18054e8764f059894533dc6f2f6a5393.png 不是一个好的估值,你也可以随时返回去重新计算,像我们将要做的那样,看看在不同的假设下这个值是如何变化的。如果认为发生盗窃时看到这些数据的概率只有 9479b14dd1434dac0e96b75744d73a07.png,你可以将这个值重新代入进行计算。贝叶斯统计让人们以一种可度量的方式产生不同的信念。因为是以量化的方式处理信念的,所以你可以重做本章所做的一切计算,看看不同的概率是否会对最终的结果产生实质性影响。

2 计算先验概率

接下来,我们需要确定家里被盗的概率。这也是本例的先验概率。先验概率非常重要,因为它允许我们使用背景信息对似然进行调整。假设前面描述的场景发生在一个荒岛上,而你是岛上唯一的居民,那么你家几乎不可能被盗(至少是被人类)。换一个场景,如果你家位于犯罪率很高的街区,那么盗窃事件就可能会经常发生。为简单起见,我们将被盗的先验概率设定为:

dacf27f4c3d5948e3750eca1ae437c92.png被盗9cef7f36b46db6a5ab76ebf35394a80f.png

请记住,如果有不同的或额外的数据,随时可以调整这个概率。

现在,我们几乎有了计算后验概率的所有条件,只差对数据进行归一化处理。在继续之前,先来看看未归一化的后验概率:

cf215e1d1b49c094f54378ed9bea3fdc.png被盗23f64ea052b054b172489d07a98a203b.png窗户玻璃碎了, 前门开着, 笔记本计算机不见了 | 被盗677a5ef7c87710e9f11db4d048d9a541.png

这个值非常小,太令人惊讶了!这是因为直觉告诉我们,根据观察到的数据,家里被盗的概率看起来非常大。这里,我们还没有分析观察到这些数据的概率。

3 归一化数据

我们的公式中还缺少所观察到的这些数据发生的概率 610b26d29c5db7072f59a9671a325d0a.png,无论家里是否被盗。在这个例子中,这是指无论出于什么原因,同时观察到“家里窗户玻璃碎了,前门开着,并且笔记本计算机不见了”的概率。现在的公式是这样的:

57522f58c8a1cd75f990acbf202f7eda.png被盗 | 窗户玻璃碎了, 前门开着, 笔记本计算机不见了43d17d953317d0cd70af2af6945bddc5.png

分子中的概率相当小,因为没有对它进行归一化处理。

我们可以在表8-1中看到当 0b1e82be002a5a202c0aa1ce9a252f7d.png 改变时,后验概率的变化情况。

表8-1 b661e1d78813b73a788c0c3ed8d78293.png 对后验概率的影响

0d6529ac078dc5f88bbc6fb3cd7cf2c1.png

后验概率

0.050

0.006

0.010

0.030

0.005

0.060

0.001

0.300

可以发现,当 f3c804f5c87eb841959d3a51db3d924c.png 减小时,后验概率就会增大。这是因为,随着观察到这些数据的概率越来越小,整个事件发生的概率在增大(见图8-2)。

dc64b2cc64f54142e8a7f3e7321d5e9d.jpeg

图8-2 随着数据发生概率的减小,后验概率会增大

思考下面这个极端的例子:你朋友成为百万富翁的“唯二”途径是中彩票或者从某个连他自己都不知道的家族成员那里继承遗产。因此,你朋友成为百万富翁的概率非常小。然而,你发现这位朋友确实成了百万富翁。那么,他中彩票的概率就变大了,因为这是他成为百万富翁仅有的两种方式之一。

当然,被盗只是你所见数据的一种可能解释,还有其他更多可能的解释。然而,如果不知道所见数据的概率,我们就无法将所有其他概率归一化。那么 00aed0f5efab96d778080dc472eb31a9.png 是多少呢?这才是最棘手的问题。

与 55c33e3ad0374a1cec30ea4f7ed236ec.png 相关的最常见的问题是,在很多现实情况下,它很难精确计算。对于公式中的其他值(虽然对这个例子来说只是猜测了一个值),我们都可以通过收集真实的数据来提供更准确的概率。对先验概率 344a16eaba1743b8217161331e5a41e4.png(被盗),似乎只能通过查看历史犯罪数据来确定你家所在的街道上,某一特定人家在某一天被盗的概率。同样,理论上我们可以调查过去的盗窃案件,得到针对所观察到的数据的一个更准确的似然。但又怎么得到 de50173ef4c266e4815b2feece2891e8.png(窗户玻璃碎了, 前门开着, 笔记本计算机不见了)的实际值呢?

与其研究所看到的数据发生的概率,不如试着计算所有其他能够解释所看到数据的可能事件的概率。由于它们的和必须等于1,因此我们可以倒过来计算 36e8f060a4ca5f8b657b6a4482e177b7.png。只是对这份特殊的数据来说,几乎有无限的可能性。

没有 d9294804680c1aa9ed47de32f20c7bdb.png 的值,我们似乎卡住了。在第6章和第7章中,我们分别计算了客服代表是男性的概率和选择不同颜色乐高凸粒的概率,当时有大量关于 c5aa24389dbd5599e431e2a371df31c7.png 的信息。这使我们可以根据观察到的情况,得出我们相信自己假设的准确概率。如果没有 038d4af609a4d6d9860cbbd1aae69d43.png,就无法求出 84e7d219e0b4dab1e313da46cfc56cf9.png(被盗 | 窗户玻璃碎了, 前门开着, 笔记本计算机不见了)的值。然而,我们并不是完全没办法。

好消息是,在某些情况下,并不需要明确知道 608c6aa8425d9b56aabbf84666c0c15e.png 的值,因为我们通常只是想对假设进行比较。具体到这个例子,我们将用另一种可能的解释与家里被盗的概率进行比较。我们可以通过观察非归一化后验分布的比值做到这一点。因为 1e3efaca2b6cf3dcac6e6f23be7632de.png 是一个常数,所以去掉它,分析结果也不会改变。

因此,在本章剩余的内容里,我们不再计算 6bb04d5ac3e5fa563fc82a0f0be64a2a.png,而是提出一个备择假设,计算它的后验概率,然后将其与原始假设的后验概率进行比较。虽然这意味着我们不能给出被盗(作为所观察到的数据的唯一可能解释)的确切概率,但我们仍然可以利用贝叶斯定理来进行推理,并分析其他的可能性。

3 考虑备择假设

现在提出另外一个假设,并将它与原来的假设进行比较。新假设包括以下3个事件。

  (1) 邻居家孩子把棒球打到了窗户上。

  (2) 你离开家时忘了锁门。

  (3) 你忘了自己带笔记本计算机去上班并把它落在了办公室。

我们用事件前面的编号来指代这些事件,并将它们统称为 c235f904466f9a5e239a29ba9a7a0a37.png,所以 5101bc14de482c136ebf55d3e32313ac.png。现在求解这些数据的似然和先验概率。

1 备择假设的似然

对似然,我们想计算的是在给定假设下所观察到的事件的概率,或者说是 079a300410425da202559534d28e4821.png。有趣的是,这个假设的似然是1:e848edecaad23b5b3bdf5f6207e6fb43.png

如果假设中的所有事件都发生了,那么你肯定会观察到窗户玻璃碎了、前门开着以及笔记本计算机不见了。

2 备择假设的先验概率

先验概率表示的是这3个事件都发生的可能性,这也意味着需要先计算出其中每个事件的概率,然后通过乘法法则来确定先验概率。在这个例子中,我们假设每个可能的结果都是条件独立的。

备择假设的第一项内容是,邻居家孩子打棒球时不小心打碎了窗户玻璃。虽然这在电影中很常见,但现实中我从未听说过这种情况,更多的情况是发生了盗窃,所以我们假设棒球打碎窗户玻璃的概率是被盗概率的一半:

4250c46f640fa46b9f0541521356cb4d.png

备择假设的第二项内容是你忘了锁门。这种情况相当普遍,所以假设它每月发生一次:

baf33f1e37bbf98e77f8a8bbdbdb2d37.png

最后,让我们来看看将笔记本计算机落在办公室的概率。虽然带着笔记本计算机去上班并将它落在办公室可能很常见,但完全忘记带着它去上班的情况不太常见。假设这种情况每年会发生一次:

b74149cf31a0403e747e00b36c73027a.png

既然已经给假设 c1d3f668cd375f65e79b8d8f7b34a000.png 中的每一个事件都赋予了概率,那么可以用乘法法则来计算先验概率了:

26438ac178443d86f088e9100e7131c6.png

正如你看到的,这3个事件同时发生的先验概率很小。现在我们需要计算出这两个假设的后验概率以进行比较。

3 备择假设的后验概率

我们知道似然 9a5101fafd8f44c95be9cb789544f6e5.png 等于1,所以如果第二个假设是真的,那么我们就一定会看到这些数据。如果没有先验概率,看起来这个新假设的后验概率要比原假设(家里被盗了)恰当得多,因为即使被盗了,我们也不太可能看到这些数据。现在我们可以看到,先验概率是如何从根本上改变非归一化的后验概率的:

2fb368144d0561896f9f40c5d6dbc42d.png

现在我们想用一个比值来比较后验信念以及假设的强度。你会发现,做这件事并不需要 a744ff6f8b7ca7fcb24ee8decefa8ff8.png

4 比较非归一化的后验概率

首先,我们需要求出两个后验概率的比值。比值能够告诉我们一个假设的可能性是另一个假设的多少倍。将原假设定义为 b9f8996186153ae89eefdf91c09bf9c1.png,这两个假设为真的概率之比如下:

5328b2e38f075a05d4e2dce4f27c04d4.png

接下来,用贝叶斯定理将其中的每一项都展开。这里将贝叶斯定理写为 1238742e2dc4a30fb03440beb0115602.png,以使下面这个公式更易于阅读:

6f1548a6e1db848d91c98bd2ff3c387f.png

请注意,分子和分母中都有 43d180970ceb807296912b602bac15a5.png,这意味着可以直接消去它,比值保持不变。这就是在比较假设时 d088e10e7e793b82e9be7833d8d3af0e.png 并不重要的原因。现在我们得到了未归一化的后验概率的比值。后验概率表示信念的坚定程度,所以这个后验概率比值告诉我们,在不知道 e7fdfba9045309efaf18312f1e7935f1.png 的情况下,0489ae0c53fe17db7c5326d7276b0102.png 对数据的解释比 f9a055f03c45d0656b7f7be71bebc3e3.png 好多少倍。消去 be1307378a6e911e08e96698b8ca6bde.png 并将实际数值代入:

ba067b65405f02fdceae3566675b04a5.png

这意味着 e8bb03435b60234a7e0a67eb7ea78492.png 对所观察数据的解释能力是 db74704396c0a8562810da64cf4b3804.png 的6570倍。换句话说,我们的分析表明,原始假设(0ff4e160e0b4930444652cd30e57840f.png)比备择假设(0b976f5f59f0ea21e73b93a7960013b0.png)更能解释所观察的数据。这也符合我们的直觉:根据观察到的场景,盗窃看上去更可能是事情的真相。

我们想用数学方式表达非归一化后验概率的这一性质,以便于进行比较。为此,需要使用如下版本的贝叶斯定理,其中符号 e2aa1f74acc030c6dd7c44cec36bb3a2.png 的意思是“成正比”:

38fe3458a8782ec21e88831086b7e0ef.png

这个公式可以理解为“后验概率,即给定数据下假设的概率,与 6e49565ab965bcb7e909c228a28bc2cc.png 的先验概率和在假设 9b555effd1422103a511e4c1801e3b29.png 下数据概率的乘积成正比”。

当想比较两个假设的概率,但 1ffe57f2a0d5955b17f454f33b2758f9.png 的值又不容易计算时,贝叶斯定理的这种形式就非常有用。虽然这无法单独得出一个有意义的假设概率值,但我们仍然能够使用贝叶斯定理的这一形式来比较假设。对假设进行比较,就意味着我们可以确切地知道一种假设对观察内容的解释要比另外一种假设的可信度强多少。

  推荐阅读

080c24d11f88f4b0954e57b457645c1f.jpeg

《趣学贝叶斯统计:橡皮鸭、乐高和星球大战中的统计学》

作者:[美] 威尔·库尔特(Will Kurt)

译者:王凌云

本书用十余个趣味十足、脑洞大开的例子,将贝叶斯统计的原理和用途娓娓道来。你将从直觉出发,自然而然地习得数学思维。读完本书,你会发现自己开始从概率角度思考每一个问题,并能坦然面对不确定性,做出更好的决策。

02

4d853f10981d5861598f66a495a47114.png

《贝叶斯的博弈:数学、思维与人工智能》

作者:黄黎原

译者:方弦

法国数学类科普书、大学数学参考及教材类图书畅销书目,在机器学习、人工智能、逻辑学和哲学等众多领域中,探索贝叶斯定理蕴藏的智慧与哲理。

贝叶斯定理一旦与算法相结合,就不再是一套枯燥的数学理论或认识论,而变成了应用广泛的知识宝库,催生了众多现代数学定理,以及令人称道的实践成果。

85dd7bbda48ba3d75e66be8474601105.png

《概率论沉思录》

作者:埃德温·汤普森·杰恩斯

译者:廖海仁

著名数学物理学家,圣路易斯华盛顿大学和斯坦福大学教授,统计力学和概率统计推断方面权谋埃德温·汤普森·杰恩斯,40年思想著作;

无数读者苦等15年的概率论神作,英文版豆瓣评分9.4高分;

概率论作为逻辑的延伸,是所有科学推断的基础。本书收集了概率统计的各种线索,将概率和统计推断融合在一起,用新的观点生动地描述了概率论在物理学、数学、经济学、化学和生物学等领域中的广泛应用,尤其是阐述了贝叶斯理论的丰富应用,弥补了传统概率论和统计学的不足,并揭开了众多悖论背后的玄机。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值