Abstract
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 family, its complexity is uniformly bounded over all ranges of parameters. This new algorithm outperforms all existing schemes. In particular, it is more efficient than the well-known double rejection scheme, which is the only algorithm with uniformly bounded complexity that we can find in the current literature. Beside the ETS family, our scheme is also flexible to be further extended for generating the GTS family, which cannot easily be done by extending the double rejection scheme. Our algorithms are straightforward to implement, and numerical experiments and tests are conducted to demonstrate the accuracy and efficiency.
- J. Abate and W. Whitt. 1992. The Fourier-series method for inverting transforms of probability distributions. Queueing Systems 10, 1 (1992), 5--87.Google ScholarDigital Library
- J. Abate and W. Whitt. 1995. Numerical inversion of Laplace transforms of probability distributions. ORSA Journal on Computing 7, 1 (1995), 36--43.Google ScholarCross Ref
- J. Abate and W. Whitt. 2006. A unified framework for numerically inverting Laplace transforms. INFORMS Journal on Computing 18, 4 (2006), 408--421.Google ScholarDigital Library
- C. Andrieu, A. Doucet, and R. Holenstein. 2010. Particle Markov chain Monte Carlo methods. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 72, 3 (2010), 269--342.Google ScholarCross Ref
- O. E. Barndorff-Nielsen and N. Shephard. 2001. Normal modified stable processes. Theory of Probability and Mathematical Statistics. Online. 1--19.Google Scholar
- O. E. Barndorff-Nielsen and N. Shephard. 2002. Econometric analysis of realized volatility and its use in estimating stochastic volatility models. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 64, 2 (2002), 253--280.Google ScholarCross Ref
- O. E. Barndorff-Nielsen and N. Shephard. 2003. Integrated OU processes and non-Gaussian OU-based stochastic volatility models. Scandinavian Journal of Statistics 30, 2 (2003), 277--295.Google ScholarCross Ref
- N. Batir. 2008. Inequalities for the gamma function. Archiv der Mathematik 91, 6 (2008), 554--563.Google ScholarCross Ref
- A. Brix. 1999. Generalized gamma measures and shot-noise Cox processes. Advances in Applied Probability 31, 4 (1999), 929--953.Google ScholarCross Ref
- P. Carr and L. Wu. 2003. The finite moment log stable process and option pricing. Journal of Finance 58, 2 (2003), 753--777.Google ScholarCross Ref
- A. Dassios, J. W. Lim, and Y. Qu. 2020. Exact simulation of truncated Lévy subordinator. ACM Transactions on Modeling and Computer Simulation 30, 3 (2020), Article 17, 26 pages.Google ScholarDigital Library
- A. Dassios, Y. Qu, and H. Zhao. 2018. Exact simulation for a class of tempered stable and related distributions. ACM Transactions on Modeling and Computer Simulation 28, 3 (2018), Article 20, 21 pages.Google ScholarDigital Library
- L. Devroye. 2009. Random variate generation for exponentially and polynomially tilted stable distributions. ACM Transactions on Modeling and Computer Simulation 19, 4 (2009), 1--20.Google ScholarDigital Library
- S. Favaro, B. Nipoti, and Y. W. Teh. 2015. Random variate generation for Laguerre-type exponentially tilted -stable distributions. Electronic Journal of Statistics 9, 1 (2015), 1230--1242.Google ScholarCross Ref
- J. D. P. Gaver. 1966. Observing stochastic processes, and approximate transform inversion. Operations Research 14, 3 (1966), 444--459.Google ScholarDigital Library
- M. Hofert. 2011a. Efficiently sampling nested Archimedean copulas. Computational Statistics & Data Analysis 55, 1 (2011), 57--70.Google ScholarDigital Library
- M. Hofert. 2011b. Sampling exponentially tilted stable distributions. ACM Transactions on Modeling and Computer Simulation 22, 1 (2011), 1--11.Google ScholarDigital Library
- P. Hougaard. 1986. Survival models for heterogeneous populations derived from stable distributions. Biometrika 73, 2 (1986), 387--396.Google ScholarCross Ref
- A. M. Law. 2015. Simulation Modeling and Analysis (5th ed.). McGraw-Hill Education, New York, NY.Google Scholar
- P. Lévy. 1925. Calcul des Probabilités. Gauthier-Villars, Paris, France.Google Scholar
- L. Li and V. Linetsky. 2013. Optimal stopping and early exercise: An eigenfunction expansion approach. Operations Research 61, 3 (2013), 625--643.Google ScholarCross Ref
- L. Li and V. Linetsky. 2014. Time-changed Ornstein-Uhlenbeck processes and their applications in commodity derivative models. Mathematical Finance 24, 2 (2014), 289--330.Google ScholarCross Ref
- L. Li and V. Linetsky. 2015. Discretely monitored first passage problems and barrier options: An eigenfunction expansion approach. Finance and Stochastics 19, 4 (2015), 941--977.Google ScholarCross Ref
- B. Mandelbrot. 1961. Stable Paretian random functions and the multiplicative variation of income. Econometrica 29, 4 (1961), 517--543.Google ScholarCross Ref
- B. Mandelbrot. 1963a. New methods in statistical economics. Journal of Political Economy 71, 5 (1963), 421--440.Google ScholarCross Ref
- B. Mandelbrot. 1963b. The variation of certain speculative prices. Journal of Business 36, 4 (1963), 394--419.Google ScholarCross Ref
- R. Mendoza-Arriaga and V. Linetsky. 2014. Time-changed CIR default intensities with two-sided mean-reverting jumps. Annals of Applied Probability 24, 2 (2014), 811--856.Google ScholarCross Ref
- R. Mendoza-Arriaga and V. Linetsky. 2016. Multivariate subordination of Markov processes with financial applications. Mathematical Finance 26, 4 (2016), 699--747.Google ScholarCross Ref
- Y. Nishiyama and K. Fukumizu. 2016. Characteristic kernels and infinitely divisible distributions. Journal of Machine Learning Research 17, 1 (2016), 6240--6267.Google ScholarDigital Library
- Y. Qu, A. Dassios, and H. Zhao. 2019. Efficient simulation of Lévy-driven point processes. Advances in Applied Probability 51, 4 (2019), 927--966.Google ScholarCross Ref
- Y. Qu, A. Dassios, and H. Zhao. 2021. Exact simulation of Ornstein-Uhlenbeck tempered stable processes. Journal of Applied Probability 58, 2 (2021), XX--XX.Google ScholarCross Ref
- H. Stehfest. 1970. Algorithm 368: Numerical inversion of Laplace transforms. Communications of the ACM 13, 1 (1970), 47--49.Google ScholarDigital Library
- V. Todorov. 2015. Jump activity estimation for pure-jump semimartingales via self-normalized statistics. Annals of Statistics 43, 4 (2015), 1831--1864.Google ScholarCross Ref
- M. Tweedie. 1984. An index which distinguishes between some important exponential families. In Statistics: Applications and New Directions: Proceedings of the Indian Statistical Institute Golden Jubilee International Conference, J. K. Ghosh and L. Roy (Eds.). Indian Statistical Institute, Calcutta, 579--604.Google Scholar
- L. Wu. 2006. Dampened power law: Reconciling the tail behavior of financial security returns. Journal of Business 79, 3 (2006), 1445--1473.Google ScholarCross Ref
- V. M. Zolotarev. 1966. On representation of stable laws by integrals. In Selected Translations in Mathematical Statistics and Probability. Vol. 6. American Mathematical Society, Providence, RI, 84--88.Google Scholar
- V. M. Zolotarev. 1986. One-Dimensional Stable Distributions. American Mathematical Society, Providence, RI.Google Scholar
Index Terms
- Random Variate Generation for Exponential and Gamma Tilted Stable Distributions
Recommendations
Random variate generation for exponentially and polynomially tilted stable distributions
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 ...
Exact Simulation for a Class of Tempered Stable and Related Distributions
In this article, we develop a new scheme of exact simulation for a class of tempered stable (TS) and other related distributions with similar Laplace transforms. We discover some interesting integral representations for the underlying density functions ...
Quantile-Based Inference for Tempered Stable Distributions
We introduce a simple, fast, and accurate way for the estimation of numerous distributions that belong to the class of tempered stable probability distributions. Estimation is based on the method of simulated quantiles (Dominicy and Veredas in J Econom ...
Comments