Convergence to stationarity in the Moran model

Journal article

Consider a population of fixed size consisting of N haploid individuals. Assume that this population evolves according to the two-allele neutral Moran model in mathematical genetics. Denote the two alleles by A1 and A2. Allow mutation from one type to another and let 0 < γ < 1 be the sum of mutation probabilities. All the information about the population is recorded by the Markov chain X = (X(t))t>0 which counts the number of individuals of type A1. In this paper we study the time taken for the population to 'reach' stationarity (in the sense of separation and total variation distances) when initially all individuals are of one type. We show that after t* = N γ-1 log N + c N the separation distance between the law of X(t*) and its stationary distribution converges to 1 - exp(-γe-γc) as N → ∞. For the total variation distance an asymptotic upper bound is obtained. The results depend on a particular duality, and couplings, between X and a genealogical process known as the lines of descent process.

Authors: Donnelly, P; Rodrigues, E

• Publication status: Published
• Journal: JOURNAL OF APPLIED PROBABILITY
• Volume: 37
• Issue: 3
• Pages: 705-717
• Publication date: 2000-09-01
• DOI: 10.1239/jap/1014842830
• ISSN: 0021-9002
• Language: English
• Keywords: coupling; rate of convergence to stationarity; genealogical processes; cut-off phenomenon