Friday, August 8, 2014

Jim Gatheral, Merrill Lynch Lecture 1: Stochastic Volatility and Local Volatility

Lecture 1: Stochastic Volatility and Local Volatility
Jim Gatheral, Merrill Lynch∗
Case Studies in Financial Modelling Course Notes, Courant Institute of Mathematical Sciences, Fall Term, 2002
Abstract In the course of the following lectures, we will study why equity options are priced as they are. In so doing, we will apply many of the techniques students will have learned in previous semesters and developsomeintuitionforthepricingofbothvanillaandexoticequity options. By considering specific examples, we will see that in pricing options, it is often as important to take into account the dynamics of underlying variables as it is to match known market prices of other claims. My hope is that these lectures will prove particularly useful to those who end up specializing in the structuring, pricing, trading and risk management of equity derivatives.
∗I am indebted to Peter Friz for carefully reading these notes, providing corrections and suggesting useful improvements.
1 Stochastic Volatility
1.1 Motivation That it might make sense to model volatility as a random variable should be clear to the most casual observer of equity markets. To be convinced, one only needs to remember the stock market crash of October 1987. Neverthe- less, given the success of the Black-Scholes model in parsimoniously describ- ing market options prices, it’s not immediately obvious what the benefits of making such a modelling choice might be. Stochastic volatility models are useful because they explain in a self- consistent way why it is that options with different strikes and expirations have different Black-Scholes implied volatilities (“implied volatilities” from now on) – the “volatility smile”. In particular, traders who use the Black- Scholes model to hedge must continuously change the volatility assumption in order to match market prices. Their hedge ratios change accordingly in an uncontrolled way. More interestingly for us, the prices of exotic options given by models based on Black-Scholes assumptions can be wildly wrong and dealers in such options are motivated to find models which can take the volatility smile into account when pricing these. From Figure 1, we see that large moves follow large moves and small moves follow small moves (so called “volatility clustering”). From Figures 2 and 3 (which shows details of the tails of the distribution), we see that the distribution of stock price returns is highly peaked and fat-tailed relative to the Normal distribution. Fat tails and the high central peak are character- istics of mixtures of distributions with different variances. This motivates us to model variance as a random variable. The volatility clustering feature implies that volatility (or variance) is auto-correlated. In the model, this is a consequence of the mean reversion of volatility 1. There is a simple economic argument which justifies the mean reversion of volatility (the same argument that is used to justify the mean reversion of interest rates). Consider the distribution of the volatility of IBM in one hundred years time say. If volatility were not mean-reverting ( i.e. if the distribution of volatility were not stable), the probability of the volatility of IBM being between 1% and 100% would be rather low. Since we believe that it is overwhelmingly likely that the volatility of IBM would in fact lie
1Note that simple jump-diffusion models do not have this property. After a jump, the stock price volatility does not change.
2
Figure 1: SPX daily log returns from 1/1/1990 to 31/12/1999
-0.06
-0.04
-0.02
0
0.02
0.04
0.06
0.08
Figure 2: Frequency distribution of SPX daily log returns from 1/1/1990 to 31/12/1999 compared with the Normal distribution
-0.06 -0.04 -0.02 0.02 0.04 0.06
10
20
30
40
50
60
in that range, we deduce that volatility must be mean-reverting. Having motivated the description of variance as a mean-reverting random variable, we are now ready to derive the valuation equation.
3
Figure 3: Tails of SPX frequency distribution
-0.06 -0.04 -0.02 0.02 0.04 0.06
0.25
0.5
0.75
1
1.25
1.5
1.75
2
1.2 Derivation of the Valuation Equation In this section, we follow Wilmott (1998) closely. We suppose that the stock price S and its variance v satisfy the following SDEs: dS(t) = µ(t)S(t)dt +qv(t)S(t)dZ1 (1) dv(t) = α(S,v,t)dt + ηβ(S,v,t)qv(t)dZ2 (2)with hdZ1 dZ2i= ρdt where µ(t) is the (deterministic) instantaneous drift of stock price returns, η is the volatility of volatility and ρ is the correlation between random stock price returns and changes in v(t). dZ1 and dZ2 are Wiener processes. The stochastic process (1) followed by the stock price is equivalent to the one assumed in the derivation of Black and Scholes (1973). This ensures that the standard time-dependent volatility version of the Black-Scholes formula (as derived in section 8.6 of Wilmott (1998) for example) may be retrieved in the limit η → 0. In practical applications, this is a key requirement of a stochastic volatility option pricing model as practitioners’ intuition for the behavior of option prices is invariably expressed within the framework of the Black-Scholes formula.
4
In the Black-Scholes case, there is only one source of randomness – the stock price, which can be hedged with stock. In the present case, random changes in volatility also need to be hedged in order to form a riskless port- folio. So we set up a portfolio Π containing the option being priced whose value we denote by V(S,v,t), a quantity −∆ of the stock and a quantity −∆1 of another asset whose value V1 depends on volatility. We have Π = V −∆S−∆1 V1 The change in this portfolio in a time dt is given by dΠ = (∂V ∂t + 1 2 vS2(t)∂2V ∂S2 + ρηvβ S(t) ∂2V ∂v∂S + 1 2 η2vβ2∂2V ∂v2)dt −∆1(∂V1 ∂t + 1 2 vS2(t)∂2V1 ∂S2 + ρηvβ S(t) ∂2V1 ∂v∂S + 1 2 η2vβ2∂2V1 ∂v2)dt +(∂V ∂S −∆1 ∂V1 ∂S −∆)dS +(∂V ∂v −∆1 ∂V1 ∂v)dv To make the portfolio instantaneously risk-free, we must choose ∂V ∂S −∆1 ∂V1 ∂S −∆ = 0 to eliminate dS terms, and ∂V ∂v −∆1 ∂V1 ∂v = 0 to eliminate dv terms. This leaves us with dΠ = (∂V ∂t + 1 2 vS2∂2V ∂S2 + ρηvβ S ∂2V ∂v∂S + 1 2 η2 vβ2∂2V ∂v2)dt −∆1(∂V1 ∂t + 1 2 vS2∂2V1 ∂S2 + ρηvβ S ∂2V1 ∂v∂S + 1 2 η2 vβ2∂2V1 ∂v2)dt = rΠdt = r(V −∆S−∆1V1)dt where we have used the fact that the return on a risk-free portfolio must equal the risk-free rate r which we will assume to be deterministic for our
5
purposes. Collecting all V terms on the left-hand side and all V1 terms on the right-hand side, we get
∂V ∂t + 1 2vS2∂2V ∂S2 + ρηvβ S ∂2V ∂v∂S + 1 2η2vβ2∂2V ∂v2 + rS∂V ∂S −rV∂V ∂v
=
∂V1 ∂t + 1 2vS2∂2V1 ∂S2 + ρηvβ S ∂2V1 ∂v∂S + 1 2η2vβ2∂2V1 ∂v2 + rS∂V1 ∂S −rV1∂V 1 ∂v The left-handside is a function of V only andthe right-handside is a function of V1 only. The only way that this can be is for both sides to be equal to some function f of the independent variables S, v and t. We deduce that
∂V ∂t
+
1 2
vS2∂2V ∂S2
+ρηvβ S
∂2V ∂v∂S
+
1 2
η2vβ2∂2V ∂v2
+rS
∂V ∂S −rV = − (α−ϕβ)
∂V ∂v (3) where, without loss of generality, we have written the arbitrary function f of S, v and t as (α−ϕβ). Conventionally, ϕ(S,v,t) is called the market price of volatility risk because it tells us how much of the expected return of V is explained by the risk (i.e. standard deviation) of v in the Capital Asset Pricing Model framework.
2 Local Volatility
2.1 History Given the computational complexity of stochastic volatility models and the extreme difficulty of fitting parameters to the current prices of vanilla op- tions, practitioners sought a simpler way of pricing exotic options consis- tently with the volatility skew. Since before Breeden and Litzenberger (1978), it was understood that the risk-neutral pdf could be derived from the market prices of European options. The breakthrough came when Dupire (1994) and Derman and Kani (1994) noted that under risk-neutrality, there was a unique diffusion process consistent with these distributions. The cor- responding unique state-dependent diffusion coefficient σL(S,t) consistent with current European option prices is known as the local volatility func- tion. It is unlikely that Dupire, Derman and Kani ever thought of local volatil- ity as representing a model of how volatilities actually evolve. Rather, it is
6
likely that they thought of local volatilities as representing some kind of average over all possible instantaneous volatilities in a stochastic volatility world (an “effective theory”). Local volatility models do not therefore really represent a separate class of models; the idea is more to make a simplify- ing assumption that allows practitioners to price exotic options consistently with the known prices of vanilla options. As if any proof had been needed, Dumas, Fleming, and Whaley (1998) performed an empirical analysis which confirmed that the dynamics of the implied volatilitysurface were not consistentwith the assumption of constant local volatilities. In section 2.5, we will show that local volatility is indeed an average over instantaneous volatilities, formalizing the intuition of those practition- ers who first introduced the concept.
2.2 A Brief Review of Dupire’s Work For a given expiration T and current stock price S0 , the collection {C (S0,K,T);K ∈(0,∞)}of undiscounted option prices of different strikes yields the risk neutral density function ϕ of the final spot ST through the relationship C (S0,K,T) =Z∞ K dST ϕ(ST,T;S0) (ST −K) Differentiate this twice with respect to K to obtain
ϕ(K,T;S0) =
∂2C ∂K2 so the Arrow-Debreu prices for each expiration may be recovered by twice differentiating the undiscounted option price with respect to K. This process will be familiar to any option trader as the construction of an (infinite size) infinitesimally tight butterfly around the strike whose maximum payoff is one. Given the distribution of final spot prices ST for each time T conditional on some starting spot price S0, Dupire shows that there is a unique risk neu- tral diffusion process which generates these distributions. That is, given the set of all European option prices, we may determine the functional form of the diffusion parameter (local volatility) of the unique risk neutral diffusion
7
process which generates these prices. Noting that the local volatility will in general be a function of the current stock price S0, we write this process as
dS S
= µ(t) dt + σ(S,t;S0) dZ
Application of Itˆo’s Lemma together with risk neutrality, gives rise to a par- tial differential equation for functions of the stock price which is a straightfor- ward generalization of Black-Scholes. In particular, the pseudo probability densities ϕ(K,T;S0) = ∂2C ∂K2 must satisfy the Fokker-Planck equation. This leads to the following equation for the undiscounted option price C in terms of the strike price K:
∂C ∂T
=
σ2K2 2
∂2C ∂K2
+ (r(T)−D(T))ÃC −K ∂C ∂K! (4)
where r(t) is the risk-free rate, D(t) is the dividend yield and C is short for C (S0,K,T). See the Appendix for a derivation of this equation. Were we to express the option price as a function of the forward price FT = S0 expnRT 0 µ(t)dto2, we would get the same expression minus the driftterm. That is ∂C ∂T = 1 2 σ2 K2 ∂2C ∂K2 where C now represents C (FT,K,T). Inverting this gives
σ2(K,T,S0) =
∂C ∂T
1 2 K2 ∂2C ∂K2
(5)
The right hand side of equation (5) can be computed from known European option prices. So, given a complete set of European option prices for all strikes and expirations, local volatilities are given uniquely by equation (5). We can view equation (5) as a definition of the local volatility function regardless of what kind of process (stochastic volatility for example) actually governs the evolution of volatility.
2From now on, µ(T) represents the risk-neutral drift of the stock price process which is the risk-free rate r(T) minus the dividend yield D(T)
8
2.3 Transforming to Black-Scholes Implied Volatility Space Market prices of options are quoted in terms of Black-Scholes implied volatil- ity σBS (K,T;S0). In other words, we may write
C (S0,K,T) = CBS (S0,K,σBS (S0,K,T),T)
It will be more convenient for us to work in terms of two dimensionless variables: the Black-Scholes implied total variance w defined by w(S0,K,T)≡σ2BS (S0,K,T) T and the log-strike y defined by y = lnµK FT¶ where FT = S0 expnRT 0 dtµ(t)o gives the forward price of the stock at time 0. In terms of these variables, the Black-Scholes formula for the future value of the option price becomes CBS (FT,y,w) = FT {N (d1)−eyN (d2)} = FT(NÃ− y √w + √w 2 !−eyNÃ− y √w − √w 2 !)(6)
and the Dupire equation (4) becomes
∂C ∂T
=
vL 2 (∂2C ∂y2 −
∂C ∂y)+ µ(T) C (7)
with vL = σ2 (S0,K,T) representing the local variance. Now, by taking derivatives of the Black-Scholes formula, we obtain ∂2CBS ∂w2 = Ã−1 8 − 1 2w + y2 2w2! ∂CBS ∂w ∂2CBS ∂y∂w = µ1 2 − y w¶ ∂CBS ∂w ∂2CBS ∂y2 − ∂CBS ∂y = 2 ∂CBS ∂w (8)
9
We may transform equation (7) into an equation in terms of implied variance by making the substitutions ∂C ∂y = ∂CBS ∂y + ∂CBS ∂w ∂w ∂y ∂2C ∂y2 = ∂2CBS ∂y2 + 2 ∂2CBS ∂y∂w ∂w ∂y + ∂2CBS ∂w2 Ã∂w ∂y!2 + ∂CBS ∂w ∂2w ∂y2 ∂C ∂T = ∂CBS ∂T + ∂CBS ∂w ∂w ∂T = ∂CBS ∂w ∂w ∂T + µ(T) CBS where the last equality follows from the fact that the only explicit depen- dence of the option price on T in equation (6) is through the forward price FT = S0 expnRT 0 dtµ(t)o . Equation (4) now becomes (cancelling µ(T) Cterms on each side) ∂CBS ∂w ∂w ∂T
=
vL 2   −
∂CBS ∂y
+
∂2CBS ∂y2 −
∂CBS ∂w
∂w ∂y
+ 2
∂2CBS ∂y∂w
∂w ∂y
+
∂2CBS ∂w2 Ã∂w ∂y!2 + ∂CBS ∂w
∂2w ∂y2  
=
vL 2
∂CBS ∂w   
2−
∂w ∂y
+ 2µ1 2 −
y w¶∂w ∂y
+Ã−1 8 −
1 2w
+
y2 2w2!Ã∂w ∂y!2 + ∂2w ∂y2  Then, taking out a factor of ∂CBS ∂w and simplifying, we get ∂w ∂T = vL   1− y w ∂w ∂y + 1 4Ã−1 4 − 1 w + y2 w2!Ã∂w ∂y!2 + 1 2 ∂2w ∂y2  Inverting this gives our final result:
vL =
∂w ∂T
1− y w
∂w ∂y + 1 4³−1 4 − 1 w + y2 w2´³∂w ∂y´2 + 1 2
∂2w ∂y2
2.4 Special Case: No Skew If the skew ∂w ∂y is zero3 , we must have
vL =
∂w ∂T
3Note that this implies that ∂ ∂KσBS (S0,K,T) is zero
10
SothelocalvarianceinthiscasereducestotheforwardBlack-Scholesimplied variance. The solution to this is of course w(T) =ZT 0 vL (t) dt 2.5 Local Variance as a Conditional Expectation of Instantaneous Variance In this section, we review the elegant derivation of Derman and Kani (1998). We assume the same stochastic process for the stock price as in equation (1) but write it in terms of the forward price Ft,T = St expnRT t dsµso. dFt,T =√vtFt,TdZ (9)
Note that dFT,T = dST. The undiscounted value of a European option with strike K expiring at time T is given by C (S0,K,T) = Eh(ST −K)+i Differentiating once with respect to K gives
∂C ∂K
=−E[θ(ST −K)] where θ(·) is the Heaviside function. Differentiating again with respect to K gives ∂2C ∂K2 = E[δ(ST −K)] where δ(·) is the Dirac δ function. Now, a formal application of Itˆo’s Lemma to the terminal payoff of the option (and using dFT,T = dST) gives the identity d(ST −K)+ = θ(ST −K) dST + 1 2 vT S2 T δ(ST −K) dT Taking conditional expectations of each side, and using the fact that Ft,T is a Martingale, we get dC = dEh(ST −K)+i= 1 2 EhvT S2 T δ(ST −K)idT 11
Also, we can write EhvTS2 T δ(ST −K)i = E[vT |ST = K] 1 2
K2 E[ δ(ST −K)]
= E[vT |ST = K]
1 2
K2 ∂2C ∂K2
Putting this together, we get
∂C ∂T
= E[vT |ST = K]
1 2
K2 ∂2C ∂K2 Comparing this with the definition of local volatility (equation (5)), we see that σ2(K,T,S0) = E[vT |ST = K] That is, local variance is the risk-neutral expectation of the instantaneous variance conditional on the final stock price ST being equal to the strike price K.
3 The Heston Model
3.1 The Model The Heston model (Heston (1993)) corresponds to choosing α(S,v(t),t) = −λ(v(t)− ¯ v) and β(S,v,t) = 1 in equations (1) and (2). These equations then become dS(t) = µ(t)S(t)dt +qv(t)S(t)dZ1 (10) and dv(t) =−λ(v(t)−¯ v)dt + ηqv(t)dZ2 (11) with hdZ1 dZ2i= ρdt where λ is the speed of reversion of v(t) to its long term mean ¯ v. The process followed by v(t) may be recognized as a version of the square root process described by Cox, Ingersoll, and Ross (1985). It is a (jump- free) special case of a so-called affine jump diffusion (AJD) which is roughly speaking a jump-diffusion process for which the drifts and covariances and jump intensities are linear in the state vector (which is{x,v}in this case with
12
x = log(S)). Duffie, Pan, and Singleton (2000) show that AJD processes are analytically tractable in general. The solution technique involves computing an “extended transform” which in the Heston case is a conventional Fourier transform. We now substitute the above values for α(S,v,t) and β(S,v,t) into the general valuation equation (equation (3)). We obtain ∂V ∂t + 1 2 vS2∂2V ∂S2 + ρηvS ∂2V ∂v∂S + 1 2 η2v∂2V ∂v2 + rS ∂V ∂S − rV = (λ(v−¯ v)−ϕ) ∂V ∂v (12) Now, to be able to use the AJD results, the market price of volatility risk also needs to be affine. Various economic arguments can be made (see for example Wiggins (1987)) that the market price of volatility risk ϕ should be proportional to the variance v. Then, let ϕ = θv for some constant θ. Now define the risk-adjusted parameters λ0 and ¯ v0 through λ0 = λ −θ, λ 0¯ v0 = λ¯ v. Substituting this into equation (12) gives ∂V ∂t + 1 2 vS2∂2V ∂S2 + ρηvS ∂2V ∂v∂S + 1 2 η2v∂2V ∂v2 + rS ∂V ∂S − rV = (λ0(v−¯ v0)) ∂V ∂v (13) Note that equation (13) is now identical to equation (12) with no explicit risk preference related parameters except that the parameters λ0 and ¯ v0 are now risk adjusted. From now on we will drop the primes on λ0 and ¯ v0 and assume that we are dealing with the risk-adjusted parameters.
3.2 The Heston Solution for European Options This section repeats the derivation of the Heston formula for the value of a European-style option first presented in Heston (1993) but with rather more detail than is provided in that paper. Before solving equation (13) with the appropriate boundary conditions, we can simplify it by making some suitable changes of variable. Let K be the strike price of the option, T be its expiry date and F(t,T) the forward price of the stock index to expiry. Then let x(t) = lnÃF(t,T) K !
13
Further, suppose that we consider only the future value to expiration C of the European option price rather than its value today and define τ = T−t. Then equation (13) simplifies to −
∂C ∂τ
+
1 2
v(t)C11 +
1 2
v(t)C1 +
1 2
η2v(t)C22 + ρηv(t)C12 −λ(v(t)−¯ v)C2 = 0 (14) where the subscripts 1 and 2 refer to differentiation with respect to x and v respectively. According to Duffie, Pan, and Singleton (2000), the solution of equation (14) has the form C(x,v,τ) =ex P1(x,v,τ)−P0(x,v,τ) (15) where the first term represents the pseudo-expectation of the final index level given that the option is in-the-money and the second term represents the pseudo-probability of exercise. Substituting the proposed solution (15) into equation (14) shows that P0 and P1 must satisfy the equation −
∂Pj ∂τ
+
1 2
v
∂Pj ∂x2 −³1 2 −j´v∂Pj ∂x
+
1 2
η2v∂2Pj ∂v2
+ρηv
∂2Pj ∂x∂v
+(a−bjv)
∂Pj ∂v
= 0 (16)
for j = 0,1 where
a=λ¯v, bj=λ−jρη subject to the terminal condition
lim τ→0
Pj(x,v,τ) = (1 if x > 0 0 if x≤0 ≡ θ(x) (17) We solve equation (16) subject to the condition (17) using a Fourier trans- form technique. To this end define the Fourier transform of Pj through ˜ P(k,v,τ) =Z∞ −∞ dxe−ikx P(x,v,τ) Then ˜ P(k,v,0) =Z∞ −∞ dxe−ikx θ(x) = 1 ik 14
The inverse transform is given by P(x,v,τ) =Z∞ −∞
dk 2π
eikx ˜ P(k,v,τ) (18)
Substituting this into equation (16) gives − ∂ ˜ Pj ∂τ − 1 2 k2 v ˜ Pj −(1 2 −j)ikv ˜ Pj
+
1 2
η2 v∂2 ˜ Pj ∂v2
+ ρηikv
∂ ˜ Pj ∂v
+ (a−bjv)
∂ ˜ Pj ∂v
= 0 (19)
Now define
α = −
k2 2 −
ik 2
+ ijk β = λ−ρηj−ρηik γ = η2 2
Then equation (19) becomes v(α ˜ Pj −β ∂ ˜ Pj ∂v
+ γ
∂2 ˜ Pj ∂v2)+a ∂ ˜ Pj ∂v −
∂ ˜ Pj ∂τ
= 0 (20)
Now substitute ˜ Pj(k,v,τ) = exp{C(k,τ)¯ v + D(k,τ)v} ˜ Pj(k,v,0) = 1 ik exp{C(k,τ)¯ v + D(k,τ)v} It follows that ∂ ˜ Pj ∂τ = (¯ v ∂C ∂τ + v ∂D ∂τ)˜ Pj ∂ ˜ Pj ∂v = D ˜ Pj ∂2 ˜ Pj ∂v2 = D2 ˜ Pj
15
Then equation (20) is satisfied if
∂C ∂τ
= λD
∂D ∂τ
= α−β D + γ D2 = γ(D−r+)(D−r−) (21)
where we define
r± =
β±√β2 −4αγ 2γ ≡
β±d η2 Integrating (21) with the terminal conditions C(k,0) = 0 and D(k,0) = 0 gives
D(k,τ) = r−
1−e−dτ 1−ge−dτ C(k,τ) = λ(r−τ − 2 η2
lnÃ1−ge−dτ 1−g !)
where we define
g ≡
r− r+ Taking the inverse transform using equation (18) and performing the com- plex integration carefully gives the final form of the pseudo-probabilities Pj in the form of an integral of a real-valued function.
Pj(x,v,τ) =
1 2
+
1 πZ∞ 0
dkRe(exp{Cj(k,τ)¯ v + Dj(k,τ)v + ikx} ik )
This integration may be performed using standard numerical methods. As a final point, it is worth noting that taking derivatives of the Heston formula with respect to x or v in order to derive risk parameters is extremely straightforward because the functions C(k,τ) and D(k,τ) are independent of x and v.
16
A Derivation of the Dupire Equation
Suppose the stock price diffuses with risk-neutral drift µ(t) and local volatil- ity σ(S,t) according to the equation: dS S = µ(t) dt + σ(S,t) dZ The undiscounted risk-neutral value C (S0,K,T) of a European option with strike K and expiration T is given by C (S0,K,T) =Z∞ K dST ϕ(ST,T;S0) (ST −K) (22) Here ϕ(ST,T;S0) is the pseudo probability density of the final spot at time T . It evolves according to the Fokker-Planck equation4 : 1 2 ∂2 ∂S2 T ³σ2S2 Tϕ´− S ∂ ∂ST (µSTϕ) = ∂ϕ ∂T Differentiating with respect to K gives ∂C ∂K = −Z∞ K dST ϕ(ST,T;S0) ∂2C ∂K2 = ϕ(K,T;S0)
Now, differentiating (22) with respect to time gives
∂C ∂T
= Z∞ K
dST ( ∂ ∂T
ϕ(ST,T;S0))(ST −K)
= Z∞ K
dST (1 2
∂2 ∂S2 T ³σ2S2 Tϕ´− ∂ ∂ST
(µSTϕ))(ST −K)
Integrating by parts twice gives: ∂C ∂T = σ2K2 2 ϕ +Z∞ K
dST µSTϕ
=
σ2K2 2
∂2C ∂K2
+ µ(T)ÃC −K ∂C ∂K!
4See Section 5 of Robert Kohn’s PDE for Finance Class Notes for a very readable account of this topic
17
which is the Dupire equation when the underlying stock has risk-neutral drift µ. That is, the forward price of the stock at time T is given by F (T) = S0 exp(ZT 0 dtµ(t))
References
Black, F., and M. Scholes, 1973, The pricing of options and corporate liabilities, Journal of Political Economy 81, 631–659. Breeden, D., and R. Litzenberger, 1978, Prices of state-contingent claims implicit in option prices, Journal of Business 51, 621–651. Cox, John C., Jonathan E. Ingersoll, and Steven A. Ross, 1985, A theory of the term structure of interest rates, Econometrica 53, 385–407. Derman, Emanuel, and Iraj Kani, 1994, Riding on a smile, Risk 7, 32–39. , 1998, Stochasticimpliedtrees: Arbitragepricingwithstochastictermand strike structure of volatility, International Journal of Theoretical and Applied Finance 1, 61–110. Duffie, Darrell, Jun Pan, and Kenneth Singleton, 2000, Transform analysis and asset pricing for affine jump diffusions, Econometrica 68, 1343–1376. Dumas, B., J. Fleming, and R. E. Whaley, 1998, Implied volatility functions: Empirical tests, The Journal of Finance 53. Dupire, Bruno, 1994, Pricing with a smile, Risk 7, 18–20. Heston, Steven L., 1993, A closed-form solution for options with stochastic volatil- ity, with application to bond and currency options, Review of Financial Studies 6, 327–343. Kohn, Robert V., 2000, PDE for Finance class notes, . Wiggins, J.B., 1987, Option values under stochastic volatility, Journal of Finan- cial Economics 15, 351–372. Wilmott, Paul, 1998, Derivatives. The Theory and Practice of Financial Engi- neering . chap. 23, pp. 299–304 (John Wiley & Sons: Chichester).

No comments:

Post a Comment