Consistently Modeling Joint Dynamics of Volatility and Underlying To Enable Effective Hedging
|
|
- Karin Harper
- 5 years ago
- Views:
Transcription
1 Consistently Modeling Joint Dynamics of Volatility and Underlying To Enable Effective Hedging Artur Sepp Bank of America Merrill Lynch, London Global Derivatives Trading & Risk Management 2013 Amsterdam April 16-18,
2 Plan of the presentation 1) Analyze the dependence between returns and volatility in conventional stochastic volatility (SV) models 2) Introduce the beta SV model by Karasinski-Sepp, Beta Stochastic Volatility Model, Risk, October ) Illustrate intuitive and robust calibration of the beta SV model to historical and implied data 4) Mix local and stochastic volatility in the beta SV model to produce different volatility regimes and equity delta 2
3 References Some theoretical and practical details for my presentation can be found in: Karasinski, P., Sepp, A. (2012) Beta Stochastic Volatility Model, Risk Magazine October, Sepp, A. (2013) Empirical Calibration and Minimum-Variance Delta Under Log-Normal Stochastic Volatility Dynamics, Working paper Sepp, A. (2013) Log-Normal Stochastic Volatility Model: Pricing of Vanilla Options and Econometric Estimation, Working paper 3
4 Empirical analysis I First, start with some empirical analysis to motivate the choice of beta SV model Data: 1) realized one month volatility of daily returns on the S&P500 index from January 1990 to March 2013 (276 observations) computed from daily returns within single month 2) the VIX index at the end of each month as a measure of implied one-month volatility of options on the S&P500 index 4
5 Empirical analysis II One month realized volatility is strongly correlated to implied volatility Left: time series of realized and implied volatilities Right: scatter plot of implied volatility versus realized volatility 80% 70% 60% Realized 1m vol VIX 80% 70% 60% 50% 50% 40% 30% VIX 40% 30% y = 0.73x R² = % 20% 10% 10% 0% 1-Feb-90 1-Dec-90 1-Oct-91 1-Aug-92 1-Jun-93 1-Apr-94 1-Feb-95 1-Dec-95 1-Oct-96 1-Aug-97 1-Jun-98 1-Apr-99 1-Feb-00 1-Dec-00 1-Oct-01 1-Aug-02 1-Jun-03 1-Apr-04 1-Feb-05 1-Dec-05 1-Oct-06 1-Aug-07 1-Jun-08 1-Apr-09 1-Feb-10 1-Dec-10 1-Oct-11 1-Aug-12 0% 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% Realized 1m volatility Observation: the level of realized volatility explains 78% of the level of implied volatility 5
6 Empirical analysis III Price returns are negatively correlated with changes in volatility Left: scatter plot of monthly changes in the VIX versus monthly returns on the S&P 500 index Right: the same for the realized volatility hly change in VIX Month y = -0.67x R² = % -25% -20% -15% -10% -5% 0% 5% 10% 15% 25% 20% 15% 10% 5% -5% -10% -15% -20% SP500 monthly return Monthly change in 1m realized vol y = -0.50x R² = % -25% -20% -15% -10% -5% 0% 5% 10% 15% 40% 30% 20% 10% -10% -20% -30% SP500 monthly return Observation: Changes in the S&P 500 index explain over 50% of variability in the implied volatility (the explanatory power is stronger for daily and weekly changes, about 80%) The impact is muted for realized volatility even though the beta coefficient is about the same in magnitude 6
7 Empirical analysis IV Conclusions for a robust SV model: 1) SV model should be consistent with the dynamics of the realized volatilities (for short-term vols) 2) SV should describe a robust dependence between the spot and implied volatility (incorporating of a local vol component, jumps) Impact on the delta-hedging: 1) leads to estimation of correct gamma P&L (expected implied vol should be above expected realized vol) 2) leads to estimation of correct vega and vanna, DvegaDspot 7
8 Conventional SV models I. The dynamics Start with analysis of dependence between returns and volatility in a typical SV model: ds(t) S(t) = V (t)dw (0) (t), S(0) = S dv (t) = a(v )dt + b(v )dw (v) (t), V (0) = V (1) with E[dW (0) (t)dw (v) (t)] = ρdt Here: V (t) - the instantaneous volatility of price returns b(v ) - volatility-of volatility that measures the overall uncertainty about the dynamic of volatility ρ - the correlation coefficient that measures linear dependence between returns in spot and changes in volatility 8
9 Conventional SV models II. Correlation Correlation is a linear measure of degree of strength between two variables Correlation can be sufficient for description of static data, but it can be useless for dynamic data Question: Given that the correlation between between volatility and returns is 0.80 and the realized spot return is 1% what is the expected change in volatility? Very practical question for risk computation The concept of correlation alone is of little help for an SV model 9
10 Volatility beta I Lets make simple transformation for SV model (1) using decomposition of Brownian W (v) (t) for volatility process: dw (v) (t) = ρdw (0) (t) + 1 ρ 2 dw (1) (t) with E[dW (0) (t)dw (1) (t)] = 0 Thus: Now where ds(t) S(t) = V (t)dw (0) (t) dv (t) = a(v )dt + b(v )ρdw (0) (t) + b(v ) 1 ρ 2 dw (1) (t) dv (t) = β(v ) ds(t) S(t) + ɛ(v )dw (1) (t) + a(v )dt β(v ) = b(v )ρ V, ɛ(v ) = b(v ) 1 ρ 2 (2) 10
11 Volatility beta II The volatility beta β(v ) is interpreted as a rate of change in the volatility given change in the spot with functional dependence on V : b(v )ρ β(v ) = V For log-normal volatility (SABR): b(v ) = ɛv so β(v ) is constant: β(v ) = ρɛ For normal volatility (Heston): b(v ) = ɛ so β(v ) is inversely proportional to V: β(v ) = ρɛ V For quadratic volatility (3/2 SV model): b(v ) = ɛv 2 so β(v ) proportional to V: β(v ) = (ρɛ)v We can measure the elasticity of volatility, α, with b(v ) = ɛv α, through its impact on the volatility beta 11
12 Volatility beta III. Empirically, the dependence between the volatility beta and the level of volatility is weak Left: scatter plot of VIX beta (regression of daily changes in the VIX versus daily returns on the S&P500 index within given month) versus average VIX within given month Right: the same scatter plot for logarithms of these variables and corresponding regression model (Heston model would imply the slope of 1, log-normal - of 0, 3/2 model - of 1) % 10% 20% 30% 40% 50% 60% 70% % -200% -150% -100% -50% 0% Monthly VIX beta onthly VIX beta) ) log ( abs (Mo y = 0.31x R² = Average monthly VIX log (Average monthly VIX) This rather implies a log-normal model for the volatility process, or that the elasticity of volatility α is slightly above one 12
13 Beta SV model with the elasticity of volatility: Let us consider the beta SV model with the elasticity of volatility: ds(t) S(t) = V (t)dw (0) (t) dv (t) = κ(θ V (t))dt + β[v (t)] α 1dS(t) S(t) + ɛ[v (t)]α dw (1) (t) Here β - the constant rate of change in volatility given change in the spot ɛ - the idiosyncratic volatility of volatility α - the elasticity of volatility κ - the mean-reversion speed θ - the long-term level of the volatility (3) 13
14 Beta SV model. Maximum likelihood estimation I Let x n denote log-return, x n = ln(s(t n )/S(t n 1 )) and v n volatility, v n = V (t n ) Apply discretization conditional on x n and y n 1 = V (t n 1 ): v n v n 1 = κ(θ v n 1 )dt n + β[v n 1 ] α 1 x n + ɛ[v n 1 ] α dt n ζ n (4) where ζ n are iid normals Apply the maximum likelihood estimator for the above model with the following specifications: α = 0 - the normal SV beta model α = 1 - the log-normal SV beta model α unrestricted - the SV beta model with the elasticity of volatility Estimation sample: Realized vols: compute monthly realized volatility over non-overlapping time intervals from April 1990 to March 2013 (sample size N = 277) Implied vols: one-month implied volatilities over weekly non-overlapping time intervals from October 2007 to March 2013 (N = 272) Four indices: S&P 500, FTSE 100, NIKKEI 225, STOXX 50 14
15 S&P500 (RV - estimation using realized vol; IV - estimation using implied vol) RV RV RV IV IV IV α β ɛ κ θ M L(Ω) Observations: 1) α is close to 1 for both realized and implied vols 2) β is larger for implied volatility 3) ɛ is larger for realized vols 4) κ is about the same for both realized and implied vols 5) θ is higher for implied vols 6) likelihood is larger for the log-normal SV than for normal SV 15
16 FTSE 100 (RV - estimation using realized vol; IV - estimation using implied vol) RV RV RV IV IV IV α β ɛ κ θ M L(Ω) Observations: 1) α is close to 1 for both realized and implied vols 2) β is larger for implied volatility 3) ɛ is larger for realized vols 4) κ is about the same for both realized and implied vols 5) θ is higher for implied vols 6) likelihood is larger for the log-normal SV than for normal SV 16
17 NIKKEI 225 (RV - estimation using realized vol; IV - estimation using implied vol) RV RV RV IV IV IV α β ɛ κ θ M L(Ω) Observations: 1) α is less than 1 for both realized and implied vols 2) β is larger for implied volatility 3) ɛ is about the same for both vols 4) κ is higher for implied vols 5) θ is about the same for both vols 6) likelihood is larger for the log-normal SV than for normal SV 17
18 STOXX 50 (RV - estimation using realized vol; IV - estimation using implied vol) RV RV RV IV IV IV α β ɛ κ θ M L(Ω) Observations: 1) α is close to 1 for both realized and implied vols 2) β is larger for implied volatility 3) ɛ is larger for realized vols 4) κ is higher for implied vols 5) θ is higher for implied vols 6) likelihood is larger for the log-normal SV than for normal SV 18
19 Conclusions The volatility process is closer to being log-normal (α 1) so that the log-normal volatility is a robust assumption The volatility is negatively proportional to changes in spot (to lesser degree for realized volatility) The volatility beta has and idiosyncratic volatility have similar magnitude among 4 indices (apart from NIKKEI 225 which typically has more convexity in implied vol), summarized below β β ɛ ɛ RV IV RV IV S&P FTSE Nikkei STOXX The beta SV model shares universal features across different underlyings and is robust for both realized and implied volatilities 19
20 Beta SV model. Option Pricing Introduce log-normal beta SV model under pricing measure: ds(t) S(t) = µ(t)dt + V (t)dw (0) (t), S(0) = S dv (t) = κ(θ V (t))dt + βv (t)dw (0) (t) + εv (t)dw (1) (t), V (0) = V (5) where E[dW (0) (t)dw (1) (t)] = 0, µ(t) is the risk-neutral drift Pricing equation for value function U(t, T, X, V ) with X = ln S(t): U t V 2 [U XX U X ] + µ(t)u X ( ε 2 + β 2) V 2 U V V + κ(θ V )U V + βv 2 U XV r(t)u = 0 (6) The beta SV model is not affine in volatility variable Nevertheless, I develop an accurate affine-like approximation 20
21 Beta SV model. Approximation I Fix expiry time T and introduce mean-adjusted volatility process Y (t): Y (t) = V (t) θ, θ E[V (T )] = θ + (V (0) θ)e κt Consider affine approximation for the moment generating function of X with second order in Y : G(t, T, X, Y ; Φ) = exp { ΦX + A (0) (t; T ) + A (1) (t; T )Y + A (2) (t; T )Y 2} where Y = V (0) θ and A (n) (T ; T ) = 0, n = 0, 1, 2 Substituting into (6) and keeping only quadratic terms in Y, yields: A (0) t + (1/2)ϑθ 2 [ 2A (2) + (A (1) ) 2] + κϖa (1) ΦβA (1) θ 2 + (1/2)θ 2 q = 0 A (1) t + ϑθ [ 2A (2) + (A (1) ) 2] + 2ϑθ 2 A (1) A (2) κa (1) + 2κϖA (2) 2ΦβA (1) θ 2ΦβA (2) θ 2 + θq = 0 A (2) t + (1/2)ϑ [ 2A (2) + (A (1) ) 2] + 4ϑθA (1) A (2) + 2ϑθ 2 (A (2) ) 2 2κA (2) ΦβA (1) 4ΦβA (2) θ + (1/2)q = 0 where q = Φ 2 + Φ, ϑ = ε 2 + β 2, ϖ = θ θ 21
22 Beta SV model. Approximation II This is system of ODE-s is solved by means of Runge-Kutta fourth order method in a fast way Thus, for pricing vanilla options under the beta SV model, we can apply the standard methods based on Fourier inversion (like in Heston and Stein-Stein SV models) In particular, the value of the call option with strike K is computed by applying Lipton s formula (Lipton (2001)): C(t, T, S, Y ) = e T t r(t )dt (e T t µ(t )dt S K [ G(t, T, x, Y ; ik 1/2) π R 0 k 2 + 1/4 where x = ln(s/k) + T t µ(t )dt The approximation formula is very accurate with differences between it and a PDE solver are less than 0.20% in terms of implied volatility It is straightforward to incorporate time-dependent model parameters (but not space-dependent local volatility) 22 ] dk )
23 Beta SV model. Implied volatility asymptotic I Now study the short-term implied volatility under the beta SV model Consider beta SV model with no mean-reversion: ds(t) S(t) = V (t)dw (0) (t) dv (t) = βv (t)dw (0) (t) + ɛv (t)dw (1) (t), V (0) = σ 0 where dw (0) dw (1) = 0 Follow idea from Andreasen-Huge (2013) and obtain the approximation for the implied log-normal volatility in the beta SV model: σ IMP (S, K) = ln(s/k) f(y), y = ln(s/k) V (0) y 1 f(y) = 0 J 1 (u)du = β ln J(y) β 2 + ɛ 2 + (β 2 + ɛ 2 )y β 2 + ɛ 2 β 2 + ɛ 2 β (1 ( J(y) = + β 2 + ɛ 2) y 2 2βy ) (7) 23
24 Beta SV model. Implied volatility asymptotic II Illustration of implied log-normal volatility computed by means PDE solver (PDE), ODE approximation (ODE) and implied vol asymptotic (IV Asymptotic) for maturity of one-month, T = 1/12 Left: V 0 = θ = 0.12, β = 1.0, κ = 2.8, ɛ = 0.5; Right: V 0 = θ = 0.12, β = 1.0, κ = 2.8, ɛ = % 25% Implied vol 21% 19% 17% 15% 13% PDE ODE IV Asymptotic Implied vol 23% 21% 19% 17% 15% PDE ODE IV Asymptotic 11% 13% 9% 11% 7% 9% 5% 80% 85% 90% 95% 100% 105% 110% Strike % 7% 80% 85% 90% 95% 100% 105% 110% Strike % Conclusion: approximation for short-term implied vol is very good for strikes in range [90%, 110%] even in the presence of mean-revertion 24
25 Beta SV model. Implied Volatility Asymptotic III Expand approximation for implied volatility (7) around k = 0, where k = ln(k/s) is log-strike: σ BSM (k) = σ βk Define skew (typically s = 5%): Thus: ( β2 σ ε2 σ 0 Skew s = 1 2s (σ IMP (+s) σ IMP ( s)) β = 2Skew s Define convexity (typically c = 5%): Convexity c = 1 c 2 (σ IMP (+c) + σ IMP ( c) 2σ IMP (0)) ) k 2 Thus: ε = 3σ IMP (0)Convexity c + 2(Skew s ) 2 25
26 Beta SV model. Implied volatility asymptotic IV As a result, volatility beta, β, can be interpreted as twice the implied volatility skew Idiosyncratic volatility of volatility, ɛ, can be interpreted as proportional to the square root of the implied convexity Time series of these implied parameters are illustrated below Oct Jan Apr Jul Oct Jan Apr Jul Oct Jan Apr Jul Oct Jan Apr Jul Oct Jan Apr Jul Oct Jan implied volatility beta implied idiosyncratic volatility of volatility Oct Jan Apr Jul Oct Jan Apr Jul Oct Jan Apr Jul Oct Jan Apr Jul Oct Jan Apr Jul Oct Jan-13 Observation: volatility beta, β, and idiosyncratic volatility of volatility, ɛ, exhibit range-bounded behavior (ranges are narrower when using longer dated implied vols to infer β and ɛ) 26
27 Summary so far I introduced the beta SV model and provided the intuition behind its key parameter - the volatility beta, β I showed that the beta SV model can adequately describe the historic dynamics of implied and realized volatilities with stable model parameters In terms of quality in fitting the implied volatility surface, the beta SV model is similar to other SV models - the model can explain the term structure of ATM volatility and longer-term skews (above 6m-1y) but it cannot reproduce steep short-term skews For short-term skews, we can introduce local volatility, jumps, a combination of both The important consideration is the impact on option delta In the second part of presentation, I concentrate on different volatility regimes and how to model them using the beta SV model 27
28 Volatility regimes and sticky rules (Derman) I 1) Sticky-strike: ( ) K σ(k; S) = σ 0 + Skew 1, σ AT M (S) = σ 0 + Skew S 0 ( ) S 1 S 0 ATM vol increase as the spot declines - typical of range-bounded markets 2) Sticky-delta: ( ) K S σ(k; S) = σ 0 + Skew, σ AT M (S) σ(s; S) = σ 0 S 0 The level of the ATM volatility does not depend on spot price -typical of stable trending markets 3) Sticky local volatility: σ(k; S) = σ 0 + Skew ( K + S S 0 2 ), σ AT M (S) = σ 0 + 2Skew ( ) S 1 S 0 ATM vol increase as the spot declines twice as much as in the sticky strike case - typical of stressed markets 28
29 Sticky rules II 40% 35% 30% S(0) = 1.00 goes down to S(1)=0.95 sticky strike sticky delta sticky local vol sigma(s0,k) 25% Volatility Sticky delta OldATM Vol 20% Implied 15% New S(1) Old S(0) 10% Spot Price Given: Skew = 1.0 and σ AT M (0) = 25.00% Spot change: down by 5% from S(0) = 1.00 to S(1) = 0.95 Sticky-strike regime: the ATM volatility moves along the original skew increasing by 5% Skew = 5% Sticky-local regime: the ATM volatility increases by 5% 2Skew = 10% and the volatility skew moves upwards Sticky-delta regime: the ATM volatility remains unchanged with the volatility skew moving downwards 29
30 Impact on option delta The key implication of the volatility rules is the impact on option delta We can show the following rule for call options: Sticky Local Sticky Strike Sticky Delta As a result, for hedging call options, one should be over-hedged (as compared to the BSM delta) in a trending market and under-hedged in a stressed market Thus, the identification of market regimes plays an important role to compute option hedges While computation of hedges is relatively easy for vanilla options and can be implemented using the BSM model, for path-dependent exotic options, we need a dynamic model consistent with different volatility regimes 30
31 Stickiness ratio I Given price return from time t n 1 to t n, X(t n ) = (S(t n ) S(t n 1 ))/S(t n 1 ), We make prediction for change in the ATM volatility: σ AT M (t n ) = σ AT M (t n 1 ) + Skew R(t n ) X(t n ) where the stickiness ratio R(t n ) indicates the rate of change in the ATM volatility predicted by skew and price return Informal definition of the stickiness ratio: R(t n ) = σ AT M(t n ) σ AT M (t n 1 ) X(t n )Skew 5% (t n 1 ) To estimate stickiness ratio, R,empirically, we apply regression model: σ AT M (t n ) σ AT M (t n 1 ) = R Skew 5% (t n 1 ) X(t n ) + ɛ n where X(t n ) is realized return for day n; ɛ n are iid normal residuals We expect that the average value of R, R, as follows: R = 1 under the sticky-strike regime R = 0 under the sticky-delta regime R = 2 under the sticky-local regime 31
32 Stickiness ratio II Empirical test is based on using market data for S&P500 (SPX) options from 9-Oct-07 to 1-Jul-12 divided into three zones crisis recovery range-bound start date 9-Oct-07 5-Mar Feb-11 end date 5-Mar Feb Jul-12 number days start SPX end SPX return % 96.76% 3.06% start ATM 1m 14.65% 45.28% 12.81% end ATM 1m 45.28% 12.81% 15.90% vol change 30.63% % 3.09% start Skew 1m % % % end Skew 1m % % % skew change 14.40% -8.20% 14.00% 32
33 Stickiness ratio III S&P500 1m ATM vol RANGE-BOUND 80% 70% RECOVERY 60% S&P500 CRISIS 1m ATM vol 50% 40% 30% % % Oct-07 Dec-07 Mar-08 Jun-08 Sep-08 Dec-08 Mar-09 Jun-09 Aug-09 Nov-09 Feb-10 May-10 Aug-10 Nov-10 Feb-11 Apr-11 Jul-11 Oct-11 Jan-12 Apr-12 Jul-12 33
34 Stickiness ratio (crisis) for 1m and 1y ATM vols Daily change in 1m ATM vol Stickeness for 1m ATM vol, crisis period Oct 07 -Mar 09 y = x R² = % -12% -10% -8% -6% -4% -2% 0% 2% 4% 6% 8% 15% 10% 5% -5% Daily change in 1y ATM vol Stickeness for 1y ATM vol, crisis period Oct 07 -Mar 09 y = x R² = % -4% -3% -2% -1% 0% 1% 2% 3% 15% 10% 5% -5% -10% -10% -15% 1m Skew * price return -15% 1y Skew * price return 34
35 Stickiness ratio (recovery) for 1m and 1y ATM vols Daily change in 1m ATM vol Stickeness for 1m ATM vol, recovery period Mar 09-Feb 11 15% 10% 5% y = x R² = % -6% -5% -4% -3% -2% -1% 0% 1% 2% 3% 4% 5% -5% Daily change in 1y ATM vol Stickeness for 1y ATM vol, recovery period Mar 09-Feb 11 y = x R² = % -2% -2% -1% -1% 0% 1% 1% 2% 15% 10% 5% -5% -10% -10% -15% 1m Skew * price return -15% 1y Skew * price return 35
36 Stickiness ratio (range) for 1m and 1y ATM vols Daily change in 1m ATM vol Stickeness for 1m ATM vol, range-bnd period Feb11-Aug12 15% 10% 5% y = x R² = % -6% -4% -2% 0% 2% 4% 6% 8% -5% Daily change in 1y ATM vol Stickeness for 1m ATM vol, range-bnd period Feb11-Aug12 15% 10% 5% y = x R² = % -2% -2% -1% -1% 0% 1% 1% 2% 2% 3% -5% -10% -10% -15% 1m Skew * price return -15% 1y Skew * price return 36
37 Stickiness ratio V. Conclusions Summary of the regression model: crisis recovery range-bound Stickiness, 1m Stickiness, 1y R 2, 1m 77% 65% 68% R 2, 1y 82% 68% 72% 1) The concept of the stickiness is statistically significant explaining about 80% of the variation in ATM volatility during crisis period and about 70% of the variation during recovery and range-bound periods 2) Stickiness ratio is stronger during crisis period, R 1.6 (closer to sticky local vol) less strong during recovery period, R 1.5 weaker during range-bound period, R 1.35 (closer to sticky-strike) 3) The volatility regime is typically neither sticky-local nor stickystrike but rather a combination of both 37
38 Stickiness ratio VII. Dynamic models A Now we consider how to model the stickiness ratio within the dynamic SV models The primary driver is change in the spot price, S/S The key in this analysis is what happens to the level of model volatility given change in the spot price The model-consistent hedge: The level of volatility changes proportional to (approximately): Skew SV Model S/S The model-inconsistent hedge: The level of volatility remains unchanged Implication for the stickiness under pure SV models: R = 2 under the model-consistent hedge R = 0 under the model-inconsistent hedge 38
39 Stickiness ratio VII. Dynamic models B How to make R = 1.5 using an SV model Under the model-consistent hedge: impossible Under the model-inconsistent hedge: mix SV with local volatility Remedy: add jumps in returns and volatility Under any spot-homogeneous jump model, R = 0 The only way to have a model-consistent hedging that fits the desired stickiness ratio is to mix stochastic volatility with jumps: the higher is the stickiness ratio, the lower is the jump premium the lower is the stickiness ratio, the higher is the jump premium Jump premium is lower during crisis periods (after a big crash or excessive market panic, the probability of a second one is lower because of realized de-leveraging and de-risking of investment portfolios, central banks interventions) Jump premium is higher during recovery and range-bound periods (renewed fear of tail events, increased leverage and risk-taking given small levels of realized volatility and related hedging) 39
40 Stickiness ratio VII. Dynamic models C The above consideration explain that the stickiness ratio is stronger during crisis period, R 1.6 (closer to sticky local vol) weaker during range-bound and recovery periods, R 1.35 (closer to sticky-strike) To model this feature within an SV model, we need to specify a proportion of the skew attributed to jumps (see my 2012 presentation on Global Derivatives) During crisis periods, the weight of jumps is about 20% During range-bound and recovery periods, the weight of jumps is about 40% For the rest of my presentation, we assume a model-inconsistent hedge and apply the beta SV model to model different volatility regimes Goal: create a dynamic model where R can be model as an input parameter 40
41 Beta SV model with CEV vol. Incorporate CEV volatility in the beta SV model (short-term analysis with no mean-reversion): ds(t) S(t) = V (t)[s(t)]β SdW (0) (t), S(0) = S 0 dv (t) = β V ds(t) S(t) + ɛv (t)dw (1) (t) = β V V (t)[s(t)] β SdW (0) (t) + ɛv (t)dw (1) (t), V (0) = σ 0 To produce volatility skew and price-vol dependence: β S is the backbone beta, β S 0 β V is the volatility beta, β V 0 Connection to SABR model (Hagan et al (2002)): V (0) = ˆα, β S = β 1, β V [S(t)] β S = νρ, ɛ = ν 1 ρ 2 Note that volatility beta, β V, measures the sensitivity of instantaneous volatility to changes in the spot independent on assumption about the local volatility of the spot Term β V [S(t)] β S, if β S < 0, increases skew in put wing 41 (8)
42 Beta SV model with CEV vol. Implied volatility To obtain approximation for the implied volatility under the beta CEV SV model, we apply formula (7) for implied vol asymptotic with β = β V and y = 1 σ 0 S β S K β S β S Limit in k = 0, k = ln(k/s): [ σ BSM (k) = S β S σ (σ 0β S + β V ) k ( σ 0 β 2 S β2 V σ ε2 σ 0 ) k 2 ] (9) 42
43 Beta SV model with CEV vol. Implied volatility I Illustration of implied log-normal volatility computed by means PDE solver (PDE), implied vol asymptotic (IV Asymptotic) and the second order expansion (2-nd order) for maturity of one-month, T = 1/12 Left: V 0 = θ = 0.12, β V = 1.0, β S = 2.0, κ = 2.8, ɛ = 0.5; Right: V 0 = θ = 0.12, β V = 1.0, β S = 2.0, κ = 2.8, ɛ = % 35% Implied vol 25% 20% PDE IV Asymptotic 2-nd order Implied vol 30% 25% 20% PDE IV Asymptotic 2-nd order 15% 15% 10% 10% 5% 80% 85% 90% 95% 100% 105% 110% Strike % 5% 80% 85% 90% 95% 100% 105% 110% Strike % Conclusion: approximation for short-term implied vol is very good for strikes in range [90%, 110%] even in the presence of mean-revertion 43
44 Interpretation of model parameters Volatility beta β V is a measure of linear dependence between daily returns and changes in the ATM volatility: σ AT M (S(t n )) σ AT M (S(t n 1 )) = β V S(t n ) S(t n 1 ) S(t n 1 ) The backbone beta β S is a measure of daily changes in the logarithm of the ATM volatility to daily returns on the stock ln [σ AT M (S(t n ))] ln [ σ AT M (S(t n 1 )) ] = β S S(t n ) S(t n 1 ) S(t n 1 ) Next I examine these regression models empirically, assuming: 1) all skew is generated by β V with β S = 0 2) all skew is generated by β S with β V = 0 44
45 Volatility beta (crisis) for 1m and 1y ATM vols Daily change in 1m ATM vol Daily change in 1m ATM vol, crisis period Oct 07 - Mar 09 15% 10% 5% y = x R² = % -15% -10% -5% 0% 5% 10% 15% -5% Daily change in 1y ATM vol Daily change in 1y ATM vol, crisis period Oct 07 - Mar 09 15% 10% 5% y = x R² = % -15% -10% -5% 0% 5% 10% 15% -5% -10% -10% Daily price return -15% Daily price return -15% 45
46 Volatility beta (recovery) for 1m and 1y ATM vols Daily change in 1m ATM vol Daily change in 1m ATM vol, recovery period Mar 09-Feb 11 15% 10% 5% y = x R² = % -6% -4% -2% 0% 2% 4% 6% 8% -5% Daily change in 1y ATM vol Daily change in 1y ATM vol, recovery period Mar 09-Feb 11 15% 10% 5% y = x R² = % -6% -4% -2% 0% 2% 4% 6% 8% -5% -10% -10% Daily price return -15% Daily price return -15% 46
47 Volatility beta (range) for 1m and 1y ATM vols Daily change in 1m ATM vol Daily change in 1m ATM vol, range-bnd period Feb11-Aug12 y = x R² = % -8% -6% -4% -2% 0% 2% 4% 6% 15% 10% 5% -5% Daily change in 1y ATM vol Daily change in 1y ATM vol, range-bnd period Feb11-Aug12 y = x R² = % -8% -6% -4% -2% 0% 2% 4% 6% 15% 10% 5% -5% -10% -10% Daily price return -15% Daily price return -15% 47
48 Volatility beta. Summary crisis recovery range-bound Volatility beta 1m Volatility beta 1y R 2 1m 76% 60% 68% R 2 1y 80% 65% 71% The volatility beta is pretty stable across different market regimes The longer term ATM volatility is less sensitive to changes in the spot Changes in the spot price explain about: 80% in changes in the ATM volatility during crisis period 60% in changes in the ATM volatility during recovery period (ATM volatility reacts slower to increases in the spot price) 70% in changes in the ATM volatility during range-bound period (jump premium start to play bigger role n recovery and range-bound periods) 48
49 Backbone beta (crisis) for 1m and 1y ATM vols Daily change in log 1m ATM vol Daily change in log 1m ATM vol, crisis period Oct07 -Mar09 30% 20% 10% y = x R² = % -15% -10% -5% 0% 5% 10% 15% -10% Daily change in 1y log ATM vol Daily change in log 1y ATM vol, crisis period Oct 07 -Mar 09 30% 20% 10% y = x R² = % -15% -10% -5% 0% 5% 10% 15% -10% -20% -20% Daily price return -30% Daily price return -30% 49
50 Backbone beta (recovery) for 1m and 1y ATM vols Daily change in log 1m ATM vol Daily change in log 1m ATM vol, recov period Mar09-Feb11 30% 20% 10% y = x R² = % -6% -4% -2% 0% 2% 4% 6% 8% -10% Daily change in log 1y ATM vol Daily change in log 1y ATM vol, recov period Mar09-Feb11 30% 20% 10% y = x R² = % -6% -4% -2% 0% 2% 4% 6% 8% -10% -20% Daily price return -30% -20% Daily price return -30% 50
51 Backbone beta (range) for 1m and 1y ATM vols Daily change in log 1m ATM vol Daily change in log 1m ATM vol, range period Feb11-Aug12 y = x R² = % -8% -6% -4% -2% 0% 2% 4% 6% 30% 20% 10% -10% Daily change in log 1y ATM vol Daily change in log 1y ATM vol, range period Feb11-Aug12 y = x R² = % -8% -6% -4% -2% 0% 2% 4% 6% 30% 20% 10% -10% -20% -20% Daily price return -30% Daily price return -30% 51
52 The backbone beta. Summary A crisis recovery range-bound Backbone beta 1m Backbone beta 1y R 2 1m 67% 54% 62% R 2 1y 77% 61% 70% The value of the backbone beta appears to be less stable across different market regimes (compared to volatility beta) Explanatory power is somewhat less (by 5-7%) for 1m ATM vols (compared to volatility beta) Similar explanatory power for 1y volatilities 52
53 The backbone beta. Summary B Change in the level of the ATM volatility implied by backbone beta β S is proportional to initial value of the ATM volatility High negative value of β S implies a big spike in volatility given a modest drop in the price - a feature of sticky local volatility model In the figure, using estimated parameters β V range-bound period, σ(0) = 20% = 1.07, β S 4.54 in Chan nge in vol level 16% 14% 12% 10% 8% 6% 4% Predicted change in volatility volatility beta_v=-1.07 backbone beta_s= % Spot return % 0% -10% -9% -7% -6% -4% -2% -1% 53
54 Model implied skew Using approximation (9) for short-term implied volatility, obtain the following approximate but accurate relationship between the model parameters and short-term implied ATM volatility, σ AT M (S), and skew Skew s : σ 0 S β S = σ AT M (S) σ 0 β S + β V = 2Skew s The first equation is known as the backbone that defines the trajectory of the ATM volatility given a change in the spot price: σ AT M (S) σ AT M (S 0 ) σ AT M (S 0 ) β S S S 0 S 0 (10) 54
55 Model implied stickiness and volatility regimes If we insist on model-inconsistent delta (change in spot with volatility level unchanged): fit backbone beta β S to reproduce specified stickiness ratio adjust β V so that the model fits the market skew Using stickiness ratio R(t n ) along with (10), we obtain that empirically: β S (t n ) = Skew s(t n 1 ) σ AT M (t n 1 ) R(t n) Thus, given an estimated value of the stickiness rate we imply β S Finally, by mixing parameters β S and β V we can produce different volatility regimes: sticky-delta with β S = 0 and β V 2Skew s sticky-local volatility with β V = 0 and β S 2Skew s /σ 0 From the empirical data we infer that, approximately, β S 60% 2Skew s /σ 0 and β V = 40% 2Skew s 55
56 Illustration of model implied stickiness T = 1/12, V 0 = 0.12, ɛ = 0.5; w SV = {1.00, 0.75, 0.50, 0.25, 0.00} β V = { 1.00, 0.75, 0.50, 0.25, 0.00}, β S = {0.00, 2.08, 4.17, 6.25, 8.33} The initial skew and ATM vol is the same for all values of w SV Left: skew change given spot return of 5% Right: corresponding ATM vol (lhs) and stickeness (rhs) 60% 19% 3.00 Implied volatility 50% 40% 30% 20% 10% Initial Skew SV=1 SV=0.75 SV=0.5 SV=0.25 SV=0 TM implied vol AT 18% 17% 16% 15% 14% 13% ATM Vol (left) Stickiness (right) Stickines 0% 75% 80% 85% 90% 95% 100% 105% 110% 115% 120% 125% K/S(0) 12% 11% Initial Skew SV=1 SV=0.75 SV=0.5 SV=0.25 SV= Conclusion: the stickiness ratio is approximately equal to twice the weight of the SV implied skew 56
57 Full beta SV model. Dynamics Pricing version of the beta SV model is specified under the pricing measure in terms of a normalized volatility process Y (t): ds(t) S(t) = µ(t)dt + (1 + Y (t))σdw (0) (t), S(0) = S dy (t) = κy (t)dt + β V (1 + Y (t))σdw (0) (t) + ε(1 + Y (t))dw (1) (t), Y (0) = 0 (11) where E[dW (0) (t)dw (1) (t)] = 0 σ is the overall level of volatility: it can be set constant, deterministic, or local stochastic volatility, σ LSV (t, S) (parametric, like CEV, or nonparametric) β V is the rate of change in the normalized SV process Y (t) to changes in the spot For constant volatility σ and θ = σ, parameters of the normalized SV beta are related as follows: Y (t) = V (t) σ 1, β V = β V σ, κ = κ, ε = ε 57
58 Beta stochastic volatility model. Calibration Parameters of SV process, β V, ε and κ are specified before calibration We calibrate the local volatility σ σ LSV (t, S), using either a parametric local volatility (CEV) or non-parametric local volatility, so that the vanilla surface is matched by construction For calibration of σ LSV (t, S) we apply the conditional expectation (Lipton 2002): σ 2 LSV (T, K)E [ (1 + Y (T )) 2 S(T ) = K ] = σ 2 LV (T, K) where σ 2 LV (T, K) is Dupire local volatility The above expectation is computed by solving the forward PDE corresponding to dynamics (11) using finite-difference methods and computing σlsv 2 (T, K) stepping forward in time (for details, see my GB presentation in 2011) Once σ LSV (t, S) is calibrated we use either backward PDE-s or MC simulation for valuation of non-vanilla options 58
59 Full beta SV model. Monte-Carlo simulation Simulate process Z(t): Z(t) = ln (1 + Y (t)), Z(0) = 0 with dynamics: ( dz(t) = κ [ ( β V ) 2 + ( ε) 2] ) e Z(t) dt + β V σdw (0) (t) + εdw (1) (t) 2 The domain of definition: Y (t) ( 1, ), Z(t) (, ) Obtain the instanteneous volatility by inversion: Y (t) = e Z(t) 1 No problem with boundary at zero that exist in some SV models 59
60 Full beta SV model. Multi-asset Dynamics I. N-asset dynamics: ds i (t) S i (t) = µ i(t)dt + (1 + Y i (t))σ i dw (0) i (t) dy i (t) = κ i Y i (t)dt + β V,i (1 + Y i (t))σ i dw (0) i (t) + ε i (1 + Y i (t))dw (1) i (t) where E[dW (0) i (t)dw (0) j (t)] = ρ (0) ij dt, {i, j} = 1,..., N, where ρ(0) ij is ith asset - jth asset correlation E[dW (0) i (t)dw (1) i (t)] = 0, i = 1,..., N E[dW (0) i (t)dw (1) j (t)] = 0, {i, j} = 1,..., N E[dW (1) i (t)dw (1) j (t)] = ρ (1) dt, {i, j} = 1,..., N, idiosyncratic volatilities are correlated (we take ρ (1) = 1 to avoid de-correlation) MC simulation: 1) Simulate N Brownian increments dw (0) i (t) with correlation matrix {ρ (0) ij }; 2) Simulate N Brownian increments dw (1) (t) with correlation matrix {ρ (1) } (only one Brownian is needed if ρ (1) = 1) i 60
61 Full beta SV model. Multi-asset Dynamics II Implied instantaneous cross asset-volatility correlation: ρ ( ) dsi (t) S i (t), dy j(t) = β V,j σ j ρ (0) ij ( β V,j σ j ) 2 + ( ε j ) 2 ρ(0) ij so the correlation is bounded from below by ρ (0) ij ( εj β V,j σ j ) 2 and declines less for large β V,j and small ε j Instantaneous volatility-volatility correlation (taking ρ (1) = 1): ρ ( dy i (t), dy j (t) ) = ρ (0) ij ρ (0) ij ( εi β V,i σ i + ( εi β V,i σ i + β V,i β V,j σ i σ j ρ (0) ij + ε i ε j ρ (1) ( β V,i ) 2 σ 2 i + ( ε i) 2 ( β V,j ) 2 σ 2 j + ( ε j) 2 ε j β V,j σ j ε j β V,j σ j ) 2 ) 2 + ρ (1) ε i ε j β V,i σ i β V,j σ j so volatility-volatility correlation is increased if ρ (1) = 1 61
62 Summary 1) Presented the beta SV model and illustrated that the model can describe very well the dynamics of both implied and realized volatilities 2) Obtained an accurate short-term asymptotic for the implied volatility in the beta SV model and showed how to express the key model parameters, the volatility beta and idiosyncratic volatility, in terms of the market implied skew and convexity 3) Derived an accurate closed form solution for pricing vanilla option in the SV beta model using Fourier inversion method 4) Presented the beta SV model with CEV local volatility to model different volatility regimes and compute appropriate option delta 5) Extended the beta SV model for multi-asset dynamics 62
63 Final words I am thankful to the members of BAML Global Quantitative Analytics The opinions and views expressed in this presentation are those of the author alone and do not necessarily reflect the views and policies of Bank of America Merrill Lynch Thank you for your attention! 63
64 References Andreasen, J. and B Huge (2013). Risk, January, Expanded forward volatility, Hagan, P., Kumar, D., Lesniewski, A., and Woodward, D. (2002), Managing smile risk, Wilmott Magazine, September, Karasinski, P and A. Sepp (2012). Beta Stochastic Volatility Model, Risk, October, ( Lipton, A. (2001). Mathematical Methods for Foreign Exchange: a Financial Engineer s Approach, World Scientific, Singapore Lipton, A. (2002). The vol smile problem, Risk, February, Sepp, A. (2012), Achieving Consistent Modeling Of VIX and Equities Derivatives, Global Derivatives conference in Barcelona 64
65 Sepp, A. (2011), Efficient Numerical PDE Methods to Solve Calibration and Pricing Problems in Local Stochastic Volatility Models, Global Derivatives conference in Paris Sepp, A., (2013), When You Hedge Discretely: Optimization of Sharpe Ratio for Delta-Hedging Strategy under Discrete Hedging and Transaction Costs, The Journal of Investment Strategies 3(1), Sepp, A., (2014), Empirical Calibration and Minimum-Variance Delta Under Log-Normal Stochastic Volatility Dynamics, Working paper, Sepp, A., (2014), Log-Normal Stochastic Volatility Model: Pricing of Vanilla Options and Econometric Estimation, Working paper,
Beta stochastic volatility model. Artur Sepp Bank of America Merrill Lynch, London
Beta stochastic volatility model Artur Sepp Bank of America Merrill Lynch, London artur.sepp@baml.com Financial Engineering Workshop Cass Business School, London October 24, 2012 1 Introduction. Quoting
More informationAdvanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives
Advanced Topics in Derivative Pricing Models Topic 4 - Variance products and volatility derivatives 4.1 Volatility trading and replication of variance swaps 4.2 Volatility swaps 4.3 Pricing of discrete
More informationPricing Barrier Options under Local Volatility
Abstract Pricing Barrier Options under Local Volatility Artur Sepp Mail: artursepp@hotmail.com, Web: www.hot.ee/seppar 16 November 2002 We study pricing under the local volatility. Our research is mainly
More information1. What is Implied Volatility?
Numerical Methods FEQA MSc Lectures, Spring Term 2 Data Modelling Module Lecture 2 Implied Volatility Professor Carol Alexander Spring Term 2 1 1. What is Implied Volatility? Implied volatility is: the
More informationPrincipal Component Analysis of the Volatility Smiles and Skews. Motivation
Principal Component Analysis of the Volatility Smiles and Skews Professor Carol Alexander Chair of Risk Management ISMA Centre University of Reading www.ismacentre.rdg.ac.uk 1 Motivation Implied volatilities
More informationRealized and implied index skews, jumps, and the failure of the minimum-variance hedging. Artur Sepp
Realized and implied index skews, jumps, and the failure of the minimum-variance hedging Artur Sepp Global Risk Analytics Bank of America Merrill Lynch, London artur.sepp@baml.com Global Derivatives Trading
More informationCalibration of SABR Stochastic Volatility Model. Copyright Changwei Xiong November last update: October 17, 2017 TABLE OF CONTENTS
Calibration of SABR Stochastic Volatility Model Copyright Changwei Xiong 2011 November 2011 last update: October 17, 2017 TABLE OF CONTENTS 1. Introduction...2 2. Asymptotic Solution by Hagan et al....2
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets (Hull chapter: 12, 13, 14) Liuren Wu ( c ) The Black-Scholes Model colorhmoptions Markets 1 / 17 The Black-Scholes-Merton (BSM) model Black and Scholes
More informationRisk managing long-dated smile risk with SABR formula
Risk managing long-dated smile risk with SABR formula Claudio Moni QuaRC, RBS November 7, 2011 Abstract In this paper 1, we show that the sensitivities to the SABR parameters can be materially wrong when
More informationMLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models
MLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models Matthew Dixon and Tao Wu 1 Illinois Institute of Technology May 19th 2017 1 https://papers.ssrn.com/sol3/papers.cfm?abstract
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets Liuren Wu ( c ) The Black-Merton-Scholes Model colorhmoptions Markets 1 / 18 The Black-Merton-Scholes-Merton (BMS) model Black and Scholes (1973) and Merton
More informationDynamic Relative Valuation
Dynamic Relative Valuation Liuren Wu, Baruch College Joint work with Peter Carr from Morgan Stanley October 15, 2013 Liuren Wu (Baruch) Dynamic Relative Valuation 10/15/2013 1 / 20 The standard approach
More informationSimple Robust Hedging with Nearby Contracts
Simple Robust Hedging with Nearby Contracts Liuren Wu and Jingyi Zhu Baruch College and University of Utah October 22, 2 at Worcester Polytechnic Institute Wu & Zhu (Baruch & Utah) Robust Hedging with
More informationExploring Volatility Derivatives: New Advances in Modelling. Bruno Dupire Bloomberg L.P. NY
Exploring Volatility Derivatives: New Advances in Modelling Bruno Dupire Bloomberg L.P. NY bdupire@bloomberg.net Global Derivatives 2005, Paris May 25, 2005 1. Volatility Products Historical Volatility
More informationInterest Rate Volatility
Interest Rate Volatility III. Working with SABR Andrew Lesniewski Baruch College and Posnania Inc First Baruch Volatility Workshop New York June 16-18, 2015 Outline Arbitrage free SABR 1 Arbitrage free
More informationStochastic Volatility (Working Draft I)
Stochastic Volatility (Working Draft I) Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu 1 Introduction When using the Black-Scholes-Merton model to price derivative
More informationEuropean option pricing under parameter uncertainty
European option pricing under parameter uncertainty Martin Jönsson (joint work with Samuel Cohen) University of Oxford Workshop on BSDEs, SPDEs and their Applications July 4, 2017 Introduction 2/29 Introduction
More informationSmile in the low moments
Smile in the low moments L. De Leo, T.-L. Dao, V. Vargas, S. Ciliberti, J.-P. Bouchaud 10 jan 2014 Outline 1 The Option Smile: statics A trading style The cumulant expansion A low-moment formula: the moneyness
More informationExtrapolation analytics for Dupire s local volatility
Extrapolation analytics for Dupire s local volatility Stefan Gerhold (joint work with P. Friz and S. De Marco) Vienna University of Technology, Austria 6ECM, July 2012 Implied vol and local vol Implied
More informationLecture 11: Stochastic Volatility Models Cont.
E4718 Spring 008: Derman: Lecture 11:Stochastic Volatility Models Cont. Page 1 of 8 Lecture 11: Stochastic Volatility Models Cont. E4718 Spring 008: Derman: Lecture 11:Stochastic Volatility Models Cont.
More informationLeverage Effect, Volatility Feedback, and Self-Exciting MarketAFA, Disruptions 1/7/ / 14
Leverage Effect, Volatility Feedback, and Self-Exciting Market Disruptions Liuren Wu, Baruch College Joint work with Peter Carr, New York University The American Finance Association meetings January 7,
More informationA Consistent Pricing Model for Index Options and Volatility Derivatives
A Consistent Pricing Model for Index Options and Volatility Derivatives 6th World Congress of the Bachelier Society Thomas Kokholm Finance Research Group Department of Business Studies Aarhus School of
More informationThe Black-Scholes Model
IEOR E4706: Foundations of Financial Engineering c 2016 by Martin Haugh The Black-Scholes Model In these notes we will use Itô s Lemma and a replicating argument to derive the famous Black-Scholes formula
More informationAn Analytical Approximation for Pricing VWAP Options
.... An Analytical Approximation for Pricing VWAP Options Hideharu Funahashi and Masaaki Kijima Graduate School of Social Sciences, Tokyo Metropolitan University September 4, 215 Kijima (TMU Pricing of
More informationCalibration Lecture 4: LSV and Model Uncertainty
Calibration Lecture 4: LSV and Model Uncertainty March 2017 Recap: Heston model Recall the Heston stochastic volatility model ds t = rs t dt + Y t S t dw 1 t, dy t = κ(θ Y t ) dt + ξ Y t dw 2 t, where
More informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationLecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing
Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing We shall go over this note quickly due to time constraints. Key concept: Ito s lemma Stock Options: A contract giving
More informationMonte Carlo Simulations
Monte Carlo Simulations Lecture 1 December 7, 2014 Outline Monte Carlo Methods Monte Carlo methods simulate the random behavior underlying the financial models Remember: When pricing you must simulate
More informationMSC FINANCIAL ENGINEERING PRICING I, AUTUMN LECTURE 9: LOCAL AND STOCHASTIC VOLATILITY RAYMOND BRUMMELHUIS DEPARTMENT EMS BIRKBECK
MSC FINANCIAL ENGINEERING PRICING I, AUTUMN 2010-2011 LECTURE 9: LOCAL AND STOCHASTIC VOLATILITY RAYMOND BRUMMELHUIS DEPARTMENT EMS BIRKBECK The only ingredient of the Black and Scholes formula which is
More informationPricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case
Pricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case Guang-Hua Lian Collaboration with Robert Elliott University of Adelaide Feb. 2, 2011 Robert Elliott,
More informationDefinition Pricing Risk management Second generation barrier options. Barrier Options. Arfima Financial Solutions
Arfima Financial Solutions Contents Definition 1 Definition 2 3 4 Contenido Definition 1 Definition 2 3 4 Definition Definition: A barrier option is an option on the underlying asset that is activated
More informationPricing of a European Call Option Under a Local Volatility Interbank Offered Rate Model
American Journal of Theoretical and Applied Statistics 2018; 7(2): 80-84 http://www.sciencepublishinggroup.com/j/ajtas doi: 10.11648/j.ajtas.20180702.14 ISSN: 2326-8999 (Print); ISSN: 2326-9006 (Online)
More informationEquity correlations implied by index options: estimation and model uncertainty analysis
1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to
More informationDevelopments in Volatility Derivatives Pricing
Developments in Volatility Derivatives Pricing Jim Gatheral Global Derivatives 2007 Paris, May 23, 2007 Motivation We would like to be able to price consistently at least 1 options on SPX 2 options on
More information2.1 Mathematical Basis: Risk-Neutral Pricing
Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t
More informationLinearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing
Linearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing Liuren Wu, Baruch College Joint work with Peter Carr and Xavier Gabaix at New York University Board of
More informationMultiscale Stochastic Volatility Models
Multiscale Stochastic Volatility Models Jean-Pierre Fouque University of California Santa Barbara 6th World Congress of the Bachelier Finance Society Toronto, June 25, 2010 Multiscale Stochastic Volatility
More informationSimple Robust Hedging with Nearby Contracts
Simple Robust Hedging with Nearby Contracts Liuren Wu and Jingyi Zhu Baruch College and University of Utah April 29, 211 Fourth Annual Triple Crown Conference Liuren Wu (Baruch) Robust Hedging with Nearby
More informationHedging Barrier Options through a Log-Normal Local Stochastic Volatility Model
22nd International Congress on Modelling and imulation, Hobart, Tasmania, Australia, 3 to 8 December 2017 mssanz.org.au/modsim2017 Hedging Barrier Options through a Log-Normal Local tochastic Volatility
More informationAn Overview of Volatility Derivatives and Recent Developments
An Overview of Volatility Derivatives and Recent Developments September 17th, 2013 Zhenyu Cui Math Club Colloquium Department of Mathematics Brooklyn College, CUNY Math Club Colloquium Volatility Derivatives
More informationA New Framework for Analyzing Volatility Risk and Premium Across Option Strikes and Expiries
A New Framework for Analyzing Volatility Risk and Premium Across Option Strikes and Expiries Liuren Wu, Baruch College Joint work with Peter Carr from Morgan Stanley Singapore Management University July
More information( ) since this is the benefit of buying the asset at the strike price rather
Review of some financial models for MAT 483 Parity and Other Option Relationships The basic parity relationship for European options with the same strike price and the same time to expiration is: C( KT
More informationHeston Stochastic Local Volatility Model
Heston Stochastic Local Volatility Model Klaus Spanderen 1 R/Finance 2016 University of Illinois, Chicago May 20-21, 2016 1 Joint work with Johannes Göttker-Schnetmann Klaus Spanderen Heston Stochastic
More informationMultiscale Stochastic Volatility Models Heston 1.5
Multiscale Stochastic Volatility Models Heston 1.5 Jean-Pierre Fouque Department of Statistics & Applied Probability University of California Santa Barbara Modeling and Managing Financial Risks Paris,
More informationLocal Volatility Dynamic Models
René Carmona Bendheim Center for Finance Department of Operations Research & Financial Engineering Princeton University Columbia November 9, 27 Contents Joint work with Sergey Nadtochyi Motivation 1 Understanding
More informationAN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL
AN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL FABIO MERCURIO BANCA IMI, MILAN http://www.fabiomercurio.it 1 Stylized facts Traders use the Black-Scholes formula to price plain-vanilla options. An
More informationLocal vs Non-local Forward Equations for Option Pricing
Local vs Non-local Forward Equations for Option Pricing Rama Cont Yu Gu Abstract When the underlying asset is a continuous martingale, call option prices solve the Dupire equation, a forward parabolic
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November
More informationSkew Hedging. Szymon Borak Matthias R. Fengler Wolfgang K. Härdle. CASE-Center for Applied Statistics and Economics Humboldt-Universität zu Berlin
Szymon Borak Matthias R. Fengler Wolfgang K. Härdle CASE-Center for Applied Statistics and Economics Humboldt-Universität zu Berlin 6 4 2.22 Motivation 1-1 Barrier options Knock-out options are financial
More informationSaddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance
Saddlepoint Approximation Methods for Pricing Financial Options on Discrete Realized Variance Yue Kuen KWOK Department of Mathematics Hong Kong University of Science and Technology Hong Kong * This is
More informationHeston Model Version 1.0.9
Heston Model Version 1.0.9 1 Introduction This plug-in implements the Heston model. Once installed the plug-in offers the possibility of using two new processes, the Heston process and the Heston time
More informationRohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)
Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,
More informationCEV Implied Volatility by VIX
CEV Implied Volatility by VIX Implied Volatility Chien-Hung Chang Dept. of Financial and Computation Mathematics, Providence University, Tiachng, Taiwan May, 21, 2015 Chang (Institute) Implied volatility
More informationNumerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps
Numerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps, Senior Quantitative Analyst Motivation: Swing Options An electricity or gas SUPPLIER needs to be capable,
More informationSensex Realized Volatility Index (REALVOL)
Sensex Realized Volatility Index (REALVOL) Introduction Volatility modelling has traditionally relied on complex econometric procedures in order to accommodate the inherent latent character of volatility.
More informationVaR Estimation under Stochastic Volatility Models
VaR Estimation under Stochastic Volatility Models Chuan-Hsiang Han Dept. of Quantitative Finance Natl. Tsing-Hua University TMS Meeting, Chia-Yi (Joint work with Wei-Han Liu) December 5, 2009 Outline Risk
More informationOptimal Hedging of Variance Derivatives. John Crosby. Centre for Economic and Financial Studies, Department of Economics, Glasgow University
Optimal Hedging of Variance Derivatives John Crosby Centre for Economic and Financial Studies, Department of Economics, Glasgow University Presentation at Baruch College, in New York, 16th November 2010
More informationShort-time-to-expiry expansion for a digital European put option under the CEV model. November 1, 2017
Short-time-to-expiry expansion for a digital European put option under the CEV model November 1, 2017 Abstract In this paper I present a short-time-to-expiry asymptotic series expansion for a digital European
More informationA Brief Introduction to Stochastic Volatility Modeling
A Brief Introduction to Stochastic Volatility Modeling Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction When using the Black-Scholes-Merton model to
More informationThe vanna-volga method for implied volatilities
CUTTING EDGE. OPTION PRICING The vanna-volga method for implied volatilities The vanna-volga method is a popular approach for constructing implied volatility curves in the options market. In this article,
More informationEmpirical Distribution Testing of Economic Scenario Generators
1/27 Empirical Distribution Testing of Economic Scenario Generators Gary Venter University of New South Wales 2/27 STATISTICAL CONCEPTUAL BACKGROUND "All models are wrong but some are useful"; George Box
More informationChapter 6 Forecasting Volatility using Stochastic Volatility Model
Chapter 6 Forecasting Volatility using Stochastic Volatility Model Chapter 6 Forecasting Volatility using SV Model In this chapter, the empirical performance of GARCH(1,1), GARCH-KF and SV models from
More informationLocally risk-minimizing vs. -hedging in stochastic vola
Locally risk-minimizing vs. -hedging in stochastic volatility models University of St. Andrews School of Economics and Finance August 29, 2007 joint work with R. Poulsen ( Kopenhagen )and K.R.Schenk-Hoppe
More information1.1 Basic Financial Derivatives: Forward Contracts and Options
Chapter 1 Preliminaries 1.1 Basic Financial Derivatives: Forward Contracts and Options A derivative is a financial instrument whose value depends on the values of other, more basic underlying variables
More informationA Closed-form Solution for Outperfomance Options with Stochastic Correlation and Stochastic Volatility
A Closed-form Solution for Outperfomance Options with Stochastic Correlation and Stochastic Volatility Jacinto Marabel Romo Email: jacinto.marabel@grupobbva.com November 2011 Abstract This article introduces
More informationMachine Learning for Quantitative Finance
Machine Learning for Quantitative Finance Fast derivative pricing Sofie Reyners Joint work with Jan De Spiegeleer, Dilip Madan and Wim Schoutens Derivative pricing is time-consuming... Vanilla option pricing
More informationValuation of Volatility Derivatives. Jim Gatheral Global Derivatives & Risk Management 2005 Paris May 24, 2005
Valuation of Volatility Derivatives Jim Gatheral Global Derivatives & Risk Management 005 Paris May 4, 005 he opinions expressed in this presentation are those of the author alone, and do not necessarily
More informationPortfolio Management Using Option Data
Portfolio Management Using Option Data Peter Christoffersen Rotman School of Management, University of Toronto, Copenhagen Business School, and CREATES, University of Aarhus 2 nd Lecture on Friday 1 Overview
More informationA SUMMARY OF OUR APPROACHES TO THE SABR MODEL
Contents 1 The need for a stochastic volatility model 1 2 Building the model 2 3 Calibrating the model 2 4 SABR in the risk process 5 A SUMMARY OF OUR APPROACHES TO THE SABR MODEL Financial Modelling Agency
More informationManaging the Newest Derivatives Risks
Managing the Newest Derivatives Risks Michel Crouhy IXIS Corporate and Investment Bank / A subsidiary of NATIXIS Derivatives 2007: New Ideas, New Instruments, New markets NYU Stern School of Business,
More informationEconomic Scenario Generator: Applications in Enterprise Risk Management. Ping Sun Executive Director, Financial Engineering Numerix LLC
Economic Scenario Generator: Applications in Enterprise Risk Management Ping Sun Executive Director, Financial Engineering Numerix LLC Numerix makes no representation or warranties in relation to information
More informationLecture 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, 2003 Abstract
More informationHedging Under Jump Diffusions with Transaction Costs. Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo
Hedging Under Jump Diffusions with Transaction Costs Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo Computational Finance Workshop, Shanghai, July 4, 2008 Overview Overview Single factor
More informationApplication of Moment Expansion Method to Option Square Root Model
Application of Moment Expansion Method to Option Square Root Model Yun Zhou Advisor: Professor Steve Heston University of Maryland May 5, 2009 1 / 19 Motivation Black-Scholes Model successfully explain
More informationSTOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL
STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL YOUNGGEUN YOO Abstract. Ito s lemma is often used in Ito calculus to find the differentials of a stochastic process that depends on time. This paper will introduce
More informationEstimation of Stochastic Volatility Models with Implied. Volatility Indices and Pricing of Straddle Option
Estimation of Stochastic Volatility Models with Implied Volatility Indices and Pricing of Straddle Option Yue Peng Steven C. J. Simon June 14, 29 Abstract Recent market turmoil has made it clear that modelling
More informationZABR -- Expansions for the Masses
ZABR -- Expansions for the Masses Preliminary Version December 011 Jesper Andreasen and Brian Huge Danse Marets, Copenhagen want.daddy@danseban.com brno@danseban.com 1 Electronic copy available at: http://ssrn.com/abstract=198076
More informationUnifying Volatility Models
The University of Reading THE BUSINESS SCHOOL FOR FINANCIAL MARKETS Unifying Volatility Models Carol Alexander (Co-authored works with E. Lazar and L. Nogueira) ISMA Centre, University of Reading Email:
More informationModeling the Implied Volatility Surface. Jim Gatheral Global Derivatives and Risk Management 2003 Barcelona May 22, 2003
Modeling the Implied Volatility Surface Jim Gatheral Global Derivatives and Risk Management 2003 Barcelona May 22, 2003 This presentation represents only the personal opinions of the author and not those
More informationThe Use of Importance Sampling to Speed Up Stochastic Volatility Simulations
The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations Stan Stilger June 6, 1 Fouque and Tullie use importance sampling for variance reduction in stochastic volatility simulations.
More information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationHedging under Model Uncertainty
Hedging under Model Uncertainty Efficient Computation of the Hedging Error using the POD 6th World Congress of the Bachelier Finance Society June, 24th 2010 M. Monoyios, T. Schröter, Oxford University
More informationVega Maps: Predicting Premium Change from Movements of the Whole Volatility Surface
Vega Maps: Predicting Premium Change from Movements of the Whole Volatility Surface Ignacio Hoyos Senior Quantitative Analyst Equity Model Validation Group Risk Methodology Santander Alberto Elices Head
More informationMarket risk measurement in practice
Lecture notes on risk management, public policy, and the financial system Allan M. Malz Columbia University 2018 Allan M. Malz Last updated: October 23, 2018 2/32 Outline Nonlinearity in market risk Market
More informationLIBOR models, multi-curve extensions, and the pricing of callable structured derivatives
Weierstrass Institute for Applied Analysis and Stochastics LIBOR models, multi-curve extensions, and the pricing of callable structured derivatives John Schoenmakers 9th Summer School in Mathematical Finance
More informationGeneralized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models
Generalized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models Bilkan Erkmen (joint work with Michael Coulon) Workshop on Stochastic Games, Equilibrium, and Applications
More informationLinear-Rational Term-Structure Models
Linear-Rational Term-Structure Models Anders Trolle (joint with Damir Filipović and Martin Larsson) Ecole Polytechnique Fédérale de Lausanne Swiss Finance Institute AMaMeF and Swissquote Conference, September
More informationOption Pricing for a Stochastic-Volatility Jump-Diffusion Model
Option Pricing for a Stochastic-Volatility Jump-Diffusion Model Guoqing Yan and Floyd B. Hanson Department of Mathematics, Statistics, and Computer Science University of Illinois at Chicago Conference
More informationEconophysics V: Credit Risk
Fakultät für Physik Econophysics V: Credit Risk Thomas Guhr XXVIII Heidelberg Physics Graduate Days, Heidelberg 2012 Outline Introduction What is credit risk? Structural model and loss distribution Numerical
More informationEmpirical Approach to the Heston Model Parameters on the Exchange Rate USD / COP
Empirical Approach to the Heston Model Parameters on the Exchange Rate USD / COP ICASQF 2016, Cartagena - Colombia C. Alexander Grajales 1 Santiago Medina 2 1 University of Antioquia, Colombia 2 Nacional
More informationDELTA HEDGING VEGA RISK?
DELTA HEDGING VEGA RISK? Stéphane CRÉPEY, Évry University stephane.crepey@univ-evry.fr QMF Conference, Sydney, December 17 2004 Figure 1: The Volatility Smile (E. Derman) Contents 1 Basics of the smile
More informationAnalytical formulas for local volatility model with stochastic. Mohammed Miri
Analytical formulas for local volatility model with stochastic rates Mohammed Miri Joint work with Eric Benhamou (Pricing Partners) and Emmanuel Gobet (Ecole Polytechnique Modeling and Managing Financial
More informationStochastic Volatility and Jump Modeling in Finance
Stochastic Volatility and Jump Modeling in Finance HPCFinance 1st kick-off meeting Elisa Nicolato Aarhus University Department of Economics and Business January 21, 2013 Elisa Nicolato (Aarhus University
More informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationLeverage Effect, Volatility Feedback, and Self-Exciting Market Disruptions 11/4/ / 24
Leverage Effect, Volatility Feedback, and Self-Exciting Market Disruptions Liuren Wu, Baruch College and Graduate Center Joint work with Peter Carr, New York University and Morgan Stanley CUNY Macroeconomics
More informationMulti-factor Stochastic Volatility Models A practical approach
Stockholm School of Economics Department of Finance - Master Thesis Spring 2009 Multi-factor Stochastic Volatility Models A practical approach Filip Andersson 20573@student.hhs.se Niklas Westermark 20653@student.hhs.se
More informationFINANCIAL MATHEMATICS WITH ADVANCED TOPICS MTHE7013A
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 2016 17 FINANCIAL MATHEMATICS WITH ADVANCED TOPICS MTHE7013A Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other
More informationChapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets
Chapter 5: Jump Processes and Incomplete Markets Jumps as One Explanation of Incomplete Markets It is easy to argue that Brownian motion paths cannot model actual stock price movements properly in reality,
More informationCounterparty Credit Risk Simulation
Counterparty Credit Risk Simulation Alex Yang FinPricing http://www.finpricing.com Summary Counterparty Credit Risk Definition Counterparty Credit Risk Measures Monte Carlo Simulation Interest Rate Curve
More informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU CBOE Conference on Derivatives and Volatility, Chicago, Nov. 10, 2017 Peter Carr (NYU) Volatility Smiles and Yield Frowns 11/10/2017 1 / 33 Interest Rates
More information