[論文レビュー] Beta-coalescents when sample size is large

Jonathan A Chetwynd-Diggle, Bjarki Eldon|arXiv (Cornell University)|Jan 13, 2026

Evolution and Genetic Dynamics被引用数 0

ひとこと要約

大規模サンプルサイズ下でのBeta-coalescentを分析し、時間を世代で測定する共alescent ratesとスケーリング領域を導出し、Nとパラメータψ(N)により有効集団サイズが変化することを示す

ABSTRACT

Sweepstakes reproduction refers to a highly skewed individual recruitment success without involving natural selection and may apply to individuals in broadcast spawning populations characterised by Type III survivorship. We consider an extension of the model of sweepstakes reproduction for a haploid panmictic population of constant size $N$; the extension also works as an alternative to the Wright-Fisher model. Our model incorporates an upper bound on the random number of potential offspring (juveniles) produced by a given individual. Depending on how the bound behaves relative to the total population size, we obtain the Kingman coalescent, an incomplete Beta-coalescent, or the (complete) Beta-coalescent. We argue that applying such an upper bound is biologically reasonable. Moreover, we estimate the error of the coalescent approximation. The error estimates reveal that convergence can be slow, and small sample size can be sufficient to invalidate convergence, for example if the stated bound is of the form $N/\log N$. We use simulations to investigate the effect of increasing sample size on the site-frequency spectrum. When the limit is a Beta-coalescent, the site frequency spectrum will be as predicted by the limiting tree even though the full coalescent tree may deviate from the limiting one. When in the domain of attraction of the Kingman coalescent the effect of increasing sample size depends on the effective population size as has been noted in the case of the Wright-Fisher model. Conditioning on the population ancestry (the random ancestral relations of the entire population at all times) may have little effect on the site-frequency spectrum for the models considered here (as evidenced by simulation results).

研究の動機と目的

• 大規模サンプルサイズ下での共alescentモデルの研究動機づけと共alescent近似の誤差源を特定する。
• k個を合併するnブロックの共alescent rate λ_{n,k}の一般形を導出し、Nとψ(N)によるスケーリングを分析する。
• ψ(N)の選択が時間スケールと誤差に与える影響を検討し、KingmanおよびBeta-coalescentsとの関連を説明する。
• 1<α<2の範囲で異なるスケーリング領域が示す有効集団サイズの範囲を特徴づける。

提案手法

• nブロックをk個合併する一般的な共alescent rateを提示: λ_{n,k}=c 1{k=2}+c′ ∫_0^1 1{0<x≤γ} x^{k-2}(1−x)^{n−k} Λ_{+}(dx) で Λ+は0で非原子的。
• ψ1,N=O(1)および ψ1,N=N/logNのとき、N世代に比例する時間スケーリングを分析する。
• ψ(N)の選択と誤差項を無視することが共alescent近似の精度に与える影響を論じる。
• 命題 CNversion0 と Lambda-convergence を用いて誤差源とKingman共alescenceとの結びつきを抽出する。
• 特定の領域で極限過程をKingman共alscenceに関連づけつつ、有効サイズの変化を強調する。
• 1<α<2の場合、有効サイズが cN から c′N^{α−1}(log N)^{2−α} の範囲で変動することを述べ、スケーリングに応じた顕著な変動を説明する。

実験結果

リサーチクエスチョン

• RQ1 大規模サンプル下でk個を合併するλ_{n,k}の形はどのようになるか。
• RQ2 異なる ψ(N) のスケーリングは共alescentの時間スケールと誤差項の大きさにどう影響するか。
• RQ3 どの条件で過程がKingmanまたはBeta-coalescentsに類似し、有効集団サイズがどのように特徴づけられるか。
• RQ4 1<α<2 のとき有効サイズの範囲はどのように現れ、対数因子はこの範囲をどう修正するか。

主な発見

• k個を合併するnブロックの共alescent rateは、γを閾値とする指示関数を含むΛ+への積分を含む一般形を持つ。
• ψ1,N=O(1)および ψ1,N=N/logN のとき、N世代に比例する新しい時間スケーリング領域が現れる。
• 共alsecent近似の誤差は ψ(N) の不適切な上界モデルと無視された誤差項から生じる。
• Kingmanの引力領域にある場合でも（ψ(N)=N/logN）、スケーリングの影響により時間は実質的にBeta-coalescentsに近い形で測定される。
• 1<α<2 の場合、有効サイズは cN から c′N^{α−1}(log N)^{2−α} の範囲にわたり、スケーリングに応じた大きな可変性を反映する。
• Propositions CNversion0 と Lambda-convergence は近似誤差の源と極限共alescents との関係を明らかにする。