[論文レビュー] A Mean Field Approximation of the Bolker-Pacala Population Model

Mariya Bessonov, Stanislav Molchanov+1 more|arXiv (Cornell University)|Feb 5, 2013

Stochastic processes and statistical mechanics参考文献 19被引用数 3

ひとこと要約

この論文は、ボルカー＝パケーラ空間的集団動態モデルを、格子上のロジスティックマルコフ連鎖としてモデル化することで、平均場近似を展開する。不変分布に関して局所中心極限定理および大偏差結果を確立し、大偏差が $ L^{-1/2} e^{-L} $ のように衰えることを示し、大規模な対称区間の境界への最初の通過時刻の漸近的性質を導出。その結果、期待到達時間がシステムサイズ $ L $ とともに指数関数的に増加することが明らかになった。解析は逐次的超幾何関数の漸近的性質に依拠し、カーツの関数法則を用いてグローバル極限定理を導出する。

ABSTRACT

We approximate the Bolker-Pacala model of population dynamics with the logistic Markov chain and analyze the latter. We find the asymptotics of the degenerated hypergeometric function and use these to prove a local CLT and large deviations result. We also state global limit theorems and obtain asymptotics for the first passage time to the boundary of a large interval.

研究の動機と目的

• 空間的ダイナミクスと競争を組み込んだ、確率的空間的ボルカー＝パケーラ集団モデルの取り扱いやすい平均場近似の構築。
• 平均場スケーリング下での不変分布および長時間挙動を理解するために、$ \mathbb{Z}^d $ 上のロジスティックマルコフ連鎖を分析すること。
• 平均場スケーリング下での不変測度に対して、局所中心極限定理および大偏差を含む極限定理を確立すること。
• 大きな対称区間の境界への平衡状態からの最初の通過時刻の漸近的表現を導出すること。
• カーツの結果を用いて、ロジスティックマルコフ連鎖に対するグローバル極限定理（関数的弱大数法則および中心極限定理）を確立すること。

提案手法

• ボルカー＝パケーラ過程を、出生、死、競争ダイナミクスを伴う $ \mathbb{Z}^d $ 上の離散時間マルコフ連鎖としてモデル化する。
• 空間的構造を有するプロセスを、状態依存のレートを持つ1次元のロジスティック出生死滅過程に、平均場近似を用いて簡略化する。
• 遷移確率を用いて不変分布を解析し、出生率と死滅率の比を含む再帰的関係式を解く。
• 不変測度に現れる逐次的超幾何関数の漸近的性質を、積分近似およびオイラー＝マクローリンの公式を用いて導出する。
• 不変測度に対して局所中心極限定理および大偏差原理を適用し、大偏差率が $ L^{-1/2} e^{-L} $ に支配されることを示す。
• ポテンシャル理論および最初の通過時刻の式を用いて、$ \mathbb{E}\tau_{n^* \to \{n_1,n_2\}} $ の漸近的性質を導出し、$ L $ に指数関数的に依存することを示す。

実験結果

リサーチクエスチョン

• RQ1 平均場スケーリング下でのロジスティックマルコフ連鎖の不変分布の漸近的挙動は何か？
• RQ2 ロジスティックマルコフ連鎖において、平衡状態からの大偏差はシステムサイズ $ L $ とともにどのようにスケーリングされるか？
• RQ3 大きな対称区間の境界への平衡状態からの期待最初の通過時刻の漸近的性質は何か？
• RQ4 既知の確率過程論の結果を用いて、ロジスティックマルコフ連鎖に対してグローバル極限定理（関数的弱大数法則および中心極限定理）を確立できるか？
• RQ5 平衡点 $ n^* $ の周りにおけるプロセスの対称性は、最初の通過時刻分布にどのように影響するか？

主な発見

• 不変分布の大偏差確率は、$ L^{-1/2} $ に比例する定数と $ e^{-L} $ の積として衰える。これは、システムサイズ $ L $ に伴い強く指数関数的に減少することを示している。
• ロジスティックマルコフ連鎖の不変分布に対して、逐次的超幾何関数の漸近的解析から導かれる収束速度を伴う局所中心極限定理が成立する。
• 平衡点 $ n^* $ から区間境界 $ n_1 $ および $ n_2 $ への期待最初の通過時刻は、$ \sqrt{L}^{-1} \exp\left( \frac{b}{\gamma}L \ln \rho_1 + \delta_1(1 - \ln \rho_1)n^* \right) $ と漸近的にスケーリングされ、$ L $ に指数関数的に依存することが示された。
• 不変測度の再帰的関係式の解 $ \psi_2(x) $ の漸近的解析により、$ \psi_2(n_2) \asymp \exp\left( \frac{1}{\omega} \int_0^{\delta_2(1 - \mu/b)} \ln(1+x) \, dx \right) $ が得られ、ここで $ \omega = \gamma/(bL) $ である。
• 対称性の仮定の下で、$ n^* $ から $ \{n_1,n_2\} $ への最初の通過時刻は、$ n_2 $ から $ n_1 $ への到達時間の期待値の半分に漸近的に一致し、式 (4.23) の主要な推定式が得られる。
• 平衡状態 $ n^* $ への期待再帰時間は $ O(\sqrt{L}) $ である一方、状態 $ n^* + \delta L $ への期待再帰時間は $ \sqrt{L} e^{O(L)} $ であり、これは大規模な偏差に対して極めて長い帰還時間を示している。