Padé approximants and exact two-locus sampling distributions
Abstract
For population genetics models with recombination, obtaining an exact, analytic sampling distribution has remained a challenging open problem for several decades. Recently, a new perspective based on asymptotic series has been introduced to make progress on this problem. Specifically, closed-form expressions have been derived for the first few terms in an asymptotic expansion of the two-locus sampling distribution when the recombination rate $\rho$ is moderate to large. In this paper, a new computational technique is developed for finding the asymptotic expansion to an arbitrary order. Computation in this new approach can be automated easily. Furthermore, it is proved here that only a finite number of terms in the asymptotic expansion is needed to recover (via the method of Padé approximants) the exact two-locus sampling distribution as an analytic function of $\rho$; this function is exact for all values of $\rho\in[0,\infty)$. It is also shown that the new computational framework presented here is flexible enough to incorporate natural selection.
- Publication:
-
arXiv e-prints
- Pub Date:
- July 2011
- DOI:
- 10.48550/arXiv.1107.3897
- arXiv:
- arXiv:1107.3897
- Bibcode:
- 2011arXiv1107.3897J
- Keywords:
-
- Mathematics - Probability;
- Quantitative Biology - Populations and Evolution
- E-Print:
- Published in at http://dx.doi.org/10.1214/11-AAP780 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)