Skip to content

贝叶斯决策理论

贝叶斯决策理论把不确定性与行动代价分开处理:先由先验与似然得到后验,再用损失函数把后验逐点映射为最优行动。

机器学习里,「贝叶斯」这个词至少承担了三项不同的工作,而它们经常被混为一谈。

第一项是贝叶斯公式。它是一条关于条件概率的恒等式,本身没有争议,也不要求任何哲学立场:

p(θ∣D)=p(D∣θ) p(θ)p(D).p(\theta\mid D)=\frac{p(D\mid\theta)\,p(\theta)}{p(D)}.

第二项是贝叶斯决策。它问的不是「参数是多少」,而是「在还不确定的情况下,应该采取什么行动」。给定后验分布与一个损失函数,最优行动就是使后验期望损失最小的那个。

第三项是贝叶斯学习。它把模型参数也看成随机变量,用后验分布而不是点估计来描述知识状态;预测时对参数做积分,而不是代入一个「最佳值」。

三者可以独立使用,也可以串成一条链:先验与似然给出后验,后验与损失给出行动,行动的结果又成为新的数据。本文沿着这条链走一遍:先补齐概率论的语言,再建立决策论的骨架,最后看这套骨架如何解释机器学习里的分类器、正则化、错误率、泛化界与不确定性。

这条链的终点是一个看起来过于理想化的量,即贝叶斯风险。它假设真实分布已知,因而在实践中永远算不出来。但正是这个「算不出来的量」充当了学习理论的坐标原点:它告诉我们,在给定问题里错误率的下限在哪里,模型的哪些部分在逼近这个下限,哪些部分只是在拟合噪声。

1. 一个决策问题的解剖

把所有监督学习问题写成同一个形式,只需要四样东西:

  • 状态 Y∈YY\in\mathcal Y:真实类别、待预测的实数值,或环境里我们关心但看不见的量;
  • 观测 X∈XX\in\mathcal X:我们拿到的特征;
  • 行动 a∈Aa\in\mathcal A:模型输出的预测、决策或控制量;
  • 损失 L:Y×A→RL:\mathcal Y\times\mathcal A\to\mathbb R:当真实状态为 yy、采取行动 aa 时付出的代价。

四者合起来称为一个决策问题。一个决策规则(decision rule)是一个映射

δ:X→A,x↦δ(x),\delta:\mathcal X\to\mathcal A,\qquad x\mapsto \delta(x),

它的风险(risk)定义为

R(δ)=E(X,Y)∼P[L(Y,δ(X))].R(\delta)=\mathbb E_{(X,Y)\sim P}\big[L(Y,\delta(X))\big].

这里期望是对联合分布 P(X,Y)P(X,Y) 取的。这个分布在实践中未知,于是我们退而求其次:用样本来估计它,或者直接在一个受限的函数类里最小化经验风险。这个「退而求其次」的过程,就是统计学习理论。

同一套语言可以覆盖很多看似不同的任务。医疗诊断中,YY 是是否患病,aa 是是否治疗,损失把漏诊与误诊的代价写成两个不同的数;垃圾邮件过滤中,YY 是邮件类别,aa 是放行或拦截;回归中,aa 是一个实数,损失通常是平方误差;强化学习中,XX 是历史状态,aa 是动作,损失是负回报。它们的共同点是:不确定性由 P(Y∣X)P(Y\mid X) 描述,偏好由 LL 描述,决策则是两者的交汇。

把联合风险按 XX 分解,可以得到一个更有用的形式:

R(δ)=∫X(∫YL(y,δ(x)) p(y∣x) dy)p(x) dx=EX[r(δ(X),X)],\begin{aligned} R(\delta) &=\int_{\mathcal X}\left(\int_{\mathcal Y}L(y,\delta(x))\,p(y\mid x)\,dy\right)p(x)\,dx\\ &=\mathbb E_X\big[r(\delta(X),X)\big], \end{aligned}

其中

r(a,x)=E[L(Y,a)∣X=x]r(a,x)=\mathbb E\big[L(Y,a)\mid X=x\big]

称为在 xx 处采取行动 aa 的条件风险。这个分解是整个理论的第一块基石:联合风险被拆成了「每个观测点上的局部风险」再对 XX 平均。于是,只要在每个 xx 上把 r(a,x)r(a,x) 压到最低,全局风险就自动最低,不同 xx 之间没有耦合。

2. 概率的语言

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

条件概率的定义是

P(A∣B)=P(A∩B)P(B),P(B)>0.P(A\mid B)=\frac{P(A\cap B)}{P(B)},\qquad P(B)>0.

它不是在「重新定义」概率,而是在缩小样本空间:一旦知道 BB 发生,所有不含 BB 的结果都被排除,原来的概率需要在剩下的世界里重新归一化。连续情形下,条件分布由条件密度 p(y∣x)=p(x,y)/p(x)p(y\mid x)=p(x,y)/p(x) 给出;更一般地,它是测度论意义下的正则条件分布,本文不展开技术细节。

由条件概率的定义立刻得到乘法公式

P(A∩B)=P(A∣B)P(B)=P(B∣A)P(A),P(A\cap B)=P(A\mid B)P(B)=P(B\mid A)P(A),

两边同除以 P(B)P(B),就是贝叶斯公式。

2.2 贝叶斯公式

把 BB 换成数据 DD,把 AA 换成参数或隐变量 θ\theta,得到

p(θ∣D)=p(D∣θ) p(θ)p(D),p(D)=∫p(D∣θ) p(θ) dθ.p(\theta\mid D)=\frac{p(D\mid\theta)\,p(\theta)}{p(D)}, \qquad p(D)=\int p(D\mid\theta)\,p(\theta)\,d\theta .

公式的四个部分各有名字:

名称记号含义
先验p(θ)p(\theta)看到数据前的信念
似然p(D∣θ)p(D\mid\theta)给定参数时数据的概率
证据p(D)p(D)与 θ\theta 无关的归一化常数,也是模型比较的关键
后验p(θ∣D)p(\theta\mid D)看到数据后的信念

下面的图把这条链画全。注意「后验」不是终点:它只是决策的输入。

贝叶斯决策流程:先验与似然合成后验,后验与损失共同决定最优行动
贝叶斯决策流程:先验与似然合成后验,后验与损失共同决定最优行动

公式右端的 p(D)p(D) 与 θ\theta 无关,所以在求后验时常常被略去,写成

p(θ∣D)∝p(D∣θ) p(θ).p(\theta\mid D)\propto p(D\mid\theta)\,p(\theta).

但这个省略只在「给定模型」内部成立。一旦要比较不同模型,p(D)p(D) 就是主角,第 6.5 节会回到这一点。

还要强调一点:贝叶斯公式本身只是概率论的一条定理,接受它不需要接受任何统计哲学。真正带有立场的是下一步——把 θ\theta 看成随机变量,并为它指定先验。

2.3 Beta–Bernoulli:一个可以手算的更新

设 θ∈[0,1]\theta\in[0,1] 是一枚硬币正面向上的概率,先验取 Beta 分布

p(θ)=θα−1(1−θ)β−1B(α,β).p(\theta)=\frac{\theta^{\alpha-1}(1-\theta)^{\beta-1}}{B(\alpha,\beta)}.

观测到 nn 次投掷中 kk 次正面,似然为

p(D∣θ)=θk(1−θ)n−k.p(D\mid\theta)=\theta^k(1-\theta)^{n-k}.

两者相乘,后验的核为

p(θ∣D)∝θα+k−1(1−θ)β+n−k−1,p(\theta\mid D)\propto \theta^{\alpha+k-1}(1-\theta)^{\beta+n-k-1},

这正是 Beta(α+k, β+n−k)\mathrm{Beta}(\alpha+k,\ \beta+n-k)。归一化常数不用算,因为 Beta 分布的形式已经把它确定了。这就是共轭先验的好处:后验与先验属于同一个分布族,更新只需要改两个参数。

后验均值是

E[θ∣D]=α+kα+β+n,\mathbb E[\theta\mid D]=\frac{\alpha+k}{\alpha+\beta+n},

把它改写成加权平均:

α+kα+β+n=α+βα+β+n⋅αα+β+nα+β+n⋅kn.\frac{\alpha+k}{\alpha+\beta+n} = \frac{\alpha+\beta}{\alpha+\beta+n}\cdot\frac{\alpha}{\alpha+\beta} + \frac{n}{\alpha+\beta+n}\cdot\frac{k}{n}.

右边第一项是先验均值,第二项是样本频率。权重由「先验有多强」(α+β\alpha+\beta)与「数据有多少」(nn)决定。数据越多,先验越被稀释。这就是最朴素的收缩(shrinkage)。

如果下一步要预测新的一次投掷,贝叶斯答案是后验预测分布

p(xnew=1∣D)=E[θ∣D]=α+kα+β+n.p(x_{\text{new}}=1\mid D)=\mathbb E[\theta\mid D]=\frac{\alpha+k}{\alpha+\beta+n}.

这里恰好等于后验均值,是因为伯努利概率是 θ\theta 的线性函数。对一般模型,后验预测分布并不等于「把后验均值代进去」,这个区别在第 6 节会反复出现。

2.4 先验是什么,不是什么

先验常被批评为「主观」。这个批评有一半是对的:先验确实是一个假设,而且不同先验会给出不同后验。但「主观」不等于「任意」。先验可以被数据稀释:在正则条件下,后验的渐近行为与先验无关(见 6.4)。先验也可以层级化:给超参数再加一层先验,让数据同时决定参数与超参数。先验还可以被检验:通过后验预测检验、敏感性分析,看结论是否依赖于某个特定选择。

更重要的是,非贝叶斯方法同样有假设。模型类、正则项、损失函数、优化算法,都是假设,只是它们不以「先验」的名字出现。贝叶斯的优点不是没有假设,而是把假设写在了明处,并且提供了一个用概率语言修改假设的机制。

第 7.8 节会给出一个更强的说法:在相当一般的条件下,所有可容许(admissible)的决策规则都可以看成某个先验(或先验序列)下的贝叶斯规则。换句话说,贝叶斯规则不是一个「流派」,而是可容许决策规则的一般形式。这个结论把「为什么要用先验」从个人偏好提升到了决策论的结构性事实。

3. 从后验到行动

3.1 条件风险与逐点最优

有了后验 p(y∣x)p(y\mid x),决策问题就变成了一个逐点的优化问题。

定义 3.1(贝叶斯决策规则) 若决策规则 δ∗\delta^* 满足

δ∗(x)∈arg⁡min⁡a∈A E[L(Y,a)∣X=x]对几乎所有的 x,\delta^*(x)\in\arg\min_{a\in\mathcal A}\ \mathbb E\big[L(Y,a)\mid X=x\big] \qquad \text{对几乎所有的 }x,

则称 δ∗\delta^* 为贝叶斯决策规则(Bayes decision rule)。

定理 3.1(逐点最优即全局最优) 对任意决策规则 δ\delta,都有

R(δ) ≥ R(δ∗)=EX[min⁡a∈Ar(a,X)].R(\delta)\ \ge\ R(\delta^*)=\mathbb E_X\Big[\min_{a\in\mathcal A}r(a,X)\Big].

证明。 由塔性质,

R(δ)=EX[E[L(Y,δ(X))∣X]]=EX[r(δ(X),X)].R(\delta)=\mathbb E_X\Big[\mathbb E\big[L(Y,\delta(X))\mid X\big]\Big] =\mathbb E_X\big[r(\delta(X),X)\big].

对每个 xx,r(δ(x),x)≥min⁡ar(a,x)r(\delta(x),x)\ge \min_a r(a,x),因此逐点取期望即得结论。□\square

证明只用了两件事:期望的塔性质,以及「逐点最小 ≤ 先对点取最小再平均」这样一个显而易见的事实。它不要求损失凸、不要求后验有密度、也不要求行动空间是有限的。这种一般性是贝叶斯决策理论最漂亮的地方之一。严格地说,若逐点下确界不可达,把 min⁡\min 换成 inf⁡\inf、把 arg⁡min⁡\arg\min 换成「取一列逼近的行动」即可,证明不变。

把上式改写成差值,可以得到一个更有解释力的恒等式:

R(δ)−R∗=EX[r(δ(X),X)−min⁡ar(a,X)⏟≥0].R(\delta)-R^* =\mathbb E_X\Big[\underbrace{r(\delta(X),X)-\min_{a}r(a,X)}_{\ge 0}\Big].

右边每一项都是非负的,称为在 xx 处的逐点遗憾(pointwise regret)。它说明:一个决策规则的次优性完全来自它在某些 xx 上没有选到局部最优行动。不同 xx 之间的错误不会相互抵消,也不会相互放大。

3.2 0-1 损失:最大后验

分类问题最常用的损失是 0-1 损失:

L(y,a)=1{y≠a}.L(y,a)=\mathbf 1\{y\ne a\}.

此时条件风险为

r(a,x)=E[1{Y≠a}∣X=x]=1−p(a∣x),r(a,x)=\mathbb E\big[\mathbf 1\{Y\ne a\}\mid X=x\big]=1-p(a\mid x),

最小化 r(a,x)r(a,x) 等价于最大化后验概率:

δ∗(x)=arg⁡max⁡y∈Y p(y∣x).\delta^*(x)=\arg\max_{y\in\mathcal Y}\ p(y\mid x).

这就是最大后验决策(maximum a posteriori decision)。二分类时,若 η(x)=p(Y=1∣x)\eta(x)=p(Y=1\mid x),规则简化为

δ∗(x)={1,η(x)≥12,0,η(x)<12,\delta^*(x)= \begin{cases} 1, & \eta(x)\ge \tfrac12,\\[2pt] 0, & \eta(x)<\tfrac12, \end{cases}

决策边界是集合 {x:η(x)=1/2}\{x:\eta(x)=1/2\}。当 η(x)=1/2\eta(x)=1/2 时两种行动同样好,实践中用固定的打破平局规则即可。

3.3 平方损失:后验均值

回归问题常用平方损失 L(y,a)=(y−a)2L(y,a)=(y-a)^2。此时

r(a,x)=E[(Y−a)2∣X=x]=Var(Y∣X=x)+(a−E[Y∣X=x])2,\begin{aligned} r(a,x) &=\mathbb E\big[(Y-a)^2\mid X=x\big]\\ &=\mathrm{Var}(Y\mid X=x)+\big(a-\mathbb E[Y\mid X=x]\big)^2, \end{aligned}

第二项在 a=E[Y∣X=x]a=\mathbb E[Y\mid X=x] 时取到最小值。因此平方损失下的贝叶斯行动是后验均值,也就是回归函数

f∗(x)=E[Y∣X=x].f^*(x)=\mathbb E[Y\mid X=x].

这个结论还可以从另一个角度理解:后验均值是 YY 在 L2L^2 意义下最好的预测,任何其他预测都会多出一项方差。

3.4 绝对损失:后验中位数

若损失为 L(y,a)=∣y−a∣L(y,a)=|y-a|,条件风险为 r(a,x)=E[∣Y−a∣∣X=x]r(a,x)=\mathbb E[|Y-a|\mid X=x]。对 aa 求次梯度,在 YY 的分布连续时它就是导数:

∂ar(a,x)=P(Y<a∣x)−P(Y>a∣x),\partial_a r(a,x)=P(Y<a\mid x)-P(Y>a\mid x),

令 00 属于次梯度,等价于

P(Y≤a∣x)≥12,P(Y≥a∣x)≥12.P(Y\le a\mid x)\ge \frac12,\qquad P(Y\ge a\mid x)\ge \frac12 .

满足这两个不等式的 aa 就是后验中位数。直观上,中位数把后验概率一分为二,左右两侧的「质量」相等,因此向左或向右移动 aa 都不会降低期望绝对误差。

平方损失对异常值敏感,因为误差被平方;绝对损失只看见误差的大小。这个差异在决策层面的表现,就是均值与中位数的差异。

3.5 非对称损失:阈值不再在 1/21/2

现实中的两类错误很少同等昂贵。设预测为正类的代价为 cFPc_{\mathrm{FP}}(假阳性),预测为负类的代价为 cFNc_{\mathrm{FN}}(假阴性),正确分类的代价为 00。此时

r(1,x)=cFP p(0∣x),r(0,x)=cFN p(1∣x).r(1,x)=c_{\mathrm{FP}}\,p(0\mid x),\qquad r(0,x)=c_{\mathrm{FN}}\,p(1\mid x).

选择正类当且仅当 r(1,x)≤r(0,x)r(1,x)\le r(0,x),即

p(1∣x) ≥ cFPcFP+cFN.p(1\mid x)\ \ge\ \frac{c_{\mathrm{FP}}}{c_{\mathrm{FP}}+c_{\mathrm{FN}}}.

阈值不再固定在 1/21/2,而是由两类代价的比值决定。若漏诊的代价是误诊的十倍,阈值就降到 1/111/11 左右,模型会倾向于把更多样本判为正类。

这个式子也说明:类别不平衡本身不是问题,损失不对称才是问题。重采样改变了训练分布,相当于改变了先验;如果最终要解决的是原分布上的决策问题,那么重采样之后还需要把阈值校正回来。

把常见损失与对应的贝叶斯行动放在一起:

损失贝叶斯行动
0-1 损失后验众数
平方损失后验均值
绝对损失后验中位数
非对称 0-1 损失阈值 cFP/(cFP+cFN)c_{\mathrm{FP}}/(c_{\mathrm{FP}}+c_{\mathrm{FN}})
带拒绝的 0-1 损失最大后验低于 1−c1-c 时拒绝

最后一行是拒绝选项(reject option):除了预测类别,模型还可以说「我不确定」。若拒绝的代价为 c∈(0,1)c\in(0,1),预测的期望代价为 1−max⁡yp(y∣x)1-\max_y p(y\mid x),拒绝的代价为 cc,因此最优规则是

拒绝⟺max⁡yp(y∣x)<1−c.\text{拒绝}\quad\Longleftrightarrow\quad \max_y p(y\mid x)<1-c.

这就是 Chow 规则。它把「不确定性」直接翻译成了一个可执行的动作,而不是一个只用来画图的分数。

3.6 贝叶斯风险

把定理 3.1 中的最优值单独命名:

R∗=inf⁡δR(δ)=EX[min⁡a∈Ar(a,X)].R^*=\inf_{\delta}R(\delta) =\mathbb E_X\Big[\min_{a\in\mathcal A}r(a,X)\Big].

R∗R^* 称为贝叶斯风险(Bayes risk),也叫贝叶斯错误率(分类情形)。它是这个问题在信息层面的下限:即使知道真实分布、知道每个 xx 上的后验,也只能做到这里。它不是某个算法的性能指标,而是问题本身的属性。

对二分类 0-1 损失,把 3.2 节的规则代回去,得到

R∗=E[min⁡{η(X),1−η(X)}].R^*=\mathbb E\big[\min\{\eta(X),1-\eta(X)\}\big].

对多分类,则是

R∗=1−E[max⁡y∈Yp(y∣X)].R^*=1-\mathbb E\Big[\max_{y\in\mathcal Y}p(y\mid X)\Big].

这两个式子在理论里出现得非常多,因为它们把「分类有多难」压缩成了一个标量。

进一步,二分类的超额风险有一个很干净的形式:

R(δ)−R∗=E[∣2η(X)−1∣ 1{δ(X)≠δ∗(X)}].R(\delta)-R^* =\mathbb E\Big[|2\eta(X)-1|\,\mathbf 1\{\delta(X)\ne\delta^*(X)\}\Big].

证明。 对固定的 xx,若 δ(x)=δ∗(x)\delta(x)=\delta^*(x),逐点遗憾为 00;若不等,则所选行动的后验概率是 max⁡{η,1−η}\max\{\eta,1-\eta\},最优行动的后验概率是 min⁡{η,1−η}\min\{\eta,1-\eta\},两者之差为

max⁡{η,1−η}−min⁡{η,1−η}=∣2η−1∣.\max\{\eta,1-\eta\}-\min\{\eta,1-\eta\}=|2\eta-1|.

对 XX 取期望即得。□\square

这个恒等式很有解释力:只有越过决策边界的样本才贡献错误,而且每个样本的权重是 ∣2η(X)−1∣|2\eta(X)-1|。后验接近 1/21/2 的样本位于边界附近,即使分错,代价也小;后验接近 00 或 11 的样本一旦分错,代价就接近 11。换句话说,贝叶斯错误率不是由边界的长度决定的,而是由边界附近的后验质量决定的。

回归情形也有类似的恒等式。设 f∗(x)=E[Y∣X=x]f^*(x)=\mathbb E[Y\mid X=x],则对任意 f^\hat f,

R(f^)−R∗=E[(f^(X)−f∗(X))2].R(\hat f)-R^*=\mathbb E\big[(\hat f(X)-f^*(X))^2\big].

证明。 把平方展开:

E[(Y−f^(X))2]−E[(Y−f∗(X))2]=E[(f^(X)−f∗(X))2]+2 E[(Y−f∗(X))(f∗(X)−f^(X))].\begin{aligned} &\mathbb E\big[(Y-\hat f(X))^2\big]-\mathbb E\big[(Y-f^*(X))^2\big]\\ &=\mathbb E\big[(\hat f(X)-f^*(X))^2\big] +2\,\mathbb E\big[(Y-f^*(X))(f^*(X)-\hat f(X))\big]. \end{aligned}

交叉项对 XX 取条件期望后,E[Y−f∗(X)∣X]=0\mathbb E[Y-f^*(X)\mid X]=0,因此消失。□\square

这两个恒等式把「模型误差」翻译成了「决策遗憾」:分类器的错误只发生在越过边界的点上,回归器的误差则是函数逼近误差的平方。贝叶斯风险就是那个无法被任何模型消除的起点。

3.7 MAP 不是贝叶斯行动

这是最常见的混淆之一。MAP 是后验的众数,是一个「点」;贝叶斯行动是使后验期望损失最小的动作,依赖损失函数。两者只有在参数空间上取 0-1 损失时才重合。

更麻烦的是,MAP 对重参数化不协变。若 ϕ=h(θ)\phi=h(\theta) 是一个非线性变换,p(θ∣D)p(\theta\mid D) 的众数经过 hh 之后,一般不是 p(ϕ∣D)p(\phi\mid D) 的众数,因为密度会多出一个 Jacobian 因子。后验中位数在单调变换下等变,后验均值在仿射变换下等变,而 MAP 连这个性质都没有。这也是为什么在贝叶斯推断里,除非有明确的 0-1 损失,否则「取众数」并不是一个默认正确的选择。

4. 贝叶斯错误率

4.1 两类情形:全变差与重叠面积

设两类先验相等,类条件密度为 p0,p1p_0,p_1。由 3.2 节,贝叶斯规则在 η(x)=1/2\eta(x)=1/2 处切换,其中

η(x)=p1(x)p0(x)+p1(x).\eta(x)=\frac{p_1(x)}{p_0(x)+p_1(x)}.

于是

pe∗=∫min⁡{12p0(x), 12p1(x)}dx=12∫min⁡{p0(x),p1(x)} dx.\begin{aligned} p_e^* &=\int \min\left\{\tfrac12 p_0(x),\ \tfrac12 p_1(x)\right\}dx\\ &=\frac12\int \min\{p_0(x),p_1(x)\}\,dx . \end{aligned}

利用恒等式

min⁡{p0,p1}=p0+p1−∣p0−p1∣2,\min\{p_0,p_1\}=\frac{p_0+p_1-|p_0-p_1|}{2},

以及 ∫(p0+p1)/2=1\int(p_0+p_1)/2=1,得到

pe∗=12(1−12∫∣p0(x)−p1(x)∣ dx)=12(1−∥P0−P1∥TV).p_e^*=\frac12\left(1-\frac12\int|p_0(x)-p_1(x)|\,dx\right) =\frac12\big(1-\|P_0-P_1\|_{\mathrm{TV}}\big).

这里 ∥⋅∥TV\|\cdot\|_{\mathrm{TV}} 是全变差距离。公式的几何意义很直接:两个类条件分布的重叠面积越大,错误率越高;当两个分布完全重合时 pe∗=1/2p_e^*=1/2,当支撑集不相交时 pe∗=0p_e^*=0。

若先验不相等,公式改为

pe∗=∫min⁡{π0p0(x), π1p1(x)} dx,p_e^*=\int \min\{\pi_0p_0(x),\ \pi_1p_1(x)\}\,dx,

它仍然是一个「重叠面积」,只是两侧的权重不同。先验越偏向某一类,另一类的判定区域就越小,错误率的构成也随之改变。

4.2 高斯情形:马氏距离与 QQ 函数

设

X∣Y=0∼N(μ0,Σ),X∣Y=1∼N(μ1,Σ),X\mid Y=0\sim\mathcal N(\mu_0,\Sigma),\qquad X\mid Y=1\sim\mathcal N(\mu_1,\Sigma),

先验相等。对数似然比为

log⁡p1(x)p0(x)=(x−μ0)⊤Σ−1(μ1−μ0)−12(μ1+μ0)⊤Σ−1(μ1−μ0).\log\frac{p_1(x)}{p_0(x)} =(x-\mu_0)^\top\Sigma^{-1}(\mu_1-\mu_0) -\frac12(\mu_1+\mu_0)^\top\Sigma^{-1}(\mu_1-\mu_0).

这是一个 xx 的线性函数。定义马氏距离的平方

Δ2=(μ1−μ0)⊤Σ−1(μ1−μ0),\Delta^2=(\mu_1-\mu_0)^\top\Sigma^{-1}(\mu_1-\mu_0),

则在 Y=1Y=1 下,对数似然比服从 N(Δ2/2, Δ2)\mathcal N(\Delta^2/2,\ \Delta^2)。贝叶斯规则在似然比为零处切换,因此

pe∗=P(N(Δ2/2,Δ2)<0)=Φ(−Δ2)=Q(Δ2),p_e^*=P\big(\mathcal N(\Delta^2/2,\Delta^2)<0\big) =\Phi\left(-\frac{\Delta}{2}\right) =Q\left(\frac{\Delta}{2}\right),

其中 Φ\Phi 是标准正态分布函数,Q=1−ΦQ=1-\Phi。

几个数值可以帮助建立直觉:

Δ\Deltape∗=Q(Δ/2)p_e^*=Q(\Delta/2)Bhattacharyya 上界 12e−Δ2/8\tfrac12e^{-\Delta^2/8}
10.30850.4412
20.15870.3033
30.06680.1623
40.02280.0677

上界来自

pe∗≤12 BC(P0,P1)=12exp⁡(−Δ28),p_e^*\le \frac12\,\mathrm{BC}(P_0,P_1) =\frac12\exp\left(-\frac{\Delta^2}{8}\right),

其中 BC(P0,P1)=∫p0p1\mathrm{BC}(P_0,P_1)=\int\sqrt{p_0p_1} 是 Bhattacharyya 系数。当 Δ\Delta 很大时,

Q(Δ2)∼2Δ2πe−Δ2/8,Q\left(\frac{\Delta}{2}\right)\sim \frac{2}{\Delta\sqrt{2\pi}}e^{-\Delta^2/8},

与上界只差一个多项式因子。指数部分完全一致,这正是 Chernoff 界想要抓住的东西:错误率衰减的速率。

4.3 用散度给错误率定界

全变差是精确的,但它通常难以估计。实践中更常用的是 Hellinger 距离与 KL 散度。它们与全变差之间有一条标准的不等式链:

∥P0−P1∥TV≤H(P0,P1),H2(P0,P1)=1−BC(P0,P1),\|P_0-P_1\|_{\mathrm{TV}}\le H(P_0,P_1), \qquad H^2(P_0,P_1)=1-\mathrm{BC}(P_0,P_1),

以及

H2(P0,P1)≤12DKL(P0∥P1).H^2(P_0,P_1)\le \frac12 D_{\mathrm{KL}}(P_0\|P_1).

把第一条与 pe∗=12(1−TV)p_e^*=\tfrac12(1-\mathrm{TV}) 结合,得到

pe∗ ≥ 12(1−DKL(P0∥P1)2).p_e^*\ \ge\ \frac12\left(1-\sqrt{\frac{D_{\mathrm{KL}}(P_0\|P_1)}{2}}\right).

这是 Le Cam 两点法的核心不等式:如果两个分布很难区分(KL 很小),那么任何分类器的错误率都不可能低。注意方向——KL 给的是错误率的下界,不是上界。KL 很大并不意味着问题容易,只意味着这个下界变得松甚至变得平凡。

上界则需要另一个量:

pe∗≤12 BC(P0,P1).p_e^*\le \frac12\,\mathrm{BC}(P_0,P_1).

这是因为 ∫min⁡{p0,p1}≤∫p0p1\int\min\{p_0,p_1\}\le\int\sqrt{p_0p_1}(算术-几何平均不等式)。Bhattacharyya 系数可以直接计算的情形(例如高斯),这个上界就能给出正确的指数衰减率。

把两个方向放在一起,可以看到散度在决策问题里的角色:它们不是「错误率」,而是错误率的证书。全变差是唯一精确的那个,但它通常最难估计;Bhattacharyya 系数给出上界,KL 给出下界,Hellinger 则通过 TV≤H\mathrm{TV}\le H 把两边连起来。关于 f-散度的那篇文章,正是把这些证书放在同一个凸函数框架下;而贝叶斯错误率则告诉我们,这些证书最终约束的是什么。

4.4 错误率住在哪里

由 3.6 节的恒等式,只有满足 δ(X)≠δ∗(X)\delta(X)\ne\delta^*(X) 的样本才贡献错误,且权重是 ∣2η(X)−1∣|2\eta(X)-1|。这意味着贝叶斯错误率不是由决策边界的「长度」或「复杂度」决定的,而是由边界附近的后验质量决定的。

高维空间里,决策边界可以非常复杂,但如果后验在边界两侧都很尖锐(η\eta 接近 00 或 11),错误率仍然可以很低。反过来,一个简单的线性边界,如果穿过一大片 η≈1/2\eta\approx 1/2 的区域,错误率就会很高。这解释了为什么「边界复杂度」本身不是泛化误差的可靠指标:真正决定错误率的是后验在边界附近的分布,而不是边界的几何形状。

4.5 回归的贝叶斯风险

平方损失下,贝叶斯风险为

R∗=E[(Y−f∗(X))2]=E[Var(Y∣X)].R^*= \mathbb E\big[(Y-f^*(X))^2\big]=\mathbb E\big[\mathrm{Var}(Y\mid X)\big].

这是数据生成过程本身的噪声,任何模型都消不掉。对一个在训练集 DD 上学到的回归器 f^D\hat f_D,有标准的偏差-方差分解:

ED[(Y−f^D(X))2]=σε2⏟贝叶斯风险+(f∗(X)−EDf^D(X))2⏟偏差2+VarD(f^D(X))⏟方差.\mathbb E_D\big[(Y-\hat f_D(X))^2\big] = \\ \underbrace{\sigma_\varepsilon^2}_{\text{贝叶斯风险}} +\underbrace{\big(f^*(X)-\mathbb E_D\hat f_D(X)\big)^2}_{\text{偏差}^2} +\underbrace{\mathrm{Var}_D\big(\hat f_D(X)\big)}_{\text{方差}}.

第一项是问题的下限,后两项是估计带来的。模型太简单,偏差大;模型太复杂,方差大。贝叶斯风险给出了这条权衡曲线的起点:无论怎么调模型,都不可能把期望误差压到它下面。

5. 贝叶斯决策如何长成算法

5.1 生成式与判别式

贝叶斯规则只依赖后验概率 η(x)=p(Y=1∣x)\eta(x)=p(Y=1\mid x)。因此有两条路可以走:直接建模 η\eta,或者分别建模类条件密度 p(x∣y)p(x\mid y) 与先验 p(y)p(y),再用贝叶斯公式合成。前者叫判别式(discriminative),后者叫生成式(generative)。

生成式的合成公式是

η(x)=π1p1(x)π0p0(x)+π1p1(x).\eta(x)=\frac{\pi_1p_1(x)}{\pi_0p_0(x)+\pi_1p_1(x)}.

只要两边的密度估计足够准,这个式子就给出贝叶斯后验。线性判别分析(LDA)假设两类共享协方差 Σ\Sigma,此时对数似然比是 xx 的线性函数:

log⁡p1(x)p0(x)=(μ1−μ0)⊤Σ−1x+常数,\log\frac{p_1(x)}{p_0(x)} =(\mu_1-\mu_0)^\top\Sigma^{-1}x+\text{常数},

决策边界是一个超平面。二次判别分析(QDA)允许两类有不同的协方差 Σ0,Σ1\Sigma_0,\Sigma_1,对数似然比出现 x⊤(Σ1−1−Σ0−1)xx^\top(\Sigma_1^{-1}-\Sigma_0^{-1})x 这样的二次项,边界变成二次曲面。朴素贝叶斯再加一层条件独立假设,把 p(x∣y)p(x\mid y) 拆成各个特征的乘积,用很大的偏差换取很小的方差。

两条路的目标是同一个 η\eta,但统计性质不同。生成式利用了 p(x∣y)p(x\mid y) 的结构,在模型正确时样本效率更高;判别式直接优化决策边界,在模型误设时通常更稳健。这不是「谁更好」的问题,而是「你愿意把假设放在哪一侧」的问题。

这里说的是 oracle 贝叶斯规则:它假设 η\eta 已知。实践中我们把估计的 η^\hat\eta 代入,得到的只是 plug-in 分类器;文献里的「贝叶斯分类器」有时指这个 oracle 规则,有时指朴素贝叶斯之类的具体算法,读的时候要分清。plug-in 分类器的风险与 R∗R^* 之间的差距,正是第 7 节要分解的估计误差。

值得一提的是,朴素贝叶斯的条件独立假设几乎总是错的,但它的分类性能常常不差。原因在于 0-1 损失只关心 η(x)\eta(x) 是否越过 1/21/2,而不要求 η\eta 的数值完全准确;只要假设带来的偏差没有把决策边界推得太远,分类结果就可以是对的。这个现象也提醒我们:概率估计的质量与决策的质量不是同一件事。

这个区别在损失函数上有一个精确的对应。0-1 损失只关心 η(x)\eta(x) 是否越过阈值,因此它对概率的微小误差不敏感;交叉熵(对数损失)则是一个严格 proper scoring rule,它的总体最优解恰好是真实的条件分布 p(y∣x)p(y\mid x):

E[−log⁡q(Y∣X)]=E[H(p(⋅∣X))]+E[KL(p(⋅∣X)∥q(⋅∣X))],\mathbb E[-\log q(Y\mid X)] =\mathbb E\big[\mathbb H(p(\cdot\mid X))\big] +\mathbb E\big[\mathrm{KL}(p(\cdot\mid X)\|q(\cdot\mid X))\big],

右边第一项与 qq 无关,第二项在 q=p(⋅∣X)q=p(\cdot\mid X) 时为零。所以,在总体意义上,用交叉熵训练模型是在要求它给出校准的概率;用 0-1 损失训练,只要求它给出正确的分类。两者目标不同,不能互相替代。

5.2 正则化就是先验

设参数化模型 p(D∣θ)=∏ip(yi∣xi,θ)p(D\mid\theta)=\prod_i p(y_i\mid x_i,\theta),噪声为高斯,方差 σ2\sigma^2。则负对数似然正比于平方误差:

−log⁡p(D∣θ)=12σ2∑i(yi−fθ(xi))2+常数.-\log p(D\mid\theta) =\frac{1}{2\sigma^2}\sum_i\big(y_i-f_\theta(x_i)\big)^2+\text{常数}.

若给参数加上高斯先验 θ∼N(0,τ2I)\theta\sim\mathcal N(0,\tau^2I),则

−log⁡p(θ)=12τ2∥θ∥2+常数.-\log p(\theta)=\frac{1}{2\tau^2}\|\theta\|^2+\text{常数}.

MAP 估计为

θ^MAP=arg⁡min⁡θ ∑i(yi−fθ(xi))2+σ2τ2∥θ∥2.\hat\theta_{\mathrm{MAP}} =\arg\min_\theta\ \sum_i\big(y_i-f_\theta(x_i)\big)^2 +\frac{\sigma^2}{\tau^2}\|\theta\|^2 .

这就是岭回归,也是神经网络里的权重衰减。若把目标写成均方误差,系数变为 σ2/(nτ2)\sigma^2/(n\tau^2)。若先验换成 Laplace 分布 p(θj)=12be−∣θj∣/bp(\theta_j)=\frac{1}{2b}e^{-|\theta_j|/b},负对数先验正比于 ∥θ∥1\|\theta\|_1,MAP 就变成 Lasso。

先验负对数先验等价惩罚
均匀先验常数无(MLE)
高斯 N(0,τ2)\mathcal N(0,\tau^2)∥θ∥2/(2τ2)\|\theta\|^2/(2\tau^2)岭回归、权重衰减
Laplace (0,b)(0,b)∥θ∥1/b\|\theta\|_1/bLasso
稀疏诱导先验非凸惩罚特征选择

这张表解释了一个常见的困惑:为什么正则化看起来像是「人为加上去的技巧」,却总是有效。它其实是在补充一个先验:参数不应该太大,或者不应该同时有很多非零分量。正则项不是外来的约束,而是模型假设的一部分,只是以惩罚的形式出现。

需要谨慎的是,深度网络里的权重衰减与高斯先验的对应只在全批量 MAP 的意义上成立。SGD 的噪声、批量大小、参数化方式(例如 BatchNorm 的缩放不变性)都会改变有效的先验。把它当成一个有用的类比,而不是一个严格的等价。

5.3 似然比、ROC 与代价

二分类的贝叶斯规则可以写成似然比检验:

δ(x)=1⟺p1(x)p0(x)≥π0 cFPπ1 cFN.\delta(x)=1 \quad\Longleftrightarrow\quad \frac{p_1(x)}{p_0(x)} \ge \frac{\pi_0\,c_{\mathrm{FP}}}{\pi_1\,c_{\mathrm{FN}}}.

左边是数据提供的证据,右边是先验与代价决定的阈值。先验偏向负类,阈值升高;漏诊更贵,阈值降低。这个分解把「数据说什么」与「我们更怕哪种错误」分开了。

把阈值从低到高扫一遍,就得到 ROC 曲线。设阈值为 tt,记真正率为 TPR(t)\mathrm{TPR}(t)、假正率为 FPR(t)\mathrm{FPR}(t),则期望代价为

C(t)=π0 cFP FPR(t)+π1 cFN (1−TPR(t)).C(t)=\pi_0\,c_{\mathrm{FP}}\,\mathrm{FPR}(t) +\pi_1\,c_{\mathrm{FN}}\,\big(1-\mathrm{TPR}(t)\big).

贝叶斯决策做的不是「选一个看起来最好的点」,而是根据先验与代价,在 ROC 上选一个使这个线性函数最小的点。等代价线的斜率是

d TPRd FPR=π0 cFPπ1 cFN.\frac{d\,\mathrm{TPR}}{d\,\mathrm{FPR}} =\frac{\pi_0\,c_{\mathrm{FP}}}{\pi_1\,c_{\mathrm{FN}}}.

这也说明 AUC 只是一个与阈值无关的排序指标。两个模型的 AUC 可以一个高一个低,但在某个特定的代价结构下,AUC 低的模型反而更优。脱离代价谈「最好的分类器」是没有意义的。

5.4 拒绝选项与选择性分类

3.5 节已经导出了 Chow 规则:当最大后验低于 1−c1-c 时拒绝。把这句话放到整个数据集上,就得到一条风险-覆盖率曲线:随着拒绝阈值升高,被拒绝的样本增多,覆盖率下降,但在被接受的样本上错误率也下降。

选择性分类的理论目标不是最大化准确率,而是最小化

选择性风险=被接受样本上的错误数被接受样本数,\text{选择性风险}=\frac{\text{被接受样本上的错误数}}{\text{被接受样本数}},

同时满足覆盖率约束。这两个量之间的权衡,仍然由损失函数决定:拒绝的代价 cc 越小,模型越倾向于拒绝;cc 越大,模型越倾向于冒险。

在医疗、自动驾驶、内容审核这些场景里,拒绝选项往往比「提高模型容量」更有效。一个知道自己不知道的模型,比一个在所有样本上都自信地犯错的模型更有用。

5.5 类别不平衡:改损失,还是改数据

类别不平衡经常被当成一个数据处理问题,但在决策论的框架里,它首先是一个先验与损失的问题。

设原分布的类先验为 π0,π1\pi_0,\pi_1,重采样后的先验为 π0′,π1′\pi'_0,\pi'_1。同一个似然比 p1(x)/p0(x)p_1(x)/p_0(x) 在两种分布下给出的后验几率不同:

porig(1∣x)porig(0∣x)=pres(1∣x)pres(0∣x)⋅π1 π0′π1′ π0.\frac{p_{\text{orig}}(1\mid x)}{p_{\text{orig}}(0\mid x)} = \frac{p_{\text{res}}(1\mid x)}{p_{\text{res}}(0\mid x)} \cdot \frac{\pi_1\,\pi'_0}{\pi'_1\,\pi_0}.

如果重采样之后不把这个因子校正回去,决策阈值就偏了。反过来,如果不想重采样,也可以直接在原分布上训练,然后按 3.5 节的公式移动阈值。两条路在理论上等价,差别只在实现难度与估计方差。

真正不能省掉的是损失函数。若两类错误的代价确实不同,那么无论怎么采样,都必须把代价写进决策规则;否则得到的只是一个在改动后的分布上「看起来平衡」的分类器,而不是原问题的最优决策。

6. 从参数到分布

6.1 三个预测器

给定训练集 DD,预测 YY 至少有三层做法:

  1. 代入式:先估计一个参数 θ^\hat\theta,再用 p(y∣x,θ^)p(y\mid x,\hat\theta) 预测;
  2. 后验预测:对参数积分, pˉ(y∣x)=∫p(y∣x,θ) p(θ∣D) dθ;\bar p(y\mid x)=\int p(y\mid x,\theta)\,p(\theta\mid D)\,d\theta;
  3. 贝叶斯行动:在预测分布上再做一次决策, a∗(x)=arg⁡min⁡a∫L(y,a) pˉ(y∣x) dy.a^*(x)=\arg\min_a\int L(y,a)\,\bar p(y\mid x)\,dy.

第二层是贝叶斯学习的核心,第三层才是完整的贝叶斯决策。很多讨论停在第二层,把「后验预测」直接当成答案;但在有明确损失的问题里,真正的答案是第三层。

后验预测是一个混合分布,不是「把后验均值代进去」。由于 p(y∣x,θ)p(y\mid x,\theta) 通常对 θ\theta 非线性,

∫p(y∣x,θ) p(θ∣D) dθ≠p(y∣x,E[θ∣D]).\int p(y\mid x,\theta)\,p(\theta\mid D)\,d\theta \ne p\big(y\mid x,\mathbb E[\theta\mid D]\big).

差别在不确定性大的区域尤其明显。设想后验在两个相距很远的参数值上各有一个峰,代入均值会得到一个位于两峰之间、几乎没有概率质量的预测;后验预测则会把两个峰的概率都保留下来。前者过度自信,后者诚实地表达了「我不知道是哪一个」。

6.2 两种不确定性

预测的不确定性可以分成两类。偶然不确定性(aleatoric uncertainty)来自数据生成过程本身的噪声:同样的 xx 可能对应不同的 yy,即使知道真实参数也无法消除。认知不确定性(epistemic uncertainty)来自知识不足:我们不知道哪个参数是对的,但随着数据增多,它会收缩。

对平方损失,这两种不确定性由全方差公式精确分解:

Var(Y∣x,D)=Eθ∣D[Var(Y∣x,θ)]⏟偶然+Varθ∣D(E[Y∣x,θ])⏟认知.\mathrm{Var}(Y\mid x,D) = \underbrace{\mathbb E_{\theta\mid D}\big[\mathrm{Var}(Y\mid x,\theta)\big]}_{\text{偶然}} + \underbrace{\mathrm{Var}_{\theta\mid D}\big(\mathbb E[Y\mid x,\theta]\big)}_{\text{认知}}.

第一项是「即使知道模型参数,预测仍然会波动的部分」;第二项是「不同参数给出的预测均值之间的差异」。在数据稀疏的区域,后验宽,第二项大;随着数据增加,后验收缩,第二项趋于零。第一项则不会消失。

这个分解也是主动学习的理论基础:我们想选的样本,是那些能最大程度降低认知不确定性的样本,而不是那些本身就带噪声的样本。

6.3 推断本身也是决策

在参数空间上定义一个损失,报一个点估计,本质上也是一次决策:

参数空间上的损失最优估计
平方损失后验均值
绝对损失后验中位数
0-1 损失后验众数(MAP)

所以贝叶斯估计不是与决策论并列的东西,而是决策论在参数空间上的应用。后验均值与 MAP 会打架,是因为它们对应不同的损失,没有谁更「正确」。若你关心的是平方误差,均值是对的;若你关心的是「最可能的那个值」,众数是对的;若你关心的是「误差不超过某个范围」,分位数才是对的。

同样的逻辑适用于区间估计。一个可信区间是一个「集合值」的行动,它的损失函数衡量区间是否覆盖真值、以及区间有多宽。置信区间与可信区间的差异当然有概率哲学的成分;但从决策论看,它们也对应不同的损失:一个关心重复抽样下的覆盖率,一个关心给定数据下的后验覆盖概率。

6.4 后验的集中

在正则条件下,后验会随着数据增多而集中到真实参数附近。Bernstein–von Mises 定理给出了更精确的结论:后验在适当的缩放与中心化之后,趋近于一个正态分布,

n (θ−θ^n)∣D ⇒ N(0, I(θ0)−1),\sqrt n\,(\theta-\hat\theta_n)\mid D \ \Rightarrow\ \mathcal N\big(0,\ I(\theta_0)^{-1}\big),

其中 I(θ0)I(\theta_0) 是 Fisher 信息。等价地,后验近似为

p(θ∣D)≈N(θ^n, 1nI(θ0)−1).p(\theta\mid D)\approx \mathcal N\left(\hat\theta_n,\ \frac{1}{n}I(\theta_0)^{-1}\right).

Doob 定理则从另一个角度保证一致性:对先验几乎所有的参数值,后验都会收敛到真值附近的任意邻域。注意「先验几乎处处」这个限定:先验可能把质量放在一个坏集合上,只是那个集合的概率为零。

这两个定理解释了为什么在大样本下,贝叶斯与频率派的答案往往一致。后验的宽度以 O(1/n)O(1/\sqrt n) 收缩,先验对形状的影响随之减弱,后验主要由似然决定。但这也划出了边界:当模型奇异(参数不可识别、网络层数过多、对称性导致 Fisher 信息退化)时,Bernstein–von Mises 失效,后验的形状本身携带信息,不能用一个高斯替代。这正是奇异学习理论要处理的情形。

6.5 边缘似然与奥卡姆剃刀

贝叶斯模型比较的核心量是边缘似然(marginal likelihood),也就是贝叶斯公式的分母:

p(D∣M)=∫p(D∣θ,M) p(θ∣M) dθ.p(D\mid M)=\int p(D\mid\theta,M)\,p(\theta\mid M)\,d\theta .

两个模型的比较用贝叶斯因子:

BF12=p(D∣M1)p(D∣M2),p(M1∣D)p(M2∣D)=BF12⋅p(M1)p(M2).\mathrm{BF}_{12}=\frac{p(D\mid M_1)}{p(D\mid M_2)}, \qquad \frac{p(M_1\mid D)}{p(M_2\mid D)} =\mathrm{BF}_{12}\cdot\frac{p(M_1)}{p(M_2)}.

边缘似然自动带有奥卡姆惩罚。用 Laplace 近似可以看出这一点。设 θ^\hat\theta 是后验众数,H=−∇2log⁡[p(D∣θ)p(θ∣M)]∣θ^H=-\nabla^2\log[p(D\mid\theta)p(\theta\mid M)]\big|_{\hat\theta},则

p(D∣M)≈p(D∣θ^) p(θ^∣M) (2π)d/2∣H∣−1/2.p(D\mid M) \approx p(D\mid\hat\theta)\,p(\hat\theta\mid M)\,(2\pi)^{d/2}|H|^{-1/2}.

在正则模型里 H≈nI(θ^)H\approx nI(\hat\theta),取对数后

log⁡p(D∣M)≈log⁡p(D∣θ^)−d2log⁡n+O(1).\log p(D\mid M) \approx \log p(D\mid\hat\theta)-\frac d2\log n+O(1).

于是

BIC=−2log⁡p(D∣θ^)+dlog⁡n\mathrm{BIC}=-2\log p(D\mid\hat\theta)+d\log n

正是 −2log⁡p(D∣M)-2\log p(D\mid M) 的主项。参数越多,惩罚越重;但更准确地说,惩罚的不是维度本身,而是参数空间中「能给出不同预测」的有效体积。一个参数如果对预测没有影响,它就不会贡献体积,也就不会被惩罚。

边缘似然也有它自己的问题:它对先验敏感。一个过分分散的先验会把概率质量摊到大量「几乎不预测数据」的区域,导致模型被过度惩罚。这就是 Lindley 悖论的一个侧面:当原假设是一个点、备择假设是一个很宽的分布时,贝叶斯因子可能强烈支持原假设,即使频率派的检验已经拒绝了它。贝叶斯方法不会自动给出「正确」的模型选择,它只是把选择的标准写清楚了。

7. 机器学习理论中的位置

7.1 风险分解

设 H\mathcal H 是模型类,RH=inf⁡h∈HR(h)R_{\mathcal H}=\inf_{h\in\mathcal H}R(h),h^\hat h 是从数据中学到的模型。则

R(h^)−R∗=(RH−R∗)⏟近似误差+(R(h^)−RH)⏟估计误差.R(\hat h)-R^* = \underbrace{\big(R_{\mathcal H}-R^*\big)}_{\text{近似误差}} + \underbrace{\big(R(\hat h)-R_{\mathcal H}\big)}_{\text{估计误差}}.

第一项由模型类决定:模型类太小,无论如何训练都逼近不了贝叶斯风险。第二项由样本量与模型复杂度决定:模型类太大,有限样本上的经验风险会偏离真实风险。

贝叶斯风险 R∗R^* 是这两项的共同起点。它既不属于模型,也不属于算法,而是问题本身的下限。学习理论的大部分内容,都可以理解为在给定 R∗R^* 的情况下,分别研究这两项能压到多低。

这也解释了为什么「更复杂的模型」有时有用、有时有害:它降低近似误差,却可能提高估计误差。深度学习成功的部分原因,是它在大规模数据上把估计误差压得足够低,同时用极宽的模型类把近似误差压得很小。这不是对偏差-方差权衡的否定,而是权衡曲线在数据量增大后的重新定位。

7.2 一致性

一个决策规则序列 δn\delta_n 称为普遍一致的(universally consistent),若对任意分布 PP,

R(δn) → P  R∗R(\delta_n)\ \xrightarrow{\ P\ }\ R^*

依概率成立。注意这里收敛的是风险,不是后验估计。一个模型可以是一致的,但给出的概率完全不校准;反过来,一个校准很好的模型也可能因为模型类太小而无法逼近贝叶斯风险。

kk 近邻是最经典的例子。Stone 定理说:若 k→∞k\to\infty 且 k/n→0k/n\to0,则 kk-NN 是普遍一致的。11-NN 更极端:它的渐近错误率不超过 2R∗2R^*。这个「二倍」的常数很说明问题——即使是最简单的局部方法,也能在渐近意义下接近贝叶斯风险,代价是收敛速度很慢,而且在高维里被维数灾难拖垮。

一致性是一个渐近性质,它不告诉你有限样本下会发生什么。一个普遍一致的估计器可能在 n=106n=10^6 时仍然表现很差;一个不一致的估计器可能在某个特定问题上表现得很好。贝叶斯风险给出的是终点,不是路线。

7.3 PAC-Bayes:把「信念」放进泛化界

PAC-Bayes 把贝叶斯式的分布对象与频率式的泛化保证放在同一个不等式里。这里的 P,QP,Q 是假设空间上的分布,不要与数据分布混淆。

定理 7.1(McAllester 界) 设先验 PP 在看到数据前固定,损失取值于 [0,1][0,1]。对任意 δ∈(0,1)\delta\in(0,1),以至少 1−δ1-\delta 的概率,对所有后验分布 QQ 同时成立:

Eh∼Q[R(h)]≤Eh∼Q[R^n(h)]+KL(Q∥P)+log⁡(2n/δ)2n,\mathbb E_{h\sim Q}\big[R(h)\big] \le \mathbb E_{h\sim Q}\big[\hat R_n(h)\big] + \sqrt{\frac{\mathrm{KL}(Q\|P)+\log(2\sqrt n/\delta)}{2n}},

其中 R^n\hat R_n 是经验风险。

这个界的美妙之处在于:PP 必须在看到数据前固定,QQ 可以是任何依赖于数据的分布。KL 项衡量「为了解释数据,你把信念从先验挪了多远」。先验越集中、后验越贴近先验,界就越紧;后验越远离先验,界就越松。

这正是奥卡姆剃刀的概率版本,只不过它作用在分布而不是单个假设上。一个在大量假设上分散的 QQ 会付出很大的 KL,即使它的经验风险很低;一个集中在少数假设上的 QQ 则可以用较小的 KL 换取较好的界。

PAC-Bayes 也有它的局限:界在深度网络里常常很松,常数项可能大到没有实际意义。但它提供了一个正确的概念框架:泛化不是「参数少」的同义词,而是「后验没有比先验复杂太多」。后来的 flatness、sharpness、信息论泛化界,都可以在这个框架里找到位置。

7.4 变分推断:当后验收敛不了

精确后验通常不可行,因为归一化常数 p(D)p(D) 是一个高维积分。变分推断用一个受限制的分布族 Q\mathcal Q 去逼近后验,目标是

q∗=arg⁡min⁡q∈Q KL(q(θ) ∥ p(θ∣D)).q^*=\arg\min_{q\in\mathcal Q}\ \mathrm{KL}\big(q(\theta)\ \|\ p(\theta\mid D)\big).

直接最小化这个 KL 仍然需要 p(D)p(D),于是利用恒等式

log⁡p(D)=Eq[log⁡p(D,θ)q(θ)]⏟ELBO(q)+KL(q(θ) ∥ p(θ∣D)),\log p(D) = \underbrace{\mathbb E_{q}\left[\log\frac{p(D,\theta)}{q(\theta)}\right]}_{\mathrm{ELBO}(q)} + \mathrm{KL}\big(q(\theta)\ \|\ p(\theta\mid D)\big),

把问题转化为最大化证据下界(ELBO):

ELBO(q)=Eq[log⁡p(D∣θ)]−KL(q(θ) ∥ p(θ)).\mathrm{ELBO}(q) = \mathbb E_q\big[\log p(D\mid\theta)\big] - \mathrm{KL}\big(q(\theta)\ \|\ p(\theta)\big).

第一项鼓励 qq 解释数据,第二项惩罚 qq 偏离先验。最大化 ELBO 等价于最小化 KL(q∥p(θ∣D))\mathrm{KL}(q\|p(\theta\mid D)),因为 log⁡p(D)\log p(D) 与 qq 无关。

KL 的方向很重要。KL(q∥p)\mathrm{KL}(q\|p) 倾向寻模(mode-seeking):只要 qq 覆盖了后验的一个峰,它就不太在意漏掉其他峰,因此常常低估方差。反方向 KL(p∥q)\mathrm{KL}(p\|q) 倾向覆盖质量(mass-covering),但它的计算需要 p(θ∣D)p(\theta\mid D) 的归一化常数。变分推断通常用前者,因为 ELBO 只需要未归一化的联合分布。

从决策论的角度看,变分推断本身也是一个决策问题:在受限的分布族里,选一个使 KL 最小的近似。这里的「损失」是 KL,行动是一个分布,而不是一个点。贝叶斯决策理论的框架在这里再次出现,只是行动空间从 A\mathcal A 变成了 Q\mathcal Q。

7.5 贝叶斯神经网络与不确定性

把参数换成网络权重 ww,后验 p(w∣D)p(w\mid D) 就是贝叶斯神经网络。精确后验不可行,实践中用 MC dropout、Laplace 近似、变分推断或深度集成来近似。

对分类问题,预测熵可以分解为

H[Y∣x,D]=Ew∣D[H[Y∣x,w]]+I(Y;w∣x,D).\mathbb H\big[Y\mid x,D\big] = \mathbb E_{w\mid D}\big[\mathbb H[Y\mid x,w]\big] + \mathbb I\big(Y;w\mid x,D\big).

第一项是偶然不确定性:即使固定权重,类别仍然有噪声。第二项是认知不确定性:不同权重给出的预测分布之间的差异。BALD 等主动学习准则,就是用它来挑选「最能减少认知不确定性」的样本。

但「贝叶斯」这个标签不保证校准。一个近似的后验可能给出过度自信的预测;深度集成不是严格的贝叶斯方法,却常常校准得更好。判断的标准应该是决策质量与校准,而不是方法的名字。一个模型如果能在高损失的区域拒绝预测,或者把不确定性传递给下游决策,它就已经在做贝叶斯决策了,哪怕它没有显式地计算后验。

7.6 序贯决策:后验采样与探索

在序贯决策里,后验不只是用来「估计」,还用来「探索」。Thompson 采样的规则很简单:每一步从后验中采样一个环境参数 θ~∼p(θ∣Dt)\tilde\theta\sim p(\theta\mid D_t),然后按这个采样出的环境做最优动作。

这个规则自动平衡了利用与探索。不确定性大的动作在采样中更可能显得最优,因此会被更频繁地尝试;随着数据增多,后验收窄,采样出的环境越来越接近真实环境,策略逐渐转向利用。它的贝叶斯遗憾有理论保证,而且实现起来只需要一个后验采样器。

从决策论的角度看,Thompson 采样把「信息」当成了决策的一部分:当前动作的收益不只是它带来的即时回报,还包括它揭示的信息。完整的贝叶斯解决方案是动态规划,在信念空间(belief space)上求解;Thompson 采样是它的一个可计算的近似。

7.7 主动学习与贝叶斯优化

主动学习问的是:下一个标签应该买哪一个?从决策论看,这是一个两级决策:先选一个样本,再在拿到标签后做预测。理想的采集函数是「期望后验风险降低」:

xnext=arg⁡min⁡x Ey∼p(y∣x,D)[R∗(D∪{(x,y)})],x_{\text{next}} = \arg\min_{x}\ \mathbb E_{y\sim p(y\mid x,D)} \left[ R^*(D\cup\{(x,y)\}) \right],

其中 R∗(D)R^*(D) 是当前数据下的贝叶斯风险。这个量通常算不出来,于是有了各种近似:熵、互信息、BALD、期望模型变化。

贝叶斯优化把目标函数本身当成随机变量,用后验选择下一个评估点。Expected Improvement 就是

EI(x)=E[max⁡{0, fmin⁡−f(x)}∣D],\mathrm{EI}(x)=\mathbb E\big[\max\{0,\ f_{\min}-f(x)\}\mid D\big],

它把「可能变好多少」与「变好的概率」乘在一起。探索与利用的权衡,在这里被写成了一个一维积分的最大化问题。

7.8 可容许性与完整类定理

最后回到 2.4 节留下的那个说法。一个决策规则 δ\delta 称为可容许的(admissible),若不存在另一个规则 δ′\delta',使得

R(δ′)≤R(δ)对所有分布 P,R(\delta')\le R(\delta)\quad\text{对所有分布 }P,

且至少对一个 PP 严格小于。换句话说,可容许就是「没有被另一个规则全面支配」。一个处处比别的规则差的决策规则,理应被淘汰;可容许性把这个直觉形式化了。

定理 7.2(Wald 完整类定理) 在损失函数凸、行动空间紧、实验族满足适当正则条件(例如被某个 σ\sigma-有限测度支配)时,每一个可容许的决策规则都是某个先验(或先验序列的极限)下的贝叶斯规则。

这个定理的意义在于:在适当条件下,它把「频率派的极大极小规则」「看起来没有先验的规则」也纳入了贝叶斯框架。先验不再是一个可选的信仰,而是可容许性这个要求所对应的对偶变量。你可以拒绝某个具体的先验,但很难拒绝「存在某个先验」这件事。

定理的条件不能省。损失不凸、行动空间不紧,或者问题是非参数的,结论可能失效;此时需要广义先验、极限贝叶斯规则,或者放弃完整类定理的结论。这些技术细节不改变核心信息:贝叶斯规则不是许多决策规则中的一种,而是所有不被支配的决策规则的一般形式。

从实践的角度看,这个定理不告诉你该选哪个先验,只告诉你:如果你想要一个不被全面支配的规则,你总可以把它写成某个先验下的贝叶斯规则。先验的选择仍然是建模的一部分,但「要不要用先验」已经不再是一个真正的选项。

8. 几个可以手算的例子

8.1 Beta–Bernoulli 的完整更新

取先验 Beta(3,3)\mathrm{Beta}(3,3),观测到 1010 次投掷中 77 次正面。后验为

Beta(3+7, 3+3)=Beta(10,6).\mathrm{Beta}(3+7,\ 3+3)=\mathrm{Beta}(10,6).

把几个估计量放在一起:

量值
先验均值3/6=0.53/6=0.5
最大似然估计7/10=0.77/10=0.7
后验均值10/16=0.62510/16=0.625
后验众数9/14≈0.6439/14\approx0.643
后验预测 p(xnew=1∣D)p(x_{\text{new}}=1\mid D)0.6250.625

后验均值被先验从 0.70.7 拉回到 0.6250.625,这就是收缩。数据越多,拉力越弱:若观测到 100100 次中 7070 次正面,后验均值变为 73/106≈0.68973/106\approx0.689,已经非常接近 0.70.7。

如果决策是「猜下一次是否为正面」,且两类错误的代价相同,那么只要后验预测大于 1/21/2 就猜正面,0.6250.625 满足条件。若两类错误的代价不同:把反面误判为正面的代价为 22,把正面误判为反面的代价为 11,则阈值变为

cFPcFP+cFN=22+1=23,\frac{c_{\mathrm{FP}}}{c_{\mathrm{FP}}+c_{\mathrm{FN}}} =\frac{2}{2+1}=\frac23,

而 0.625<2/30.625<2/3,最优行动变成猜反面。同一个后验,换一个损失,结论就变了。这正是贝叶斯决策理论与「最大后验」的分水岭:后验只有一个,行动却取决于代价。

8.2 高斯先验与收缩

设观测模型 yi∼N(θ,σ2)y_i\sim\mathcal N(\theta,\sigma^2),先验 θ∼N(μ0,τ2)\theta\sim\mathcal N(\mu_0,\tau^2),样本均值为 xˉ\bar x,样本量为 nn。后验仍是高斯:

θ∣D∼N(μn,τn2),\theta\mid D\sim\mathcal N(\mu_n,\tau_n^2),

其中

μn=μ0/τ2+nxˉ/σ21/τ2+n/σ2,1τn2=1τ2+nσ2.\mu_n= \frac{\mu_0/\tau^2+n\bar x/\sigma^2}{1/\tau^2+n/\sigma^2}, \qquad \frac{1}{\tau_n^2}=\frac{1}{\tau^2}+\frac{n}{\sigma^2}.

把 μn\mu_n 写成加权平均:

μn=σ2/nτ2+σ2/n μ0+τ2τ2+σ2/n xˉ.\mu_n = \frac{\sigma^2/n}{\tau^2+\sigma^2/n}\,\mu_0 + \frac{\tau^2}{\tau^2+\sigma^2/n}\,\bar x .

第二项前面的系数就是收缩因子。取 μ0=0\mu_0=0、τ=1\tau=1、σ=1\sigma=1、xˉ=1\bar x=1:

nn后验均值 μn\mu_n后验方差 τn2\tau_n^2
11/2=0.51/2=0.51/2=0.51/2=0.5
44/5=0.84/5=0.81/5=0.21/5=0.2
2525/26≈0.96225/26\approx0.9621/26≈0.0381/26\approx0.038
100100/101≈0.990100/101\approx0.9901/101≈0.00991/101\approx0.0099

样本越多,后验越靠近样本均值,方差越小。先验的作用不是「永远拉住估计」,而是在小样本时提供一个稳定的起点,随着证据积累自动退场。

8.3 医疗筛查:先验、似然与代价

设某疾病的患病率为 1%1\%,检测的灵敏度为 99%99\%,特异度为 95%95\%。一个人检测阳性,患病概率是多少?

用贝叶斯公式:

P(D∣+)=P(+∣D)P(D)P(+∣D)P(D)+P(+∣¬D)P(¬D)=0.99×0.010.99×0.01+0.05×0.99=0.00990.0099+0.0495≈0.167.\begin{aligned} P(D\mid +) &=\frac{P(+\mid D)P(D)}{P(+\mid D)P(D)+P(+\mid \neg D)P(\neg D)}\\ &=\frac{0.99\times0.01}{0.99\times0.01+0.05\times0.99}\\ &=\frac{0.0099}{0.0099+0.0495} \approx0.167 . \end{aligned}

尽管检测看起来很准,阳性者真正患病的概率只有 16.7%16.7\%,因为健康人占了绝大多数,5%5\% 的假阳性率乘以 99%99\% 的基础人数,产生了一个比真阳性更大的绝对数量。这就是基础率谬误:把 P(+∣D)P(+\mid D) 与 P(D∣+)P(D\mid +) 混为一谈。

患病率变化时,同一个检测给出的后验完全不同:

患病率P(D∣+)P(D\mid +)
0.5%0.5\%0.0900.090
1%1\%0.1670.167
5%5\%0.5100.510
10%10\%0.6880.688

现在加入代价。设漏诊的代价是误诊的 1010 倍,即 cFN=10c_{\mathrm{FN}}=10、cFP=1c_{\mathrm{FP}}=1,则治疗阈值是

cFPcFP+cFN=111≈0.091.\frac{c_{\mathrm{FP}}}{c_{\mathrm{FP}}+c_{\mathrm{FN}}} =\frac{1}{11}\approx0.091.

在患病率 1%1\% 时,P(D∣+)≈0.167>0.091P(D\mid +)\approx0.167>0.091,应当治疗;在患病率 0.5%0.5\% 时,P(D∣+)≈0.090P(D\mid +)\approx0.090,已经落到阈值以下,治疗与不治疗的期望代价几乎相同。同一个阳性结果,在不同的先验与代价下导向不同的行动,这不是检测本身变了,而是决策问题变了。

这个例子也说明了为什么「准确率」不是一个好的汇报指标。检测的准确率是

0.99×0.01+0.95×0.99≈0.9504,0.99\times0.01+0.95\times0.99\approx0.9504,

看起来很高,但它几乎完全由健康人贡献;在真正需要决策的那一小部分人身上,后验概率才是关键。

8.4 高斯两类问题

设两类共享方差 σ2=1\sigma^2=1,均值分别为 00 与 22,先验相等。则 Δ=2\Delta=2,贝叶斯错误率为

pe∗=Q(1)≈0.1587.p_e^*=Q(1)\approx0.1587.

若把均值差扩大到 44,Δ=4\Delta=4,错误率降为 Q(2)≈0.0228Q(2)\approx0.0228。错误率随马氏距离以 e−Δ2/8e^{-\Delta^2/8} 的速率衰减,但永远不为零,因为两个高斯的支撑集在整个实轴上重叠。

Bhattacharyya 上界给出

pe∗≤12e−Δ2/8,p_e^*\le\frac12e^{-\Delta^2/8},

在 Δ=2\Delta=2 时为 0.3030.303,在 Δ=4\Delta=4 时为 0.06770.0677。上界不紧,但它的指数部分与真实错误率一致;这正是我们在 4.2 节强调的「速率」而非「常数」。

结语

把全文的线索压缩成一行:

p(θ)⏟先验×p(D∣θ)⏟似然 ⟶ p(θ∣D)⏟后验 ⟶ arg⁡min⁡a E[L(Y,a)∣D]⏟行动 ⟶ R∗⏟下限.\underbrace{p(\theta)}_{\text{先验}} \times \underbrace{p(D\mid\theta)}_{\text{似然}} \ \\ \longrightarrow\ \underbrace{p(\theta\mid D)}_{\text{后验}} \ \\ \longrightarrow\ \underbrace{\arg\min_a\ \mathbb E[L(Y,a)\mid D]}_{\text{行动}} \ \\ \longrightarrow\ \underbrace{R^*}_{\text{下限}}.

这条链的每一段都可以独立展开,也可以独立质疑。它有两个层次的角色。

它是规范性框架。它把「假设」「证据」「偏好」「行动」四件事分开:先验与似然负责描述世界,损失负责描述我们想要什么,决策规则负责把两者接起来。很多机器学习里的争论,本质上是在争论这四个部分中的某一个,却常常被包装成「模型好不好」的单一问题。它也是理论基准。贝叶斯风险给出了任何算法都无法突破的下限,超额风险分解把模型误差拆成近似误差与估计误差,贝叶斯错误率把分布差异翻译成可分性。没有这个基准,「泛化」就只是一个没有坐标的词。

最后一点也许最重要:贝叶斯决策理论并不要求你真的相信某个先验。它要求你把不确定性写下来,把代价写下来,然后接受由此推出的结论。剩下的,是计算问题。