2. 概率基础与贝叶斯公式

贝叶斯参数估计的全部计算,都建立在三个概率工具之上:条件概率、全概率公式和贝叶斯公式。它们不是三条孤立的公式,而是一条层层递进的推导链——条件概率定义了"在已知信息下如何重新分配概率",全概率公式把未知原因的影响汇总成可观测结果的总概率,贝叶斯公式则把这条链反向使用,从观测结果反推原因的可信度。理解这条链的方向与含义,比记住公式本身更重要。

条件概率:信息如何改变概率

设事件 A、B 满足 P(B) > 0,则在 B 发生的条件下 A 的条件概率定义为

P(A | B) = P(A ∩ B) / P(B)

直观理解:把样本空间从"全部可能"收缩到"B 已发生"这一子集,再在这个子集内看 A 占多大比例。分母 P(B) 起归一化作用,保证条件概率在 B 上仍然是一个合法的概率分布。条件概率不改变事件本身,只改变我们衡量它的参照范围。例如某疾病在人群中的患病率为 1%,检测灵敏度(患病者检出阳性)为 95%,特异度(健康者检出阴性)为 90%。记 D 为患病、+ 为阳性,则 P(+|D)=0.95,P(−|D^c)=0.90,即假阳性率 P(+|D^c)=0.10。这些数字都是条件概率,它们描述的是"给定状态后检测的表现",而不是"检测阳性后患病的概率"——后者才是贝叶斯公式要回答的问题。

全概率公式:把原因汇总为结果

若事件 B₁, B₂, …, Bₙ 构成样本空间的一个划分(两两互斥且并为全集),且每个 P(Bᵢ) > 0,则对任意事件 A:

P(A) = Σᵢ P(A | Bᵢ) · P(Bᵢ)

它的作用是"按原因分解结果":A 的总概率等于各条成因路径上"先到 Bᵢ、再由 Bᵢ 到 A"的联合概率之和。回到检测例子,阳性结果的总概率为

P(+) = P(+|D)·P(D) + P(+|D^c)·P(D^c)
     = 0.95 × 0.01 + 0.10 × 0.99
     = 0.0095 + 0.099 = 0.1085

也就是说,在患病率 1% 的人群中随机检测一人,阳性率约为 10.85%。注意其中 0.099 来自健康人的假阳性,远大于真阳性的 0.0095——这个结构差异正是下一节"阳性并不等于大概率患病"的根源。

贝叶斯公式:从结果反推原因

由条件概率的对称性 P(A∩B) = P(A|B)P(B) = P(B|A)P(A),解出 P(B|A) 即得贝叶斯公式:

P(Bᵢ | A) = P(A | Bᵢ) · P(Bᵢ) / P(A)
          = P(A | Bᵢ) · P(Bᵢ) / Σⱼ P(A | Bⱼ) · P(Bⱼ)

公式中每一项都有明确的名字与角色:P(Bᵢ) 是先验概率,代表观测前的信念;P(A|Bᵢ) 是似然,代表该原因下出现观测的可能性;分母 P(A) 是证据的边际概率,由全概率公式算出,只起归一化作用;P(Bᵢ|A) 是后验概率,即观测后更新过的信念。用检测例子计算:

P(D | +) = 0.95 × 0.01 / 0.1085 ≈ 0.0876

阳性者真正患病的概率只有约 8.76%,远低于多数人的直觉。原因不在检测不准,而在先验极低:健康人基数大,10% 的假阳性率乘以 0.99 的先验,贡献了绝大部分阳性样本。这说明后验由先验与似然共同决定,任何一方被忽略都会导致误判。若对高风险人群(先验 20%)做同一检测,则 P(D|+) = 0.95×0.20 / (0.95×0.20 + 0.10×0.80) ≈ 0.704,结论截然不同。

条件概率 → 全概率 → 贝叶斯:一条可逆的推导链 先验 P(Bᵢ) 观测前的信念 P(D)=0.01,P(Dᶜ)=0.99 似然 P(A|Bᵢ) 给定原因下观测的概率 P(+|D)=0.95,P(+|Dᶜ)=0.10 联合概率 P(A|Bᵢ)·P(Bᵢ) 0.0095 与 0.099 证据 P(A) 全概率公式求和 P(+)=0.1085 后验 P(Bᵢ|A) ≈ 0.0876 正向(全概率):先验 × 似然 → 联合 → 求和得证据 P(A),回答"结果总体多常见"。 反向(贝叶斯):联合 ÷ 证据 → 后验,回答"看到结果后,某原因有多可信"。 示例数值为教学假设:患病率 1%、灵敏度 95%、特异度 90%。

连续情形:从求和到积分,从分布到分布

参数估计中,未知量 θ 通常是连续的,观测 x 也是连续的。此时把离散公式中的求和换成积分、概率换成密度,形式完全平行:

先验密度:π(θ)
似然:    p(x | θ)
边际密度:m(x) = ∫ p(x | θ) · π(θ) dθ
后验密度:π(θ | x) = p(x | θ) · π(θ) / m(x)

由于 m(x) 不含 θ,它只是归一化常数,实际计算中常写作 π(θ|x) ∝ p(x|θ)·π(θ),即"后验正比于似然乘先验"。这是贝叶斯计算最常用的工作形式。举一个可手算的共轭示例:设 θ 为某产品的次品率,先验取 Beta(2, 8)(均值 0.2,体现"次品率偏低"的既有认识);抽检 n=20 件,发现 k=3 件次品,似然为二项分布。Beta 先验与二项似然共轭,后验仍为 Beta 分布,参数直接相加:

先验 Beta(a, b) = Beta(2, 8)
观测:n = 20,k = 3
后验 Beta(a + k, b + n − k) = Beta(5, 25)
后验均值 = 5 / (5 + 25) ≈ 0.1667
后验标准差 = sqrt(5×25 / (30²×31)) ≈ 0.0669

先验均值 0.2 被数据拉低到约 0.167,同时后验分布比先验更集中(标准差由约 0.126 降至约 0.067),体现了"数据越多、后验越尖锐"的一般规律。若换用均匀先验 Beta(1,1),后验为 Beta(4,18),均值约 0.182,说明先验的选择会实质影响小样本下的结论,这正是贝叶斯分析必须交代先验来源的原因。

判断该用哪条公式,可以看问题方向:已知原因求结果的总概率,用全概率公式;已知结果反推原因的可信度,用贝叶斯公式;只想在某个已知条件下重新衡量概率,用条件概率定义即可。三者共用同一个乘法规则,方向不同而已。
常见错误有三类:一是把 P(+|D) 与 P(D|+) 混为一谈,这是典型的"条件方向颠倒";二是计算后验时漏掉分母中的其他成因项,导致概率不归一;三是在连续情形下把密度当作概率直接比较,实际上密度值可以大于 1,只有积分才有概率含义。