Skip to content

定义Definition

分布恒等与接近性测试

施工中 · 内容逐步修订

Distribution identity testing · Distribution closeness testing · Discrete distribution goodness-of-fit testing

在有限离散域的 sample access 下,分别测试未知分布是否等于已知参考,或两份未知分布是否彼此接近。

形式陈述 ​

两个不同输入协议 ​

固定有限域 [n] 与 n≥2,距离阈值取 0<ε≤1。Identity testing 给定显式已知分布 q∈Δn 和未知SAMP oracle p,要区分

p=q与dTV(p,q)≥ε.

Closeness testing 则同时获得两个未知 oracle p,q,要区分 p=q 与 dTV(p,q)≥ε。第二个问题不能免费读取参考概率,需要从两边都抽样估计差异,因此通常更难。

Total variation在有限域上为

dTV(p,q)=12∑i=1n|pi−qi|.

两项任务都采用 exact yes 与 far no 的 gap,不要求在 0<dTV<ε 时作固定判断。若 yes 侧允许正距离,那是 tolerant 版本,样本复杂度会改变。

Identity 的最优量级 ​

以下最优量级统一取 0<ε≤1/4、成功率 2/3 与最坏样本预算。在标准 IID sample access 下,最坏已知参考分布的 identity testing 样本复杂度为

Θ(nε2).

上界可用按参考质量适当分组、截断或重加权的 Pearson/χ2 型统计量;不能把任意参考分布直接交给一个未经方差分析的平方和。若 qi=0 却抽到元素 i,可直接拒绝。对其余正参考质量,Cauchy–Schwarz 不等式说明

∑i:qi>0(pi−qi)2qi≥(∑i:qi>0|pi−qi|)2.

当 p 也支持在这些坐标上时,右端就是 4dTV(p,q)2。这解释聚合平方偏差的信号来源;达到最优界还须控制统计量的方差。[1]

分析中常取 λ>0,先独立抽取总样本数 N∼Pois(λ),再从 p 抽 N 次。此时各坐标计数精确独立,并满足 Ni∼Pois(λpi)。对任意非负整数 ki、K=∑iki,由Poisson 分布与条件计数概率直接得到

Pr[N1=k1,…,Nn=kn]=e−λλKK!K!∏iki!∏ipiki=∏ie−λpi(λpi)kiki!.

其中零强度坐标按恒为零的计数解释。若至少两个坐标具有正质量,固定正整数总样本数时,计数和受约束,这些坐标计数便不独立。Poisson 化也不自动给硬预算:可先把错误以常数次重复降到 1/6,将放大后总样本数的期望记为 Λ=O(λ),再在总数超过 ⌈12Λ⌉ 时停止并任意作答。Markov 不等式将额外失败界为 1/12,总错误至多 1/4,而样本数已有 O(λ) 硬上限。

直接估计每个 pi 到加性 O(ε/n) 会需要远多于 n 样本。最优 tester 聚合许多微小偏差,只解决是否相等,不重建整个概率向量。

下界构造一族对均匀分布作随机正负微扰的 alternatives。样本太少时,碰撞与计数分布在 null/alternative 下接近,任何统计量都无法稳定区分;这与“某个坐标差很多”的容易实例不同。

Closeness 的最优量级 ​

两份未知分布的 closeness testing 最优样本数为

Θ(max{n2/3ε4/3,nε2})

从每个 oracle 各取该量级样本。统计量要同时消除两边采样噪声,并处理重元素与轻元素:重元素可较准确估计,轻元素数量多、单项稀疏,需要碰撞型聚合。

Closeness 不能通过先用少量样本“学出 q”再直接调用已知参考的恒等测试,就自动达到最优量级。大域上的经验参考会遗漏大量轻元素;即便得到接近 q 的 q^,p=q 也不意味着 p=q^,未满足恒等测试的 exact-yes 前提。Closeness 的 n2/3 项正反映两边都未知时的额外困难。

直觉

恒等测试有一份可以直接查询概率值的参照,因此样本只需揭示未知分布偏离参照的总体信号。接近性测试的两边都带采样噪声;算法既要判断 observed difference 是否真实,又不能逐个学出大域上的轻元素,所以额外出现 n2/3/ε4/3 这一困难区间。

最优统计量的共同思路是聚合,而不是搜寻一个显眼坏坐标。TV 距离可能由大量 O(ε/n) 的小偏差共同组成,单项都淹没在抽样波动里;去偏的平方差或碰撞统计把这些微弱变化累积成可检测信号,同时控制 null 情形下的方差。

例子与边界

一个可计算的 TV 实例 ​

令已知参考为四点均匀分布

q=(1/4,1/4,1/4,1/4),

未知分布为

p=(1/2,1/6,1/6,1/6).

第一坐标多出 1/4,其余三坐标各少 1/12,所以

dTV(p,q)=12(14+3⋅112)=14.

样本序列若频繁出现元素 1 会推动统计量拒绝,但某一小批样本也可能碰巧接近均匀。测试保证来自整个计数统计的尾界,不是“看见两次 1 就判不等”之类固定阈值。

成功率与样本分配 ​

对 0<δ≤1/3,把成功率从 2/3 提高到 1−δ,可用独立重复与多数表决,通常乘 O(log⁡(1/δ)) 样本。更精细算法可能直接调阈值,但失败率依赖仍要报告。

Closeness 若两边样本预算不同,应写成 (mp,mq) 的 trade-off;“总样本 m”不能默认意味着每边各 m。Known-reference identity 不需要从 q 抽样,读取显式 qi 的计算时间则是另一资源。

失败边界 ​

上述最优式针对最坏 n 点分布和标准 SAMP。若 q 有已知小支持、概率质量下界或特殊结构,instance-optimal identity complexity 可更低;conditional/evaluation oracle 也会改变量级。

TV far 不表示存在某个坐标差至少 ε。差异可分散到全部 n 个坐标,每项只有 O(ε/n),所以逐点阈值扫描不是 sample-efficient 算法。

最后,检验 p=q 与估计 dTV(p,q) 到加性误差是不同任务。后者要在整个距离区间输出数值,不能从 gap tester 的一次接受 bit 直接恢复。

推论与应用

这两项基准任务把分布测试中的三种信息来源区分开来:显式概率表、普通样本和更强的条件或求值 oracle。更换其中任何一项都可能改变最优量级,因而调用结论时必须同时带上域大小、距离阈值、成功率和每个 oracle 的样本预算。

它们也是判断“测试是否比学习便宜”的典型例子。Identity 只需 Θ(n/ε2) 个样本就能否定最坏参考分布,而以 TV 误差学习一般 n 点分布需要线性依赖于 n;任务只要求一个 gap 判定 bit,算法便可利用全局统计信号绕过完整重建。

分布测试模型的完整均匀性下界给出成对扰动、奇数域调整、混合似然比二阶矩及四点联合表;从TV口径推到最坏样本预算的Ω(n/ε2),并明确不以该下界替代匹配上界。

参考资料
  • Tugkan Batu, Lance Fortnow, Ronitt Rubinfeld, Warren D. Smith, and Patrick White, “Testing That Distributions Are Close,” FOCS, 2000, pp. 259–269.
  • Siu-On Chan, Ilias Diakonikolas, Gregory Valiant, and Paul Valiant, “Optimal Algorithms for Testing Closeness of Discrete Distributions,” SODA, 2014, pp. 1193–1203,Theorem 1;原文用 ℓ1 距离,换成 TV 只改变常数。
  • [1] Gregory Valiant and Paul Valiant, “An Automatic Inequality Prover and Instance Optimal Identity Testing,” Theorem 1 及其最坏 n 点分布推论。
  • Clément Canonne, “A Survey on Distribution Testing: Your Data Is Big. But Is It Blue?” Theory of Computing Graduate Surveys 9, 2020, pp. 1–100.
关系图谱14 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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