Monte Carlo Methods and Path-Generation Techniques for Pricing Multi-Asset Path-Dependent Options

We consider the problem of pricing path-dependent options on a basket of underlying assets using simulations. As an example we develop our studies using Asian options. Asian options are derivative contracts in which the underlying variable is the average price of given assets sampled over a period of time. Due to this structure, Asian options display a lower volatility and are therefore cheaper than their standard European counterparts. This paper is a survey of some recent enhancements to improve efficiency when pricing Asian options by Monte Carlo simulation in the Black-Scholes model. We analyze the dynamics with constant and time-dependent volatilities of the underlying asset returns. We present a comparison between the precision of the standard Monte Carlo method (MC) and the stratified Latin Hypercube Sampling (LHS). In particular, we discuss the use of low-discrepancy sequences, also known as Quasi-Monte Carlo method (QMC), and a randomized version of these sequences, known as Randomized Quasi Monte Carlo (RQMC). The latter has proven to be a useful variance reduction technique for both problems of up to 20 dimensions and for very high dimensions. Moreover, we present and test a new path generation approach based on a Kronecker product approximation (KPA) in the case of time-dependent volatilities. KPA proves to be a fast generation technique and reduces the computational cost of the simulation procedure.

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

[2]  C. Loan,et al.  Approximation with Kronecker Products , 1992 .

[3]  Mark Broadie,et al.  Recent advances in simulation for security pricing , 1995, WSC '95.

[4]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

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

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

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

[8]  Art B. Owen,et al.  Monte Carlo extension of quasi-Monte Carlo , 1998, 1998 Winter Simulation Conference. Proceedings (Cat. No.98CH36274).

[9]  C. Loan The ubiquitous Kronecker product , 2000 .

[10]  Fred Espen Benth,et al.  Valuation of Asian Basket Options with Quasi-Monte Carlo Techniques and Singular Value Decomposition , 2001 .

[11]  J. Vecer A new PDE approach for pricing arithmetic average Asian options , 2001 .

[12]  Art B. Owen,et al.  Variance and discrepancy with alternative scramblings , 2002 .

[13]  Julian V. Noble,et al.  The full Monte , 2002, Comput. Sci. Eng..

[14]  Ken Seng Tan,et al.  Enhanced quasi-Monte Carlo methods with dimension reduction , 2002, Proceedings of the Winter Simulation Conference.

[15]  Art B. Owen,et al.  Quasi-Monte Carlo Sampling by , 2003, SIGGRAPH 2003.

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

[17]  J HickernellF,et al.  Implementing Scrambled Digital Nets , 2003 .

[18]  W. Stewart,et al.  The Kronecker product and stochastic automata networks , 2004 .

[19]  R. Rebonato Volatility and correlation : the perfect hedger and the fox , 2004 .

[20]  Ken Seng Tan,et al.  Minimizing Effective Dimension Using Linear Transformation , 2004 .

[21]  Harald Niederreiter,et al.  Monte Carlo and quasi-Monte Carlo methods 2004 , 2006 .

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