DATA NOTEBOOK阅读地图
9个主要章节
4张图片与图示
1组对照与手算表

先确认变量、单位与假设,再核对手算、图表和解释范围。

章节速览9 节 · 展开选择

分类阈值讨论给定分数以后采取什么行动,概率校准检查报出的概率是否可信。本篇再向前一步:观察到一条证据以后,原有概率应怎样更新?“发现了异常信号”与“对象一定异常”之间,还隔着基率和证据在不同状态下的发生概率。

我们只讨论两个互斥且穷尽的状态,以及离散事件。所有数字均为教学合成参数,不对应真实检测设备、生产线或业务收益。重点是把证据出现的概率、状态本身的概率,以及需要采取的行动分开,并检查第二条证据是否真的带来新信息。

先认清分母信号出现后,要在信号组内重新计数

先列联合数量,再归一化;不能把“异常时常见”直接倒读成“出现后必然异常”。

再检查新信息复制一次结果,不是又做一次观察

连续更新使用带历史条件的似然比;简化乘法需要额外的条件独立假设。

先给状态与观察命名

设H表示“物品具有目标缺陷”,¬H表示没有该缺陷;E表示筛查出现标记。这里的标记只是一项观察,H才是我们要判断的状态。P(E|H)把讨论范围限制在缺陷物品中,P(H|E)则把范围限制在被标记物品中,两者分母不同。

假定某一明确来源总体中,P(H)=0.02;在H成立时,标记率P(E|H)=0.9;在¬H成立时,标记率P(E|¬H)=0.05。先验指看到本次E以前的概率,不意味着它必须来自主观猜测;它也可以来自适用于当前对象的历史总体模型。反过来,历史比例也不会自动适用于来源已改变的对象。

π = P(H)|先验
观察本次标记之前,对目标状态的概率。必须说明对象来自什么总体,不能把训练集的人为类别配比直接当部署基率。
a、b|两类标记率
a=P(E|H),b=P(E|¬H)。它们描述同一项观察在两个状态下的概率,不是两个后验概率。
P(H|E)|后验
已经知道标记出现后,对H的条件概率。以上概率均无单位;数量表则以件为单位。

一万件物品,把分母直接写出来

把上述概率按10000件展开,得到200件缺陷和9800件非缺陷。缺陷中180件被标记,非缺陷中490件被标记;因此一共有670件带标记,其中只有180件属于H。这个数量表是给定模型的等比例展开,不是实际抽样结果,也不附带真实抽样误差。

真实状态出现标记E,件没有标记¬E,件合计,件
缺陷H18020200
非缺陷¬H49093109800
合计670933010000
一万件中缺陷200件,非缺陷9800件;两组分别标记180和490件,因此标记组缺陷比例为180除以670
图1|原创合成数量账本。框与箭头表示分类和筛选关系,面积不代表比例。查看原图

所以P(H|E)=180/670=18/67≈26.865672%。它显著高于原来的2%,但远低于90%。90%回答的是“缺陷物品有多少会被标记”,26.87%回答的是“已标记物品有多少缺陷”,不是同一个问题。

这也解释了为什么“标记组里非缺陷较多”并不意味着标记毫无信息。先验只有2%,观察后升到约26.87%,证据确实改变了相对支持程度;它是否足以触发人工复核或直接报废,还需要决策成本。不要把“没有达到确定性”写成“没有提供信息”。

从联合概率推导归一化公式

条件概率定义给出P(H∩E)=P(E|H)P(H),也等于P(H|E)P(E)。同时,H与¬H把总体分成两部分,P(E)=P(E|H)P(H)+P(E|¬H)P(¬H)。代入即可得到下式;条件概率与全概率的基础见文献1。

公式B1|二状态贝叶斯更新
P(H|E) = aπ / [aπ + b(1−π)]

H与¬H互斥且穷尽,a=P(E|H)、b=P(E|¬H),要求分母P(E)>0。各项均无单位。先算联合概率aπ与b(1−π),再在观察到E的范围内归一化。

分母不是一个可以随意省去的“比例常数”。在只比较两个状态的相对大小时,未归一化权重有用;要报告概率,就必须让两种状态的后验之和为1。只写aπ=0.018,得到的是总体中“缺陷且标记”的联合概率,不是标记以后的缺陷概率。

没有标记也能更新。用1−a=0.1与1−b=0.95代入,得到P(H|¬E)=20/9330=2/933≈0.214362%。这不是零风险,也不等于“有90%把握没有缺陷”。正负证据要分别写出事件,不要让“准确率”一个词承担所有条件。

自检1:180件缺陷被标记,为什么后验不除以200?

除以200得到的是缺陷组的标记率P(E|H)。后验已经知道E发生,分母应是全部670件被标记物品;因此是180/670。先说明已知条件,再选择分母。

赔率让证据的作用更容易追踪

概率p对应赔率p/(1−p),反向转换为赔率/(1+赔率)。这里的“赔率”只是概率比,不涉及收益或投注。将H与¬H的贝叶斯表达相除,公共分母P(E)消去,得到先验赔率乘似然比;赔率形式见文献2,以下数值由本篇独立展开。

公式B2|赔率乘似然比
O(H) = P(H)/P(¬H)
LR(E) = P(E|H)/P(E|¬H)
O(H|E) = O(H) × LR(E)

此有限正赔率写法要求0<P(H)<1、P(E|¬H)>0,并且E可观察。比值均无单位。若似然分母为0,应回到B1检查联合概率,不能盲算0×∞。

本例先验赔率为1/49,标记的似然比为0.9/0.05=18,后验赔率便是18/49。再转回概率:18/(18+49)=18/67。似然比18的意思是:这个标记在H下出现的概率,是在¬H下的18倍;它没有直接说H的后验是18倍,更没有说后验是18。

若LR=1,观察对这两个状态的相对权重没有改变;LR>1倾向H,0<LR<1倾向¬H。但似然比强弱与最终概率高低仍是两回事,因为起点可能不同。先验为0.2时,同一个LR=18得到后验9/11≈81.818182%,与先验0.02时差异很大。

同一条证据,在不同基率下意味着什么

保持a=0.9、b=0.05,只改变π。π依次为0.01、0.02、0.10、0.20时,后验分别为2/13、18/67、2/3、9/11。这不是说现实中标记率总能跨总体保持;它是一次把其他条件固定的敏感性分析。

固定标记率0.9和0.05,先验1%、2%、10%、20%对应后验约15.38%、26.87%、66.67%、81.82%
图2|原创模型对照,柱高按后验概率线性绘制。只改变先验,不表示真实数据集之间只存在基率差异。查看原图

若希望标记后概率至少0.8,在这组固定标记率下,需要先验赔率至少4/18=2/9,即π至少2/11≈18.181818%。这只是反解模型条件,不是给系统选择目标概率或业务阈值的建议。

现实中总体改变还可能影响a和b,例如物品来源、缺陷类型和标记过程一起变化。因此迁移时应分别检查先验与条件分布,而不是只换一个基率数字。参数来自有限样本时,也要考虑它们的估计不确定性;把估计值代入公式,并不会消除误差。

连续观察:第二个乘数必须带历史条件

再观察F时,新的信息是“已经知道E以后,F又发生了”。准确的第二个似然比应是P(F|H,E)/P(F|¬H,E)。只有在H和¬H各自内部,F与E均满足相应条件独立,才能分别去掉E。

公式B3|保留历史的两次更新
O(H|E,F) = O(H) × LR(E)
× P(F|H,E) / P(F|¬H,E)

要求所用条件事件有正概率,并使比值分母为正。若要把第二项写成P(F|H)/P(F|¬H),需在两个状态下分别满足条件独立;仅说“做了两次”或“总体上不相关”都不够。

在一个明确假定两次标记条件独立、且两类标记率相同的新模型中,两次都标记可用18×18:后验赔率324/49,概率324/373≈86.863271%。这里第二次是一次新的随机观察,不能把这一结果套在复制日志或重发同一条消息上。

反例非常简单:令F就是E的副本。知道E已经发生以后,P(F|H,E)=P(F|¬H,E)=1,第二个似然比是1。后验仍为18/67,而非324/373。即使两个数据字段名字不同,只要它们来自同一个原始观测,也可能包含重复信息。条件独立的用途与边界可参看文献3。

自检2:相同标记出现两行,是否就应把似然比平方?

先检查两行是否来自独立的新观察。若第二行只是第一行的副本,给定第一行后它必然相同,第二次似然比为1。只有对应的条件独立假设成立时,才可使用两个相同的无历史似然比。

零概率、错误模型与行动边界

若a=b=0,E在当前模型中不可能出现,B1分母为0;此时看到E意味着模型、记录或事件定义至少有一处需要复核,不能硬填一个后验。若先验把H设为0,在有限非零证据下也不会自动把H恢复成正概率;这提示我们,结构性排除某种状态是一项很强的建模决定。

后验0.2687也不会自行决定报废。若行动是“人工复核”,与行动是“销毁物品”,误判成本可能完全不同。先冻结概率模型与核对口径,再通过阈值与决策成本讨论行动,避免用行为结果反过来证明概率计算正确。

自检3:似然比18,能否直接把2%乘18当作后验?

不能。乘法作用于赔率1/49,得到18/49后还需转回概率18/67。概率直接乘18得到36%,遗漏了观察组的重新归一化。

自检4:模型说观察不可能发生,实际却看到了,怎样处理?

检查事件定义、数据记录和模型支持范围。分母为0时,当前离散条件概率公式没有定义;不能把计算失败隐藏成0%、50%或100%。修订模型应有新的依据,而非仅为了让程序返回一个数字。

一份可复核的更新记录

实际记录至少保留对象来源、H与¬H的定义、先验、观察事件、两类观察概率以及参数版本。连续证据还要记录原始来源和依赖关系。把联合数量与归一化结果并列,通常比只留下一个后验小数更容易发现分母错误。

本文内部用精确有理数联合表核对直接归一化、赔率与顺序更新,并单独构造复制证据的依赖模型。核对的是给定离散概率模型中的算术恒等式,不是现实标记过程的有效性,也没有测量真实分类收益。

参考文献与使用位置

  1. Jeremy Orloff、Jonathan Bloom,MIT 18.05:Reading 3,Conditional Probability, Independence and Bayes' Theorem,2022。 官方讲义。用于条件概率、全概率及基率的基础关系;本篇物品数量与敏感性算例独立设置。返回B1。
  2. John Lipor,Portland State University,EE520:Session 2,Bayes' Rule and Sequential Updating。 官方课程讲义。核对赔率形式与连续观察的条件独立假设,不复制其搜索算例。返回B2。
  3. MIT 6.034:Probabilistic Inference Notes,2010。 官方概率推断笔记。用于检查把多条观察拆乘所需的独立假设;复制证据反例由本文直接构造。返回B3。