Coalescence and sampling distributions for Feller diffusions
Abstract
Consider the diffusion process defined by the forward equation $u_t(t, x) = \tfrac{1}{2}\{x u(t, x)\}_{xx} - \alpha \{x u(t, x)\}_{x}$ for $t, x \ge 0$ and $-\infty < \alpha < \infty$, with an initial condition $u(0, x) = \delta(x - x_0)$. This equation was introduced and solved by Feller to model the growth of a population of independently reproducing individuals. We explore important coalescent processes related to Feller's solution. For any $\alpha$ and $x_0 > 0$ we calculate the distribution of the random variable $A_n(s; t)$, defined as the finite number of ancestors at a time $s$ in the past of a sample of size $n$ taken from the infinite population of a Feller diffusion at a time $t$ since since its initiation. In a subcritical diffusion we find the distribution of population and sample coalescent trees from time $t$ back, conditional on non-extinction as $t \to \infty$. In a supercritical diffusion we construct a coalescent tree which has a single founder and derive the distribution of coalescent times.
- Publication:
-
arXiv e-prints
- Pub Date:
- October 2022
- DOI:
- 10.48550/arXiv.2210.12894
- arXiv:
- arXiv:2210.12894
- Bibcode:
- 2022arXiv221012894B
- Keywords:
-
- Mathematics - Probability;
- Quantitative Biology - Populations and Evolution
- E-Print:
- 32 pages, 6 figures. New diagrams have been added and some rewording to sections 1 to 6. Section 7 of the original manuscript contained an error which has necessitated rewriting the original sections 7 to 9 as a new and more straightforward Section 7 which contains new results as a theorem and a corollary