arxiv: v2 [q-fin.pr] 12 Oct 2013

Transcription

1 Lower Bound Approxmaton to Basket Opton Values for Local Volatlty Jump-Dffuson Models Guopng Xu and Harry Zheng arxv: v2 [q-fn.pr 12 Oct 213 Abstract. In ths paper we derve an easly computed approxmaton to European basket call prces for a local volatlty jump-dffuson model. We apply the asymptotc expanson method to fnd the approxmate value of the lower bound of European basket call prces. If the local volatlty functon s tme ndependent then there s a closed-form expresson for the approxmaton. Numercal tests show that the suggested approxmaton s fast and accurate n comparson wth the Monte Carlo and other approxmaton methods n the lterature. Keywords. Basket optons valuaton, local volatlty jump-dffuson model, lower bound approxmaton, second order asymptotc expanson. 1 Introducton It s n general dffcult to value basket optons due to lack of analytc characterzaton of the dstrbuton of the underlyng basket asset prce process. Research s manly focused on developng fast and accurate approxmaton methods and fndng tght lower and upper bounds for basket opton values. In the Black-Scholes settng Curran (1994 and Rogers and Sh (1995 derve a lower bound for Asan optons by the condtonng random varable and Jensen s nequalty. Deelstra et al. (24 obtan the bounds for basket optons wth the comonotoncty approach. In affne Lévy models one can derve lower bounds numercally for arthmetc Asan optons based on the characterstc functon and the method developed n Duffe et al. (2, see Albrecher et al. (28. Deelstra et al. (21 provde a good overvew of the recent development n fndng and computng the bounds. It s known that Rogers and Sh s lower bound s generally tght and s one of the most accurate approxmatons to basket call opton values. The lower bound can be calculated exactly n the Black-Scholes framework. Xu and Zheng (29 show that the lower bound can also be calculated exactly n a specal jump-dffuson model wth constant volatlty and two types of Posson jumps (systematc and dosyncratc jumps. The usefulness of Rogers and Sh s lower bound depends crucally on one s ablty of fndng some hghly correlated random varables to the basket value and computng the condtonal expectaton exactly. It s dffcult to extend Rogers and Sh s lower bound to more general models such as local volatlty models due to lack of explctly known dstrbutons for models wth non-constant volatltes, see Albrecher et al. (28. To the best of our knowledge, Rogers and Sh s lower bound for models wth local volatltes has not been dscussed n the lterature. In ths paper we am to fnd a good approxmaton to Rogers and Sh s lower bound for a local volatlty jump-dffuson model and use t to approxmate European basket call opton values. We frst apply the second order asymptotc expanson (see Benhamou et al. (29 to approxmate the basket asset value, then choose a normal varable and a Posson varable as Rsk Analytcs, Ctgroup Centre, Canary Wharf, London E14 5LB, UK. Emal: guopng.xu@ct.com Department of Mathematcs, Imperal College, London SW7 2AZ, UK. Emal: h.zheng@mperal.ac.uk 1

2 condtonng varables whch are hghly correlated to the basket asset value, and fnally apply the condtonal expectaton results of multple Wener-Itô ntegrals (see Kuntomo and Takahash (21 to approxmate Rogers and Sh s lower bound. The man contrbuton of the paper s the dervaton of an approxmaton to Rogers and Sh s lower bound for local volatlty jump-dffuson models. We suggest an easly mplemented algorthm to compute the lower bound approxmaton. If the local volatlty functon s tme ndependent then there s a closed-form expresson for the approxmaton. Numercal tests show that the lower bound approxmaton s fast and accurate n most cases n comparson wth the Monte Carlo method, the partal-exact approxmaton (Xu and Zheng (29 and the asymptotc expanson approxmaton (Xu and Zheng (21. The paper s organzed as follows. Secton 2 formulates the basket local volatlty jumpdffuson model and explans some known methods n prcng European basket call optons. Secton 3 apples the second order asymptotc expanson to derve an easly computed approxmaton to Roger and Sh s lower bound. Secton 4 compares the numercal performance of the lower bound approxmaton wth other methods n the lterature. Secton 5 concludes. 2 Model Assume (Ω,P,F,F t s a fltered rsk-neutral probablty space and F t s the augmented natural fltraton generated by correlated Brownan motons W 1,...,W n wth correlaton matrx (ρ j and a Posson process N wth ntensty λ. Assume Brownan motons W 1,...,W n and Posson process N are ndependent to each other. Assume the rskfree nterest rate s zero for the sake of clarty (otherwse one may smply consder the dscounted asset prces and the dscounted strke prce. Assume the portfolo s composed of n assets wth prces S 1,...,S n that are martngales satsfyng the followng stochastc dfferental equaton ds (t = σ (t,s (tdw (t+h S (t dm(t, (1 where h > 1 are constant jump szes (percentage changes of the values at the jump tme of N, σ (t,s are the local volatlty functons, and M(t = N(t λt s the compensated Posson process. The basket value at tme T s gven by S(T = w S (T, =1 where w are postve constant weghts. The European basket call opton prce at tme s gven by C(T,K = E[(S(T K +. Almost all works n the lterature on Asan or basket optons valuaton assume the underlyng asset prces follow lognormal processes, whch corresponds to σ (t,s = σ S wth σ beng postve constant and h = for all n (1. Snce Asan optons are smlar to basket optons, we do not dfferentate these two types of optons, even though some technques are orgnally developed for Asan optons. We now revew some well-known approaches n approxmaton and error bound estmaton for the Black-Scholes case. Levy (1992 approxmates the basket value S(T wth a lognormal varable whch matches the frst two moments of S(T and derves an approxmate closed-form prcng formula for C(T, K. The result s good when maturty T and volatltes σ are relatvely small. The performance deterorates as T or σ ncreases. Curran (1994 ntroduces the dea of condtonng varables. 2

3 Assume Λ s a random varable whch has strong correlaton wth S(T and satsfes that S(T K whenever Λ d Λ for some constant d Λ. The basket opton prce can be decomposed as E[(S(T K + = E[(S(T K1 [Λ dλ +E[(S(T K + 1 [Λ<dΛ. (2 Curran (1994 chooses the geometrc average of the asset prces as the condtonng varable Λ (a lognormal varable and fnds the closed-form expresson for the the frst term n (2 and uses the lognormal varable and the condtonal moment matchng to fnd the approxmate value of the second term n (2. Rogers and Sh (1995 use the condtonng varable Λ to show that the lower and upper bounds of the basket opton prce C(T,K are gven by LB = E [ (E[S(T Λ K + and UB = LB E[ 1 var(s(t Λ1 2 [Λ<dΛ E[1 [Λ<dΛ 1 2. These bounds may be computed analytcally for specfc condtonng varables Λ. Numercal tests show that Rogers and Sh s lower bound s normally very close to the true (Monte Carlo value. Xu and Zheng (29 extend the Black-Scholes model of Curran (1994 to the jump-dffuson model (1 wth local volatlty functons σ (t,s = σ S for all. (The jump part of the model n Xu and Zheng (29 s more general than (1 wth both common and ndvdual jumps. Asset prces S (T at tme T have closed form expressons ( S (T = S (exp 1 2 σ2 T h λt +σ W (T+N(Tln(1+h, (3 where S ( are asset prces at tme, for = 1,...,n. We can wrte the basket prce S(T as S(T = a exp(σ W (T+N(Tln(1+h, =1 where a = w S (exp(( 1 2 σ2 h λt, = 1,...,n. Snce exp(x 1+x for all x and a (weghts w nonnegatve the basket prce S(T satsfes S(T a (1+σ W (T+N(Tln(1+h = c+mn(t+σw, =1 where c = n =1 a, m = n =1 a ln(1+h, σ 2 = n,j=1 a a j ρ j σ σ j T, N(T s a Posson varable wth parameter λt, andw = 1 σ n =1 a σ W (T s a standardnormal varable and s ndependent of N(T. Note that the above nequalty holds only when all weghts w are nonnegatve, whch excludes spread optons or basket optons wth negatve weghts. If we choose a condtonng varable by Λ = (N(T,W then S(T K whenever c+mn(t+ σw K. t s easy to fnd E[S(T Λ = (k,y = =1 ( 1 a exp 2 (σ2 T R+R 2 y +kln(1+h, (4 where R = 1 σ n j=1 a jρ j σ σ j T. The lower bound can be computed from the expresson LB = p k (E[S(T Λ = (k,y K + dφ(y, k= 3

4 where p k = P(N(T = k = exp( λt(λt k /k! and Φ s the cumulatve dstrbuton functon of a standard normal varable. Xu and Zheng (29 suggest an approxmaton, called the partal exact approxmaton (PEA, to the basket call opton value C(T, K, gven by C A (T,K := E[(S(T K + 1 [c+mn(t+σw K 3 + q E[(E[S(T Λ+α K + 1 [c+mn(t+σw<k, =1 where q 1 = 1/6, q 2 = 2/3, q 3 = 1/6, and α 1 = 3ǫ, α 2 =, α 3 = 3ǫ, ǫ s some constant dependng on the condtonng varable Λ and the condtonal second moment of S(T gven Λ. The lower bound plays a domnant role n the approxmaton wth a weght 2/3, the other two parts wth a weght 1/6 each are the second moment adjustment to the lower bound. Xu and Zheng (29 show that LB C A (T,K UB, whch provdes the error bounds for the approxmate basket opton value. The PEA method s fast and accurate n comparson wth some other well known numercal schemes for basket optons, see Xu and Zheng (29 for detals. A key condton for the PEA method to work s that one must know the closed form expressons of asset prces S (T as n (3 for = 1,...,n. Ths s mpossble for general local volatlty functons. The PEA method can be easly modfed to nclude the possblty of defaults of some underlyng names n the portfolo. For example, assume that company n defaults whenever there s a jump event, then we may set h n = 1, whch mples S n (T = f N(T = k 1. The formula for condtonal expectaton (4 s stll vald but we need to dfferentate the cases of no jump (k = and at least one jump (k 1. In the former we can use the same σ 2 and c as above wth the conventon exp(ln = 1. In the latter, due to S n (T = we need to remove terms wth ndex n n all expressons, e.g., m = n 1 =1 a ln(1+h and σ 2 = n 1,j=1 a a j ρ j σ σ j T. XuandZheng(21 dscussthebasketoptonsprcngforasset model(1wththesamejump szes h = h for = 1,...,n. Snce the PEA method cannot be appled for general local volatlty functons σ (t,s, Xu and Zheng(21 frstreducethe dmensonalty of the problem by drectly workng on the portfolo asset value S(t whch satsfes the followng stochastc volatlty jump dffuson model ds(t = V(tdW(t+hS(t dm(t wth the ntal prce S( = n =1 w S (, wherew s a standard Brownan moton, ndependent of M, and V(t 2 := w w j σ (t,s (tσ j (t,s j (tρ j.,j=1 The European basket call opton prce C(T, K at tme satsfes a partal ntegral dfferental equaton (PIDE C T (T,K = λhkc K (T,K+ 1 ( K 2 σ(t,k2 C KK (T,K+λ(h+1 C(T, h+1 C(T,K wth the ntal condton C(,K = (S( K + and the local varance functon σ(t,k 2 = E[V(T 2 S(T = K. The man dffculty s to compute the condtonal expectaton as asset prces S (T have no closed-form expressons. Followng the frst order asymptotc expanson methods of Kuntomo 4

5 and Takahash (21 and Benhamou et al. (29, Xu and Zheng (21 show that the unknown local varance functon σ(t,k 2 can be approxmated by σ(t,k 2 a(t+b(t(k S(, (5 where a(t = n,j=1 w w j ρ j p p j, b(t = ( n,j=1 1 σ w c 2 w j ρ j p p q j p C + q j p j C j, p = σ (T,S (, q = S ( σ (T,S (, C = [ n j=1 w j ρ j ( T σ (t,s (σ j (t,s j (dt, and σc 2 = n =1 w C j. One can fnd the approxmate basket opton prce by solvng the PIDE wth the mplct-explct fnte dfference method. Ths approxmaton, called the asymptotc expanson approxmaton (AEA, s less accurate than the PEA when the local volatlty functon s the Black-Scholes type, but s capable of dealng wth general local volatlty functons. The man lmtatons of the AEA method are that the approxmate varance functon n (5 s not always postve, whch can cause numercal errors n solvng the PIDE, and that t requres the same jump sze for all assets, whch s unrealstc for a general asset portfolo. 3 Lower Bound Approxmaton The asymptotc expanson for basket optons prcng n a general dffuson model s dscussed n Kuntomo and Takahash (21 n whch the valuaton of condtonal expectatons s a necessary step to obtan the characterstc functon of the basket value. In ths paper we use the same method to expand the parameterzed asset prce processes to the second order, see Benhamou et al. (29, and apply the condtonal expectaton results of multple Wener-Itô ntegrals drectly to approxmate Rogers and Sh s lower bound for the basket call opton values n a jump-dffuson model (1. For ǫ [,1 defne ds ǫ (t = ǫσ (t,s ǫ (tdw (t+ǫh S ǫ (t dm(t (6 wth ntal condton S ǫ ( = S (. Note that S 1 (T = S (T. Defne S (k (t := k S ǫ (t ǫ k ǫ= and σ (k (t := k σ (t,s ǫ (t (S ǫ k ǫ= fork =,1,... Itsobvous thats ( (t = S (t = S ( andσ ( (t = σ (t,s ( forforall t and = 1,...,n. In partcular, f local volatlty functons σ are tme ndependent then functons are all constants for = 1,...,n, whch may smplfy consderably some computatons. An mportant case s the constant elastcty of varance (CEV model, where σ (t,s = α S β for = 1,...,n. To smplfy the notatons n subsequent dscussons we denote by σ (k σ (k (t = j=1 w j σ (k (tσ ( j (tρ j, k =,1. The second order asymptotc expanson around ǫ = for S ǫ (t s Expand (6 to the second order, we have S ǫ (T S ( (T+S (1 (Tǫ+ 1 2 S(2 (Tǫ 2. ds (1 (t = σ ( (tdw (t+h S (dm(t, ds (2 (t = 2σ (1 (ts (1 (tdw (t+2h S (1 (t dm(t, 5

6 wth the ntal condtons S (1 ( = S (2 ( =. Therefore, S (1 (t = λh S (t+ S (2 (t = 2λh t t σ ( S (1 (sds+2 (sdw (s+h S (N(t (7 t σ (1 (ss (1 (sdw (s+2h t for t T. Lettng ǫ = 1 we may approxmate the basket value S(T by S (1 (s dn(s (8 S(T S A (T := S(+S (1 (T+ S(2 (T, (9 2 where S (j (T := n =1 w S (j (T for j = 1,2. Snce there are no closed-form expressons for S (T and S(T, t s dffcult to compute the lower bound. If we approxmate S(T by S A (T, defned n (9, we may be able to compute the condtonal expectaton E[S A (T Λ for some condtonng varable Λ and then to approxmate the lower bound. We therefore propose that LB LBA := E[(E[S A (T Λ K +. (1 Ths approxmaton s called the lower bound approxmaton (LBA of the basket call opton prce. Note that the LBA s only an approxmaton to Rogers and Sh s lower bound and s therefore possble to have values greater than the basket call opton values. The next step s to choose the condtonng varable for the approxmaton. In the Black- Scholes framework and for European arthmetc average optons, t s wthout excepton to choose the geometrc average of asset prces as the condtonng varable. Ths choce s natural as the geometrc average of lognormal varables s agan a lognormal varable and provdes good nformaton on the arthmetc average. It s not clear what one should choose for general jump-dffuson models as there are no closed form solutons for underlyng asset prces. Wth the nsght from Xu and Zheng (29 we choose Λ = (N(T, (T for the model, where (T = n T σ( j (tdw j (t s a normal varable wth mean and varance v 2 = n =1 w and (9 we have LBA = [( p k k= j=1 w j T σ( (tdt, and N(T s a Posson varable wth parameter λt. From (1 S(+E[S (1 (T Λ = (k,vx+e[ S(2 (T Λ = (k,vx K 2 + dφ(x, wherep k = exp( λt(λt k /k!. SnceS (j (T = n =1 w S (j (T, weonlyneedtofnde[s (j (T Λ = (k,vx for j = 1,2 and = 1,...,n. We frst calculate E[S (1 (T Λ = (k,vx. E[S (1 (T Λ = (k,vx = λh S (T +E[h S (N(T N(T = k+e = λh S (T +h S (k + 1 v ( T σ ( (tdt x. σ ( (tdw (t (T = vx The valuaton of E[ S(2 (T 2 Λ = (k,vx s more nvolved. Snce S (2 (T s the sum of three terms n (8 we may fnd the condtonal expectaton of each term by substtutng S (1 (t n (7 6 (11

7 nto the ntegrands and then computng three condtonal expectatons. We now derve them one by one. The frst term can be wrtten as ( λh E S (1 (tdt Λ = (k,vx = A 1 +A 2 +A 3, where A 1 = ( λh E ( λh S (tdt Λ = (k,vx = 1 2 S (λ 2 h 2 T2 t A 2 = ( λh E ( σ ( (sdw (sdt Λ = (k,vx = ( λh E (T tσ ( (tdw (t (T = vx = ( λh 1 ( T (T t σ ( (tdt x v A 3 = ( λh E (h S (N(tdt Λ = (k,vx T = ( λs (h 2 E[N(t N(T = kdt = ( λs (h 2 (kt 2. Here we have used the fact that (N(t N(T = k s a bnomal varable wth k ndependent -1 trals and probablty λt λt = t kt T of takng 1, whch mples E[N(t N(T = k = T. The second term can be wrtten as E where σ (1 (ts (1 (tdw (t Λ = (k,vx = B 1 +B 2 +B 3, B 1 = E σ (1 (t( λh S (tdw (t Λ = (k,vx ( T = ( λh S ( t σ (1 1 (tdt v x ( t B 2 = E σ (1 (t σ ( (sdw (s = 1 ( T ( t v 2 σ ( (sds σ (1 (tdt (x 2 1 B 3 = E σ (1 (t(h S (N(tdW (t Λ = (k,vx = h S (E σ (1 (t kt T dw (t (T = vx = h S ( k ( T t σ (1 1 (tdt T v x. dw (t Λ = (k,vx The computaton of B 2 s dscussed n Kuntomo and Takahash (21, Lemma A.1. The thrd term can be wrtten as h E S (1 (t dn(t Λ = (k,vx = C 1 +C 2 +C 3, 7

8 where ( C 1 = h E λh S (t dn(t Λ = (k,vx = 1 2 S (λh 2 Tk ( t C 2 = h E σ ( (sdw (s dn(t Λ = (k,vx T ( t = h E[ σ ( k (sdw (s dt (T = vx T ( k T = h (T t σ ( 1 (tdt T v x ( C 3 = h E h S (N(t dn(t Λ = (k,vx = h 2 S (E [N(T 1 l= = h 2 S ( k2 k. 2 l N(T = k Substtutng the frst and second order condtonal expectatons nto (11 we get the lower bound approxmaton as ( + LBA = p k cx 2 +a 1 (kx+a (k dφ(x, (12 where k= c = 1 ( T ( t v 2 w =1 a (k = S(+(K λt a 1 (k = v + 1 v σ ( (sds σ (1 (tdt w h S (+ 1 2 =1 w h ( k ( T T λ =1 w S (h 2 ( (k λt 2 k c K =1 ( (T t σ ( (t+s (t σ (1 (t dt Note that n computng the lower bound approxmaton (12 there s no need to do numercal ntegraton. Ths s because for every fxed k the ntegrand s the postve part of a quadratc functon whch contans only three cases: no root, one root, and two roots. Therefore the computaton of the lower bound approxmaton LBA s easy and fast. Ths s one key advantage over the other methods such as smulaton or fnte dfference for the PIDE.. 4 Numercal Results In ths secton we conduct some numercal tests for European basket calls wth underlyng asset prce processes satsfyng (1 to test the performance of the lower bound approxmaton. The Monte Carlo smulaton provdes the benchmark results. The control varate technque s adopted to reduce the standard devatons. The control varate used s the frst order asymptotc expanson of the basket prce S ( T, whch results n a closed form prcng formula, see Xu and Zheng (21 for detals. The followng data are used n all numercal tests: the number of assets n the basket 8

9 n = 4, the portfolo weghts of each asset w =.25 for = 1,...,n, the correlaton coeffcents of Brownan motons ρ j =.3 for,j = 1,...,n, the ntal asset prces S ( = 1 for = 1,...,n, the exercse prce K = 1, and the jump ntensty λ =.3. The number of smulatons for each test case s 1, and the Posson summaton s approxmated by truncatng the nfnte seres after the frst 1 terms. Tables 1 dsplays the European basket call opton values and ther mpled volatltes wth four dfferent methods: the Monte Carlo (MC, the partal exact approxmaton (PEA, the asymptotc expanson (AEA and the lower bound approxmaton (LBA. The local volatlty functon s σ(t,s =.2S. We perform numercal tests for three constant jump szes h = e η 1 wth η =.25,.125 and.625. It s clear that all three approxmaton methods perform well wth relatve errors less than 1%. We have done other tests wth the same data except the volatlty functon beng changed to σ(t,s =.5S. The overall relatve error s.6% for the PEA, 4% for the AEA and 1.7% for the LBA. Tables 1 shows that the performance of the PEA s the best whle the LBA performs better than the AEA. The results suggest that the PEA s the best approxmaton method for the European basket call opton valuaton when local volatlty functons are of the Black-Scholes type. On the other hand, the AEA and LBA are much more flexble and can handle general local volatlty functons. The LBA has the addtonal advantage of havng closed-form solutons when the local volatlty functons are tme ndependent and can deal wth dfferent jump szes to common jumps. The AEA s slow n solvng the PIDE wth the fnte dfference method and requres the same jump sze to common jumps for all assets. It s also nterestng to note that the mpled volatltes are essentally same for T = 1 and T = 3 wth all other parameters beng the same. It s not yet clear to us the reason of nsenstvty of the mpled volatlty wth respect to the tme to maturty. Tables 2 dsplays the European basket call opton values and ther mpled volatltes wth the Monte Carlo (MC and the lower bound approxmaton (LBA. Ths test s to compare the results wth dfferent moneyness (K/S from deep n the money.7 to deep out of the money 1.3. The other parameters are also beng set near the extreme case to see the numercal approxmaton performance, ncludng the jump ntensty λ = 4 (frequent jumps, correlaton coeffcent of Brownan motons ρ =.9 (hghly correlated, and local volatlty functon σ(t, S =.5S (hgh volatlty. Jump szes of four assets are also dfferent to reflect the heterogeneous settng wth h 1 =, h 2 =.1, h 2 =.3, h 3 =.5, whch mples that jump has no mpact to asset 1, small postve mpact to assets 2 and 3, and bg negatve mpact to asset 4. We compute the prces and ther mpled volatltes. The average relatve errors are about 4 percent, whch s heavly nfluenced by errors n computng out-of-money call opton prces wth long maturty (T = 2. Ths s perhaps not surprsng as the asymptotc expanson method works well only when all coeffcents are very small whleas the parameters n ths test are chosen to be large. Table 3 dsplays the numercal results wth the MC, AEA and LBA and wth dfferent maturtes (T = 1,3 and local volatlty functons (σ(s = αs β wth α =.2,.5 and β = 1,.8,.5. The jump sze s h = e.25 1 =.2212 for all assets. The last row dsplays the average errors of the AEA and LBA. Table 3 shows that the performance of the LBA s excellent wth the average relatve error.4%. The overall performance of the LBA s better than that of the AEA, especally when the local volatlty functon s σ(s =.5S. Matlab s used for computaton. The LBA only takes a few seconds for each case, much faster than the AEA and MC. We have done other tests wth the same data as n Table 3 except dfferent ntensty rates (λ = and 1. The overall performance s essentally the same as that of Table 3. We can say wth reasonable confdence that the LBA suggested n ths paper works well. Tables 4 dsplays the numercal results wth the MC and LBA and wth dfferent maturtes T and local volatlty functons σ(t,s = αs β. The jump szes for four assets are h 1 =, h 2 =.3, h 3 =.3, and h 4 =, whch mples that the jump event has no mpact to asset prces 1 and 4, 9

10 postve mpact to asset prce 2 and negatve mpact to asset prce 3. Snce The PEA requres the Black-Scholes settng whle the AEA requres the same jump sze for all assets, nether method can be used n ths test. It s clear that the LBA performs well n comparson wth the MC wth the overall relatve error less than 1%. 5 Conclusons In ths paper we suggest to approxmate the basket call opton value of a local volatlty jumpdffuson model wth Rogers and Sh s lower bound whch s approxmated wth the second order asymptotc expanson method. The lower bound approxmaton (LBA s easly computed and, If the local volatlty functon s tme ndependent (e.g., the CEV model, has a closed-form expresson. We compare the numercal performance of the LBA wth other methods usng dfferent parameters and show that the LBA s fast and accurate n most cases. Acknowledgment. The authors are grateful to the anonymous referees for ther constructve comments and suggestons whch have helped to mprove the paper. References [1 Albrecher, H., Mayer, P., and Schoutens W. (28. General lower bounds for arthmetc Asan opton prces. Appled Mathematcal Fnance, 15, [2 Benhamou, E., Gobet, E. and Mr. M. (29. Smart expanson and fast calbraton for jump dffusons, Fnance and stochastcs, 4, [3 Curran,M. (1994. Valung Asan and portfolo optons by condtonng on the geometrc mean prce, Management Scence, 4, [4 Deelstra, G., Lnev, J., and Vanmaele, M. (24. Prcng of arthmetc basket optons by condtonng, Insurance: Mathematcs and Economcs, 34, [5 Deelstra, G., Dhaene, J., and Vanmaele, M. (21. An overvew of comonotoncty and ts applcatons n fnance and nsurance, Advanced Mathematcal Methods for Fnance, Nunno, G., and Oksendal, B. (Eds, Sprnger, Germany. [6 Duffe, D., Pan, J. and Sngleton, K. (2. Transform analyss and asset prcng for affne jump-dffusons, Econometrca, 68, [7 Kuntomo, N. and Takahash, A.(21. The asymptotc expanson approach to the valuaton of nterest rate contngent clams, Mathematcal Fnance, 11, [8 Levy, E. (1992. Prcng European average rate currency optons. Journal of Internatonal Money and Fnance, 11, [9 Rogers, L.C.G., and Sh, Z. (1995. The value of an Asan opton, Journal of Appled Probablty, 32, [1 Xu, G. and Zheng, H. (29. Approxmate basket optons valuaton for a jump-dffuson model, Insurance: Mathematcs and Economcs, 45, [11 Xu, G. and Zheng, H. (21. Basket optons valuaton for a local volatlty jump-dffuson model wth the asymptotc expanson method, Insurance: Mathematcs and Economcs, 47,

11 λ η T MC (stdev PEA (IV% AEA (IV% LBA (IV% ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( ( (15.2 Average Rel Err % Table 1: The comparson of European basket call opton prces wth the Monte Carlo (MC, the partal exact approxmaton (PEA, the asymptotc expanson approxmaton (AEA and the lower bound approxmaton (LBA. The table dsplays the results wth dfferent jump ntenstes λ, jump szes h = e η 1, maturtes T, and local volatlty functon σ (t,s =.2S. The numbers nsde brackets n the MC columns are the standard devatons and those n the PEA, AEA and LBA columns are the mpled volatltes (IV computed wth the Black-Scholes call opton formula. T Prce Imp Vol % Moneyness % MC (stdev LBA MC LBA ( ( ( ( ( ( ( ( ( ( Average Rel Err % Table 2: The comparson of European basket call opton prces wth the Monte Carlo (MC and the lower bound approxmaton (LBA. The table dsplays the prces and ther mpled volatltes wth dfferent maturtes T and moneyness K/S. The local volatlty functon s σ (t,s =.5S, thecorrelaton coeffcent ofbrownanmotons s ρ =.9, thejumpntensty λ = 4, andjumpszes h 1 =,h 2 =.1,h 3 =.3,h 4 =.5. The numbers nsde brackets n the MC columns are the standard devatons. The last row s the average relatve errors between results of MC and LBA. 11

12 T α β MC (stdev AEA LBA ( ( , ( ( ( ( ( ( ( ( ( ( Average Rel Err %.7.4 Table 3: The comparson of European basket call opton prces wth the MC, AEA and LBA. The table dsplays results wth dfferent maturtes T, local volatlty functons σ (t,s = αs β and jump szes h = e.25 1 = T α β MC (stdev LBA ( ( ( ( ( ( ( ( ( ( ( ( Average Rel Err %.3 Table 4: The comparson of European basket call opton prces wth the MC and LBA. The table dsplays results wth dfferent maturtes T, local volatlty functons σ (t,s = αs β, and jump szes h 1 =, h 2 =.3, h 3 =.3, and h 4 =. 12