Transition probabilities for general birth–death processes with applications in ecology, genetics, and evolution

A birth–death process is a continuous-time Markov chain that counts the number of particles in a system over time. In the general process with n current particles, a new particle is born with instantaneous rate λn and a particle dies with instantaneous rate μn. Currently no robust and efficient method exists to evaluate the finite-time transition probabilities in a general birth–death process with arbitrary birth and death rates. In this paper, we first revisit the theory of continued fractions to obtain expressions for the Laplace transforms of these transition probabilities and make explicit an important derivation connecting transition probabilities and continued fractions. We then develop an efficient algorithm for computing these probabilities that analyzes the error associated with approximations in the method. We demonstrate that this error-controlled method agrees with known solutions and outperforms previous approaches to computing these probabilities. Finally, we apply our novel method to several important problems in ecology, evolution, and genetics.

[1]  G. Yule,et al.  A Mathematical Theory of Evolution, Based on the Conclusions of Dr. J. C. Willis, F.R.S. , 1925 .

[2]  W. Leighton,et al.  Numerical Continued Fractions , 1942 .

[3]  D. Kendall On the Generalized "Birth-and-Death" Process , 1948 .

[4]  William Feller,et al.  An Introduction to Probability Theory and Its Applications , 1951 .

[5]  Karl P. Schmidt,et al.  Principles of Animal Ecology , 1950 .

[6]  F. A. Pitelka,et al.  PRINCIPLES OF ANIMAL ECOLOGY , 1951 .

[7]  S. Karlin,et al.  The differential equations of birth-and-death processes, and the Stieltjes moment problem , 1957 .

[8]  Samuel Karlin,et al.  The classification of birth and death processes , 1957 .

[9]  P. A. P. Moran,et al.  Random processes in genetics , 1958, Mathematical Proceedings of the Cambridge Philosophical Society.

[10]  Samuel Karlin,et al.  Many server queueing processes with Poisson input and exponential service times , 1958 .

[11]  S. Karlin,et al.  LINEAR GROWTH, BIRTH AND DEATH PROCESSES , 1958 .

[12]  Samuel Karlin,et al.  On a Genetics Model of Moran , 1962, Mathematical Proceedings of the Cambridge Philosophical Society.

[13]  G. Blanch,et al.  Numerical Evaluation of Continued Fractions , 1964 .

[14]  Norman T. J. Bailey The Elements of Stochastic Processes with Applications to the Natural Sciences , 1964 .

[15]  David Levin,et al.  Development of non-linear transformations for improving convergence of sequences , 1972 .

[16]  J. A. Murhy,et al.  Some Properties of Continued Fractions with Applications in Markov Processes , 1975 .

[17]  W J Lentz,et al.  Generating bessel functions in mie scattering calculations using continued fractions. , 1976, Applied optics.

[18]  W. Grassmann Transient solutions in Markovian queues : An algorithm for finding them and determining their waiting-time distributions , 1977 .

[19]  Winfried K. Grassmann Transient solutions in markovian queueing systems , 1977, Comput. Oper. Res..

[20]  C. J-F,et al.  THE COALESCENT , 1980 .

[21]  J. Kingman On the genealogy of large populations , 1982 .

[22]  J. Kingman On the genealogy of large populations , 1982, Journal of Applied Probability.

[23]  B. Roehner,et al.  Application of Stieltjes theory for S-fractions to birth and death processes , 1983, Advances in Applied Probability.

[24]  Peter Donnelly,et al.  The transient behaviour of the Moran model in population genetics , 1984, Mathematical Proceedings of the Cambridge Philosophical Society.

[25]  A. R. Barnett,et al.  Coulomb and Bessel functions of complex arguments and order , 1986 .

[26]  William H. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[27]  W. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[28]  M. Ismail,et al.  Linear birth and death models and associated Laguerre and Meixner polynomials , 1988 .

[29]  Virginia Held,et al.  Birth and Death , 1989, Ethics.

[30]  Eric Renshaw Modelling biological populations in space and time , 1990 .

[31]  Ward Whitt,et al.  Numerical inversion of probability generating functions , 1992, Oper. Res. Lett..

[32]  Ward Whitt,et al.  The Fourier-series method for inverting transforms of probability distributions , 1992, Queueing Syst. Theory Appl..

[33]  Sri Gopal Mohanty,et al.  On the transient behavior of a finite birth-death process with an application , 1993, Comput. Oper. Res..

[34]  M. Burrows,et al.  Modelling biological populations in space and time: Cambridge studies in mathematical biology: 11 , 1993 .

[35]  William B. Jones,et al.  A survey of truncation error analysis for Padé and continued fraction approximants , 1993 .

[36]  Haakon Waadeland,et al.  Continued fractions with applications , 1994 .

[37]  Ward Whitt,et al.  Numerical Inversion of Laplace Transforms of Probability Distributions , 1995, INFORMS J. Comput..

[38]  Stephen M. Krone,et al.  Ancestral Processes with Selection , 1997, Theoretical population biology.

[39]  Ward Whitt,et al.  Computing Laplace Transforms for Numerical Inversion Via Continued Fractions , 1999, INFORMS J. Comput..

[40]  Fabrice Guillemin,et al.  Excursions of birth and death processes, orthogonal polynomials, and continued fractions , 1999, Journal of Applied Probability.

[41]  P. Flajolet,et al.  The formal theory of birth-and-death processes, lattice path combinatorics and continued fractions , 2000, Advances in Applied Probability.

[42]  H. Wall,et al.  Analytic Theory of Continued Fractions , 2000 .

[43]  Brian Dennis,et al.  Allee effects in stochastic populations , 2002 .

[44]  Michael Mederer,et al.  Transient solutions of Markov processes and generalized continued fractions , 2003 .

[45]  P. Rousseeuw,et al.  Wiley Series in Probability and Mathematical Statistics , 2005 .

[46]  J. Felsenstein,et al.  An evolutionary model for maximum likelihood alignment of DNA sequences , 1991, Journal of Molecular Evolution.

[47]  R. Sudhesh,et al.  A formula for the coefficients of orthogonal polynomials from the three-term recurrence relations , 2006, Appl. Math. Lett..

[48]  Eugene V. Koonin,et al.  Biological applications of the theory of birth-and-death processes , 2005, Briefings Bioinform..

[49]  P. R. Parthasarathy,et al.  Exact transient solution of a state-dependent birth-death process , 2006 .

[50]  Annie A. M. Cuyt,et al.  Handbook of Continued Fractions for Special Functions , 2008 .

[51]  P. R. Parthasarathy,et al.  ZIPF AND LERCH LIMIT OF BIRTH AND DEATH PROCESSES , 2009, Probability in the Engineering and Informational Sciences.

[52]  Kenneth Dixon,et al.  Introduction to Stochastic Modeling , 2011 .