[論文レビュー] 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)}$ の漸近的挙動を調査し、その2次のモーメントの正確な漸近的表現を導出する。モーメント計算および $\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

研究の動機と目的

• Beta$(2-\alpha,\alpha)$-コアレスケントにおいて、$1<\alpha<2$ の下で全外部枝長 $L_{ext}^{(n)}$ の漸近的挙動を分析すること。
• 特に2次の中心モーメントを含む、$L_{ext}^{(n)}$ の高次モーメントの正確な漸近的表現を導出すること。
• 複数の衝突を伴うコアレスケント過程における外部枝長の揺らぎのスケーリング挙動を理解すること。
• このような過程における中心極限定理のモーメント法の限界を、正確な分散スケーリングの提供によって解消すること。

提案手法

• 著者らは、$\Lambda$ が Beta$(2-\alpha,\alpha)$ 分布である $\Lambda$-コアレスケントの文脈において、全外部枝長 $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 Beta$(2-\alpha,\alpha)$-コアレスケントにおいて、$1 < \alpha < 2$ の下で、全外部枝長 $L_{ext}^{(n)}$ の2次のモーメントの漸近的挙動は何か？
• 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$ に対する制約は何か？

主な発見

• 全外部枝長の2次の中心モーメントは、$\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$ のとき分散が発散し、$\alpha = 5/3$ で安定化し、$\alpha < 5/3$ のとき収束することを示し、コアレスケント構造における相転移を反映する。
• この結果は、$\alpha \in (1,2)$ の下でコアレスケント機構が重い尾を持つことにより、モーメント法が正しく収束速度を捉えられていないことが原因であることを確認する。
• 解析により、収束のためには $r < \alpha$ が必須であることが判明し、$r \geq \alpha$ のときバインドが失敗すること（$\omega = \alpha$ に対応）が示される。
• 外部枝長を $n^{2-\alpha}\mathbb{E}[T]$ のまわりに中心化するために、$n^{\alpha-1}T^{(n)}$ の極限分布が非退化な確率変数 $T$ に収束することは不可欠である。
• ガンマ関数の恒等式および $\Gamma(n+\bar{r})$ の漸近的近似を用いることで、モーメント展開における組合せ的および測度論的項を精密に制御できる。