[論文レビュー] The total external branch length of Beta-coalescents

Götz Kersting, Iulia Stanciu+1 more|arXiv (Cornell University)|Dec 25, 2012

Stochastic processes and statistical mechanics参考文献 14被引用数 10

ひとこと要約

本稿は、$α \in (1,2)$ であるベータコalescentsにおける全外部枝長の漸近的分布を導出し、そのフラクチュエーションが $τ_n^{2-\alpha}$ のスケールで変動することを示している。ここで $τ_n$ は全共役イベントの総数を表す。全木の長さとは異なり、$α_0 = (1+\sqrt{5})/2$ における相転移が見られないため、$1<\alpha<2$ の全範囲で外部枝長のフラクチュエーションは $τ_n^{2-\alpha}$ と一貫して関連している。

ABSTRACT

For $1

研究の動機と目的

• 1<\alpha<2 であるベータコalescentsにおける外部枝長の合計の漸近的挙動を理解すること。
• サンプルサイズ $n \to \infty$ の下での外部枝長のフラクチュエーションスケーリングを特徴づけること。
• 全木の長さが $\alpha_0 = (1+\sqrt{5})/2$ で相転移を示すのとは対照的に、外部枝長のフラクチュエーション挙動を比較すること。
• $1<\alpha<2$ の全パラメータ領域において、外部枝長のフラクチュエーションが $\tau_n^{2-\alpha}$ に従うことを確立すること。

提案手法

• 著者たちは、各共役イベントからの寄与を分解することで全外部枝長 $\ell_n$ を分析し、時刻 $T_k$ における生存する外部枝の数を追跡するプロセス $Y_k$ を用いる。
• 共役プロセス $X_k$ および確率的時間 $\tau_n$ に依存する $\ell_n$ の正規化版 $\mathfrak{l}_n'$ を導入する。
• モーメントの上限と条件付き期待値を用いて、$\ell_n$ とその近似 $\mathfrak{l}_n'$ 間の誤差を制御し、$\mathfrak{l}_n' - \ell_n = O_P(n^{3/2 - \alpha + \varepsilon})$ を示す。
• 共役間の待ち時間が i.i.d. の指数分布確率変数 $W_k$ でモデル化され、$X_k$ および $Y_k$ での条件付き分布を用いて $\ell_n$ の分布的極限を導出する。
• チェビシェフの不等式および支配収束定理を用いて誤差項を抑え、$\ell_n$ と $\gamma^{1-\alpha}\tau_n^{2-\alpha}$ の漸近的同値性を確立する。
• 安定な部分統計量の性質とモーメント推定を用いて、フラクチュエーション項 $G_{k_n} = Y_{k_n} - \mathbb{E}[Y_{k_n} \mid X]$ が $o_P(n)$ であることを示し、正規化外部枝数の確率的収束を保証する。

実験結果

リサーチクエスチョン

• RQ1 1<\alpha<2 であるベータコalescentsにおける全外部枝長のフラクチュエーションは、漸近的にどのように振る舞うか？
• RQ2 外部枝長のフラクチュエーションは、$1<\alpha<2$ の範囲で、全木の長さと同様に振る舞うか？
• RQ3 全共役イベント数 $\tau_n$ に関して、外部枝長のフラクチュエーションの正確なスケーリングは何か？
• RQ4 1<\alpha<2 の場合、$n \to \infty$ のとき正規化外部枝長 $\ell_n / n$ はどのように振る舞うか？
• RQ5 $\tau_n$、すなわち共役時間プロセスを用いて、$\ell_n$ の漸近的分布を特徴づけることができるか？

主な発見

• すべての $1<\alpha<2$ に対して、全外部枝長 $\ell_n$ は $\tau_n^{2-\alpha}$ のスケールでフラクチュエートし、フラクチュエーション挙動に相転移は見られない。
• 正規化外部枝長は、任意の $\varepsilon > 0$ に対して $\ell_n = \gamma^{1-\alpha}\tau_n^{2-\alpha} + o_P(n^{3/2 - \alpha + \varepsilon})$ を満たす。
• フラクチュエーション項 $\mathfrak{l}_n' - \ell_n$ は $O_P(n^{3/2 - \alpha + \varepsilon})$ であり、主項 $\tau_n^{2-\alpha}$ に比べて無視できる。
• $\ell_n$ の漸近的分布は、ランダム時間 $\tau_n$ によって完全に特徴づけられ、外部枝長は $\gamma^{1-\alpha}\tau_n^{2-\alpha}$ に比例してスケーリングする。
• $k_n \sim c\tau_n$ のとき、$\frac{Y_{k_n}}{n} \to (1 - c)^\alpha$ が確率的に収束することを確認し、中間時間における外部枝数の漸近的挙動を裏付ける。
• 条件付き期待値の誤差 $G_{k_n} = Y_{k_n} - \mathbb{E}[Y_{k_n} \mid X]$ は $O_P(\sqrt{n})$ であり、これは $o_P(n)$ であるため、漸近的近似の有効性が保証される。