English

Coalescence and sampling distributions for Feller diffusions

Probability 2023-09-13 v2 Populations and Evolution

Abstract

Consider the diffusion process defined by the forward equation ut(t,x)=12{xu(t,x)}xxα{xu(t,x)}xu_t(t, x) = \tfrac{1}{2}\{x u(t, x)\}_{xx} - \alpha \{x u(t, x)\}_{x} for t,x0t, x \ge 0 and <α<-\infty < \alpha < \infty, with an initial condition u(0,x)=δ(xx0)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 x0>0x_0 > 0 we calculate the distribution of the random variable An(s;t)A_n(s; t), defined as the finite number of ancestors at a time ss in the past of a sample of size nn taken from the infinite population of a Feller diffusion at a time tt since since its initiation. In a subcritical diffusion we find the distribution of population and sample coalescent trees from time tt back, conditional on non-extinction as tt \to \infty. In a supercritical diffusion we construct a coalescent tree which has a single founder and derive the distribution of coalescent times.

Keywords

Cite

@article{arxiv.2210.12894,
  title  = {Coalescence and sampling distributions for Feller diffusions},
  author = {Conrad J. Burden and Robert C. Griffiths},
  journal= {arXiv preprint arXiv:2210.12894},
  year   = {2023}
}

Comments

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

R2 v1 2026-06-28T04:18:49.737Z