[論文レビュー] 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

ひとこと要約

この論文は、一般の摂動フレームワークへスティン法を拡張し、既知の簡単な基本分布からの解を活用することで、複雑な分布（例えば、複合ポアソン分布、正規分布、ポアソン過程など）の正確な近似を可能にしている。主な貢献は、コルモゴロフ、ワサーライン、全 Variation 距離の各指標において、近似誤差にきびしいバインディングを保証する強固な摂動理論である。

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.

研究の動機と目的

• スティン法における摂動法を、ポアソン近似にとどまらず、正規分布、複合ポアソン分布、二項分布、ポアソン過程近似へ一般化すること。
• 摂動された分布のスティン方程式の解が、基本分布からの良い解析的性質をどのように継承するかの条件を確立すること。
• コルモゴロフ、ワサーライン、全 Variation などのさまざまな指標における近似誤差をバインディングする統一的なフレームワークを提供すること。
• 非標準的分布ではしばしば取り扱いが困難な、解のノルム（例：‖gh‖G）に対する解析的バインディングの制限を克服すること。
• 実数値の確率変数に限定されない一般の距離空間および確率的要素へも適用可能な方法を提示すること。

提案手法

• 目的分布 π₁ に対応するスティン作用素 A₁ を、理解が容易な基本分布 π₀ に関連する「よい」作用素 A₀ の小さな摂動とみなす摂動フレームワークを提案する。
• 基本分布のスティン方程式 A₀gₕ = h − π₀(h) の解 gₕ を用い、ネウマン級数展開により A₁gₕ = h − π₁(h) の近似解を構成する。
• 摂動方程式の解ノルム ‖gₕ‖G に対するバインディングを、基本解ノルムと摂動の大きさ ε の関数として導出し、重要な不等式 (1.4) が成立することを保証する。
• 適切なテスト関数クラスを用いて、正規分布、ポアソン分布、複合ポアソン分布、二項分布、ポアソン過程などさまざまな分布と、コルモゴロフ、ワサーライン、全 Variation などのさまざまな指標にこの手法を適用する。
• 生成子法、交換可能な対、テイラー展開などの技術を用い、摂動項に対する初期不等式 (1.1) を導出する。
• 微分方程式解析および標準正規分布への変換を用いて、特定のケース（例：局所的およびリプシッツ連続なテスト関数を用いた正規近似）における解ノルムの明示的バインディングを提供する。

実験結果

リサーチクエスチョン

• RQ1 スティン法における摂動法は、ポアソン近似にとどまらず、正規分布、複合ポアソン分布、二項分布、ポアソン過程分布へ一般化可能か？
• RQ2 摂動された分布のスティン方程式の解が、基本分布の解から良い解析的バインディングをどのように継承できるかの条件は何か？
• RQ3 摂動フレームワークは、コルモゴロフ、ワサーライン、全 Variation などの異なる確率的距離にどのように適応可能か？
• RQ4 基本測度に完全な台の条件を要求しない状況で、摂動されたスティン作用素が有効な近似バインディングを提供する十分条件は何か？
• RQ5 この方法は、実数値の確率変数に限定されず、一般の距離空間上の確率的要素に対しても適用可能か？

主な発見

• 摂動法は、正規分布、複合ポアソン分布、二項分布、ポアソン過程など、さまざまな分布に、さまざまな指標において成功裏に拡張可能である。
• h = 1(−∞,z] に対する正規近似では、解ノルムが ‖y‖∞ ≤ (1/4)√(2π/(1−ψ)) で抑えられ、微分係数は ‖y′‖∞ ≤ 1 で抑えられる。
• 有界な h に対しては、解ノルムが ‖y‖∞ ≤ √(2π/(1−ψ))‖h‖∞ で抑えられ、微分係数ノルムは ‖y′‖∞ ≤ 4‖h‖∞ で抑えられる。
• リプシッツ連続な h に対しては、解ノルムが ‖y‖∞ ≤ (2/(1−ψ))‖h′‖∞ で抑えられ、2階微分係数ノルムは ‖y′′‖∞ ≤ (2/√(1−ψ))‖h′‖∞ で抑えられる。
• 複合ポアソン近似において、コルモゴロフ距離の明示的誤差バインディングが得られ、誤差項は ∑l≥2 l(l−1)|λ₂ₗ| ≤ ∑ᵢ≥ₙ₊₁ pᵢ²/(1−2pᵢ)² で制御される。
• このフレームワークにより、π₀ に対する台の条件を緩和でき、π₀ が標本空間全体にわたって定義されていなくても応用可能となる。