形式陈述
在Bayesian 模型 理路 Bayesian 统计模型 Bayesian statistical model 以参数先验和条件抽样模型共同定义参数与数据的联合分布,并由观测条件化得到后验。 中,观测 X = x 后参数后验为 Π ( d θ ∣ x ) 。这是统计决策与风险 理路 估计量、决策规则与风险 Estimator and decision rule · Statistical risk 从观测到行动的可实施规则,以及在逐参数、Bayes 与最坏情形量词下的期望损失。 在给定数据后的条件版本:对行动 a ∈ A 与损失 L ( θ , a ) ,后验风险定义为
ρ ( a ∣ x ) = ∫ Θ L ( θ , a ) Π ( d θ ∣ x ) , 只要积分有定义。Bayes 行动是
δ B ( x ) ∈ argmin a ∈ A ρ ( a ∣ x ) . 若可测选择存在,则 x ↦ δ B ( x ) 是 Bayes 规则。全局 Bayes 风险通过迭代期望 理路 全期望公式与全方差公式 Law of total expectation · Law of total variance · Iterated expectation 借助条件信息分解总体均值,并把总波动拆成组内与组间两部分。 连接到逐数据的后验风险:
r ( Π , δ ) = E [ L ( Θ , δ ( X ) ) ] = E X [ ρ ( δ ( X ) ∣ X ) ] . 逐个 x 最小化后验风险因此也最小化先验—抽样联合分布下的 Bayes 风险。对本库采用的非负损失,可由Tonelli 定理 理路 Tonelli 定理 Tonelli's theorem 非负可测函数的二重积分与两种迭代积分相等,允许共同取无穷。 处理扩展积分;若允许带符号损失,则需另有可积性条件。无论哪种情形,都须保证选择出的行动规则可测。
直觉
三种典型损失
实值参数、平方损失 L ( θ , a ) = ( a − θ ) 2 下,若后验二阶矩有限,
ρ ( a ∣ x ) = Var ( Θ ∣ x ) + ( a − E [ Θ ∣ x ] ) 2 , 故唯一 Bayes 行动为后验均值。
绝对损失 | a − θ | 下,后验一阶绝对矩有限时,任何后验中位数都是 Bayes 行动;若中位数区间不唯一,Bayes 行动也不唯一。对离散行动和 0 –1 损失,选择后验概率最大的状态,即 MAP 行动。
这些结论由损失决定。后验均值、后验中位数和 MAP 不是三个可以随意互换的“Bayesian 点估计”。
例子与边界
取 a , b > 0 ,以 Beta( a , b ) 为先验。观察整数 n ≥ 1 次 Bernoulli 试验,其中成功 s 次(0 ≤ s ≤ n ),则
θ ∣ x ∼ Beta ( a + s , b + n − s ) . 平方损失下 Bayes 估计为
δ B ( x ) = a + s a + b + n . 它可写为
a + b a + b + n a a + b + n a + b + n s n , 即先验均值与样本比例的加权平均。取 a = b = 1 , n = 8 , s = 6 ,结果为 7 / 10 ,而样本比例为 3 / 4 。
若目标行动是预测下一次 0 或 1 且错判损失相同,则后验预测成功概率超过 1 / 2 时选 1 。这与估计连续参数的平方损失问题不同。
若低估损失权重为 c − > 0 、高估权重为 c + > 0 ,损失写成
L ( θ , a ) = c − ( θ − a ) + + c + ( a − θ ) + . 后验一阶绝对矩有限时,最小化后验风险得到分位数 q = c − / ( c − + c + ) 。当漏报代价更高时 q 上移,行动更保守地覆盖大参数;这是损失改变最优规则的可追踪机制,不需要修改后验本身。
边界与失败情形
Bayes 规则不总是 MAP。连续参数下“精确猜中”的 0 –1 损失几乎处处为一,不能简单用密度众数解释;MAP 常对应小邻域损失的极限,并且受参数化影响。
后验风险是给定模型和当前数据的条件期望,不等于频率学点风险 R ( θ , δ ) 。Bayes 规则平均意义好,不保证对每个固定 θ 都优于其他规则。
已知正态噪声的Tweedie公式 理路 Tweedie公式与正态后验矩 Tweedie's formula for normal means · Gaussian posterior score identity · Tweedie后验均值公式 从正态混合边缘密度的导数恢复后验均值和方差,证明后验均值单调,并区分固定先验公式、密度代入和重新拟合的导数。 用观测边缘密度的斜率和曲率恢复后验均值、方差;对任意正常先验,正态似然都保证每个固定观测处的后验矩有限。若先从同批数据估计先验,再把它固定来套公式,得到的是拟合模型的条件量。它不自动等于这套经验估计规则的重复抽样风险;重新拟合时的输入导数也不同于固定先验下的导数。
若后验二阶矩无穷,则对每个有限实数行动 a ,平方损失风险都无穷;即使均值存在,也不足以恢复有限风险比较。另一方面,行动空间非闭或损失不下半连续时,风险下确界可能不被达到。argmin 可以作为空集定义,但选取 δ B ( x ) ∈ argmin 前必须证明它非空。
不当先验若后验虽 proper,某些损失下的后验矩仍可能发散。代数上得到后验核不保证每种 Bayes 行动可定义。
若后验二阶矩有限且行动受约束,如 a ≥ 0 ,平方损失 Bayes 行动是后验均值在非空闭凸行动集上的投影。非凸行动空间可能产生多个全局极小点,数值优化返回其中一个时要保留其决策等价性。
推论与应用
后验风险还能比较行动的不确定代价,而不只输出点估计。非对称损失把高估和低估赋予不同权重,Bayes 行动变为相应后验分位数。
在适当条件下,Bayes 规则可为 admissible;某些 least favorable prior 的 Bayes 规则还达到minimax 风险 理路 极小极大风险 Minimax risk 在模型族最坏参数上评价算法风险,再在所有允许算法中寻找最优值。 。但这些结论需要全参数风险分析,不能从“使用了先验”直接推出。
若两个行动 a 1 , a 2 具有相同的有限最小后验风险 r ∗ ,以概率 λ 选择前者、以概率 1 − λ 选择后者,风险仍为 λ r ∗ + ( 1 − λ ) r ∗ = r ∗ 。这只用期望对随机化分布的线性性,不要求损失对行动线性。直接输出平均行动 λ a 1 + ( 1 − λ ) a 2 则是另一回事:它首先须合法,若风险对行动凸,才保证其风险不高于两者的平均。
唯一性可由凸行动集上后验风险的严格凸性保证,并不要求后验非退化。例如后验集中在一个点时,平方损失仍有唯一极小点。软件在并列极小点中固定返回第一个,只是确定了一条选择规则。
多分类行动可用损失矩阵 C i j 表示真实状态 i 却选择 j 的代价。给定后验概率 q i ( x ) ,行动 j 的风险是 ∑ i C i j q i ( x ) ,逐列计算后取最小者。若所有错分代价为同一个正常数、正确代价为零,这就退化为选择最大后验概率类别。
参考资料
James O. Berger, Statistical Decision Theory and Bayesian Analysis , 2nd ed., Springer, 1985,Ch. 4。
Christian P. Robert, The Bayesian Choice , 2nd ed., Springer, 2007,Ch. 2。
Andrew Gelman et al., Bayesian Data Analysis , 3rd ed., CRC Press, 2013,Ch. 2。