Dynamics of the time to the most recent common ancestor in a large branching population
نویسندگان
چکیده
If we follow an asexually reproducing population through time, then the amount of time that has passed since the most recent common ancestor (MRCA) of all current individuals lived will change as time progresses. The resulting “MRCA age” process has been studied previously when the population has a constant large size and evolves via the diffusion limit of standard Wright-Fisher dynamics. For any population model, the sample paths of the MRCA age process are made up of periods of linear upwards drift with slope +1 punctuated by downwards jumps. We build other Markov processes that have such paths from Poisson point processes on R++ × R++ with intensity measures of the form λ ⊗ μ, where λ is Lebesgue measure and μ (the “family lifetime measure”) is an arbitrary absolutely continuous measure satisfying μ((0,∞)) = ∞ and μ((x,∞)) < ∞ for all x > 0. Special cases of this construction describe the time evolution of the MRCA age in (1+β)-stable continuous state branching processes conditioned on non-extinction – a particular case of which, β = 1, is Feller’s continuous state branching process conditioned on non-extinction. As well as the continuous time process, we also consider the discrete time Markov chain that records the value of the continuous process just before and after its successive jumps. We find transition probabilities for both the continuous and discrete time processes, determine when these processes are transient and recurrent, and compute stationary distributions when they exist. Moreover, we introduce a new family of Markov processes that stands in a relation with respect to the general (1+ β)-stable continuous state branching process and its conditioned version that is similar to the one between the family of Bessel-squared diffusions and the unconditioned and conditioned Feller continuous state branching process.
منابع مشابه
Smaller population size at the MRCA time for stationary branching processes
Abstract. We present an elementary model of random size varying population given by a stationary continuous state branching process. For this model we compute the joint distribution of: the time to the most recent common ancestor, the size of the current population and the size of the population just before the most recent common ancestor (MRCA). In particular we show a natural mild bottleneck ...
متن کاملA population model with non-neutral mutations using branching processes with immigration
We consider a stationary continuous model of random size population with non-neutral mutations using a continuous state branching process with non-homogeneous immigration. We assume the type (or mutation) of the immigrants is random given by a constant mutation rate measure. We determine some genealogical properties of this process such as: distribution of the time to the most recent common anc...
متن کاملThe Time Back to the Most Recent Common Ancestor in Exchangeable Population Models
A class of haploid population models with non-overlapping generations and exchangeable offspring distribution is considered. Based on an analysis of the discrete ancestral process, we present solutions, algorithms and strong upper bounds for the expected time back to the most recent common ancestor. New insights into the asymptotical behavior of the expected time back to the most recent common ...
متن کاملOn the genealogy and coalescence times of Bienaym\'e-Galton-Watson branching processes
Coalescence processes have received a lot of attention in the context of conditional branching processes with fixed population size and nonoverlapping generations. Here we focus on similar problems in the context of the standard unconditional Bienaymé-Galton-Watson branching processes, either (sub)-critical or supercritical. Using an analytical tool, we derive the structure of some counting asp...
متن کاملProbability, Population Genetics and Evolution
Branching process in random environment (BPRE) are a generalization of GaltonWatson processes where the reproduction law is picked randomly in an i.i.d manner in each generation. In the supercritical case, the population Zn either explodes (with a positive probability) or becomes extincted. We are interested in the event when the process Zn takes positive but bounded values for large times n. W...
متن کامل