MANAGEMENT SCIENCE informs ® Vol. 52, No. 12, December 2006, pp. 1930–1944 doi 10.1287/mnsc.1060.0575 issn 0025-1909 eissn 1526-5501 06 5212 1930 © 2006 INFORMS
Efficient Monte Carlo and Quasi–Monte Carlo Option Pricing Under the Variance Gamma Model
Athanassios N. Avramidis, Pierre L’Ecuyer Département d’Informatique et de Recherche Opérationnelle, Université de Montréal, C.P. 6128, Succ. Centre-Ville, Montréal, H3C 3J7, Canada {[email protected], [email protected]}
e develop and study efficient Monte Carlo algorithms for pricing path-dependent options with the vari- Wance gamma model. The key ingredient is difference-of-gamma bridge sampling, based on the represen- tation of a variance gamma process as the difference of two increasing gamma processes. For typical payoffs, we obtain a pair of estimators (named low and high) with expectations that (1) are monotone along any such bridge sampler, and (2) contain the continuous-time price. These estimators provide pathwise bounds on unbi- ased estimators that would be more expensive (infinitely expensive in some situations) to compute. By using these bounds with extrapolation techniques, we obtain significant efficiency improvements by work reduction. We then combine the gamma bridge sampling with randomized quasi–Monte Carlo to reduce the variance and thus further improve the efficiency. We illustrate the large efficiency improvements on numerical examples for Asian, lookback, and barrier options. Key words: option pricing; efficiency improvement; extrapolation; quasi–Monte Carlo; variance gamma History: Accepted by Wallace J. Hopp, stochastic models and simulation; received October 19, 2004. This paper was with the authors 7 months for 2 revisions.
1. Introduction on the “more detailed” (or continuous-time) sample Madan and Seneta (1990) and Madan et al. (1998) path and approximate the corresponding payoff by introduced the variance gamma (VG) model with extrapolation, with explicit bounds on the approx- application to modeling asset returns and option pric- imation error. We simulate the process paths over ing. The VG process is a Brownian motion with a given time interval by successive refinements via a random time change that follows a stationary bridge sampling truncated after a small number of gamma process. This model permits more flexibility steps and combine this with extrapolation and ran- in modeling skewness and kurtosis relative to Brow- domized quasi–Monte Carlo sampling. nian motion. We developed closed-form solutions for Quasi–Monte Carlo (QMC) methods have attracted European option pricing and provided empirical evi- a lot of interest in finance over the past 10 years dence that this model gives a better fit to market or so (Caflisch and Moskowitz 1995, Moskowitz and option prices than the classical Black-Scholes model. Caflisch 1996, Acworth et al. 1997, Åkesson and One approach to option pricing under VG is based on Lehoczy 2000, Owen 1998, Glasserman 2004, L’Ecuyer characteristic function representations and employs 2004a). Randomized QMC (RQMC) methods replace efficient techniques such as the fast Fourier trans- the vectors of independent uniform random num- form (Carr and Madan 1999); however, it does not bers that drive the simulation by sets or sequences appear to be sufficiently well developed for path- of points that cover the space more uniformly than dependent options. Monte Carlo (MC) methods offer random points, so the variance is reduced, and yet a simple way of estimating option values in that the points are randomized in a way that provides an context (Glasserman 2004). However, these methods unbiased estimator of the integral together with a con- can be very time consuming if a precise estimator is fidence interval. needed. In applications where RQMC really improves on Our aim in this paper is to develop and study ordinary MC, the integral has typically low effec- strategies that improve the efficiency of MC estima- tive dimension in some sense (L’Ecuyer and Lemieux tors in the context of pricing path-dependent options 2002). Roughly, this means that the variance of the under the VG model. Our approach exploits the key integrand depends (almost) only on a small number property that a VG process can be written as the dif- of uniforms, or is a sum of functions that depends ference of two (increasing) gamma processes. By sim- (almost) only on a small number of uniforms. Caflisch ulating a path of these two gamma processes over a and Moskowitz (1995) and Moskowitz and Caflisch coarse time discretization, we can compute bounds (1996) introduced a bridge sampling method for
1930 Avramidis and L’Ecuyer: Efficient Monte Carlo and Quasi–Monte Carlo Option Pricing Management Science 52(12), pp. 1930–1944, © 2006 INFORMS 1931 generating the path of a Brownian motion at discrete until the gap between the two estimators becomes points over a finite time interval: First generate the zero; this procedure yields an unbiased estimator value at the end, then the value in the middle con- whose expected work requirement is often consid- ditional on that at the end, then the values at one erably smaller than full-dimensional path sampling. quarter and at three quarters conditional on the previ- Ribeiro and Webber (2003) also develop and study ous ones, and so on, halving each interval recursively. empirically a method for “correcting for simulation This method rewrites the integrand so as to concen- bias” for lookback and barrier options when estimat- trate the variance to the first few coordinates, and this ing continuous-time prices for models driven by Lévy makes QMC and RQMC much more effective. processes; this method heuristic does not yield error Path-dependent option pricing under the VG model bounds and risks increasing the bias. via QMC with bridge sampling was recently pro- We compare empirically the variance and efficiency posed and studied empirically by Ribeiro and Webber of MC and RQMC for three option pricing exam- (2004) and Avramidis et al. (2003). The algorithms ples on a single asset, namely an arithmetic average proposed there are similar in spirit and structure Asian option, a lookback option, and a barrier option. to the bridge sampling algorithm mentioned for the We experiment with different types of QMC point Brownian motion. Both are based on the premise that sets from those used by Ribeiro and Webber (2004) the bridge sampling methodology improves the effec- and randomize our QMC point sets to obtain unbi- tiveness of QMC by concentrating the variance on the ased estimators of both the mean and variance (which first few random numbers (i.e., by reducing the effec- Ribeiro and Webber do not have). Our bridge sam- tive dimension of the integrand). pling algorithms combined with RQMC improve sim- Our contributions are a theoretical and empirical ulation efficiency by large factors in these examples. analysis of a particular bridge sampling technique The remainder of the paper is organized as fol- named difference-of-gammas bridge sampling (DGBS) lows. Section 2, reviews the VG model and option (Avramidis et al. 2003). We exploit this in two ways: pricing under this model. In §3, we introduce sam- (i) to approximate continuous-time prices, or prices of pling algorithms for the gamma and VG processes. options whose payoffs are based on a large number In §4, we develop and study truncated DGBS, exam- of observation times, by simulating the process only ine special cases, and suggest a few estimators for over a limited number of observation times, and (ii) to each case. In §5, we compare the estimators (in terms improve the effectiveness of RQMC. of variance and efficiency) on a numerical example, The representation of a VG process that underlies with and without truncation and extrapolation, for DGBS yields a pair of estimators (named low and MC and RQMC. In §6, we briefly review Lévy pro- high) whose expectations are monotone along any cesses, generalize the approach to Lévy processes of sequence of paths with increasing resolution. Under finite variation, and state the practical difficulties that monotonicity assumptions on the payoff that hold for arise beyond the VG case. many frequently traded options, if bridge sampling is stopped at a coarser resolution than the one on which 2. Background and Notation the payoff is based, then the low and high estima- tors have negative and positive bias, respectively. By 2.1. The Variance Gamma Process and extrapolating estimators based on different levels of Option Pricing resolution, one can reduce the bias of these inexpen- We provide the relevant background for the VG model sive estimators and improve their efficiency relative to developed by Madan and Seneta (1990), Madan and estimators that compute the exact payoff. (As usual, Milne (1991), and Madan et al. (1998), with emphasis we define efficiency as one over the product of the on pricing path-dependent options via MC sampling mean square error by the expected computing time of methods. The notation X ∼ 2 means that X has the estimator.) the normal distribution with mean and variance 2. We focus on prototypical Asian, lookback, and bar- Gamma denotes the gamma distribution with rier options. For Asian options, we obtain the con- mean and variance 2. vergence rate of the expected gap between the low Let B = Bt = Bt t ≥ 0 be a Brownian and high estimators and develop an estimator, based motion with drift parameter and variance parame- on Richardson’s extrapolation, whose bias has a bet- ter . That is, B0 = 0, the process B has independent ter convergence rate. For the lookback case, we exhibit increments, and Bt + − Bt ∼ 2, for all the low and high estimators and provide empirical t ≥ 0 and >0. Let G = Gt = t t ≥ 0 be evidence that the extrapolation method can improve a gamma process independent of B, with drift >0 both convergence of the bias and efficiency. For and volatility >0. That is, G0 = 0, the process G the barrier case, we exhibit low and high estima- has independent increments, and Gt + − Gt ∼ tors and develop bridge sampling that continues Gamma2/ / for t ≥ 0 and >0. Avramidis and L’Ecuyer: Efficient Monte Carlo and Quasi–Monte Carlo Option Pricing 1932 Management Science 52(12), pp. 1930–1944, © 2006 INFORMS
A VG process with parameters is defined 3. Bridge Sampling for Gamma and by Variance Gamma Processes X = Xt=Xt =BGt 1 t≥0 (1) 3.1. Sampling the Brownian Motion In simple terms, the VG process is obtained by sub- A Brownian motion B with parameters can be jecting the Brownian motion to a random time change sampled at arbitrary discrete times 0 = t0 ESt = S0 expr − qt (3) sampled at discrete times 0 = t0