Simulating random variables using moment-generating functions and the saddlepoint approximation

When we are given only a transform such as the moment-generating function of a distribution, it is rare that we can efficiently simulate random variables. Possible approaches such as the inverse transform using numerical inversion of the transform are computationally very expensive. However, the saddlepoint approximation is known to be exact for the Normal, Gamma, and inverse Gaussian distribution and remarkably accurate for a large number of others. We explore the efficient use of the saddlepoint approximation for simulating distributions and provide three examples of the accuracy of these simulations.

[1]  Mark Broadie,et al.  Exact Simulation of Stochastic Volatility and Other Affine Jump Diffusion Processes , 2006, Oper. Res..

[2]  P. Carr,et al.  Saddlepoint methods for option pricing , 2009 .

[3]  S. Rice,et al.  Saddle point approximation for the distribution of the sum of independent random variables , 1980, Advances in Applied Probability.

[4]  D. Cox,et al.  Asymptotic techniques for use in statistics , 1989 .

[5]  N. L. Johnson,et al.  Continuous Univariate Distributions. , 1995 .

[6]  T. Alderweireld,et al.  A Theory for the Term Structure of Interest Rates , 2004, cond-mat/0405293.

[7]  H. Daniels Exact saddlepoint approximations , 1980 .

[8]  Hidetoshi Murakami,et al.  Saddlepoint Approximations to the Limiting Distribution of the Modified Anderson–Darling Test Statistic , 2009, Commun. Stat. Simul. Comput..

[9]  Paul Glasserman,et al.  Monte Carlo Methods in Financial Engineering , 2003 .

[10]  J. Hammersley,et al.  Monte Carlo Methods , 1965 .

[11]  S. Heston A Closed-Form Solution for Options with Stochastic Volatility with Applications to Bond and Currency Options , 1993 .

[12]  L. Devroye Non-Uniform Random Variate Generation , 1986 .

[13]  L. Devroye SIMULATING BESSEL RANDOM VARIABLES , 2002 .

[14]  Brian D. Ripley,et al.  Stochastic Simulation , 2005 .

[15]  Paul Glasserman,et al.  Gamma expansion of the Heston stochastic volatility model , 2008, Finance Stochastics.

[16]  J. Kalbfleisch,et al.  On the Bessel Distribution and Related Problems , 2000 .

[17]  Ronald W. Butler,et al.  Saddlepoint Approximations to the CDF of Some Statistics with Nonnormal Limit Distributions , 1993 .

[18]  M. S. Ridout,et al.  Generating random numbers from a distribution specified by its Laplace transform , 2009, Stat. Comput..

[19]  Xiong Lin,et al.  Inverse transform method for simulating levy processes and discrete Asian options pricing , 2011, Proceedings of the 2011 Winter Simulation Conference (WSC).