Numerical Analysis of Monte Carlo Evaluation of Greeks by Finite Differences

An error analysis of approximation of deltas (derivatives of the solution to the Cauchy problem for parabolic equations) by finite differences is given, taking into account that the value of the hedging portfolio itself (the solution of the problem) is evaluated using weak-sense numerical integration of the corresponding system of stochastic differential equations together with the Monte Carlo technique. It is shown that finite differences are effective when the method of dependent realizations is exploited in the Monte Carlo simulations. Evaluation of other Greeks is also considered. Results of numerical experiments, including those with Heston stochastic volatility model, are presented.

[1]  P. Glynn.  Optimization of stochastic systems via simulation , 1989, WSC '89.

[2]  Nigel J. Newton.  Numerical Methods in Finance: Continuous-Time Monte Carlo Methods and Variance Reduction , 1997 .

[3]  G. Milstein,et al.  Numerical construction of a hedging strategy against the multi-asset European claim , 1999 .

[4]  M. Musiela,et al.  Martingale Methods in Financial Modelling , 2002 .

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

[6]  Pierre-Louis Lions,et al.  Applications of Malliavin calculus to Monte Carlo methods in finance , 1999, Finance Stochastics.

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

[8]  W. Wagner.  Monte carlo evaluation of functionals of solutions of stochastic differential equations. variance reduction and numerical examples , 1988 .

[9]  M. V. Tretyakov,et al.  Stochastic Numerics for Mathematical Physics , 2004, Scientific Computation.

[10]  Pierre L'Ecuyer,et al.  On the Convergence Rates of IPA and FDC Derivative Estimators , 1990, Oper. Res..

[11]  P. Glasserman,et al.  Estimating security price derivatives using simulation , 1996 .

[12]  P. Glasserman,et al.  Some Guidelines and Guarantees for Common Random Numbers , 1992 .

[13]  Uwe Wystup,et al.  Efficient computation of option price sensitivities using homogeneity and other tricks , 2005 .

[14]  D. Talay,et al.  The law of the Euler scheme for stochastic differential equations , 1996 .

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

[16]  John Schoenmakers,et al.  Transition density estimation for stochastic differential equations via forward-reverse representations , 2004 .

[17]  D. Talay,et al.  Expansion of the global error for numerical schemes solving stochastic differential equations , 1990 .

[18]  P. Glynn,et al.  Efficient Monte Carlo Simulation of Security Prices , 1995 .

[19]  A. Friedman.  Partial Differential Equations of Parabolic Type , 1983 .