[論文レビュー] Perfect simulation of spatial point processes using dominated coupling from the past with application to a multiscale area-interaction point process

Graeme K. Ambler, Bernard W. Silverman|ArXiv.org|Mar 15, 2009

Point processes and geometric inequalities参考文献 16被引用数 9

ひとこと要約

本稿は、局所的に安定した空間点過程の完全シミュレーションのための、過去からの支配的カップリング（CFTP）アルゴリズムを提示する。特に、純粋に引力的でも斥力的でもないモデルに対して特に有効である。また、複数スケールでのクラスタリングをモデル化可能な、新しいマルチスケール領域相互作用過程を導入し、その応用をレッドウッドの種子lingsデータに適用することで、正確なシミュレーションとパラメータ推定を実現した。

ABSTRACT

We consider perfect simulation algorithms for locally stable point processes based on dominated coupling from the past. A version of the algorithm is developed which is feasible for processes which are neither purely attractive nor purely repulsive. Such processes include multiscale area-interaction processes, which are capable of modelling point patterns whose clustering structure varies across scales. We prove correctness of the algorithm and existence of these processes. An application to the redwood seedlings data is discussed.

研究の動機と目的

• 純粋に引力的でも斥力的でもない点過程に対して、計算コストを低減した実用的な完全シミュレーションアルゴリズムの開発。既存のCFTP手法の限界を克服することを目的とする。
• 複数の空間スケールでクラスタリング行動を捉えることが可能な、新しいマルチスケール領域相互作用点過程モデルの導入。
• 完全シミュレーション技術を用いて、複雑な点過程モデルの正確なシミュレーションとパラメータ推定を可能にすること。
• 実世界のデータ、特にマルチスケール相互作用パターンを示すレッドウッドの種子lingsデータに対して、本手法の有効性を示すこと。

提案手法

• 本稿は、非極端なモデルに対して計算量をO(2^n)からO(n)に削減する支配的CFTPアルゴリズムを提案し、局所的に安定した点過程の実用的シミュレーションを可能にする。
• この手法は確率的単調性に基づくカップリングに依存し、Papangelou条件付き強度の均一な上界λを用いて収束を保証する。
• 新しいマルチスケール領域相互作用過程は、複数の相互作用半径R₁, R₂とパラメータγ₁, γ₂を含む密度関数によって定義され、小スケールで斥力的、中スケールで引力的となる挙動を可能にする。
• 支配的過程を用いたリジェクトサンプリングを用いることで、バーンインや収束チェックを必要とせず、ターゲット点過程分布からの正確なシミュレーションが実現される。
• パラメータ推定は擬似尤度最大化により実施され、K関数およびL関数の検査を通じてモデルの適合度をガイドする。
• 複数の相互作用項を含む一般形の強度測度を用いることで、高階のマルチスケールモデルへの拡張が可能となる。

実験結果

リサーチクエスチョン

• RQ1 純粋に引力的でも斥力的でもない点過程に対して、計算コストを低減した完全シミュレーションアルゴリズムを開発できるか？
• RQ2 マルチスケール領域相互作用過程モデルは、複数の空間スケールにわたる複雑なクラスタリングパターンを捉えることができるか？
• RQ3 支配的CFTPを用いて、このような複雑なモデルの正確なシミュレーションをどのように達成できるか？
• RQ4 提案されたモデルは、レッドウッドの種子lingsデータのようにマルチスケール相互作用を示す実データにどの程度適合するか？

主な発見

• 提案された支配的CFTPアルゴリズムは、1反復あたりO(n)の計算量を達成し、純粋に引力的または斥力的でないモデルに対しても実用的である。これに対して、従来のO(2^n)手法とは対照的である。
• マルチスケール領域相互作用過程は、レッドウッドの種子lingsデータをよく再現し、小スケールでの斥力的挙動、中スケールでの引力的挙動、大スケールでの再びの斥力的挙動を捉えている。
• モデルのパラメータ推定値はλ = 0.118、γ₁ = 2000、γ₂ = 10⁻²⁰⁰、R₁ = 0.07、R₂ = 0.013であり、非常に小さなスケールでの強い斥力が示唆されている。
• 19個のモデル実現を用いたシミュレーション結果から、L関数およびT関数がレッドウッドの種子lingsデータの実測関数とよく一致しており、良好な適合が確認された。
• モデルのマルチスケール相互作用の捉え具合は妥当であるが、R₁のわずかな調整により適合度が向上する可能性があり、スケールパラメータに敏感であることが示唆された。
• 一般形の強度測度により、3スケール以上への拡張が可能であることが示されており、今後の複雑な空間構造のモデリングに向けた基盤を提供している。