Abstract
In this paper, we establish an upper bound for time to convergence to stationarity for the discrete time infinite alleles Moran model. If M is the population size and μ is the mutation rate, this bound gives a cutoff time of log(M μ)/μ generations. The stationary distribution for this process in the case of sampling without replacement is the Ewens sampling formula. We show that the bound for the total variation distance from the generation t distribution to the Ewens sampling formula is well approximated by one of the extreme value distributions, namely, a standard Gumbel distribution. Beginning with the card shuffling examples of Aldous and Diaconis and extending the ideas of Donnelly and Rodrigues for the two allele model, this model adds to the list of Markov chains that show evidence for the cutoff phenomenon. Because of the broad use of infinite alleles models, this cutoff sets the time scale of applicability for statistical tests based on the Ewens sampling formula and other tests of neutrality in a number of population genetic studies.
Original language | English (US) |
---|---|
Pages (from-to) | 189-206 |
Number of pages | 18 |
Journal | Journal of mathematical biology |
Volume | 60 |
Issue number | 2 |
DOIs | |
State | Published - Sep 2010 |
Keywords
- Cutoff phenomena
- Ewens sampling formula
- Extreme value distribution
- Hoppe's urn
- Infinite alleles Moran model
- Lines of descent
- Markov chains
ASJC Scopus subject areas
- Modeling and Simulation
- Agricultural and Biological Sciences (miscellaneous)
- Applied Mathematics