Skip to content

定义Definition

Bayes 估计量与后验风险

施工中 · 内容逐步修订

Bayes estimator · Posterior expected loss

给定观测后最小化后验期望损失的行动规则及其条件风险。

形式陈述 ​

在Bayesian 模型中,观测 X=x 后参数后验为 Π(dθ∣x)。这是统计决策与风险在给定数据后的条件版本:对行动 a∈A 与损失 L(θ,a),后验风险定义为

ρ(a∣x)=∫ΘL(θ,a)Π(dθ∣x),

只要积分有定义。Bayes 行动是

δB(x)∈argmina∈Aρ(a∣x).

若可测选择存在,则 x↦δB(x) 是 Bayes 规则。全局 Bayes 风险通过迭代期望连接到逐数据的后验风险:

r(Π,δ)=E[L(Θ,δ(X))]=EX[ρ(δ(X)∣X)].

逐个 x 最小化后验风险因此也最小化先验—抽样联合分布下的 Bayes 风险。对本库采用的非负损失,可由Tonelli 定理处理扩展积分;若允许带符号损失,则需另有可积性条件。无论哪种情形,都须保证选择出的行动规则可测。

直觉

三种典型损失 ​

实值参数、平方损失 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+sa+b+n.

它可写为

a+ba+b+naa+b+na+b+nsn,

即先验均值与样本比例的加权平均。取 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公式用观测边缘密度的斜率和曲率恢复后验均值、方差;对任意正常先验,正态似然都保证每个固定观测处的后验矩有限。若先从同批数据估计先验,再把它固定来套公式,得到的是拟合模型的条件量。它不自动等于这套经验估计规则的重复抽样风险;重新拟合时的输入导数也不同于固定先验下的导数。

若后验二阶矩无穷,则对每个有限实数行动 a,平方损失风险都无穷;即使均值存在,也不足以恢复有限风险比较。另一方面,行动空间非闭或损失不下半连续时,风险下确界可能不被达到。argmin 可以作为空集定义,但选取 δB(x)∈argmin 前必须证明它非空。

不当先验若后验虽 proper,某些损失下的后验矩仍可能发散。代数上得到后验核不保证每种 Bayes 行动可定义。

若后验二阶矩有限且行动受约束,如 a≥0,平方损失 Bayes 行动是后验均值在非空闭凸行动集上的投影。非凸行动空间可能产生多个全局极小点,数值优化返回其中一个时要保留其决策等价性。

推论与应用

后验风险还能比较行动的不确定代价,而不只输出点估计。非对称损失把高估和低估赋予不同权重,Bayes 行动变为相应后验分位数。

在适当条件下,Bayes 规则可为 admissible;某些 least favorable prior 的 Bayes 规则还达到minimax 风险。但这些结论需要全参数风险分析,不能从“使用了先验”直接推出。

若两个行动 a1,a2 具有相同的有限最小后验风险 r∗,以概率 λ 选择前者、以概率 1−λ 选择后者,风险仍为 λr∗+(1−λ)r∗=r∗。这只用期望对随机化分布的线性性,不要求损失对行动线性。直接输出平均行动 λa1+(1−λ)a2 则是另一回事:它首先须合法,若风险对行动凸,才保证其风险不高于两者的平均。

唯一性可由凸行动集上后验风险的严格凸性保证,并不要求后验非退化。例如后验集中在一个点时,平方损失仍有唯一极小点。软件在并列极小点中固定返回第一个,只是确定了一条选择规则。

多分类行动可用损失矩阵 Cij 表示真实状态 i 却选择 j 的代价。给定后验概率 qi(x),行动 j 的风险是 ∑iCijqi(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。
关系图谱20 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系