Coalescence and sampling distributions for Feller diffusions

Abstract

Consider the diffusion process defined by the forward equation ut(t, x) = 12\x u(t, x)\xx - α \x u(t, x)\x for t, x 0 and -∞ < α < ∞, with an initial condition u(0, x) = δ(x - x0). 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 α and x0 > 0 we calculate the distribution of the random variable An(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 ∞. In a supercritical diffusion we construct a coalescent tree which has a single founder and derive the distribution of coalescent times.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…