[論文レビュー] Multiple projection MCMC algorithms on submanifolds

Tony Lelièvre, Gabriel Stoltz+1 more|arXiv (Cornell University)|Mar 20, 2020

Markov Chains and Monte Carlo Methods参考文献 48被引用数 7

ひとこと要約

本稿では、制約によって定義される部分多様体上の確率分布をサンプリングするための、複数の射影を許容する新しいマルコフ連鎖モンテカルロ（MCMC）アルゴリズムを提案する。従来の手法を一般化し、射影段階で複数の解が得られる場合にも対応できる。主な貢献は、複数の射影が存在する場合でも正しくサンプリングを保証する厳密な可逆性条件の確立であり、これにより効率性が向上し、拒否率を増加させることなく非局所的移動が可能になる。

ABSTRACT

We propose new Markov Chain Monte Carlo algorithms to sample probability distributions on submanifolds, which generalize previous methods by allowing the use of set-valued maps in the proposal step of the MCMC algorithms. The motivation for this generalization is that the numerical solvers used to project proposed moves to the submanifold of interest may find several solutions. We show that the new algorithms indeed sample the target probability measure correctly, thanks to some carefully enforced reversibility property. We demonstrate the interest of the new MCMC algorithms on illustrative numerical examples.

研究の動機と目的

• 制約によって定義される部分多様体上の確率分布を正しくサンプリングするMCMCアルゴリズムの開発、特に射影段階で複数の解が得られる場合を想定する。
• 提案ステップに集合値関数を組み込むことで、既存の部分多様体上でのMCMC手法を一般化し、非局所的移動を可能にする。
• 複数の射影を考慮する一般化された可逆性条件を課すことにより、詳細平衡を保証し、正しく目標測度がサンプリングされることを確保する。
• RATTLE積分子と複数解を伴う射影処理を組み合わせることで、制約付きMCMCにおける幾何的およびシンプレクティックな視点を統合する。
• 数値例を通じて新アルゴリズムの実用的利点を示し、拒否率の低下と混合性能の向上を確認する。

提案手法

• 提案移動が複数の射影のいずれかを通じて逆にできるかどうかを検証する一般化された逆射影チェックを導入し、詳細平衡を保証する。
• 位相空間における制約付きハミルトニアン力学のシミュレーションにおいて、シンプレクティック構造を保持するRATTLE型積分子を用いる。
• 提案ステップでは、制約多様体への複数の有効な射影をもたらす写像を用い、これを集合値関数としてモデル化する。
• 受理確率は、利用可能な射影の数とその可逆性を考慮した修正されたメトロポリス・ハスティングス則により計算する。
• 複数の解が存在する場合でも、適切に構築された重み関数により遷移核の対称性を保証することで、詳細平衡を維持する。
• 状態空間および位相空間の両方の定式化に適用され、後者では運動量変数と運動量反転を組み込み、詳細平衡を達成する。

実験結果

リサーチクエスチョン

• RQ1 射影段階で複数の解が得られる場合でも、正しさを損なわずに部分多様体上でのMCMCアルゴリズムを一般化できるか？
• RQ2 射影写像が多価関数である場合に詳細平衡をどのように保てるか。また、正しい目標測度がサンプリングされるための条件は何か？
• RQ3 複数の射影を許容することで、制約付き多様体上でのMCMCサンプラーの拒否率と混合効率にどのような影響が生じるか？
• RQ4 複数の射影を考慮する場合、幾何的およびシンプレクティックな制約付きMCMCの定式化はどのように関連するか？
• RQ5 すべての可能な射影が計算可能である場合、逆射影チェックを単純化または削除できるか。どのような条件下で可能か？

主な発見

• 提案されたMCMCアルゴリズムは、射影段階で複数の解が得られる場合でも、一般化された可逆性条件を課すことにより、部分多様体上の目標確率測度を正しくサンプリングする。
• すべての可能な射影が計算可能である場合、逆射影チェックは不要となり、射影の数を数えるだけで十分になるため、アルゴリズムが簡略化される。
• 状態空間および位相空間の両方の定式化において詳細平衡が達成され、位相空間バージョンでは運動量反転までを考慮した可逆性が保証される。
• 数値実験により、新アルゴリズムは拒否率を低下させ、単一射影法よりも大きな非局所的移動を可能にする。
• 理論的解析により、[46]の幾何的アプローチと[31]のシンプレクティックアプローチの深い関係が明らかになり、複数射影フレームワークの下で両者の視点が統合される。
• 可逆性条件が受理ルールを通じて保たれていれば、射影計算における数値誤差に対してもアルゴリズムは頑健である。