[論文レビュー] A Central Limit Theorem for Fleming-Viot Particle Systems with Hard Killing

Bernard Delyon, Frédéric Cérou+2 more|arXiv (Cornell University)|Sep 20, 2017

Stochastic processes and statistical mechanics参考文献 15被引用数 6

ひとこと要約

本稿は、最小限の仮定の下で、硬い殺滅を伴うフレミング＝ビオの粒子系に対して中心極限定理（CLT）を確立し、固定時刻 T における生存確率および条件付き分布の推定子の漸近正規性を証明する。この手法は、マルティンゲール分解、二次変動解析、および時間に一様な L² 評価に依拠しており、境界で吸収される楕円型拡散過程を含む、先行研究ではカバーされていなかったケースを拡張する。

ABSTRACT

Fleming-Viot type particle systems represent a classical way to approximate the distribution of a Markov process with killing, given that it is still alive at a final deterministic time. In this context, each particle evolves independently according to the law of the underlying Markov process until its killing, and then branches instantaneously on another randomly chosen particle. While the consistency of this algorithm in the large population limit has been recently studied in several articles, our purpose here is to prove Central Limit Theorems under very general assumptions. For this, we only suppose that the particle system does not explode in finite time, and that the jump and killing times have atomless distributions. In particular, this includes the case of elliptic diffusions with hard killing.

研究の動機と目的

• 硬い殺滅を伴うフレミング＝ビオ粒子系に対して中心極限定理（CLT）を確立し、従来のソフト殺滅や有界な carré du champ 条件を要する結果を超えて拡張すること。
• 一般条件の下で、固定時刻 T における生存確率 $ p_T $ および条件付き分布 $ \eta_T $ の推定子の漸近的挙動を分析すること。
• 正規化された粒子系の摂動の分布収束を証明し、一般設定における漸近分散を同定すること。
• 従来、有界な殺滅強度が存在しないため、カバーされなかった領域の境界に達した際に殺滅される楕円型拡散過程のケースを、ジャンプおよび殺滅時刻の原子をもたない分布を仮定することで取り扱うこと。
• 基礎となる過程に最小限の正則性条件しか課さない条件下で、稀な事象の推定における逐次モンテカルロ法の厳密な基礎を提供するため、CLT を証明すること。

提案手法

• 粒子系の経験測度および生存確率推定子のマルティンゲール分解を用いて、信号とノイズを分離する。
• 粒子過程 $ p_t^N $ に対して時間に一様な $ \mathbb{L}^2 $-推定を用い、時間区間 $[0,T]$ 全体にわたる安定性を保証する。
• 二次変動推定および $ p_t^N = (1 - 1/N)^{N\mathcal{N}_t} $ を含む部分積分公式を適用し、$ \mathcal{N}_t $ は分岐回数を数える。
• 二次変動過程の収束を用いて漸近分散の表現を導出し、殺滅およびジャンプ時刻の原子をもたない分布の仮定を活用する。
• 停止時刻の議論と独立構造を用いて、ジャンプの同時発生がほとんど確実に起こらないことを証明し、粒子系が適切に定義されることを保証する。
• 正規化された摂動過程に対するマルティンゲール CLT を適用し、ガウス型極限への収束を保証する条件を満たす。

実験結果

リサーチクエスチョン

• RQ1 硬い殺滅を伴う過程のフレミング＝ビオ粒子系が、中心極限定理を満たすために必要な最小限の仮定は何か？
• RQ2 粒子数 $ N \to \infty $ のとき、推定子 $ \hat{p}_T^N $ および $ \hat{\eta}_T^N $ の分布はどのように漸近的に振る舞い、その漸近分散の形は何か？
• RQ3 殺滅強度が有界でない（硬い殺滅）ため、従来の結果ではカバーされていなかった、境界に達した際に殺滅される楕円型拡散過程に対しても CLT を確立できるか？
• RQ4 殺滅およびジャンプ時刻の分布が原子をもたないという性質が、分岐事象の非同時発生を保証し、CLT の有効性を確保する上で果たす役割は何か？
• RQ5 粒子系の漸近分散は、基礎となるマルコフ過程および時間変換過程 $ Q^{T-t}(\varphi) $ を用いてどのように特徴付けられるか？

主な発見

• 本稿は、粒子系が有限時間内で爆発しないこと、および殺滅・ジャンプ時刻の分布が原子をもたないことという最小限の仮定の下で、硬い殺滅を伴うフレミング＝ビオ粒子系に対して中心極限定理を証明する。
• 正規化された $ p_T^N $ の摂動の極限分布は、時間変換過程 $ Q^{T-t}(\varphi) $ に依存する明確に定義された漸近分散を持つガウス分布である。ここで $ Q^{T-t}(\varphi) $ はテスト関数の条件付き期待値を表す。
• 漸近分散は、粒子系の経験測度の分解におけるマルティンゲール成分の二次変動の極限として特徴付けられる。
• 境界で吸収される楕円型拡散過程（硬い殺滅）に対しても CLT が成立し、従来の研究でソフト殺滅や有界な carré du champ 条件を要していたため、カバーされなかったケースを含む。
• 時間に一様な $ \mathbb{L}^2 $-推定を $ p_t^N $ に対して確立し、時間区間 $[0,T]$ 全体にわたる摂動の制御に不可欠であることを示した。
• 付録の積分-by-parts 公式 (4.3)–(4.5) を用いることで、$ p_t^N $ を含む確率積分の精密近似が可能となり、誤差項が $ O(1/N) $ で有界になる。これは漸近的解析において重要である。