Abstract
We develop exact random variate generators for the polynomially and exponentially tilted unilateral stable distributions. The algorithms, which generalize Kanter's method, are uniformly fast over all choices of the tilting and stable parameters. The key to the solution is a new distribution which we call Zolotarev's distribution. We also present a novel double rejection method that is useful whenever densities have an integral representation involving an auxiliary variable.
- Aalen, O. O. 1992. Modelling heterogeneity in survival analysis by the compound Poisson distribution. Annals Appl. Probability 2, 951--972.Google ScholarCross Ref
- Abramowitz, M. and Stegun, I. A. 1970. Handbook of Mathematical Tables. Dover Publications, New York, NY.Google Scholar
- Ahrens, J. H. and Dieter, U. 1982. Generating gamma variates by a modified rejection technique. Comm. ACM 25, 47--54. Google ScholarDigital Library
- Barndorff-Nielsen, O. E. and Shepard, N. 2001. Normal modified stable processes. Theor. Probability Math. Statist. 65, 1--19.Google Scholar
- Best, D. J. 1983. A note on gamma variate generators with shape parameter less than unity. Computing 30, 185--188.Google ScholarCross Ref
- Brix, A. 1999. Generalized Gamma measures and shot-noise Cox processes. Advances Appl. Probability 31, 929--953.Google ScholarCross Ref
- Bucklew, J. A. 1990. Large Deviation Techniques in Decision, Simulation and Estimation. John Wiley, New York.Google Scholar
- Bucklew, J. A. 2004. Introduction to Rare Event Simulation. Springer-Verlag, Berlin, Germany. Google ScholarDigital Library
- Cerquetti, A. 2007. On a Gibbs characterization of normalized generalized Gamma processes. arXiv:0707.3408v1, 1705--1711.Google Scholar
- Chambers, J. M., Mallows, C. L., and Stuck, B. W. 1976. A method for simulating stable random variables. J. Amer. Statist. Assoc. 71, 340--344.Google ScholarCross Ref
- Cheng, R. C. H. 1977. The generation of gamma variables with non-integral shape parameter. Appl. Statist. 26, 71--75.Google ScholarCross Ref
- Cheng, R. C. H. and Feast, G. M. 1979. Some simple gamma variate generators. Appl. Statist. 28, 290--295.Google ScholarCross Ref
- Devroye, L. 1984. A simple algorithm for generating random variates with a log-concave density. Computing. 33, 247--257.Google ScholarCross Ref
- Devroye, L. 1986. Non-Uniform Random Variate Generation. Springer-Verlag, Berlin, Germany.Google Scholar
- W. Hörmann. 1994. A universal generator for discrete log-concave distributions. Computing 52, 89--96.Google ScholarCross Ref
- Hörmann, W., Leydold, J., and Derflinger, G. 2004. Automatic Nonuniform Random Variate Generation. Springer-Verlag, Berlin.Google Scholar
- Hougaard, P. 1986. Survival models for hetereneous populations derived from stable distributions. Biometrika 73, 387--396.Google ScholarCross Ref
- James, L. F. 2006. Gamma tilting calculus for GGC and Dirichlet means with applications to Linnik processes and occupation time laws for randomly skewed Bessel processes and bridges. arXiv:math/0610218v3.Google Scholar
- James, L. F. and Yor, M. 2007. Tilted stable subordinators, gamma time changes and occupation time of rays by Bessel spiders. arXiv:math/0701049v1.Google Scholar
- James, L. F. 2007. Lamperti type laws: Positive Linnik, Bessel bridge occupation and Mittag-Leffler functions. arXiv:0708, vol. 0618v1.Google Scholar
- Kanter, M. 1975. Stable densities under change of scale and total variation inequalities. Ann. Prob. 3, 697--707.Google ScholarCross Ref
- Le Minh, D. 1988. Generating gamma variates. ACM Trans. Math. Softw. 14, 261--266. Google ScholarDigital Library
- Leydold, J. 2000. A note on transformed density rejection. Computing, 65, 187--192. Google ScholarDigital Library
- J. Leydold. 2001. A simple universal generator for continuous and discrete univariate T-concave distributions. ACM Trans. Math. Softw. 27, 66--82. Google ScholarDigital Library
- G. Marsaglia. 1977. The squeeze method for generating gamma variates. Comput. Math. Appli. 3, 321--325.Google ScholarCross Ref
- Marsaglia, G. and Tsang, W. W. 2001. A simple method for generating gamma variables. ACM Trans. Math. Softw. 24, 341--350. Google ScholarDigital Library
- Perman, M., Pitman, J., and Yor, M. 1992. Size-biased sampling of Poisson point processes and excursions. Prob. Theor. Relat. Fields 92, 21--39.Google ScholarCross Ref
- Pitman, J. 2003. Poisson Kingman partitions. In Science and Statistics: A Festschrift for Terry Speed. D. R. Goldstein, Ed. Lecture Notes Monograph Series, vol. 30, Institute of Mathematical Statistics, Hayward, CA, 1--34.Google Scholar
- Rosinski, J. 2007a. Tempering stable processes. Stochastic Process. Appl. 117, 677--707.Google ScholarCross Ref
- Rosinski, J. 2007b. Simulation of Lévy processes. Tech. rep., Department of Mathematics, University of Tennessee.Google Scholar
- Sadowsky, J. S. and Bucklew, J. A. 1989. Large deviations theory techniques in Monte Carlo simulation. In Proceedings of the Winter Simulation Conference, E. A. MacNair, K. J. Musselman, and P. Heidelberger, Eds. IEEE, Los Alamitos, CA, 505--513. Google ScholarDigital Library
- Sato, K. 1999. Levy Processes and Infinitely Divisible Distributions. Cambridge University Press, Cambridge, UK.Google Scholar
- Schmeiser, B. W. and Lal, R. 1980. Squeeze methods for generating gamma variates. J. Amer. Stat. Ass. 75, 679--682.Google ScholarCross Ref
- Shohat, J. A. and Tamarkin, J. D. 1943. The Problem of Moments. Math. Surv. 1, American Mathematical Society, New York.Google Scholar
- Stoyanov, J. 2000. Krein condition in probabilistic moment problems. Bernoulli. 6, 939--949.Google ScholarCross Ref
- von Neumann, J. 1963. Various techniques used in connection with random digits. Collected Works, vol. 5, Pergamon Press, 768--770.Google Scholar
- Zolotarev, V. M. 1959. On analytic properties of stable distribution laws. Select. Trans. Math. Stat. Prob. 1, 207--211.Google Scholar
- Zolotarev, V. M. 1966. On the representation of stable laws by integrals. Select. Trans. Math. Stat. Prob. 6, 84--88.Google Scholar
- Zolotarev, V. M. 1981. Integral transformations of distributions and estimates of parameters of multidimensional spherically symmetric stable laws. In Contributions to Probability. Academic Press, 283--305.Google Scholar
- Zolotarev, V. M. 1986. One-Dimensional Stable Distributions. American Mathematical Society, Providence, RI.Google Scholar
Index Terms
- Random variate generation for exponentially and polynomially tilted stable distributions
Recommendations
Random Variate Generation for Exponential and Gamma Tilted Stable Distributions
We develop a new efficient simulation scheme for sampling two families of tilted stable distributions: exponential tilted stable (ETS) and gamma tilted stable (GTS) distributions. Our scheme is based on two-dimensional single rejection. For the ETS ...
Sampling Exponentially Tilted Stable Distributions
Several algorithms for sampling exponentially tilted positive stable distributions have recently been suggested. Three of them are known as exact methods, that is, neither do they rely on approximations nor on numerically critical procedures. One of ...
Simulating Perpetuities
AbstractA perpetuity is a random variable that can be represented as $$1 + W_1 + W_1 W_2 + W_1 W_2 W_3 + \cdot \cdot \cdot ,$$, where the W i's are i.i.d. random variables. We study exact random variate generation for perpetuities and discuss the ...
Comments