Skip to content

算法Algorithm

自归一化重要性采样

施工中 · 内容逐步修订

Self-normalized importance sampling · SNIS

用未归一化目标与提议的权重和估计未知归一化常数,并以归一化权重形成目标期望的比率估计量。

形式陈述 ​

设可测空间上的目标只以未归一化密度 γ:X→[0,∞) 给出,

0<Z=∫Xγ(x)dx<∞,π(x)=γ(x)Z.

提议密度 q 必须在 γ(x)>0 的区域为正;分母要估计整个 Z,即使 f(x)=0 也不能漏掉该处目标质量。从 q 独立抽取 X1,…,XN,定义原始权重和归一化权重

Wi=γ(Xi)q(Xi),W~i=Wi∑j=1NWj.

重要性采样的两个样本均值分别估计 ZEπf 与 Z,自归一化估计量取其比值:

I^NSN=∑i=1NW~if(Xi)=N−1∑iWif(Xi)N−1∑iWi.

上述归一化权重与比率先在 ∑iWi>0 的事件上定义。为使大数律与极限分布陈述中的统计量在每个有限 N 都有定义,可在真实总权重为零的事件上给比率指定任意有限值,例如零;这是数学上的补充约定,实现仍应报告这一批次失败。

算法输入是 γ 的逐点评值、可采样且可求密度的 q、函数 f 与样本数;输出是比率估计和权重诊断。γ 乘任何正常数,所有 W~i 都不变,因此不需要知道 Z。实现时通常对 log⁡Wi 减去最大值后指数化,避免上下溢;这只改变共同尺度,不改变结果。若真实权重中至少有一个为正、但浮点指数化把它们全变成零,这是数值下溢,应回到对数权重和稳定归一化。若所有真实权重确实为零,则这一批样本没有携带目标质量,算法应报告失败并检查采样覆盖;不能返回均匀权重掩盖问题,也不能仅凭一个这样的批次便断言提议的理论支撑必然错误。

若 EqW=Z>0、Eq|Wf|<∞,则大数律给出平均权重几乎必然趋于 Z,因此分母几乎必然最终为正。具体地,令 p0=q{W=0}<1,独立性给出 P(∑iWi=0)=p0N→0。再对分子、分母应用大数律并取比率,得到 I^NSN→Eπf 几乎必然成立,零总权重事件上的上述有限补充值不影响结论。有限 N 时,比率通常有偏;即使原始分子、分母分别无偏,随机变量之比也不是比值的无偏估计。若

Eq[W2(f−I)2]<∞,I=Eπf,

令 Yi=Wi(f(Xi)−I),则 EqYi=0。在正分母事件上,N(I^NSN−I)=(N−1/2∑iYi)/(N−1∑iWi)。若 Varq(Yi)>0,一元中心极限定理处理分子,大数律使分母趋于 Z>0,再由Slutsky 定理取商。若中心化方差为零,Yi=0 几乎必然,得到相应的退化零极限;真实零权重批次的概率趋零,补充定义同样不改变极限。因此

N(I^NSN−I)→dN(0,Varq(W(f−I))Z2).

因此未缩放估计量的常规大样本方差近似为 Z−2Varq(W(f−I))/N;要进一步声称真实有限样本方差按这个比例收敛,还需控制二阶矩的极限交换。有限二阶加权矩本身已比只要求一致性的可积条件更强。

直觉

原始权重只有相对大小有意义。分母把这一批样本携带的总“目标质量”定为一,分子再按质量份额平均函数值。未知的 Z 同时出现在目标概率和原始权重尺度中,于比值内消去。

归一化并没有创造信息。若一个样本拿走 99% 的总权重,结果几乎只是那个点的函数值;权重和仍等于一,只说明会计平衡,不说明覆盖充分。比率结构还把分子与分母的波动耦合起来:这种相关性常降低尺度误差,却带来有限样本偏差。

自归一化权重的单位质量记账
例子与边界

取三点未归一化后验

γ=(1,3,6),Z=10,f(a,b,c)=(0,1,2),

以及 q=(0.5,0.25,0.25)。若一轮样本为 (a,b,c),原始权重为 (2,12,24),总和 38,归一化权重为 (1/19,6/19,12/19)。于是

I^3SN=1190+6191+12192=3019≈1.579,

而真值为 1.5。同一轮若错误地用 N−1∑Wifi 当作后验期望,会得到 20;它实际估计的是 ZI=15。N=1 时自归一化估计恒等于 f(X1),其期望为 Eqf=0.75,直接展示了有限样本偏差。

支撑完整也不保证有限渐近方差。令目标为 (0,1) 上均匀分布,提议 q(x)=2x,函数 f(x)=x。除零测集端点外支撑条件成立,W=1/(2x) 且 EqW=1;但在 x↓0 时,

Eq[W2(f−1/2)2]=∫01(x−1/2)22xdx=∞.

因此估计量仍一致,却没有通常的有限方差中心极限。若 q 在目标有正质量的整段区域为零,归一化只能在已访问区域内重新分配,极限也会错;“权重之和为一”无法发现这类 support mismatch。

推论与应用

自归一化权重可直接形成后验均值、分位数的加权经验分布和重加权预测。要估计多个 f,样本与权重可复用,但每个函数都有自己的二阶矩风险;某个平滑均值稳定,不代表尾部概率也稳定。

权重有效样本量 1/∑iW~i2 把集中程度压成一个数,适合触发重采样或比较提议。它不是误差定理:相同权重可对应遗漏一个模态的整批样本。序贯 Monte Carlo 正是反复执行增量加权和归一化,并在权重退化后重采样;可靠实现仍需保留对数权重、归一化常数增量和祖先信息,以便检查数值稳定性与路径退化。

有限总体调查中的Hájek 比率估计也用权重和归一化分子。它可以消除与估计总体规模共同变化的波动,但同样一般有有限样本偏差;其方差由实际抽样设计决定,不能直接套用独立重要性样本的表达。

参考资料
  • Christian P. Robert and George Casella, Monte Carlo Statistical Methods, 2nd ed., Springer, 2004, §3.3, importance sampling with unknown normalizing constants.
  • Art B. Owen, Monte Carlo Theory, Methods and Examples, 2013, §9.2, self-normalized importance sampling.
  • A. B. Kong, J. S. Liu, and W. H. Wong, “Sequential Imputations and Bayesian Missing Data Problems,” Journal of the American Statistical Association 89(425), 1994, pp. 278–288.
关系图谱7 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

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

上位 / 更一般

下位 / 直接特例

暂未标注直接特例。

类型化关系

使用的工具

被这些条目使用