[論文レビュー] Moments of the frequency spectrum of a splitting tree with neutral Poissonian mutations

Nicolas Champagnat, Benoît Henry|arXiv (Cornell University)|Sep 22, 2015

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

ひとこと要約

本稿では、中立的ポisson的突然変異を伴う超臨界分裂木におけるアレル頻度スぺクトルの階乗モーメントの明示的な再帰的公式を導出する。新たな共析点過程の構成と一般化されたキャンベルの公式を用いることで、適切な指数的時間スケーリングの下で頻度スぺクトルのほとんど確実収束を確立し、一般分岐過程理論に依存しない新しい、素朴な証明を与える。

ABSTRACT

We consider a branching population where individuals live and reproduce independently. Their lifetimes are i.i.d. and they give birth at a constant rate b. The genealogical tree spanned by this process is called a splitting tree, and the population counting process is a homogeneous, binary Crump-Mode-Jagers process. We suppose that mutations affect individuals independently at a constant rate $\\\ heta$ during their lifetimes, under the infinite-alleles assumption: each new mutation gives a new type, called allele, to his carrier. We study the allele frequency spectrum which is the numbers A(k, t) of types represented by k alive individuals in the population at time t. Thanks to a new construction of the coalescent point process describing the genealogy of individuals in the splitting tree, we are able to compute recursively all joint factorial moments of (A(k, t)) k$\\ge$1. These moments allow us to give an elementary proof of the almost sure convergence of the frequency spectrum in a supercritical splitting tree.

研究の動機と目的

• 中立的ポisson的突然変異を伴う分裂木におけるアレル頻度スぺクトルの階乗モーメントの明示的再帰的公式を導出すること。
• 局所的独立構造を用いたモーメント計算を可能にする、共析点過程（CPP）の新規構成を開発すること。
• 一般分岐過程理論に依存せずに、超臨界状態における頻度スぺクトルのほとんど確実収束を素朴に証明すること。
• 高次漸近近似（例：中心極限定理型結果）に役立つモーメント公式を確立すること。

提案手法

• 共析点過程を独立同分布の確率変数列としての新規表現を導入し、再帰的モーメント計算を可能にする。
• 局所的独立性の下で、確率測度に関する積分の期待値に対する一般化キャンベルの公式を構築する。
• CPPの同時モーメント構造を用いて、頻度スぺクトル $A(k,t)$ の階乗モーメントを再帰的に計算する。
• モーメント公式を応用し、$e^{-\alpha t} A(k,t)$ が $\mathcal{E}/\psi'(\alpha) \cdot c_k$ に比例する極限にほとんど確実に収束することを証明する。ここで $\mathcal{E}$ は指数分布に従う確率変数である。
• 独立のLévy過程を用いた輪郭過程の構成により、残存寿命を特徴付け、CPPにおける独立性を確立する。
• Lévy過程のスケーリング関数とラプラス指数を用いて、集団および突然変異過程の長期的挙動を分析する。

実験結果

リサーチクエスチョン

• RQ1 分裂木におけるポisson的突然変異を伴うアレル頻度スぺクトル $A(k,t)$ の同時階乗モーメントは、どのように再帰的に計算できるか？
• RQ2 共析点過程は、時刻 $t$ における集団の系統的構造をどのように特徴付けるか？
• RQ3 一般分岐過程理論を用いずに、頻度スぺクトルのほとんど確実収束をモーメント法により証明できるか？
• RQ4 超臨界状態における $e^{-\alpha t} A(k,t)$ の正確な極限分布は何か？また、それは集団のマラトゥーシアンパラメータとどのように関係するか？
• RQ5 時刻 $t$ における個体の残存寿命は、共析点過程の独立構造にどのように寄与するか？

主な発見

• 本稿では、時刻 $t$ に $k$ 個の個体が生存するアレルの数 $A(k,t)$ の階乗モーメントについて、低次のモーメントを用いた明示的な再帰的公式を導出する。
• ほとんど確実収束 $e^{-\alpha t} A(k,t) \to \frac{\mathcal{E}}{\psi'(\alpha)} c_k$ が確立され、ここで $\mathcal{E}$ は指数分布に従う確率変数であり、$c_k$ はモデルパラメータに依存する計算可能な定数である。
• 収束速度が指数的であることが示され、$\mathbb{E}_{t}[e^{-2\alpha t}(c_k N_t - A(k,t))^2] \sim \mathcal{C} e^{-\gamma t}$ と表され、ここで $\gamma = \theta$ は亜臨界／臨界状態、$\gamma = 2\alpha - \theta$ は超臨界状態である。
• 時刻 $t$ における個体の残存寿命が互いに独立であり、特に最初の残存寿命が他のものと独立であることが、独立Lévy過程を用いた輪郭過程の構成に基づいて示される。
• 新規の共析点過程構成により、確率的特徴量の理論に依存せずに、モーメント公式の直接的導出が可能になる。
• 本手法は、頻度スぺクトルの収束誤差を定量的に評価することで、高次漸近結果（例：関数中心極限定理）の導出に基盤を提供する。