[論文レビュー] A Sharp First Order Analysis of Feynman-Kac Particle Models

Hock Peng Chan, Pierre Del Moral+1 more|arXiv (Cornell University)|Nov 14, 2014

Cold Atom Physics and Bose-Einstein Condensates参考文献 20被引用数 4

ひとこと要約

本稿は、経路空間上の有界でない関数に対しても成り立つ、$\mathcal{O}(n/N)$ 階の鋭い非漸近的バイアス推定を提供する、前向きおよび後向きフェインマン＝カツの粒子モデルに対する新しい1次漸近解析を導入する。テンソル積の粒子密度プロファイルにおける、非漸近的上界および下界の混沌の伝播推定の初の結果を確立し、核および最小化定数の1次展開を用いた粒子ギブス・サンプラーへの応用を含む。

ABSTRACT

This article provides a new theory for the analysis of forward and backward particle approximations of Feynman-Kac models. Such formulae are found in a wide variety of applications and their numerical (particle) approximation are required due to their intractability. Under mild assumptions, we provide sharp and non-asymptotic first order expansions of these particle methods, potentially on path space and for possibly unbounded functions. These expansions allows one to consider upper and lower bound bias type estimates for a given time horizon $n$ and particle number $N$; these non-asymptotic estimates are of order $\mathcal{O}(n/N)$. Our approach is extended to tensor products of particle density profiles, leading to new sharp and non-asymptotic propagation of chaos estimates. The resulting upper and lower bound propagation of chaos estimates seems to be the first result of this kind for mean field particle models. As a by-product of our results, we also provide some analysis of the particle Gibbs sampler, providing first order expansions of the kernel and minorization estimates.

研究の動機と目的

• フェインマン＝カツモデルの前向きおよび後向き粒子近似を分析するための新しい理論的枠組みを構築すること。
• 経路空間および有界でない関数に対しても有効な、非漸近的かつ鋭い1次展開を、粒子近似から導出すること。
• テンソル積の粒子密度プロファイルにおける、非漸近的上界および下界の混沌の伝播推定の初の結果を確立すること。
• マルコフ核および最小化定数の1次展開を用いて、粒子ギブス・サンプラーを分析すること。
• 従来の理論を拡張し、有界でない関数に対して成り立ち、かつ下界を含む結果を提供すること。

提案手法

• 前向きおよび後向きの設定において、粒子近似の1次展開を導出するために、摂動半群技法を用いる。
• 複数の粒子の相互作用を分析するために、経験的および組み合わせ的テンソル積解析を用いる。特に、経験的測度の文脈で検討する。
• 粒子の依存関係および選択メカニズムをモデル化するために、融合度および第二種スターリング数を用いる。
• マルチノミアル融合遷移のバインドと遷移作用素の一様バインドを用いて、主要な推定値を導出する。
• 時間経過に伴う経験的測度のテンソル積およびその相互作用を用いた、粒子系の新規な分解を活用する。
• 理論的結果は、命題B.4の繰り返し適用および粒子遷移カーネルの作用素ノルムの推定を用いて証明される。

実験結果

リサーチクエスチョン

• RQ1 経路空間上でのフェインマン＝カツモデルの前向きおよび後向き粒子近似に対して、鋭い非漸近的1次展開を導出可能か？
• RQ2 粒子近似において、有界でない関数に対して、非漸近的上界および下界のバイアス推定 $\mathcal{O}(n/N)$ が成り立つか？
• RQ3 テンソル積の粒子密度プロファイルに対して、上界および下界を併せ持つ混沌の伝播推定を確立可能か？
• RQ4 この1次分析が、粒子ギブス・サンプラーの収束特性に与える影響は何か？
• RQ5 新規な混沌の伝播推定は、有界でない関数および最小限の仮定のもとで有効か？

主な発見

• 本稿は、経路空間上の有界でない関数に対しても、前向きおよび後向き粒子近似の非漸近的バイアス推定 $\mathcal{O}(n/N)$ を確立する。
• テンソル積の粒子密度プロファイルに対する、非漸近的上界および下界の混沌の伝播推定の初の結果を提供する。これも $\mathcal{O}(n/N)$ 階である。
• 関数 $f = g^{\otimes q}$ に対して、弱い仮定のもとで、推定 $Q^{(q)}_{p,n}(g^{\times q}) \leq Q^{\times q}_{p,n}(g^{\times q}) + \rho^{q+1}eq^2(n-p)/N\cdot Q_{p,n}(g^q)$ が成り立つ。
• 分析により、粒子ギブス核の1次展開が得られ、最小化推定および明示的な $\mathcal{O}(n/N)$ 依存性を有する収束速度の分析が可能になる。
• 結果は最小限の仮定（A1）–（A2）のもとで得られ、$N \geq eq^2n\rho$ を満たすことで安定性および収束性が保証される。
• 混沌の伝播推定が鋭いことが示され、ある有限な $a(\rho)$ を用いて、$(1-a(\rho)n/N)\eta^{(q)}_n(f) \leq \mathbb{E}[m(\xi_n)^{\odot q}(f)] \leq \eta^{(q)}_n(f)(1+a(\rho)n/N)$ の形のバインドが成り立つ。