A Comparison of Some Monte Carlo and Quasi Monte Carlo Techniques for Option Pricing

This article compares the performance of ordinary Monte Carlo and quasi Monte Carlo methods in valuing moderate-and high-dimensional options. The dimensionality of the problems arises either from the number of time steps along a single path or from the number of underlying assets. We compare ordinary Monte Carlo with and without antithetic variates against Sobol’, Faure, and Generalized Faure sequences and three constructions of a discretely sampled Brownian path. We test the standard random walk construction with all methods, a Brownian bridge construction proposed by Caflisch and Morokoff with Sobol’ points and an alternative construction based on principal components analysis also with Sobol’ points. We find that the quasi Monte Carlo methods outperform ordinary Monte Carlo; the Brownian bridge construction generally outperforms the standard construction; and the principal components construction generally outperforms the Brownian bridge construction and is more widely applicable. We interpret both the Brownian bridge and principal components constructions in terms of orthogonal expansions of Brownian motion and note an optimality property of the principal components construction.

[1]  Samuel Karlin,et al.  A First Course on Stochastic Processes , 1968 .

[2]  David M. Kreps,et al.  Martingales and arbitrage in multiperiod securities markets , 1979 .

[3]  J. Harrison,et al.  Martingales and stochastic integrals in the theory of continuous trading , 1981 .

[4]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[5]  Bennett L. Fox,et al.  Algorithm 647: Implementation and Relative Efficiency of Quasirandom Sequence Generators , 1986, TOMS.

[6]  William H. Press,et al.  Numerical Recipes in FORTRAN - The Art of Scientific Computing, 2nd Edition , 1987 .

[7]  Paul Bratley,et al.  Algorithm 659: Implementing Sobol's quasirandom sequence generator , 1988, TOMS.

[8]  J. Hull Options, Futures, and Other Derivatives , 1989 .

[9]  R. Adler An introduction to continuity, extrema, and related topics for general Gaussian processes , 1990 .

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

[11]  D. Duffie Dynamic Asset Pricing Theory , 1992 .

[12]  A. Weron,et al.  Simulation and Chaotic Behavior of Alpha-stable Stochastic Processes , 1993 .

[13]  Joseph F. Traub,et al.  Faster Valuation of Financial Derivatives , 1995 .

[14]  S. Tezuka Uniform Random Numbers: Theory and Practice , 1995 .

[15]  Shu Tezuka,et al.  Uniform Random Numbers , 1995 .

[16]  R. Caflisch,et al.  Modified Monte Carlo Methods Using Quasi-Random Sequences , 1995 .

[17]  H. Niederreiter,et al.  Monte Carlo and Quasi-Monte Carlo Methods in Scientific Computing , 1995 .

[18]  P. Glasserman,et al.  Monte Carlo methods for security pricing , 1997 .

[19]  Mark Broadie,et al.  Connecting discrete and continuous path-dependent options , 1999, Finance Stochastics.