[論文レビュー] Perfect Simulation of Determinantal Point Processes

Laurent Decreusefond, Ian Flint+1 more|arXiv (Cornell University)|Nov 5, 2013

Random Matrices and Applications参考文献 14被引用数 5

ひとこと要約

本稿は、反復法の収束診断やバーンインを必要とせず、正確なサンプリングを可能にする、決定的ポイントプロセス（DPP）の完全なシミュレーション手法を提示する。この手法は、コピュリング・フォム・ザ・パスト（CFTP）を用い、理論的なコalescence時間の上界を $ O\big(\int_E J(x,x)\,d\nu(x) \log \int_E J(x,x)\,d\nu(x)\big) $ として確立した。この上界は、Rのspatstatライブラリを用いたガウス、マテルン、コーシー核を用いた定常DPPモデルのシミュレーションにより検証された。

ABSTRACT

Determinantal point processes (DPP) serve as a practicable modeling for many applications of repulsive point processes. A known approach for simulation was proposed in \cite{Hough(2006)}, which generate the desired distribution point wise through rejection sampling. Unfortunately, the size of rejection could be very large. In this paper, we investigate the application of perfect simulation via coupling from the past (CFTP) on DPP. We give a general framework for perfect simulation on DPP model. It is shown that the limiting sequence of the time-to-coalescence of the coupling is bounded by $K|Λ|\log K|Λ|$. An application is given to the stationary models in DPP.

研究の動機と目的

• 従来の拒否サンプリング手法の高い拒否率を回避する、DPPの完全なシミュレーションフレームワークの開発。
• コピュリング・フォム・ザ・パスト（CFTP）手法をDPPに適用し、均衡分布からの正確なサンプリングを可能にすること。
• 一般DPPモデルに対するCFTPアルゴリズムのコalescence時間の理論的上限を導出すること。
• 一般的な共分散関数（ガウス、マテルン、コーシー）を用いた定常DPPモデルに対して、この手法の実証的検証を行うこと。
• 理論的ペア相関関数およびサイズ分布との比較を通じて、シミュレーションの正確性を評価すること。

提案手法

• CFTPを用いてDPPを正確にシミュレートし、マルコフ連鎖の収束問題を回避する。
• CFTPに必要な確率的支配を定義するために、パパネロウ条件付き強度を用いる。
• コalescence時間の上界を $ O\big(\int_E J(x,x)\,d\nu(x) \log \int_E J(x,x)\,d\nu(x)\big) $ として導出する。ここで $ J(x,x) $ はパパネロウ強度である。
• スペクトル表現が利用できない場合の数値的取り扱いのため、核 $ \mathbb{K}(x,y) $ にフーリエ展開近似を適用する。
• $ \sum_{-N\leq i,j\leq N} \lambda_{i,j}^\Delta \approx \rho $ となるように $ N $ を選択し、強度の99.9％をカバーする。
• 統計的分析および可視化にはRのspatstatライブラリを用い、離散フーリエ変換によるPMF計算およびペア相関関数推定を実施する。

実験結果

リサーチクエスチョン

• RQ1 CFTP手法はDPPに適応可能であり、拒否サンプリングを回避する完全なシミュレーションを達成できるか？
• RQ2 DPPに対するCFTPアルゴリズムの理論的コalescence時間の上界は何か？
• RQ3 定常DPPモデルにおいて、コalescence時間は強度および核構造にどのように依存するか？
• RQ4 切断フーリエ展開は、核をどれほど正確に表現でき、DPPの統計的性質を保つことができるか？
• RQ5 シミュレートされたペア相関関数は、ガウス、マテルン、コーシーDPPモデルにおいて理論的期待とどの程度一致するか？

主な発見

• DPPに対するCFTPアルゴリズムのコalescence時間は、$ O\big(\int_E J(x,x)\,d\nu(x) \log \int_E J(x,x)\,d\nu(x)\big) $ で上界が保証され、実行時間に対する理論的保証が得られた。
• 実証的シミュレーションでは、実際の停止時間が理論的上界よりも著しく低く、$ \rho = 50 $、$ \alpha = 0.04 $ の場合に $ H \approx 57.5 $ であった。
• DPPのシミュレートされたサイズ分布は、理論的ポisson-二項分布とよく一致しており、切断フーリエ近似の正確性が裏付けられた。
• ガウス、マテルン、コーシーの3つのモデルすべてについて、シミュレートされたペア相関関数は理論的 $ g(r) $ とよく一致しており、特に $ r > 0 $ において良好であった。ただし、$ r \to 0 $ 付近ではフーリエ切断による高調波損失のため、わずかな差異が生じた。
• 本手法は、一般的な共分散関数を用いた定常DPPを、信頼性の高い収束性と統計的正確性を伴ってシミュレートできた。
• CFTPに基づくシミュレーションフレームワークは、拒否サンプリングのバイアスや非効率性を回避し、ターゲットDPP分布からの正確なサンプルを提供した。