[論文レビュー] Non-asymptotic Error Bounds for Sequential MCMC and Stability of Feynman-Kac Propagators

Nikolaus Schweizer|arXiv (Cornell University)|Apr 11, 2012

Markov Chains and Monte Carlo Methods参考文献 10被引用数 20

ひとこと要約

本稿は、フェインマン＝ハイズプロパゲーターの安定性と関連付けることで、逐次MCMCの非漸近的誤差バウンドを確立する。スペクトルギャップ、超有界性、有界相対密度を用いて明示的な誤差バウンドを導出し、高次元の積測度において次元に依存する計算コストがO(d³)として多項式的に増加することを示す。

ABSTRACT

We provide a generic way of deducing non-asymptotic error bounds for Sequential MCMC methods from suitable stability properties of Feynman-Kac propagators. We show how to derive this type of stability from mixing conditions for the MCMC dynamics, namely, spectral gaps and hyperboundedness, and from upper bounds on the relative densities in the sequence of distributions.

研究の動機と目的

• 高次元および非コンパクト状態空間に適用可能な、逐次MCMCアルゴリズムの非漸近的誤差バウンドを提供すること。
• 誤差バウンドとフェインマン＝ハイズプロパゲーターの安定性特性を結びつける汎用的フレームワークを確立すること。
• これらの安定性特性が成立する条件を導出し、特にスペクトルギャップおよび有界相対密度を特定すること。
• R^d 上の積測度の文脈において、誤差バウンドの次元依存性を分析すること。
• 有限状態空間や連続時間設定を超えて、離散時間の安定性および誤差解析を拡張すること。

提案手法

• フェインマン＝ハイズプロパゲーターの安定性を用いて、非正規化粒子近似の非漸近的平均二乗誤差バウンドを導出する。
• MCMCカーネルのスペクトルギャップおよび超有界性条件から導かれる、プロパゲーターのLp安定性バウンドを用いる。
• 連続する分布間の相対密度に対する一様上界を適用して、誤差伝搬を制御する。
• 超有界性の仮定のもとでL2p–Lpバウンドを確立し、これにより定理3.2の主誤差バウンドへの応用を可能にする。
• d次元積空間における分布の列を、相対密度が有界となるようにd−1個の中間測度を挿入することで再インデックス化する。
• MCMCダイナミクスが積測度において対数ソボレフ不等式を満たす場合、誤差バウンドの定数が次元に依存しないことを示す。

実験結果

リサーチクエスチョン

• RQ1 非漸近的誤差バウンドは、基礎となるフェインマン＝ハイズプロパゲーターの安定性からどのように導出可能か？
• RQ2 MCMCダイナミクスにどのような条件が課されると、逐次MCMCにおける誤差制御に十分な安定性が保証されるか？
• RQ3 高次元積測度において、計算コストは次元dに対してどのようにスケーリングするか？
• RQ4 適切な混合性および密度条件のもとで、誤差バウンドを明示的かつ次元に依存しない形にできるか？
• RQ5 相対密度バウンドおよびスペクトルギャップの性質は、レベル間での誤差増大を制御するために果たす役割は何か？

主な発見

• 定理3.2は、非正規化粒子近似の平均二乗積分誤差に対する汎用的非漸近的誤差バウンドを提供し、主項が中心極限定理からの漸近分散と一致する。
• スペクトルギャップおよび相対密度バウンドを用いて、フェインマン＝ハイズプロパゲーターのLp安定性が確立されるが、L2p–Lpバウンドの確立には追加の超有界性仮定が必要である。
• d次元積測度において、計算コストはO(d³)としてスケーリングする。これは、レベル数の増加（要因d）、1ステップあたりのMCMCコスト増加（要因d）、およびレベル数に比例する粒子数増加（要因d）に起因する。
• MCMCダイナミクスが対数ソボレフ不等式を満たす場合、誤差バウンドの定数は次元に依存しない。これは、積測度において既知の安定性を示す。
• 本フレームワークは一般状態空間および離散時間ダイナミクスに適用可能であり、相対密度の下界の一様性を仮定する必要がない。
• 結果は連続時間および有限状態空間の解析を実用的な離散時間アルゴリズムへと一般化し、非コンパクトかつ高次元設定においても有効な明示的バウンドを提供する。