[论文解读] The fixation line in the $Λ$-coalescent

Olivier Hénard|arXiv (Cornell University)|Jul 2, 2013

Stochastic processes and statistical mechanics参考文献 24被引用 9

一句话总结

本文引入了固定线（fixation line）——在 $\Lambda$-君主模型中追踪种群最近共同祖先谱系的马尔可夫过程。它建立了与块计数过程的对偶性，并证明了在 $n\to\infty$ 时，$\operatorname{Beta}(2-\alpha,\alpha)$-君主模型中 $n$-君主的最终共祖事件中块数的分布收敛，且显式计算了极限随机变量的生成函数。

ABSTRACT

We define a Markov process in a forward population model with backward genealogy given by the $Λ$-coalescent. This Markov process, called the fixation line, is related to the block counting process through its hitting times. Two applications are discussed. The probability that the $n$-coalescent is deeper than the $(n-1)$-coalescent is studied. The distribution of the number of blocks in the last coalescence of the $n$-$\operatorname {Beta}(2-α,α)$-coalescent is proved to converge as $n ightarrow\infty$, and the generating function of the limiting random variable is computed.

研究动机与目标

• 定义并分析作为 $\Lambda$-君主模型中块计数过程对偶的固定线马尔可夫过程。
• 研究 $n$-君主的深度，特别是 $n$-君主比 $(n-1)$-君主更深的概率。
• 研究当 $n\to\infty$ 时，$n$-君主最终共祖事件中涉及的块数的极限分布。
• 计算 $\operatorname{Beta}(2-\alpha,\alpha)$-君主模型中极限随机变量的生成函数。

提出的方法

• 通过自然耦合方法对 $n$-君主进行耦合，构造向前时间的马尔可夫过程，从而引入固定线。
• 利用固定线与块计数过程之间的对偶性，分析首次 hitting 时间与共祖动态。
• 使用强度为 $dt\Lambda(dx)x^{-2}$ 的泊松点测度构造，以建模共祖事件。
• 运用内维乌分支过程推导 $n$-君主深度 $\tau_1^n$ 的渐近行为。
• 通过分析 $\alpha_1^n$（即固定线首次 hitting 时间）的增长率，建立分布收敛性。
• 使用组合恒等式（包括范德蒙德恒等式）证明对偶方程中被积函数的等价性。

实验结果

研究问题

• RQ1 当 $n\to\infty$ 时，$n$-君主最终共祖事件中涉及的块数的分布是什么？
• RQ2 $n$-君主的深度 $\tau_1^n$ 的渐近行为如何？
• RQ3 $\operatorname{Beta}(2-\alpha,\alpha)$-君主模型中，最后共祖事件的块数的极限生成函数是什么？
• RQ4 固定线如何通过 duality 和 hitting 时间与块计数过程关联？
• RQ5 在何种条件下，$n$-君主的深度几乎必然有界或无界？

主要发现

• $n$-君主的深度 $\tau_1^n$ 满足 $\tau_1^n - \log\log n \Rightarrow -\log \mathbf{e}$，当 $n\to\infty$ 时，其中 $\mathbf{e}$ 是速率为 1 的指数分布随机变量。
• 极限随机变量 $-\log \mathbf{e}$ 服从 Gumbel 分布，其累积分布函数为 $\mathbb{P}(-\log \mathbf{e} \leq x) = e^{-e^{-x}}$。
• 在 $\operatorname{Beta}(2-\alpha,\alpha)$-君主模型中，$n$-君主最终共祖事件中块数的分布收敛。
• 在 $\operatorname{Beta}(2-\alpha,\alpha)$-君主模型中，极限随机变量的生成函数被显式计算。
• 固定线过程被证明是向前时间的马尔可夫过程，与块计数过程对偶，且 hitting 时间对应于共祖事件。
• 证明依赖于 $\alpha_1^n$ 的 a.s. 增长率，以及利用内维乌分支过程推导极限分布。