Abstract
A universal generator for integer-valued square-integrable random variables is introduced. The generator relies on a rejection technique based on a generalization of the inversion formula for integer-valued random variables. This approach allows to create a dominating probability function, whose evaluation solely involves two integrals depending on the characteristic function of the random variable to be generated. The proposal gives rise to a simple algorithm which may be implemented in a few code lines and which may show good performance when the classical families of distributions—such as the Poisson and the Binomial—are considered. In addition, applications to the Poisson-Tweedie and the Luria-Delbrück distributions are provided.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.References
Aalen, O.O.: Modelling heterogeneity in survival analysis by the compound Poisson distribution. Ann. Appl. Probab. 2, 951–972 (1992)
Ahrens, J.H., Dieter, U.: A convenient sampling method with bounded computation times for Poisson distributions. In: Nelson, P.R., Dudewicz, E.J., Öztürk, A., van der Meulen, E.C. (eds.) The Frontiers of Statistical Computation, Simulation and Modeling, pp. 137–149. American Sciences Press, Columbus (1991)
Devroye, L.: The computer generation of random variables with a given characteristic function. Comput. Math. Appl. 7, 547–552 (1981)
Devroye, L.: Non-uniform Random Variate Generation. Springer, New York (1986)
Devroye, L.: A simple generator for discrete log-concave distributions. Computing 39, 87–91 (1987)
Devroye, L.: On random variate generation when only moments or Fourier coefficients are known. Math. Comput. Simul. 31, 71–89 (1989)
Devroye, L.: Algorithms for generating discrete random variables with a given generating function or a given moment sequence. SIAM J. Sci. Stat. Comput. 12, 107–126 (1991)
Devroye, L.: A triptych of discrete distributions related to the stable law. Stat. Probab. Lett. 18, 349–351 (1993)
Devroye, L.: Random variate generation for exponentially and polynomially tilted stable distributions. ACM Trans. Model. Comput. Simul. 19, 18 (2009)
Devroye, L.: A note on generating random variables with log-concave densities. Stat. Probab. Lett. 82, 1035–1039 (2012)
El-Shaarawi, A.H., Zhu, R., Joe, H.: Modelling species abundance using the Poisson–Tweedie family. Environmetrics 22, 152–164 (2009)
Feller, W.: An Introduction to Probability Theory and Its Applications, vol. II, 2nd edn. Wiley, New York (1971)
Hörmann, W.: A universal generator for discrete log-concave distributions. Computing 52, 89–96 (1994)
Hörmann, W., Leydold, J., Derflinger, G.: Automatic Nonuniform Random Variate Generation. Springer, Berlin (2004)
Hougaard, P., Lee, M.T., Whitmore, G.A.: Analysis of overdispersed count data by mixtures of Poisson variables and Poisson processes. Biometrics 53, 1225–1238 (1997)
Johnson, N.L., Kemp, A.W., Kotz, S.: Univariate Discrete Distributions, 3rd edn. Wiley, New York (2005)
Lea, D.E., Coulson, C.A.: The distribution of the numbers of mutants in bacterial populations. J. Genet. 49, 264–285 (1949)
Leydold, J.: A simple universal generator for continuous and discrete univariate T-concave distributions. ACM Trans. Math. Softw. 27, 66–82 (2001)
Leydold, J.: Short universal generators via generalized ratio-of-uniforms method. Math. Comput. 72, 1453–1471 (2003)
Mandelbrot, B.: 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 (1974)
Marcheselli, M., Baccini, A., Barabesi, L.: Parameter estimation for the discrete stable family. Commun. Stat., Theory Methods 37, 815–830 (2008)
Ridout, M.S.: Generating random numbers from a distribution specified by its Laplace transform. Stat. Comput. 19, 439–540 (2009)
Stadlober, E.: Binomial random variate generation: a method based on ratio of uniforms. Am. J. Math. Manag. Sci. 9, 1–20 (1989)
Stadlober, E.: The ratio of uniforms approach for generating discrete random variates. J. Comput. Appl. Math. 31, 181–189 (1990)
Zheng, Q.: Progress of a half century in the study of the Luria Delbrück distribution. Math. Biosci. 162, 1–32 (1999)
Zheng, Q.: New algorithms for Luria–Delbrück fluctuation analysis. Math. Biosci. 196, 198–214 (2005)
Zheng, Q.: The Luria-Delbrück distribution. Chance 23, 15–18 (2010)
Zhu, R., Joe, H.: Modelling heavy-tailed count data using a generalized Poisson-inverse Gaussian family. Stat. Probab. Lett. 79, 1695–1703 (2009)
Acknowledgements
The authors would like to thank the two anonymous reviewers for their valuable comments and suggestions which have truly improved the early version of the paper.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Barabesi, L., Pratelli, L. A note on a universal random variate generator for integer-valued random variables. Stat Comput 24, 589–596 (2014). https://doi.org/10.1007/s11222-013-9390-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11222-013-9390-8