Electronic Journal of Statistics

Statistics for the Luria-Delbrück distribution

Agnès Hamon and Bernard Ycart

Full-text: Open access


The Luria-Delbrück distribution is a classical model of mutations in cell kinetics. It is obtained as a limit when the probability of mutation tends to zero and the number of divisions to infinity. It can be interpreted as a compound Poisson distribution (for the number of mutations) of exponential mixtures (for the developing time of mutant clones) of geometric distributions (for the number of cells produced by a mutant clone in a given time). The probabilistic interpretation, and a rigourous proof of convergence in the general case, are deduced from classical results on Bellman-Harris branching processes. The two parameters of the Luria-Delbrück distribution are the expected number of mutations, which is the parameter of interest, and the relative fitness of normal cells compared to mutants, which is the heavy tail exponent. Both can be simultaneously estimated by the maximum likehood method. However, the computation becomes numerically unstable when the maximal value of the sample is large, which occurs frequently due to the heavy tail property. Based on the empirical probability generating function, robust estimators are proposed and their asymptotic variance is given. They are comparable in precision to maximum likelihood estimators, with a much broader range of calculability, a better numerical stability, and a negligible computing time.

Article information

Electron. J. Statist., Volume 6 (2012), 1251-1272.

First available in Project Euclid: 26 July 2012

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Primary: 92D25: Population dynamics (general)
Secondary: 60J28: Applications of continuous-time Markov processes on discrete state spaces

Luria-Delbrück distribution fluctuation analysis Bellman-Harris branching process probability generating function estimator


Hamon, Agnès; Ycart, Bernard. Statistics for the Luria-Delbrück distribution. Electron. J. Statist. 6 (2012), 1251--1272. doi:10.1214/12-EJS711. https://projecteuclid.org/euclid.ejs/1343310297

Export citation


  • Anderson, T. W. (2003)., An introduction to multivariate statistical analysis, 3rd ed. Wiley, New-York.
  • Angerer, W. P. (2010). Proliferation model dependence in fluctuation analysis: the neutral case., J. Math. Biol. 61 55–93.
  • Armitage, P. (1952). The statistical theory of bacterial populations subject to mutation., J. R. Statist. Soc. B 14 1–40.
  • Athreya, K. B. and Ney, P. E. (1972)., Branching processes. Springer-Verlag, Berlin.
  • Bartlett, M. S. (1978)., An introduction to stochastic processes, with special reference to methods and applications, 3rd ed. Cambridge University Press.
  • Boe, L., Tolker-Nielsen, T., Eegholm, K. M., Spliid, H. and Vrang, A. (1994). Fluctuation analysis of mutations to nalidixic acid resistance in, Escherichia Coli. J. Bacteriol. 176 2781–2787.
  • Dewanji, A., Luebeck, E. G. and Moolgavkar, S. H. (2005). A generalized Luria-Delbrück model., Math. Biosci. 197 140–152.
  • Eldar, Y. C. (2006). Uniformly improving the Cramér-Rao bound and maximum likelihood estimation., IEEE Trans. Signal Proc. 54 2943–2956.
  • Foster, P. L. (2006). Methods for determining spontaneous mutation rates., Methods Enzymol. 409 195–213.
  • Harris, T. E. (1963)., The theory of branching processes. Springer-Verlag, Berlin.
  • Jaeger, G. and Sarkar, S. (1995). On the distribution of bacterial mutants: the effects of differential fitness of mutants and non-mutants., Genetica 96 217–223.
  • Jean, L. W., Suchorolski, M. T., Jeon, J. and Luebeck, E. G. (2010). Multiscale estimation of cell kinetics., Comput. Math. Meth. Med. 11 239–254.
  • Jones, M. E. (1994). Luria-Delbrück fluctuation experiments; accounting simultaneously for plating efficiency and differential growth rate., J. Theo. Biol. 166 355–363.
  • Jones, M. E., Wheldrake, J. and Rogers, A. (1993). Luria-Delbrück fluctuation analysis: estimating the Poisson parameter in a compound Poisson distribution., Comput. Biol. Med. 23 525–534.
  • Kendall, D. G. (1952). Les processus stochastiques de croissance en biologie., Ann. IHP 13 43–108.
  • Kendall, D. G. (1966). Branching processes since 1873., J. London Math. Soc. 41 385–406.
  • Koch, A. L. (1982). Mutation and growth rates from Luria-Delbrück fluctuation tests., Mutat. Res. 95 129–143.
  • Kuczek, T. (1982). Almost sure limit results for the supercritical Bellman-Harris process., J. Appl. Probab 19 668–674.
  • Lea, D. E. and Coulson, C. A. (1949). The distribution of the number of mutants in bacterial populations., J. Genetics 49 264–285.
  • Lehmann, E. L. and Casella, G. (1999)., Theory of point estimation, 2nd ed. Springer, New York.
  • Luria, D. E. and Delbrück, M. (1943). Mutations of bacteria from virus sensitivity to virus resistance., Genetics 28 491–511.
  • Ma, W. T., v. H. Sandri, G. and Sarkar, S. (1992). Analysis of the Luria-Delbrück distribution using discrete convolution powers., J. Appl. Probab. 29 255–267.
  • Mandelbrot, B. (1974). A population birth-and-mutation process, I: explicit distributions for the number of mutants in an old culture of bacteria., J. Appl. Probab. 11 437–444.
  • Marcheselli, M., Baccini, A. and Barabesi, L. (2008). Parameter estimation for the discrete stable family., Commun. Statist. Theory Methods 37 815–830.
  • Möhle, M. (2005). Convergence results for compound Poisson distributions and applications to the standard Luria-Delbrück distribution., J. Appl. Probab. 42 620–631.
  • Neuts, M. F. and Resnick, S. I. (1971). On the times of birth in a linear birth process., Natural Resource Modeling 12 473–475.
  • Oprea, M. and Kepler, T. B. (2001). Improved inference of mutation rates II: generalization of the Luria-Delbrück distribution for realistic cell-cycle time distributions., Theor. Pop. Biol. 59 49–59.
  • Pakes, A. G. (1993). Remarks on the Luria-Delbrück distribution., J. Appl. Probab. 30 991–994.
  • R Development Core Team (2008)., R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria.
  • Reed, W. J. (2006). A brief introduction to scale free networks., Natural Resource Modeling 19 3–14.
  • Rémillard, B. and Theodorescu, R. (2000). Inference based on the empirical probability generating function for mixtures of Poisson distributions., Statist. Decisions 18 349–366.
  • Rosche, W. A. and Foster, P. L. (2000). Determining mutation rates in bacterial populations., Methods 20 1–17.
  • Wilcox, R. (2012)., Introduction to robust estimation and hypothesis testing, 3rd ed. Elsevier, Amsterdam.
  • Wu, X., Strome, E. D., Meng, Q., Hastings, P. J., Plon, S. E. and Kimmel, M. (2009). A robust estimator of mutation rates., Mut. Res. 661 101–109.
  • Yule, G. U. (1925). A mathematical theory of evolution, based on the conclusions of Dr. J.C. Willis, F.R.S., Phil. Trans. Roy. Soc. London Ser. B 213 21–87.
  • Zheng, Q. (1999). Progress of a half century in the study of the Luria-Delbrück distribution., Math. Biosc. 162 1–32.
  • Zheng, Q. (2002). Statistical and algorithmic methods for fluctuation analysis with SALVADOR as an implementation., Math. Biosci. 176 237–252.
  • Zheng, Q. (2005). New algorithms for Luria-Delbrück fluctuation analysis., Math. Biosci. 196 198–214.
  • Zheng, Q. (2008). On Bartlett’s formulation of the Luria-Delbrück mutation model., Math. Biosci. 215 48–54.
  • Zheng, Q. (2010). The Luria-Delbrück distribution: early statistical thinking about evolution., Chance 23 15–18.

Supplemental materials

  • Supplementary material: R script. Set of R functions for statistical computation with Luria-Delbrück distributions, including random sample simulation, estimation of parameters with ML and GF methods, asymptotic variance matrices, confidence intervals and p-values for hypothesis testing.