We present a new numerical method to price vanilla options quickly in time-changed Brownian motion models. The method is based on rational function approximations of the Black-Scholes formula. Detailed numerical results are given for a number of widely used models. In particular, we use the variance-gamma model, the CGMY model and the Heston model without correlation to illustrate our results. Comparison to the standard fast Fourier transform method with respect to accuracy and speed appears to favour the newly developed method in the cases considered. We present error estimates for the option prices. Additionally, we use this method to derive a procedure to compute, for a given set of arbitrage-free European call option prices, the corresponding Black-Scholes implied volatility surface. To achieve this, rational function approximations of the inverse of the Black-Scholes formula are used. We are thus able to work out implied volatilities more efficiently than one can by the use of other common methods. Error estimates are presented for a wide range of parameters.
Abstract. We develop a new Monte Carlo variance reduction method to estimate the expectation of two commonly encountered path-dependent functionals: first-passage times and occupation times of sets. The method is based on a recursive approximation of the first-passage time probability and expected occupation time of sets of a Lévy bridge process that relies in part on a randomisation of the time parameter. We establish this recursion for general Lévy processes and derive its explicit form for mixed-exponential jump-diffusions, a dense subclass (in the sense of weak approximation) of Lévy processes, which includes Brownian motion with drift, Kou's double-exponential model and hyper-exponential jump-diffusion models. We present a highly accurate numerical realisation and derive error estimates. By way of illustration the method is applied to the valuation of range accruals and barrier options under exponential Lévy models and Bates-type stochastic volatility models with exponential jumps. Compared with standard Monte Carlo methods, we find that the method is significantly more efficient.
We develop a new Monte Carlo variance reduction method to estimate the expectation of two commonly encountered path-dependent functionals: first-passage times and occupation times of sets. The method is based on a recursive approximation of the first-passage time probability and expected occupation time of sets of a Levy bridge process that relies in part on a randomisation of the time parameter. We establish this recursion for general Levy processes and derive its explicit form for mixed-exponential jump-diffusions, a dense subclass (in the sense of weak approximation) of Levy processes, which includes Brownian motion with drift, Kou's double-exponential model, and hyperexponential jump-diffusion models. We present a highly accurate numerical realisation and derive error estimates. By way of illustration the method is applied to the valuation of range accruals and barrier options under exponential Levy models and Bates-type stochastic volatility models with exponential jumps. Compared with standard Monte Carlo methods, we find that the method is significantly more efficient.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.