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 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.

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

[2]  Luc Devroye,et al.  Random variate generation for exponentially and polynomially tilted stable distributions , 2009, TOMC.

[3]  Marius Hofert,et al.  Efficiently sampling nested Archimedean copulas , 2011, Comput. Stat. Data Anal..

[4]  P. Hougaard Survival models for heterogeneous populations derived from stable distributions , 1986 .

[5]  N. Shephard,et al.  Normal Modified Stable Processes , 2001 .

[6]  Yan Qu,et al.  Exact Simulation for a Class of Tempered Stable and Related Distributions , 2018, ACM Trans. Model. Comput. Simul..

[7]  Harald Stehfest,et al.  Remark on algorithm 368: Numerical inversion of Laplace transforms , 1970, CACM.

[8]  Neil Shephard,et al.  Integrated OU Processes and Non‐Gaussian OU‐based Stochastic Volatility Models , 2003 .

[9]  Yan Qu,et al.  Exact Simulation of a Truncated Lévy Subordinator , 2020, ACM Trans. Model. Comput. Simul..

[10]  A. Dassios,et al.  Exact simulation of Ornstein–Uhlenbeck tempered stable processes , 2020, Journal of Applied Probability.

[11]  Angelos Dassios,et al.  Efficient simulation of Lévy-driven point processes , 2019, Advances in Applied Probability.

[12]  D. P. Gaver,et al.  Observing Stochastic Processes, and Approximate Transform Inversion , 1966, Oper. Res..

[13]  Liuren Wu,et al.  Dampened Power Law: Reconciling the Tail Behavior of Financial Security Returns , 2004 .

[14]  V. Todorov Jump activity estimation for pure-jump semimartingales via self-normalized statistics , 2015, 1508.04216.

[15]  N. Shephard,et al.  Econometric analysis of realized volatility and its use in estimating stochastic volatility models , 2002 .

[16]  Ward Whitt,et al.  A Unified Framework for Numerically Inverting Laplace Transforms , 2006, INFORMS J. Comput..

[17]  Vadim Linetsky,et al.  Discretely monitored first passage problems and barrier options: an eigenfunction expansion approach , 2014, Finance Stochastics.

[18]  V. Zolotarev One-dimensional stable distributions , 1986 .

[19]  P. Carr,et al.  The Finite Moment Log Stable Process and Option Pricing , 2003 .

[20]  Vadim Linetsky,et al.  TIME‐CHANGED ORNSTEIN–UHLENBECK PROCESSES AND THEIR APPLICATIONS IN COMMODITY DERIVATIVE MODELS , 2012, 1204.3679.

[21]  A. Brix Generalized Gamma measures and shot-noise Cox processes , 1999, Advances in Applied Probability.

[22]  B. Mandelbrot New Methods in Statistical Economics , 1963, Journal of Political Economy.

[23]  S. Favaro,et al.  Random variate generation for Laguerre-type exponentially tilted α-stable distributions , 2015 .

[24]  V. Linetsky,et al.  MULTIVARIATE SUBORDINATION OF MARKOV PROCESSES WITH FINANCIAL APPLICATIONS , 2014 .

[25]  B. Mandelbrot STABLE PARETIAN RANDOM FUNCTIONS AND THE MULTIPLICATIVE VARIATION OF INCOME , 1961 .

[26]  V. Linetsky,et al.  Time-Changed CIR Default Intensities with Two-Sided Mean-Reverting Jumps , 2014, 1403.5402.

[27]  Marius Hofert,et al.  Sampling Exponentially Tilted Stable Distributions , 2011, TOMC.

[28]  B. Mandlebrot The Variation of Certain Speculative Prices , 1963 .

[29]  Kenji Fukumizu,et al.  Characteristic Kernels and Infinitely Divisible Distributions , 2014, J. Mach. Learn. Res..

[30]  A. Doucet,et al.  Particle Markov chain Monte Carlo methods , 2010 .

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

[32]  H. Alzer Inequalities for the gamma function , 1999 .