1. XenForo 1.5.14 中文版——支持中文搜索!现已发布!查看详情
  2. Xenforo 爱好者讨论群:215909318 XenForo专区

有哪些贝叶斯推理入门的书籍?

本帖由 漂亮的石头2022-07-10 发布。版面名称:知乎日报

  1. 漂亮的石头

    漂亮的石头 版主 管理成员

    注册:
    2012-02-10
    帖子:
    486,020
    赞:
    46
    [​IMG] 猴子,著有畅销书《数据分析思维》科普中国专家 阅读原文

    把我下面写的内容看懂就可以了,有些书上反而太复杂,一般人儿看不懂。

    贝叶斯定理太有用了,不管是在投资领域,还是机器学习,或是日常生活中几乎都在用到它。

    例如,生命科学家用贝叶斯定理研究基因是如何被控制的;教育学家意识到,学生的学习过程其实就是贝叶斯法则的运用;基金经理用贝叶斯法则找到投资策略;谷歌用贝叶斯定理改进搜索功能,帮助用户过滤垃圾邮件;无人驾驶汽车接收车顶传感器收集到的路况和交通数据,运用贝叶斯定理更新从地图上获得的信息;人工智能、机器翻译中大量用到贝叶斯定理...

    我将从以下 4 个角度来科普贝叶斯定理及其背后的思维:

    1.贝叶斯定理有什么用?

    2.什么是贝叶斯定理?

    3.贝叶斯定理的应用案例

    4.生活中的贝叶斯思维



    1.贝叶斯定理有什么用?

    英国数学家托马斯·贝叶斯(Thomas Bayes)在 1763 年发表的一篇论文中,首先提出了这个定理。而这篇论文是在他死后才由他的一位朋友发表出来的。

    (ps:贝叶斯定理其实就是下面图片中的概率公式,这里先不讲这个公式,而是重点关注它的使用价值,因为只有理解了它的应用意义,你才会更有兴趣去学习它。)

    [​IMG]

    在这篇论文中,他为了解决一个“逆概率”问题,而提出了贝叶斯定理。

    在贝叶斯写这篇文章之前,人们已经能够计算“正向概率”。什么是正向概率呢?举个例子,杜蕾斯举办了一个抽奖,抽奖桶里有 10 个球,其中 2 个白球,8 个黑球,抽到白球就算你中奖。你伸手进去随便摸出 1 颗球,摸出是中奖球的概率是多大。

    [​IMG]

    根据频率概率的计算公式,你可以轻松的知道中奖的概率=中奖球数(2 个白球)/ 球总数(2 个白球 +8 个黑球)=2/10

    如果还不懂怎么算出来的,可以看我之前写的科普概率的回答:猴子:如何理解条件概率?

    而贝叶斯在他的文章中是为了解决一个“逆概率”的问题。比如上面的例子我们并不知道抽奖桶里有什么,而是摸出一个球,通过观察这个球的颜色,来预测这个桶里里白色球和黑色球的比例。

    [​IMG]

    这个预测其实就可以用贝叶斯定理来做。贝叶斯当时的论文只是对“逆概率”这个问题的求解尝试,这哥们当时并不清楚这里面这里面包含着的深刻思想。

    然而后来,贝叶斯定理席卷了概率论,并将应用延伸到各个领域。可以说,所有需要作出概率预测的地方都可以见到贝叶斯定理的影子,特别地,贝叶斯是机器学习的核心方法之一。

    为什么贝叶斯定理在现实生活中这么有用呢?

    这是因为现实生活中的问题,大部分都是像上面的“逆概率”问题。因为生活中绝大多数决策面临的信息都是不全的,我们手中只有有限的信息。既然无法得到全面的信息,我们就只能在信息有限的情况下,尽可能做出一个好的预测。

    比如天气预报说,明天降雨的概率是 30%,这是什么意思呢?

    我们无法像计算频率概率那样,重复地把明天过上 100 次,然后计算出大约有 30 次会下雨(下雨的天数 / 总天数)

    而是只能利用有限的信息(过去天气的测量数据),用贝叶斯定理来预测出明天下雨的概率是多少。

    同样的,在现实世界中,我们每个人都需要预测。想要深入分析未来、思考是否买股票、政策给自己带来哪些机遇、提出新产品构想,或者只是计划一周的饭菜。

    贝叶斯定理就是为了解决这些问题而诞生的,它可以根据过去的数据来预测出未来事情发生概率。

    贝叶斯定理的思考方式为我们提供了有效的方法来帮助我们做决策,以便更好地预测未来的商业、金融、以及日常生活。

    总结下第 1 部分:贝叶斯定理有什么用?

    在有限的信息下,能够帮助我们预测出概率。

    所有需要作出概率预测的地方都可以见到贝叶斯定理的影子,特别地,贝叶斯是机器学习的核心方法之一。例如垃圾邮件过滤,中文分词,艾滋病检查,肝癌检查等。

    2.什么是贝叶斯定理?

    贝叶斯定理长这样:

    [​IMG]

    到这来,你可能会说:猴子,说人话,我一看到公式就头大啊。

    其实,我和你一样,不喜欢公式。我们还是从一个例子开始聊起。

    我的朋友小鹿说,他的女神每次看到他的时候都冲他笑,他现在想知道女神是不是喜欢他呢?

    谁让我学过统计概率知识呢,下面我们一起用贝叶斯帮小鹿预测下女神喜欢他的概率有多大,这样小鹿就可以根据概率的大小来决定是否要表白女神。

    首先,我分析了给定的已知信息和未知信息:

    1)要求解的问题:女神喜欢你,记为 A 事件

    2)已知条件:女神经常冲你笑,记为 B 事件

    所以,P(A|B)表示女神经常冲你笑这个事件(B)发生后,女神喜欢你(A)的概率。

    [​IMG]

    从公式来看,我们需要知道这么 3 个事情:

    1)先验概率

    我们把 P(A)称为"先验概率"(Prior probability),也就是在不知道 B 事件的前提下,我们对 A 事件概率的一个主观判断。

    对应这个例子里就是在不知道女神经常对你笑的前提下,来主观判断出女神喜欢一个人的概率。这里我们假设是 50%,也就是不喜欢你,可能不喜欢你的概率都是一半。

    2)可能性函数

    P(B|A)/P(B)称为"可能性函数"(Likelyhood),这是一个调整因子,也就是新信息 B 带来的调整,作用是将先验概率(之前的主观判断)调整到更接近真实概率。

    可能性函数你可以理解为新信息过来后,对先验概率的一个调整。比如我们刚开始看到“人工智能”这个信息,你有自己的理解(先验概率 - 主观判断),但是当你学习了一些数据分析,或者看了些这方面的书后(新的信息),然后你根据掌握的最新信息优化了自己之前的理解(可能性函数 - 调整因子),最后重新理解了“人工智能”这个信息(后验概率)

    如果"可能性函数"P(B|A)/P(B)>1,意味着"先验概率"被增强,事件 A 的发生的可能性变大;

    如果"可能性函数"=1,意味着 B 事件无助于判断事件 A 的可能性;

    如果"可能性函数"<1,意味着"先验概率"被削弱,事件 A 的可能性变小。

    还是刚才的例子,根据女神经常冲你笑这个新的信息,我调查走访了女神的闺蜜,最后发现女神平日比较高冷,很少对人笑,也就是对你有好感的可能性比较大(可能性函数>1)。所以我估计出"可能性函数"P(B|A)/P(B)=1.5(具体如何估计,省去 1 万字,后面会有更详细科学的例子)

    3)后验概率

    P(A|B)称为"后验概率"(Posterior probability),即在 B 事件发生之后,我们对 A 事件概率的重新评估。这个例子里就是在女神冲你笑后,对女神喜欢你的概率重新预测。

    带入贝叶斯公式计算出 P(A|B)=P(A)* P(B|A)/P(B)=50% *1.5=75%

    因此,女神经常冲你笑,喜欢上你的概率是 75%。这说明,女神经常冲你笑这个新信息的推断能力很强,将 50%的"先验概率"一下子提高到了 75%的"后验概率"。

    [​IMG]

    在得到概率值后,小鹿自信满满的发了下面的表白微博:

    [​IMG]

    稍后,果然收到了女神的回复。预测成功。

    [​IMG]

    现在我们再看一遍贝叶斯公式,你现在就能明白这个公式背后的关键思想了:

    我们先根据以往的经验预估一个"先验概率"P(A),然后加入新的信息(实验结果 B),这样有了新的信息后,我们对事件 A 的预测就更加准确。

    [​IMG]

    因此,贝叶斯定理可以理解成下面的式子:

    后验概率(新信息出现后的 A 概率) = 先验概率(A 概率) x 可能性函数(新信息带来的调整)

    贝叶斯的底层思想就是:

    如果我能掌握一个事情的全部信息,我当然能计算出一个客观概率(古典概率)。

    可是生活中绝大多数决策面临的信息都是不全的,我们手中只有有限的信息。既然无法得到全面的信息,我们就在信息有限的情况下,尽可能做出一个好的预测。也就是,在主观判断的基础上,你可以先估计一个值(先验概率),然后根据观察的新信息不断修正(可能性函数)。

    如果用图形表示就是这样的:

    [​IMG]

    其实阿尔法狗也是这么战胜人类的,简单来说,阿尔法狗会在下每一步棋的时候,都可以计算自己赢棋的最大概率,就是说在每走一步之后,他都可以完全客观冷静的更新自己的概率值,完全不受其他环境影响。

    3.贝叶斯定理的应用案例

    前面我们介绍了贝叶斯定理公式,及其背后的思想。现在我们来举个应用案例,你会更加熟悉这个牛瓣的工具。

    为了后面的案例计算,我们需要先补充下面这个知识。

    1.全概率公式

    这个公式的作用是计算贝叶斯定理中的 P(B)。

    假定样本空间 S,由两个事件 A 与 A'组成的和。例如下图中,红色部分是事件 A,绿色部分是事件 A',它们共同构成了样本空间 S。

    [​IMG]

    这时候来了个事件 B,如下图:

    [​IMG]

    全概率公式:

    [​IMG]

    它的含义是,如果 A 和 A'构成一个问题的全部(全部的样本空间),那么事件 B 的概率,就等于 A 和 A'的概率分别乘以 B 对这两个事件的条件概率之和。

    看到这么复杂的公式,记不住没关系,因为我也记不住,下面用的时候翻到这里来看下就可以了。

    案例 1:贝叶斯定理在做判断上的应用

    有两个一模一样的碗,1 号碗里有 30 个巧克力和 10 个水果糖,2 号碗里有 20 个巧克力和 20 个水果糖。

    [​IMG]

    然后把碗盖住。随机选择一个碗,从里面摸出一个巧克力。

    问题:这颗巧克力来自 1 号碗的概率是多少?

    好了,下面我就用套路来解决这个问题,到最后我会给出这个套路。

    第 1 步,分解问题

    1)要求解的问题:取出的巧克力,来自 1 号碗的概率是多少?

    来自 1 号碗记为事件 A1,来自 2 号碗记为事件 A2

    取出的是巧克力,记为事件 B,

    那么要求的问题就是 P(A1|B),也就是取出的是巧克力(B),来自 1 号碗(A1)的概率

    2)已知信息:

    1 号碗里有 30 个巧克力和 10 个水果糖

    2 号碗里有 20 个巧克力和 20 个水果糖

    取出的是巧克力

    第 2 步,应用贝叶斯定理

    [​IMG]

    1)求先验概率

    由于两个碗是一样的,所以在得到新信息(取出是巧克力之前),这两个碗被选中的概率相同,因此 P(A1)=P(A2)=0.5,(其中 A1 表示来自 1 号碗,A2 表示来自 2 号碗)

    这个概率就是"先验概率",即没有做实验之前,来自一号碗、二号碗的概率都是 0.5。

    2)求可能性函数

    P(B|A1)/P(B)

    其中,P(B|A1)表示从 1 号碗中(A1)取出是巧克力(B)的概率。

    因为 1 号碗里有 30 个巧克力和 10 个水果糖,所以 P(B|A1)=巧克力数(30)/(糖果总数 30+10)=75%

    现在贝叶斯公式里只剩 P(B)了,只有求出 P(B)就可以得到答案。

    根据全概率公式,可以用下图求得 P(B):

    [​IMG]

    图中 P(B|A1)是 1 号碗中巧克力的概率,我们根据前面的已知条件,很容易求出。

    同样的,P(B|A2)是 2 号碗中巧克力的概率,也很容易求出(图中已给出)。

    而 P(A1)=P(A2)=0.5

    将这些数值带入公式中就是小学生也可以算出来的事情了。最后 P(B)=62.5%

    所以,可能性函数 P(B|A1)/P(B)=75%/62.5%=1.2。

    可能性函数>1.表示新信息 B 对事情 A1 的可能性增强了。

    3)带入贝叶斯公式求后验概率

    将上述计算结果,带入贝叶斯定理,即可算出 P(A1|B)=60%

    [​IMG]

    这个例子中我们需要关注的是约束条件:抓出的是巧克力。如果没有这个约束条件在,来自一号碗这件事的概率就是 50%了,因为巧克力的分布不均把概率从 50%提升到 60%。

    现在,我总结下刚才的贝叶斯定理应用的套路,你就更清楚了,会发现像小学生做应用题一样简单:

    第 1 步. 分解问题

    简单来说就像做应用题的感觉,先列出解决这个问题所需要的一些条件,然后记清楚哪些是已知的,哪些是未知的。

    1)要求解的问题是什么?

    识别出哪个是贝叶斯中的事件 A(一般是想要知道的问题),哪个是事件 B(一般是新的信息,或者实验结果)

    2)已知条件是什么?

    第 2 步.应用贝叶斯定理

    第 3 步,求贝叶斯公式中的 2 个指标

    1)求先验概率

    2)求可能性函数

    3)带入贝叶斯公式求后验概率

    [​IMG]

    案例 2:贝叶斯定理在医疗行业的应用

    每一个医学检测,都存在假阳性率和假阴性率。假阳性,就是没病,但是检测结果显示有病。假阴性正好相反,有病但是检测结果正常。

    即使检测准确率是 99%,如果医生完全依赖检测结果,也会误诊。也就是说假阳性的情况,根据检测结果显示有病,但是你实际并没有得病。

    举个更具体的例子,因为艾滋病潜伏期很长,所以即便感染了也可能在很长的一段时间,身体没有任何感觉,所以艾滋病检测的假阳性会导致被测人非常大的心理压力。

    你可能会觉得,检测准确率都 99%了,误测几乎可以忽略不计了吧?所以你觉得这人肯定没有患艾滋病了对不对?

    让我们用贝叶斯定理算一下,就会发现你的直觉是错误的。

    假设某种疾病的发病率是 0.001,即 1000 人中会有 1 个人得病。现在有一种试剂可以检验患者是否得病,它的准确率是 0.99,即在患者确实得病的情况下,它有 99%的可能呈现阳性。它的误报率是 5%,即在患者没有得病的情况下,它有 5%的可能呈现阳性。

    现在有一个病人的检验结果为阳性,请问他确实得病的可能性有多大?

    好了,我知道你面对这一大推信息又头大了,我也是。但是我们不是有贝叶斯模板套路嘛,下面开始。

    第 1 步,分解问题

    1)要求解的问题:病人的检验结果为阳性,他确实得病的概率有多大?

    病人的检验结果为阳性(新的信息)记为事件 B,他得病记为事件 A,

    那么要求的问题就是 P(A|B),也就是病人的检验结果为阳性(B),他确实得病的概率(A)

    2)已知信息

    这种疾病的发病率是 0.001,即 P(A)=0.001

    试剂可以检验患者是否得病,准确率是 0.99,也就是在患者确实得病的情况下(A),它有 99%的可能呈现阳性(B),所以 P(B|A)=0.99

    试剂的误报率是 5%,即在患者没有得病的情况下,它有 5%的可能呈现阳性。得病我们记为事件 A,那么没有得病就是事件 A 的反面,记为 A',所以这句话就可以表示为 P(B|A')=5%

    2.应用贝叶斯定理

    [​IMG]

    1)求先验概率

    疾病的发病率是 0.001,即 P(A)=0.001

    2)求可能性函数

    P(B|A)/P(B)

    其中,P(B|A)表示在患者确实得病的情况下(A),试剂呈现阳性的概率,从前面的已知条件中我们已经知道 P(B|A)=0.99

    现在只有求出 P(B)就可以得到答案。根据全概率公式,可以用下图求得 P(B)=0.05094

    [​IMG]

    所以可能性函数 P(B|A)/P(B)=0.99/0.05094=19.4346

    3)带入贝叶斯公式求后验概率

    我们得到了一个惊人的结果,P(A|B)等于 1.94%。

    也就是说,筛查的准确率都到了 99%了,通过体检结果有病(阳性)确实得病的概率也只有 1.94%

    [​IMG]

    你可能会说,再也不相信那些吹的天花乱坠的技术了,说好了筛查准确率那么高,结果筛查的结果对于确诊疾病一点用都没有,这还要医学技术干什么?

    没错,这就是贝叶斯分析告诉我们的。我们拿艾滋病来说,由于发艾滋病实在是小概率事件,所以当我们对一大群人做艾滋病筛查时,虽说准确率有 99%,但仍然会有相当一部分人因为误测而被诊断为艾滋病,这一部分人在人群中的数目甚至比真正艾滋病患者的数目还要高。

    你肯定要问了,那该怎样纠正测量带来这么高的误诊呢?

    造成这么不靠谱的误诊的原因,是无差别地给一大群人做筛查,而不论测量准确率有多高,因为正常人的数目远大于实际的患者,所以误测造成的干扰就非常大了。

    根据贝叶斯定理,我们知道提高先验概率,可以有效的提高后验概率。

    所以解决的办法倒也很简单,就是先锁定可疑的人群,比如 10000 人中检查出现问题的那 10 个人,再独立重复检测一次。因为正常人连续两次体检都出现误测的概率极低,这时筛选出真正患者的准确率就很高了,这也是为什么许多疾病的检测,往往还要送交独立机构多次检查的原因。

    这也是为什么艾滋病检测第一次呈阳性的人,还需要做第二次检测,第二次依然是阳性的还需要送交国家实验室做第三次检测。

    在《医学的真相》这本书里举了个例子,假设检测艾滋病毒,对于每一个呈阳性的检测结果,只有 50%的概率能证明这位患者确实感染了病毒。但是如果医生具备先验知识,先筛选出一些高风险的病人,然后再让这些病人进行艾滋病检查,检查的准确率就能提升到 95%。

    案例 4:贝叶斯垃圾邮件过滤器

    垃圾邮件是一种令人头痛的问题,困扰着所有的互联网用户。全球垃圾邮件的高峰出现在 2006 年,那时候所有邮件中 90%都是垃圾,2015 年 6 月份全球垃圾邮件的比例数字首次降低到 50%以下。

    最初的垃圾邮件过滤是靠静态关键词加一些判断条件来过滤,效果不好,漏网之鱼多,冤枉的也不少。

    2002 年,Paul Graham 提出使用"贝叶斯推断"过滤垃圾邮件。他说,这样做的效果,好得不可思议。1000 封垃圾邮件可以过滤掉 995 封,且没有一个误判。

    因为典型的垃圾邮件词汇在垃圾邮件中会以更高的频率出现,所以在做贝叶斯公式计算时,肯定会被识别出来。之后用最高频的 15 个垃圾词汇做联合概率计算,联合概率的结果超过 90%将说明它是垃圾邮件。

    用贝叶斯过滤器可以识别很多改写过的垃圾邮件,而且错判率非常低。甚至不要求对初始值有多么精确,精度会在随后计算中逐渐逼近真实情况。

    (ps:如果留言想详细了解这个知识的很多,我后面会专门写文章来回答大家)

    4.生活中的贝叶斯思维

    贝叶斯定理与人脑的工作机制很像,这也是为什么它能成为机器学习的基础。

    如果你仔细观察小孩学习新东西的这个能力,会发现,很多东西根本就是看一遍就会。比如我 3 岁的外甥,看了我做俯卧撑的动作,也做了一次这个动作,虽然动作不标准,但也是有模有样。

    同样的,我告诉他一个新单词,他一开始并不知道这个词是什么意思,但是他可以根据当时的情景,先来个猜测(先验概率 / 主观判断)。一有机会,他就会在不同的场合说出这个词,然后观察你的反应。如果我告诉他用对了,他就会进一步记住这个词的意思,如果我告诉他用错了,他就会进行相应调整。(可能性函数 / 调整因子)。经过这样反复的猜测、试探、调整主观判断,就是贝叶斯定理思维的过程。

    同样的,我们成人也在用贝叶斯思维来做出决策。比如,你和女神在聊天的时候,如果对方说出“虽然”两个字,你大概就会猜测,对方后面九成的可能性会说出“但是”。我们的大脑看起来就好像是天生在用贝叶斯定理,即根据生活的经历有了主观判断(先验概率),然后根据搜集新的信息来修正(可能性函),最后做出高概率的预测(后验概率)。

    其实这个过程,就是下图的大脑决策过程:

    [​IMG]

    所以,在生活中涉及到预测的事情,用贝叶斯的思维可以提高预测的概率。你可以分 3 个步骤来预测:

    1.分解问题

    简单来说就像小学生做应用题的感觉,先列出要解决的问题是什么?已知条件有哪些?

    2. 给出主观判断

    不是瞎猜,而是根据自己的经历和学识来给出一个主观判断。

    3.搜集新的信息,优化主观判断

    持续关于你要解决问题相关信息的最新动态,然后用获取到的新信息来不断调整第 2 步的主观判断。如果新信息符合这个主观判断,你就提高主观判断的可信度,如果不符合,你就降低主观判断的可信度。

    比如我们刚开始看到“人工智能是否造成人类失业”这个信息,你有自己的理解(主观判断),但是当你学习了一些数据分析,或者看了些这方面的最新研究进展(新的信息),然后你根据掌握的最新信息优化了自己之前的理解(调整因子),最后重新理解了“人工智能”这个信息(后验概率)。这也就是胡适说的“大胆假设,小心求证”。

    概率的基础知识补充:

    概率(Probability)的本质是什么?

    如何理解条件概率?

    如何更好地理解「选择比努力更重要」?

    参考资料:

    YouTube 英文视频《Thomas Bayes: Probability for Success》

    YouTube 英文视频《Everything You Ever Wanted to Know About Bayes' Theorem But Were Afraid To Ask.》

    贝叶斯垃圾邮件过滤器:http://www.paulgraham.com/spam.html

    贝叶斯垃圾邮件过滤 Wiki:https://en.wikipedia.org/wiki/Naive_Bayes_spam_filtering

    贝叶斯推断及其互联网应用(一)

    《联邦党人文集》背后的统计学幽灵

    猴子数据分析

    阅读原文
     
正在加载...