形式陈述
设一维目标密度为 π ( x ) = γ ( x ) / Z ,其中 0 < Z < ∞ 且可逐点评值 γ ( x ) ≥ 0 。在扩展空间定义联合密度
p ( x , u ) = 1 Z 1 { 0 < u < γ ( x ) } . 对 u 积分得到 π ( x ) 。两个条件分布 公理库 条件分布 Conditional distribution · Regular conditional distribution 给定观测值后随机量的概率律,以及它与原联合分布相容的核表示。 分别是
U ∣ X = x ∼ Uniform ( 0 , γ ( x ) ) , X ∣ U = u ∼ Uniform ( S u ) , 其中水平切片 S u = { x : γ ( x ) > u } ,第二式要求该切片有有限正测度。交替精确抽这两个条件分布就是扩展目标上的 Gibbs 核,边缘 X 因而保持 π ;归一化常数 Z 从未参与计算。这里的均匀分布 公理库 均匀分布 Uniform distribution 在有限集合或有限正测度区域上按计数或体积等比例分配概率的分布。 是在切片的参考测度上归一化,不是“随便在附近挑一点”。
一维中通常无法直接枚举 S u ,可用宽度 w > 0 的 stepping-out 与 shrinkage。给定旧点 x 0 ,先抽 u ∈ ( 0 , γ ( x 0 ) ) ;随机放置包含 x 0 的宽 w 区间,若最多允许 m − 1 次扩展,先均匀抽整数 J ∈ { 0 , … , m − 1 } ,令 K = m − 1 − J ,分别作为左、右扩展额度;每侧以宽度 w 扩展,直到端点落到切片外或该侧额度用尽。不能给两侧各自固定同样次数而直接沿用此证明。随后在区间内均匀提议 x ′ :若 γ ( x ′ ) > u 就接受;否则按 x ′ 位于 x 0 左右收缩对应端点并重抽。正确的随机初始偏移、扩展与收缩规则共同保证核不偏向旧点;把它们替换成任意优化搜索需要重新证明不变性。
作为MCMC 公理库 Markov 链 Monte Carlo Markov chain Monte Carlo · MCMC 构造以目标分布为不变分布的 Markov 转移,并用一条遍历轨道的时间平均估计目标期望。 ,切片核还需在目标支持上不可约并满足返性。扩展目标的不变性本身不给有限步独立性;X t 仍是相关轨迹,轨道均值的偏差、一致性和中心极限定理依赖与其他 MCMC 相同层次的假设。多维逐坐标切片还会继承坐标几何问题。
直觉
密度高的点之所以更常出现,是因为它们上方有更高的竖直柱,可容纳更多高度 u 。先在当前柱内选择高度,再在该高度的横截面上均匀移动,长期看每个 x 获得的扩展面积正比于 γ ( x ) 。
步长 w 不直接规定最终移动尺度,只决定寻找切片边界的起始尺子。w 太小会多次向外扩展,太大则在收缩中反复击中切片外;两者主要改变成本。然而若人为限制扩展次数,参数也会限制可达距离,可能从效率问题变成遍历问题。
图片加载失败 切片采样的水平集更新
例子与边界
对三点未归一化目标 γ ( a , b , c ) = ( 1 , 3 , 6 ) ,从 X t = b 开始。若本轮在 ( 0 , 3 ) 抽得 u = 2.4 ,则
S 2.4 = { b , c } . 在这个两点切片上均匀抽样,若离散均匀数落在后半区,就取 X t + 1 = c 。下一轮高度可在 ( 0 , 6 ) 变化:当 u > 3 时切片只有 { c } ,链留在高密度点;当 u < 1 时三点全可访问。各点被包含的高度长度分别是 1 , 3 , 6 ,所以边缘概率正比于目标权重,而不是正比于切片内点数。从 b 转到 c 的概率可直接积分:高度在 ( 0 , 1 ) 的概率为 1 / 3 ,此时选 c 的概率为 1 / 3 ;高度在 ( 1 , 3 ) 的概率为 2 / 3 ,此时选 c 的概率为 1 / 2 ,合计 4 / 9 。反向从 c 到 b 的概率为 2 / 9 ,因此 π ( b ) P ( b , c ) = ( 3 / 10 ) ( 4 / 9 ) = ( 6 / 10 ) ( 2 / 9 ) = π ( c ) P ( c , b ) ,与不变目标一致。
断裂支持需要结合搜索规则判断。例如目标只在 [ 0 , 1 ] ∪ [ 3 , 4 ] 上为正,且采用最大区间长度 m w < 2 的有限扩展规则,则任何候选区间都无法同时接触两个分量,链永远不能跨越空隙。仅凭“遇到切片外端点就停止”不能推出这一结论:较宽初始区间或较大扩展步可能直接跳过空隙。两个分量的任意混合仍可不变,但从单一分量出发的链不可约。必须使用能覆盖多个分量的全局区间、额外跨模态更新或分别估计分量质量;延长运行无法修复零转移概率。
无界或无限测度切片也破坏朴素条件抽样。例如不正规“目标” γ ( x ) = 1 定义在整条实线上时,S u = R 没有均匀概率分布,这同时暴露原目标 Z = ∞ 。对正规但极重尾目标,切片有限却可能极宽,使扩展成本与跳跃尺度具有重尾;应设置可报告的预算上限,而不是静默截断并假装仍从完整切片抽样。
推论与应用
切片构造把密度值转成随机等高集,适合归一化常数未知、局部尺度随位置变化而又能廉价评值的目标。它自动调整单步跨度的说法只在边界搜索正确且切片几何可处理时成立;多峰、强相关或高维薄壳仍可能让水平集合难以探索。
离散变量可直接枚举切片,单调或对数凹一维条件则可用解析边界减少评值。若把切片步骤嵌进 Gibbs 更新,每个坐标的扩展变量都属于相应条件目标,整体正确性来自条件核的复合。调参应同时报告每个有效样本的密度评值数、扩展次数和跨模态行为,不能只比较“接受率”,因为合规的 shrinkage 最终总会接受。
参考资料
Radford M. Neal, “Slice Sampling,” The Annals of Statistics 31(3), 2003, pp. 705–767,§§3–4、Figures 3 and 5,作者全文 。
Brian D. Ripley, APTS 2010–11: Computer-Intensive Statistics , University of Oxford 课程讲义,2011,p.28 “Slice sampling”,讲义 。