[论文解读] Asympotic behavior of the total length of external branches for Beta-coalescents

Jean‐Stéphane Dhersin, Linglong Yuan|arXiv (Cornell University)|Feb 27, 2012

Stochastic processes and statistical mechanics参考文献 31被引用 8

一句话总结

本文研究了在 $1 < \alpha < 2$ 条件下，Beta$(2-\alpha,\alpha)$-君士坦丁的总外部分支长度 $L_{ext}^{(n)}$ 的渐近行为，推导出其二阶矩的精确渐近表达式。通过矩计算和 $\Lambda$-君士坦丁的性质，作者证明了 $\lim_{n\to\infty} n^{3\alpha-5} \mathbb{E}[(L_{ext}^{(n)} - n^{2-\alpha}\mathbb{E}[T])^2] = \frac{((\alpha-1)\Gamma(\alpha+1))^2 \Gamma(4-\alpha)}{(3-\alpha)\Gamma(4-2\alpha)}$，从而精确刻画了在重尾祖先过程下外部分支长度的波动特性。

ABSTRACT

We consider a $Λ$-coalescent and we study the asymptotic behavior of the total length $L^{(n)}_{ext}$ of the external branches of the associated $n$-coalescent. For Kingman coalescent, i.e. $Λ=δ_0$, the result is well known and is useful, together with the total length $L^{(n)}$, for Fu and Li's test of neutrality of mutations% under the infinite sites model asumption . For a large family of measures $Λ$, including Beta$(2-α,α)$ with $0

研究动机与目标

• 分析在 $1<\alpha<2$ 条件下，Beta$(2-\alpha,\alpha)$-君士坦丁中总外部分支长度 $L_{ext}^{(n)}$ 的渐近行为。
• 推导 $L_{ext}^{(n)}$ 的高阶矩（特别是二阶中心矩）的精确渐近表达式。
• 理解在允许多重共祖的君士坦丁过程中，外部分支长度波动的缩放行为。
• 通过提供精确的方差缩放表达式，弥补矩方法在这些过程中中心极限定理中无法捕捉正确收敛速度的局限性。

提出的方法

• 作者在 $\Lambda$-君士坦丁框架下，基于 $\Lambda$ 为 Beta$(2-\alpha,\alpha)$ 分布的背景，对总外部分支长度 $L_{ext}^{(n)}$ 进行基于矩的分析。
• 利用已知的弱收敛关系 $n^{\alpha-1}T^{(n)} \xrightarrow{d} T$，其中 $T^{(n)}$ 是随机选取的外部分支长度，$T$ 是极限的随机变量。
• 分析涉及对存活谱系数的尾概率和比值期望的控制，利用测度 $\nu(dx) = x^{-2}\Lambda(dx)$ 的性质。
• 关键估计通过分部积分和伽马函数渐近展开得出，特别是针对 $\rho^{(-\bar{r})}(x)$，其捕捉共祖测度的尾部行为。
• 证明依赖于一个关键引理，即对 $0 \leq r < \alpha$，有 $\mathbb{E}[(n/(n-X_1^{(n)}))^{r} \mathbf{1}_{X_1^{(n)} \geq ns}]$ 的有界性，这对矩收敛至关重要。
• 使用涉及 $ (a+b)^\beta $ 的凸性不等式，控制矩展开中的高阶项，确保在所需缩放下实现收敛。

实验结果

研究问题

• RQ1 在 $1 < \alpha < 2$ 条件下，Beta$(2-\alpha,\alpha)$-君士坦丁中总外部分支长度 $L_{ext}^{(n)}$ 的二阶矩的渐近行为是什么？
• RQ2 在该类君士坦丁过程中，$L_{ext}^{(n)}$ 的方差如何随种群大小 $n$ 缩放？
• RQ3 为何矩方法在这些过程的中心极限定理中无法捕捉正确的收敛速度？
• RQ4 当 $n \to \infty$ 时，$n^{3\alpha-5} \mathbb{E}[(L_{ext}^{(n)} - n^{2-\alpha}\mathbb{E}[T])^2]$ 的精确极限表达式是什么？
• RQ5 在该君士坦丁框架中，矩估计有效性的时刻阶数 $r$ 需要满足什么约束？

主要发现

• 总外部分支长度的二阶中心矩满足 $\lim_{n\to\infty} n^{3\alpha-5} \mathbb{E}[(L_{ext}^{(n)} - n^{2-\alpha}\mathbb{E}[T])^2] = \frac{((\alpha-1)\Gamma(\alpha+1})^2 \Gamma(4-\alpha)}{(3-\alpha)\Gamma(4-2\alpha)}$，提供了精确的渐近方差表达式。
• 缩放指数 $3\alpha - 5$ 表明，当 $\alpha > 5/3$ 时，$L_{ext}^{(n)}$ 的方差发散，当 $\alpha = 5/3$ 时趋于稳定，当 $\alpha < 5/3$ 时收敛，反映了共祖结构中的相变现象。
• 该结果证实，矩方法无法捕捉中心极限定理中正确收敛速度的原因在于 $\alpha \in (1,2)$ 时共祖机制的重尾性质。
• 分析表明，临界矩阶数 $r < \alpha$ 是收敛所必需的，这由当 $r \geq \alpha$ 时界不成立（对应 $\omega = \alpha$）所证实。
• $n^{\alpha-1}T^{(n)}$ 对非退化随机变量 $T$ 的极限分布是将外部分支长度围绕 $n^{2-\alpha}\mathbb{E}[T]$ 中心化并实现正确归一化的关键。
• 利用伽马函数恒等式和 $\Gamma(n+\bar{r})$ 的渐近近似，能够精确控制矩展开中组合与测度论项。