[论文解读] On Stein's method and perturbations

Andrew D. Barbour, V. Čekanavičius+1 more|ArXiv.org|Feb 1, 2007

Random Matrices and Applications参考文献 13被引用 5

一句话总结

本文将 Stein 方法扩展至一般的扰动框架，通过利用已知的简单基分布解，实现了对复杂数分布（如复合泊松分布、正态分布和泊松过程）的精确近似。核心贡献是提出了一套稳健的扰动理论，确保在 Kolmogorov、Wasserstein 和总变差度量下近似误差的紧致界。

ABSTRACT

Stein's (1972) method is a very general tool for assessing the quality of approximation of the distribution of a random element by another, often simpler, distribution. In applications of Stein's method, one needs to establish a Stein identity for the approximating distribution, solve the Stein equation and estimate the behaviour of the solutions in terms of the metrics under study. For some Stein equations, solutions with good properties are known; for others, this is not the case. Barbour and Xia (1999) introduced a perturbation method for Poisson approximation, in which Stein identities for a large class of compound Poisson and translated Poisson distributions are viewed as perturbations of a Poisson distribution. In this paper, it is shown that the method can be extended to very general settings, including perturbations of normal, Poisson, compound Poisson, binomial and Poisson process approximations in terms of various metrics such as the Kolmogorov, Wasserstein and total variation metrics. Examples are provided to illustrate how the general perturbation method can be applied.

研究动机与目标

• 将 Stein 方法中的扰动方法从泊松近似推广至正态分布、复合泊松分布、二项分布和泊松过程近似。
• 建立扰动分布的 Stein 方程解继承基分布解良好解析性质的条件。
• 提供一个统一的框架，用于在多种度量（包括 Kolmogorov、Wasserstein 和总变差度量）下界近似误差。
• 克服解析解范数（如 ‖gh‖G）在非标准分布下通常难以计算的局限性。
• 证明该方法可应用于一般度量空间和随机元，而不仅限于实值随机变量。

提出的方法

• 提出一种扰动框架，将目标分布 π₁ 的 Stein 算子 A₁ 视为与一个‘良好’基算子 A₀ 的小扰动，其中 A₀ 对应于一个理解充分的分布 π₀。
• 利用基分布的 Stein 方程解 gₕ（满足 A₀gₕ = h − π₀(h)）通过 Neumann 级数展开构造 A₁gₕ = h − π₁(h) 的近似解。
• 推导出扰动方程解范数 ‖gₕ‖G 的界，其表达式以基解范数和扰动大小 ε 表示，确保关键不等式 (1.4) 成立。
• 针对各种分布（正态分布、泊松分布、复合泊松分布、二项分布、泊松过程）和度量（Kolmogorov、Wasserstein、总变差度量），采用适当的测试函数类应用该方法。
• 采用生成元方法、交换对方法和泰勒展开等技术，推导扰动项的初始不等式 (1.1)。
• 通过微分方程分析和标准化正态形式变换，对特定情形（如局部有界和利普希茨测试函数的正态近似）提供解范数的显式界。

实验结果

研究问题

• RQ1 Stein 方法中的扰动方法能否从泊松近似推广至正态分布、复合泊松分布、二项分布和泊松过程分布？
• RQ2 在何种条件下，扰动分布的 Stein 方程解可继承基分布解的良好解析界？
• RQ3 如何将该扰动框架适配以处理不同概率度量，如 Kolmogorov、Wasserstein 和总变差度量？
• RQ4 扰动 Stein 算子在不依赖基测度全支撑条件的情况下，其有效近似界成立的充分条件是什么？
• RQ5 该方法能否应用于一般度量空间上的随机元，而不仅限于实值随机变量？

主要发现

• 该扰动方法在多种度量下成功扩展至一般分布，包括正态分布、复合泊松分布、二项分布和泊松过程。
• 对于正态近似中 h = 1(−∞,z] 的情形，解范数满足 ‖y‖∞ ≤ (1/4)√(2π/(1−ψ))，其导数满足 ‖y′‖∞ ≤ 1。
• 对于有界 h，解范数满足 ‖y‖∞ ≤ √(2π/(1−ψ))‖h‖∞，其导数范数满足 ‖y′‖∞ ≤ 4‖h‖∞。
• 对于利普希茨 h，解范数满足 ‖y‖∞ ≤ (2/(1−ψ))‖h′‖∞，其二阶导数范数满足 ‖y′′‖∞ ≤ (2/√(1−ψ))‖h′‖∞。
• 该方法为复合泊松近似提供了 Kolmogorov 距离下的显式误差界，误差项受 ∑l≥2 l(l−1)|λ₂ₗ| ≤ ∑ᵢ≥ₙ₊₁ pᵢ²/(1−2pᵢ)² 控制。
• 该框架可放宽对 π₀ 的支撑条件，使方法适用于 π₀ 并未完全支撑于样本空间的分布。