Please use this identifier to cite or link to this item:
Lord, Roger
Koekkoek, Remmert
van Dijk, Dick
Year of Publication: 
Series/Report no.: 
Tinbergen Institute Discussion Paper 06-046/4
When using an Euler discretisation to simulate a mean-reverting square root process, one runs into the problem that while the process itself is guaranteed to be nonnegative, the discretisation is not. Although an exact and efficient simulation algorithm exists for this process, at present this is not the case for the Heston stochastic volatility model, where the variance is modelled as a square root process. Consequently, when using an Euler discretisation, one must carefully think about how to fix negative variances. Our contribution is threefold. Firstly, we unify all Euler fixes into a single general framework. Secondly, we introduce the new full truncation scheme, tailored to minimise the upward bias found when pricing European options. Thirdly and finally, we numerically compare all Euler fixes to a recent quasi-second order scheme of Kahl and Jäckel and the exact scheme of Broadie and Kaya. The choice of fix is found to be extremely important. The full truncation scheme by far outperforms all biased schemes in terms of bias, root-mean-squared error, and hence should be the preferred discretisation method for simulation of the Heston model and extensions thereof.
Stochastic volatility
square root process
strong convergence
weak convergence
boundary behaviour
Document Type: 
Working Paper

Files in This Item:
394.21 kB

Items in EconStor are protected by copyright, with all rights reserved, unless otherwise indicated.