[论文解读] Quasi-stationary distribution for multi-dimensional birth and death processes conditioned to survival of all coordinates

Nicolas Champagnat, Denis Villemonais|arXiv (Cornell University)|Aug 13, 2015

Probability and Risk Models参考文献 18被引用 5

一句话总结

该论文为多物种出生死亡过程在所有物种存活条件下的唯一拟平稳分布建立了指数收敛性，即使吸收率无界。通过新颖的李雅普诺夫函数方法，证明了在强种内竞争或与维度相关的中性竞争条件下，对初始条件具有统一的指数混合性。

ABSTRACT

This article studies the quasi-stationary behaviour of multidimensional birth and death processes, modeling the interaction between several species, absorbed when one of the coordinates hits 0. We study models where the absorption rate is not uniformly bounded, contrary to most of the previous works. To handle this natural situation, we develop original Lyapunov function arguments that might apply in other situations with unbounded killing rates. We obtain the exponential convergence in total variation of the conditional distributions to a unique stationary distribution, uniformly with respect to the initial distribution. Our results cover general birth and death models with stronger intra-specific than inter-specific competition, and cases with neutral competition with explicit conditions on the dimension of the process.

研究动机与目标

• 分析当任一物种灭绝时即发生灭绝的多维出生死亡过程的拟平稳行为。
• 解决在生态系系统中常见的无界吸收（杀灭）率模型缺乏结果的问题。
• 建立在给定非吸收条件下，过程的条件分布以指数速度收敛到唯一拟平稳分布的条件。
• 将现有理论扩展至无界杀灭率模型，特别是具有强种内竞争或中性竞争的模型。
• 提供一个通用框架，利用李雅普诺夫函数适用于更广泛的具有无界吸收的相互作用粒子系统。

提出的方法

• 提出一种基于新李雅普诺夫函数的方法，以处理多维出生死亡过程中的无界吸收率。
• 利用过程的生成元 $ L $ 分析测试函数的无穷小演化，整合出生、死亡和竞争项。
• 对出生率 $ b_j(n) $、死亡率 $ d_j(n) $ 和竞争系数 $ c_{jk}(n) $ 施加条件，其中竞争通过 $ \left(\sum_k c_{jk}(n)n_k\right)^\gamma $ 影响死亡率。
• 应用 $ Q $-过程框架对过程进行非吸收条件化，确保极限平稳测度的存在性。
• 通过基于李雅普诺夫函数 $ V_\varepsilon(n) = \sum_{j=1}^r n_j^{-\varepsilon} $（其中 $ \varepsilon \in (0,1) $）的耦合论证，建立总变差范数下的统一指数收敛性。
• 将该方法适配至扩展情形，如突发灭绝事件（灾难）和多重出生，证明在矩条件和速率条件下仍可实现收敛。

实验结果

研究问题

• RQ1 在什么条件下，给定非吸收时，多物种出生死亡过程的条件分布会以指数速度收敛到唯一的拟平稳分布？
• RQ2 拟平稳分布理论能否扩展至具有无界吸收率的模型，特别是在具有强竞争的生态系系统中？
• RQ3 种内竞争强度（$ \gamma > 0 $）如何影响拟平稳分布的存在性与唯一性？
• RQ4 在何种维度 $ r $ 和竞争结构条件下，中性竞争模型可实现指数收敛？
• RQ5 李雅普诺夫函数方法能否适配至具有灾难或多重出生的模型？何种矩条件或速率条件可保证统一混合性？

主要发现

• 论文证明了在总变差范数下，条件分布对唯一拟平稳分布实现统一指数收敛，且满足界 $ \|\mathbb{P}_\mu(X_t \in \cdot \mid t < \tau_\partial) - \alpha\|_{TV} \leq C e^{-\lambda t} $，对所有初始测度 $ \mu $ 成立。
• 该收敛性在一般条件下成立，包括强种内竞争（$ \gamma > 0 $）和中性竞争，后者具有显式的维度相关条件。
• 该方法依赖于一种新颖的李雅普诺夫函数论证，可处理无界杀灭率，克服了先前工作中的一项关键局限。
• 对于具有灾难的一维过程，若出生率亚线性增长、死亡率二次增长，且灾难率被限制为死亡率的某一比例，则可实现指数收敛。
• 结果可推广至具有灾难和多重出生的多维过程，前提是灾难率的增长速度慢于竞争驱动的死亡率，且平均后代数量有统一有界。
• 在所提条件下，唯一拟平稳分布的存在性得到保证，且相应的 $ Q $-过程是指数遍历的。