Quasi-Monte Carlo Method for Infinitely Divisible Random Vectors via Series Representations

An infinitely divisible random vector without Gaussian component admits representations of shot noise series. Due to possible slow convergence of the series, they have not been investigated as a device for Monte Carlo simulation. In this paper, we investigate the structure of shot noise series representations from a simulation point of view and discuss the effectiveness of quasi-Monte Carlo methods applied to series representations. The structure of series representations in nature tends to decrease their effective dimension and thus increase the efficiency of quasi-Monte Carlo methods, thanks to the greater uniformity of low-discrepancy sequence in the lower dimension. We illustrate the effectiveness of our approach through numerical results of moment and tail probability estimations for stable and gamma random variables.

[1]  D. Applebaum.  Stable non-Gaussian random processes , 1995, The Mathematical Gazette.

[2]  L. Bondesson.  On simulation from infinitely divisible distributions , 1982, Advances in Applied Probability.

[3]  I. Sobol.  Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates , 2001 .

[4]  Wolfgang Hörmann,et al.  Continuous random variate generation by fast numerical inversion , 2003, TOMC.

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

[6]  Reiichiro Kawai.  An importance sampling method based on the density transformation of Lévy processes , 2006, Monte Carlo Methods Appl..

[7]  Harald Niederreiter,et al.  Random number generation and Quasi-Monte Carlo methods , 1992, CBMS-NSF regional conference series in applied mathematics.

[8]  Kai-Tai Fang,et al.  The effective dimension and quasi-Monte Carlo integration , 2003, J. Complex..

[9]  S. Asmussen,et al.  Approximations of small jumps of Lévy processes with a view towards simulation , 2001, Journal of Applied Probability.

[10]  A. Owen,et al.  Valuation of mortgage-backed securities using Brownian bridges to reduce effective dimension , 1997 .

[11]  Jan Rosiński,et al.  On Series Representations of Infinitely Divisible Random Vectors , 1990 .

[12]  Ken Seng Tan,et al.  Dimension reduction approach to simulating exotic options in a Meixner Levy market , 2009 .

[13]  Art B. Owen,et al.  Latin supercube sampling for very high-dimensional simulations , 1998, TOMC.

[14]  Ken Seng Tan,et al.  An Accelerating Quasi-Monte Carlo Method for Option Pricing Under the Generalized Hyperbolic L[e-acute]vy Process , 2009, SIAM J. Sci. Comput..

[15]  W. Vervaat.  On a stochastic difference equation and a representation of non–negative infinitely divisible random variables , 1979, Advances in Applied Probability.

[16]  T. Ferguson,et al.  A Representation of Independent Increment Processes without Gaussian Components , 1972 .

[17]  C. Houdré,et al.  On layered stable processes , 2005, math/0503742.

[18]  A. Owen.  Randomly Permuted (t,m,s)-Nets and (t, s)-Sequences , 1995 .

[19]  M. Yor,et al.  The Fine Structure of Asset Retums : An Empirical Investigation ' , 2006 .

[20]  M. Yor,et al.  Representing the CGMY and Meixner Lévy processes as time changed Brownian motions , 2008 .

[21]  R. Eiichiro.  An importance sampling method based on the density transformation of Lévy processes , 2008 .

[22]  Robert L. Brown.  An Accelerating Quasi-Monte Carlo Method for Option Pricing Under the Generalized Hyperbolic Lévy Process , 2008 .

[23]  Ken Seng Tan,et al.  A general dimension reduction technique for derivative pricing , 2006 .

[24]  J. Rosínski.  Tempering stable processes , 2007 .

[25]  C. Houdré,et al.  On fractional tempered stable motion , 2005, math/0503741.

[26]  J. Rosínski.  Series Representations of Lévy Processes from the Perspective of Point Processes , 2001 .

[27]  R. Kawai,et al.  On finite truncation of infinite shot noise series representation of tempered stable laws , 2011 .