Short-Time Asymptotic Methods in Financial Mathematics

Size: px
Start display at page:

Download "Short-Time Asymptotic Methods in Financial Mathematics"

Transcription

1 Short-Time Asymptotic Methods in Financial Mathematics José E. Figueroa-López Department of Mathematics Washington University in St. Louis Probability and Mathematical Finance Seminar Department of Mathematical Sciences Carnegie Mellon University February 20, 2017 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 1 / 39

2 Outline 1 Short-time asymptotics of Option Prices Short-time ATM Skew Asymptotics A Calibration Case-Study 2 High-Frequency Based Estimation Methods Multipower Variations and Truncated Realized Variations via Expected number of jump misclassifications via conditional Mean Square Error (cmse) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 2 / 39

3 Short-time asymptotics of Option Prices Outline 1 Short-time asymptotics of Option Prices Short-time ATM Skew Asymptotics A Calibration Case-Study 2 High-Frequency Based Estimation Methods Multipower Variations and Truncated Realized Variations via Expected number of jump misclassifications via conditional Mean Square Error (cmse) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 3 / 39

4 Short-time asymptotics of Option Prices Short-time ATM Skew Asymptotics Outline 1 Short-time asymptotics of Option Prices Short-time ATM Skew Asymptotics A Calibration Case-Study 2 High-Frequency Based Estimation Methods Multipower Variations and Truncated Realized Variations via Expected number of jump misclassifications via conditional Mean Square Error (cmse) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 4 / 39

5 Short-time asymptotics of Option Prices Short-time ATM Skew Asymptotics Short-time ATM Skew Asymptotics Theorem (F-L & Ólafsson, F&S 2016) For Y (1, 2), the implied volatility ˆσ (κ, t) of a vanilla option with log-moneyness κ and time-to-maturity t is such that d 1 t 1 2, if pure jump model with C + C, ˆσ (κ, t) d 1 t Y, if pure jump model with C + = C, κ κ=0 d 1 t 1 Y 2, if mixed model with C + C, d 1, if mixed model with C + = C. Furthermore, in the asymmetric cases, the sign of d 1 is the same as that of C + C. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 5 / 39

6 Short-time asymptotics of Option Prices A Calibration Case-Study Outline 1 Short-time asymptotics of Option Prices Short-time ATM Skew Asymptotics A Calibration Case-Study 2 High-Frequency Based Estimation Methods Multipower Variations and Truncated Realized Variations via Expected number of jump misclassifications via conditional Mean Square Error (cmse) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 6 / 39

7 Short-time asymptotics of Option Prices A Calibration Case-Study Calibration of S&P 500 Option Prices Objectives: To assess the plausibility of the divergence of the skew to as a power of time-to-maturity. Calibrate the parameter Y, which measures the activity of small jumps in the price process. Data: Daily closing bid and ask prices for S&P 500 index options, across all strikes K, and maturities t from Jan. 2 to Jun. 30, J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 7 / 39

8 Short-time asymptotics of Option Prices A Calibration Case-Study Implied Volatility Smile and ATM Skew (a) 0.14 (b) ATM skew -6-8 Implied vol Time-to-maturity (t) log-moneyness (κ) Figure: (a) The graphs of the ATM implied volatility skew are consistent with divergence to as a power law. (b) The implied volatility smiles on Jan. 15, 2014, corresponding to maturities ranging from (8 days) to 0.25 (3 months), show greater steepness as time-to-maturity decreases, which is consistent with the presence of jumps. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 8 / 39

9 Short-time asymptotics of Option Prices A Calibration Case-Study Calibration of the Index Of Jump Activity Y (a) 1-day regression 21-day regression (b) Y pj 1.5 Y cc day regression 21-day regression 1 Jan Feb Mar Apr May Jun Day 1 Jan Feb Mar Apr May Jun Day Figure: The Y -estimates implied by daily and monthly regressions of on log t assuming a pure-jump model (panel (a)), and a mixed κ=0 log ˆσ(κ,t) κ model (panel (b)). The first and third quantiles of the regressions R 2 are and 0.995, respectively. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 9 / 39

10 Short-time asymptotics of Option Prices A Calibration Case-Study ATM Implied Volatility (a) (b) VIX ATM vol. ATM implied vol Volatility Time-to-maturity (t) 0.06 Jan Feb Mar Apr May Jun Day Figure: (a) ATM implied volatility as a function of time-to-maturity in years, for each business day in Jan Jun, (b) ATM implied volatility for the shortest outstanding maturity compared to the VIX index. These suggest a mixed model rather than a pure jump model since ˆσ(t, 0) σ 0 > 0. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 10 / 39

11 Short-time asymptotics of Option Prices Outline 1 Short-time asymptotics of Option Prices Short-time ATM Skew Asymptotics A Calibration Case-Study 2 High-Frequency Based Estimation Methods Multipower Variations and Truncated Realized Variations via Expected number of jump misclassifications via conditional Mean Square Error (cmse) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 11 / 39

12 Short-time asymptotics of Option Prices 1 Besides the skew, the convexity of the implied volatility, κ=0, is very important in trading and, thus, its asymptotic 2ˆσ(κ,t) κ 2 behavior is of great interest. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 12 / 39

13 Short-time asymptotics of Option Prices 1 Besides the skew, the convexity of the implied volatility, κ=0, is very important in trading and, thus, its asymptotic 2ˆσ(κ,t) κ 2 behavior is of great interest. 2 Small-time asymptotics of options written on Leveraged Exchange-Traded Fund (LETF). LETF is a managed portfolio, which seeks to multiply the instantaneous returns of a reference Exchange-Traded Fund (ETF). J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 12 / 39

14 Short-time asymptotics of Option Prices 1 Besides the skew, the convexity of the implied volatility, κ=0, is very important in trading and, thus, its asymptotic 2ˆσ(κ,t) κ 2 behavior is of great interest. 2 Small-time asymptotics of options written on Leveraged Exchange-Traded Fund (LETF). LETF is a managed portfolio, which seeks to multiply the instantaneous returns of a reference Exchange-Traded Fund (ETF). 3 There is no known results for American options. It is expected that the early exercise feature won t affect the leading order terms of the asymptotics but they should contribute to the high-order terms. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 12 / 39

15 Outline 1 Short-time asymptotics of Option Prices Short-time ATM Skew Asymptotics A Calibration Case-Study 2 High-Frequency Based Estimation Methods Multipower Variations and Truncated Realized Variations via Expected number of jump misclassifications via conditional Mean Square Error (cmse) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 13 / 39

16 Merton Log-Normal Model Consider the following model for the log-return process X t = log St S 0 of a financial asset: X t = at + σw t + ζ i, i.i.d. ζ i N(µ jmp, σjmp), 2 {τ i } i 1 Poisson(λ) i:τ i t Goal: Estimate the volatility σ based on a discrete record X t1,..., X tn Log-Normal Merton Model Stock Price Price process Continuous process Times of jumps Time in years (252 days) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 14 / 39

17 Multipower Variations and Truncated Realized Variations Outline 1 Short-time asymptotics of Option Prices Short-time ATM Skew Asymptotics A Calibration Case-Study 2 High-Frequency Based Estimation Methods Multipower Variations and Truncated Realized Variations via Expected number of jump misclassifications via conditional Mean Square Error (cmse) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 15 / 39

18 Multipower Variations and Truncated Realized Variations Two main classes of estimators Precursor. Realized Quadratic Variation: QV [X] n := 1 T n 1 ( Xti+1 X ) 2 t i i=0 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 16 / 39

19 Multipower Variations and Truncated Realized Variations Two main classes of estimators Precursor. Realized Quadratic Variation: QV [X] n := 1 T n 1 ( Xti+1 X ) 2 n t i σ 2 + i=0 j:τ j T ζ 2 j J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 16 / 39

20 Multipower Variations and Truncated Realized Variations Two main classes of estimators Precursor. Realized Quadratic Variation: QV [X] n := 1 T n 1 ( Xti+1 X ) 2 n t i σ 2 + i=0 j:τ j T 1 Bipower Realized Variations (Barndorff-Nielsen and Shephard): n 1 BPV [X] n := Xti+1 X ti Xti+2 X ti+1, i=0 ζ 2 j J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 16 / 39

21 Multipower Variations and Truncated Realized Variations Two main classes of estimators Precursor. Realized Quadratic Variation: QV [X] n := 1 T n 1 ( Xti+1 X ) 2 n t i σ 2 + i=0 j:τ j T 1 Bipower Realized Variations (Barndorff-Nielsen and Shephard): n 1 BPV [X] n := Xti+1 X ti Xti+2 X ti+1, i=0 2 Truncated Realized Variations (Mancini): n 1 ( TRV n [X](ε) := Xti+1 X ) 2 t i 1 { Xti+1 X ti i=0 ε ζ 2 j }, (ε [0, )). J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 16 / 39

22 Multipower Variations and Truncated Realized Variations Calibration or Tuning of the Estimator Problem: How do you choose the threshold parameter ε? J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 17 / 39

23 Multipower Variations and Truncated Realized Variations Calibration or Tuning of the Estimator Problem: How do you choose the threshold parameter ε? Log-Normal Merton Model Performace of Truncated Realized Variations Stock Price Price process Continuous process Times of jumps Truncated Realized Variation (TRV) TRV True Volatility (0.4) Realized Variation (0.51) Bipower Variation (0.42) Time in years (252 days) Truncation level (epsilon) Figure: (left) 5-min Merton observations with σ = 0.4, σ jmp = 3(h), µ jmp = 0, λ = 200; (right) TRV estimates for all the truncation levels. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 17 / 39

24 Multipower Variations and Truncated Realized Variations Calibration or Tuning of the Estimator Problem: How do you choose the threshold parameter ε? Log-Normal Merton Model Performace of Truncated Realized Variations Stock Price Price process Continuous process Times of jumps Truncated Realized Variation (TRV) TRV True Volatility (0.2) Realized Variation (0.4) Bipower Variation (0.26) Cont. Realized Variation (0.19) Time in years (252 days) Truncation level (epsilon) Figure: (left) 5 minute Merton observations with σ = 0.2, σ jmp = 1.5(h), µ jmp = 0, λ = 1000; (right) TRV performance wrt the truncation level J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 18 / 39

25 Multipower Variations and Truncated Realized Variations Popular truncations ε Literature consists of mostly ad hoc" selection methods for ε, aimed to satisfy sufficient conditions for the consistency and asymptotic normality of the associated estimators. The most popular is the so-called Power Threshold: ε Pwr α,ω := α h ω, for α > 0 and ω (0, 1/2). The rule of thumb says to use a value of ω close to 1/2 (typically, 0.495). J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 19 / 39

26 Outline 1 Short-time asymptotics of Option Prices Short-time ATM Skew Asymptotics A Calibration Case-Study 2 High-Frequency Based Estimation Methods Multipower Variations and Truncated Realized Variations via Expected number of jump misclassifications via conditional Mean Square Error (cmse) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 20 / 39

27 General Approach 1 Fix a suitable and sensible metric for estimation error; say, [ ( ) ] 1 2 MSE(ε) = E T TRV n(ε) σ 2 2 Show the existence of the optimal threshold ε n minimizing the error function; 3 Characterize the asymptotic behavior ε n (when n ) to infer qualitative information such as rate of convergence on n and dependence on the underlying parameters of the model (σ, σ J, λ) Devise a plug-in type implementation of ε by estimating those parameters (if possible). J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 21 / 39

28 Via Expected number of jump misclassifications J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 22 / 39

29 Via Expected number of jump misclassifications 1 Estimation Error: (F-L & Nisen, SPA 2013) 2 Notation: Loss n (ε) := E n i=1 N t := # of jumps by time t n i X := X t i X ti 1 (1 [ ni X >ε, ni N=0] + 1 [ ni X ε, ni N 0] ). n i N := N t i N ti 1 = # of jumps during (t i 1, t i ] J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 22 / 39

30 Existence and Infill Asymptotic Characterization Theorem (F-L & Nisen, SPA 2013)) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 23 / 39

31 Existence and Infill Asymptotic Characterization Theorem (F-L & Nisen, SPA 2013)) 1 For n large enough, the loss function Loss n (ε) is convex and, moreover, enjoys a unique global minimum ε n. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 23 / 39

32 Existence and Infill Asymptotic Characterization Theorem (F-L & Nisen, SPA 2013)) 1 For n large enough, the loss function Loss n (ε) is convex and, moreover, enjoys a unique global minimum ε n. 2 As n, the optimal threshold sequence (ε n) n is such that ( ) 1 ε n = 3σ 2 h n log + h.o.t., where h n := t i t i 1 is the time step. h n J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 23 / 39

33 Existence and Infill Asymptotic Characterization Theorem (F-L & Nisen, SPA 2013)) 1 For n large enough, the loss function Loss n (ε) is convex and, moreover, enjoys a unique global minimum ε n. 2 As n, the optimal threshold sequence (ε n) n is such that ( ) 1 ε n = 3σ 2 h n log + h.o.t., where h n := t i t i 1 is the time step. Remark: The leading order term is proportional to the Lévy s modulus of continuity of the Brownian motion!. h n J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 23 / 39

34 A Feasible Implementation based on ε n (i) Get a rough" estimate of σ 2 via, e.g., realized QV: ˆσ n,0 2 := 1 n X ti X ti 1 2 T i=1 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 24 / 39

35 A Feasible Implementation based on ε n (i) Get a rough" estimate of σ 2 via, e.g., realized QV: ˆσ n,0 2 := 1 n X ti X ti 1 2 T i=1 (ii) Use ˆσ n,0 2 to estimate the optimal threshold ˆε n,0 := (3 ˆσ 2 n,0 h n log(1/h n ) ) 1/2 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 24 / 39

36 A Feasible Implementation based on ε n (i) Get a rough" estimate of σ 2 via, e.g., realized QV: ˆσ n,0 2 := 1 n X ti X ti 1 2 T i=1 (ii) Use ˆσ n,0 2 to estimate the optimal threshold ˆε n,0 := (3 ˆσ 2 n,0 h n log(1/h n ) ) 1/2 (iii) Refine ˆσ n,0 2 using thresholding, ˆσ n,1 2 = 1 n X ti X ti T [ X ti X ti 1 ˆε n,0 i=1 ] J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 24 / 39

37 A Feasible Implementation based on ε n (i) Get a rough" estimate of σ 2 via, e.g., realized QV: ˆσ n,0 2 := 1 n X ti X ti 1 2 T i=1 (ii) Use ˆσ n,0 2 to estimate the optimal threshold ˆε n,0 := (3 ˆσ 2 n,0 h n log(1/h n ) ) 1/2 (iii) Refine ˆσ n,0 2 using thresholding, ˆσ n,1 2 = 1 n X ti X ti T [ X ti X ti 1 ˆε n,0 i=1 (iv) Iterate Steps (ii) and (iii): ˆσ n,0 2 ˆε n,0 ˆσ2 n,1 ˆε n,1 ˆσ2 n,2 ˆσ2 n, ] J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 24 / 39

38 Illustration I Log-Normal Merton Model Performace of Truncated Realized Variations Stock Price Truncated Realized Variation (TRV) TRV True Volatility (0.4) Realized Variation (0.46) Bipower Variation (0.41) Cont. Realized Variation (0.402) Time in years (252 days) Truncation level (epsilon) Figure: (left) Merton Model with σ = 0.4, σ jmp = 3(h), µ jmp = 0, λ = 200; (right) TRV performance wrt the truncation level. Red dot is ˆσ n,1 = 0.409, while purple dot is the limiting estimator ˆσ n, = J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 25 / 39

39 Illustration II Log-Normal Merton Model Performace of Truncated Realized Variations Stock Price Truncated Realized Variation (TRV) TRV True Volatility (0.2) Realized Variation (0.33) Bipower Variation (0.24) Cont. Realized Variation (0.195) Time in years (252 days) Truncation level (epsilon) Figure: (left) Merton Model with σ = 0.2, σ jmp = 1.5(h), µ jmp = 0, λ = 1000; (right) TRV performance wrt the truncation level. Red dot is ˆσ n,1 = 0.336, while purple dot is the limiting estimator ˆσ n, = J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 26 / 39

40 Via conditional Mean Square Error (cmse) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 27 / 39

41 Via conditional Mean Square Error (cmse) 1 We now propose a second approach in which we aim to control ( T 2 MSE c (ε) := E TRV n (ε) σsds) 2 0 σ, J. This is in the more general semimartingale setting: where J is a pure-jump process. t X t = X 0 + σ s dw s + J t 0 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 27 / 39

42 Via conditional Mean Square Error (cmse) 1 We now propose a second approach in which we aim to control ( T 2 MSE c (ε) := E TRV n (ε) σsds) 2 0 σ, J. This is in the more general semimartingale setting: where J is a pure-jump process. 2 Key Assumptions: t X t = X 0 + σ s dw s + J t 0 σ t > 0, t, and σ and J are independent of W. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 27 / 39

43 Key Relationships b i (ε) := E [ ( i X) 2 1 { i X ε} σ, J ] J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 28 / 39

44 Key Relationships b i (ε) := E [ ( i X) 2 1 { i X ε} σ, J ] Then, with the notation m i := n i J = J ti J ti 1, σi 2 := ti t i 1 σ 2 s ds, J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 28 / 39

45 Key Relationships Then, with the notation b i (ε) := E [ ( i X) 2 1 { i X ε} σ, J ] we have m i := n i J = J ti J ti 1, σ 2 i := b i (ε) = σ ( i 2π e (ε mi ) 2 ti t i 1 σ 2 s ds, 2σ i 2 (ε + m i ) + e (ε+m i ) 2 ) 2σ i 2 (ε m i ) ε+m i + (mi 2 + σi 2 σ i ) ε m i σ i e x 2 /2 dx 2π J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 28 / 39

46 Key Relationships Then, with the notation b i (ε) := E [ ( i X) 2 1 { i X ε} σ, J ] we have m i := n i J = J ti J ti 1, σ 2 i := b i (ε) = σ ( i 2π e (ε mi ) 2 ti t i 1 σ 2 s ds, 2σ i 2 (ε + m i ) + e (ε+m i ) 2 ) 2σ i 2 (ε m i ) ε+m i + (mi 2 + σi 2 σ i ) ε m i σ i e x 2 /2 dx 2π db i (ε) dε = J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 28 / 39

47 Key Relationships Then, with the notation b i (ε) := E [ ( i X) 2 1 { i X ε} σ, J ] we have m i := n i J = J ti J ti 1, σ 2 i := b i (ε) = σ ( i 2π e (ε mi ) 2 ti t i 1 σ 2 s ds, 2σ i 2 (ε + m i ) + e (ε+m i ) 2 ) 2σ i 2 (ε m i ) db i (ε) dε ε+m i + (mi 2 + σi 2 σ i ) ε m i σ i e x 2 /2 dx 2π = ε 2 a i (ε), with a i (ε) := e (ε m i ) 2 2σ 2 i + e (ε+m i ) σ i 2π 2σ 2 i 2 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 28 / 39

48 Conditional Mean Square Error (MSE c ) Theorem (F-L & Mancini, 2017) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 29 / 39

49 Conditional Mean Square Error (MSE c ) Theorem (F-L & Mancini, 2017) 1 MSE c (ε) := E [(TRV n (ε) IV ) 2 ] σ, J is such that d MSE c (ε) = ε 2 G(ε), dε where G(ε) := ( i a i(ε) ε j i b j(ε) 2 ) j σ2 j. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 29 / 39

50 Conditional Mean Square Error (MSE c ) Theorem (F-L & Mancini, 2017) 1 MSE c (ε) := E [(TRV n (ε) IV ) 2 ] σ, J is such that d MSE c (ε) = ε 2 G(ε), dε where G(ε) := ( i a i(ε) ε j i b j(ε) 2 ) j σ2 j. 2 Furthermore, there exists an optimal threshold ε n MSE c (ε) and is such that G(ε n ) = 0. that minimizes J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 29 / 39

51 Asymptotics: FA process with constant variance Theorem (F-L & Mancini, 2017) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 30 / 39

52 Asymptotics: FA process with constant variance Theorem (F-L & Mancini, 2017) 1 Suppose that σ t σ is constant and J is a finite jump activity process (with or without drift; not necessarily Lévy): N t X t = σw t + j=1 ζ j J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 30 / 39

53 Asymptotics: FA process with constant variance Theorem (F-L & Mancini, 2017) 1 Suppose that σ t σ is constant and J is a finite jump activity process (with or without drift; not necessarily Lévy): N t X t = σw t + j=1 2 Then, as n, the optimal threshold ε n cmse is such that ε n ( 1 2σ 2 h n log h n ζ j that minimizes the ), h n 0. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 30 / 39

54 A Feasible Implementation based on ε n (i) Get a rough" estimate of σ 2 via, e.g., the QV: σ n,0 2 := 1 n X ti X ti 1 2 T i=1 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 31 / 39

55 A Feasible Implementation based on ε n (i) Get a rough" estimate of σ 2 via, e.g., the QV: σ n,0 2 := 1 n X ti X ti 1 2 T i=1 (ii) Use σ n,0 2 to estimate the optimal threshold ˆε n,0 := (2 σ 2 n,0 h n log(1/h n ) ) 1/2 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 31 / 39

56 A Feasible Implementation based on ε n (i) Get a rough" estimate of σ 2 via, e.g., the QV: σ n,0 2 := 1 n X ti X ti 1 2 T i=1 (ii) Use σ n,0 2 to estimate the optimal threshold ˆε n,0 := (2 σ 2 n,0 h n log(1/h n ) ) 1/2 (iii) Refine σ n,0 2 using thresholding, σ n,1 2 = 1 n X ti X ti T [ ] X ti X ti 1 ˆε n,0 i=1 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 31 / 39

57 A Feasible Implementation based on ε n (i) Get a rough" estimate of σ 2 via, e.g., the QV: σ n,0 2 := 1 n X ti X ti 1 2 T i=1 (ii) Use σ n,0 2 to estimate the optimal threshold ˆε n,0 := (2 σ 2 n,0 h n log(1/h n ) ) 1/2 (iii) Refine σ n,0 2 using thresholding, σ n,1 2 = 1 n X ti X ti T [ ] X ti X ti 1 ˆε n,0 i=1 (iv) Iterate Steps (ii) and (iii): σ 2 n,0 ˆε n,0 σ2 n,1 ˆε n,1 σ2 n,2 σ2 n, J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 31 / 39

58 Illustration II. Continued... Log-Normal Merton Model Performace of Truncated Realized Variations Stock Price Truncated Realized Variation (TRV) TRV True Volatility (0.2) Realized Variation (0.33) Bipower Variation (0.24) Cont. Realized Variation (0.195) Time in years (252 days) Truncation level (epsilon) Figure: (left) Merton Model with λ = Red dot is ˆσ n,1 = 0.336, while purple dot is the limiting ˆσ n,k = Orange square is σ n,1 = 0.225, while brown square is the limiting estimator σ n, = J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 32 / 39

59 Monte Carlo Simulations Estimator ˆσ std(ˆσ) Loss ε N RQV ˆσ n, ˆσ n, σ n, σ n, BPV Table: Estimation of the volatility σ = 0.2 for a log-normal Merton model based on simulations of 5-minute observations over a 1 month time horizon. The jump parameters are λ = 1000, σ Jmp = 1.5 h and µ Jmp = 0. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 33 / 39

60 Outline 1 Short-time asymptotics of Option Prices Short-time ATM Skew Asymptotics A Calibration Case-Study 2 High-Frequency Based Estimation Methods Multipower Variations and Truncated Realized Variations via Expected number of jump misclassifications via conditional Mean Square Error (cmse) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 34 / 39

61 1 In principle, we can apply the proposed method for varying volatility t σ t by localization; i.e., applying it to periods where σ is approximately constant. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 35 / 39

62 1 In principle, we can apply the proposed method for varying volatility t σ t by localization; i.e., applying it to periods where σ is approximately constant. 2 However, is it possible to analyze the asymptotic behavior of MSE c = ε 2 G(ε) in terms of certain estimable summary measures? Say, G(ε) G 0 (ε, m 1,..., m n, σ, σ), where σ := inf t T σ t and σ := sup t T σ t J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 35 / 39

63 3 As it turns, for a Lévy process J and constant σ, the expected [ MSE(ε) := E (TRV n (ε) IV ) 2] is such that d dε MSE(ε) = nε2 E[a 1 (ε)] (ε 2 + 2(n 1)E[b 1 (ε)] 2nh n σ 2) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 36 / 39

64 3 As it turns, for a Lévy process J and constant σ, the expected [ MSE(ε) := E (TRV n (ε) IV ) 2] is such that d dε MSE(ε) = nε2 E[a 1 (ε)] (ε 2 + 2(n 1)E[b 1 (ε)] 2nh n σ 2) Therefore, there exists a unique minimum point ε n solution of the equation ε 2 + 2(n 1)E[b 1 (ε)] 2nh n σ 2 = 0. which is a J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 36 / 39

65 3 Furthermore, J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 37 / 39

66 3 Furthermore, As expected, in the finite jump activity case, it turns out ( ) 1 2σ 2 h n log ε n h n J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 37 / 39

67 3 Furthermore, As expected, in the finite jump activity case, it turns out ( ) 1 2σ 2 h n log ε n But, surprisingly, if J is a Y -stable Lévy process (IA), ( ) 1 (2 Y )σ 2 h n log ε n Thus, the higher the jump activity is, the lower the optimal threshold has to be to discard the higher noise represented by the jumps. h n h n J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 37 / 39

68 3 Furthermore, As expected, in the finite jump activity case, it turns out ( ) 1 2σ 2 h n log ε n But, surprisingly, if J is a Y -stable Lévy process (IA), ( ) 1 (2 Y )σ 2 h n log ε n Thus, the higher the jump activity is, the lower the optimal threshold has to be to discard the higher noise represented by the jumps. Does this phenomenon holds for the minimizer ε n h n h n of the cmse? J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 37 / 39

69 3 Furthermore, As expected, in the finite jump activity case, it turns out ( ) 1 2σ 2 h n log ε n But, surprisingly, if J is a Y -stable Lévy process (IA), ( ) 1 (2 Y )σ 2 h n log ε n Thus, the higher the jump activity is, the lower the optimal threshold has to be to discard the higher noise represented by the jumps. Does this phenomenon holds for the minimizer ε n of the cmse? Can we generalize it to Lévy processes with stable like jumps? h n h n J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 37 / 39

70 Further Reading I J.E. Figueroa-López, R. Gong, and C. Houdré, High-order short-time expansions for ATM option prices of exponential Lévy models. Mathematical Finance, 26(3), , J.E. Figueroa-López and S. Ólafsson, Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility. Finance & Stochastics 20(1), , J.E. Figueroa-López and S. Ólafsson, Short-time asymptotics for the implied volatility skew under a stochastic volatility model with Lévy jumps. Finance & Stochastics 20(4), , J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 38 / 39

71 Further Reading II J.E. Figueroa-López & J. Nisen. Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models. Stochastic Processes and their Applications 123(7), , J.E. Figueroa-López & C. Mancini. Optimum thresholding using mean and conditional mean square error. Available at J.E. Figueroa-López, C. Li, & J. Nisen. Optimal iterative threshold-kernel estimation of jump diffusion processes. In preparation, J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics CCF Seminar 39 / 39

Optimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error

Optimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error Optimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error José E. Figueroa-López Department of Mathematics Washington University in St. Louis Spring Central Sectional Meeting

More information

Asymptotic Methods in Financial Mathematics

Asymptotic Methods in Financial Mathematics Asymptotic Methods in Financial Mathematics José E. Figueroa-López 1 1 Department of Mathematics Washington University in St. Louis Statistics Seminar Washington University in St. Louis February 17, 2017

More information

Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models

Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models José E. Figueroa-López 1 1 Department of Statistics Purdue University University of Missouri-Kansas City Department of Mathematics

More information

Short-Time Asymptotic Methods In Financial Mathematics

Short-Time Asymptotic Methods In Financial Mathematics Short-Time Asymptotic Methods In Financial Mathematics José E. Figueroa-López Department of Mathematics Washington University in St. Louis Department of Applied Mathematics, Illinois Institute of Technology

More information

Short-Time Asymptotic Methods In Financial Mathematics

Short-Time Asymptotic Methods In Financial Mathematics Short-Time Asymptotic Methods In Financial Mathematics José E. Figueroa-López Department of Mathematics and Statistics Washington University in St. Louis School Of Mathematics, UMN March 14, 2019 Based

More information

Optimally Thresholded Realized Power Variations for Stochastic Volatility Models with Jumps

Optimally Thresholded Realized Power Variations for Stochastic Volatility Models with Jumps Optimally Thresholded Realized Power Variations for Stochastic Volatility Models with Jumps José E. Figueroa-López 1 1 Department of Mathematics Washington University ISI 2015: 60th World Statistics Conference

More information

A Simulation Study of Bipower and Thresholded Realized Variations for High-Frequency Data

A Simulation Study of Bipower and Thresholded Realized Variations for High-Frequency Data Washington University in St. Louis Washington University Open Scholarship Arts & Sciences Electronic Theses and Dissertations Arts & Sciences Spring 5-18-2018 A Simulation Study of Bipower and Thresholded

More information

Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models

Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models José E. Figueroa-López 1 1 Department of Statistics Purdue University High Dimensional Probability VII Institut d Études Scientifiques

More information

Volatility. Roberto Renò. 2 March 2010 / Scuola Normale Superiore. Dipartimento di Economia Politica Università di Siena

Volatility. Roberto Renò. 2 March 2010 / Scuola Normale Superiore. Dipartimento di Economia Politica Università di Siena Dipartimento di Economia Politica Università di Siena 2 March 2010 / Scuola Normale Superiore What is? The definition of volatility may vary wildly around the idea of the standard deviation of price movements

More information

Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models

Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models José E. Figueroa-López 1 1 Department of Statistics Purdue University Financial Statistics Stevanovich Center The University

More information

Limit 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 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 information

Short-time asymptotics for ATM option prices under tempered stable processes

Short-time asymptotics for ATM option prices under tempered stable processes Short-time asymptotics for ATM option prices under tempered stable processes José E. Figueroa-López 1 1 Department of Statistics Purdue University Probability Seminar Purdue University Oct. 30, 2012 Joint

More information

Economics 883: The Basic Diffusive Model, Jumps, Variance Measures. George Tauchen. Economics 883FS Spring 2015

Economics 883: The Basic Diffusive Model, Jumps, Variance Measures. George Tauchen. Economics 883FS Spring 2015 Economics 883: The Basic Diffusive Model, Jumps, Variance Measures George Tauchen Economics 883FS Spring 2015 Main Points 1. The Continuous Time Model, Theory and Simulation 2. Observed Data, Plotting

More information

Near-expiration behavior of implied volatility for exponential Lévy models

Near-expiration behavior of implied volatility for exponential Lévy models Near-expiration behavior of implied volatility for exponential Lévy models José E. Figueroa-López 1 1 Department of Statistics Purdue University Financial Mathematics Seminar The Stevanovich Center for

More information

Optimal Placement of a Small Order Under a Diffusive Limit Order Book (LOB) Model

Optimal Placement of a Small Order Under a Diffusive Limit Order Book (LOB) Model Optimal Placement of a Small Order Under a Diffusive Limit Order Book (LOB) Model José E. Figueroa-López Department of Mathematics Washington University in St. Louis INFORMS National Meeting Houston, TX

More information

Economics 883: The Basic Diffusive Model, Jumps, Variance Measures, and Noise Corrections. George Tauchen. Economics 883FS Spring 2014

Economics 883: The Basic Diffusive Model, Jumps, Variance Measures, and Noise Corrections. George Tauchen. Economics 883FS Spring 2014 Economics 883: The Basic Diffusive Model, Jumps, Variance Measures, and Noise Corrections George Tauchen Economics 883FS Spring 2014 Main Points 1. The Continuous Time Model, Theory and Simulation 2. Observed

More information

Ultra High Frequency Volatility Estimation with Market Microstructure Noise. Yacine Aït-Sahalia. Per A. Mykland. Lan Zhang

Ultra High Frequency Volatility Estimation with Market Microstructure Noise. Yacine Aït-Sahalia. Per A. Mykland. Lan Zhang Ultra High Frequency Volatility Estimation with Market Microstructure Noise Yacine Aït-Sahalia Princeton University Per A. Mykland The University of Chicago Lan Zhang Carnegie-Mellon University 1. Introduction

More information

Optimal Kernel Estimation of Spot Volatility of SDE

Optimal Kernel Estimation of Spot Volatility of SDE Optimal Kernel Estimation of Spot Volatility of SDE José E. Figueroa-López Department of Mathematics Washington University in St. Louis figueroa@math.wustl.edu (Joint work with Cheng Li from Purdue U.)

More information

Equity correlations implied by index options: estimation and model uncertainty analysis

Equity 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 information

I Preliminary Material 1

I Preliminary Material 1 Contents Preface Notation xvii xxiii I Preliminary Material 1 1 From Diffusions to Semimartingales 3 1.1 Diffusions.......................... 5 1.1.1 The Brownian Motion............... 5 1.1.2 Stochastic

More information

Modeling the extremes of temperature time series. Debbie J. Dupuis Department of Decision Sciences HEC Montréal

Modeling the extremes of temperature time series. Debbie J. Dupuis Department of Decision Sciences HEC Montréal Modeling the extremes of temperature time series Debbie J. Dupuis Department of Decision Sciences HEC Montréal Outline Fig. 1: S&P 500. Daily negative returns (losses), Realized Variance (RV) and Jump

More information

Estimation methods for Levy based models of asset prices

Estimation methods for Levy based models of asset prices Estimation methods for Levy based models of asset prices José Enrique Figueroa-López Financial mathematics seminar Department of Statistics and Applied Probability UCSB October, 26 Abstract Stock prices

More information

Asymptotic Theory for Renewal Based High-Frequency Volatility Estimation

Asymptotic Theory for Renewal Based High-Frequency Volatility Estimation Asymptotic Theory for Renewal Based High-Frequency Volatility Estimation Yifan Li 1,2 Ingmar Nolte 1 Sandra Nolte 1 1 Lancaster University 2 University of Manchester 4th Konstanz - Lancaster Workshop on

More information

A Consistent Pricing Model for Index Options and Volatility Derivatives

A 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 information

Valuation of Volatility Derivatives. Jim Gatheral Global Derivatives & Risk Management 2005 Paris May 24, 2005

Valuation 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 information

European option pricing under parameter uncertainty

European 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 information

Mgr. Jakub Petrásek 1. May 4, 2009

Mgr. Jakub Petrásek 1. May 4, 2009 Dissertation Report - First Steps Petrásek 1 2 1 Department of Probability and Mathematical Statistics, Charles University email:petrasek@karlin.mff.cuni.cz 2 RSJ Invest a.s., Department of Probability

More information

Lecture on advanced volatility models

Lecture on advanced volatility models FMS161/MASM18 Financial Statistics Stochastic Volatility (SV) Let r t be a stochastic process. The log returns (observed) are given by (Taylor, 1982) r t = exp(v t /2)z t. The volatility V t is a hidden

More information

Estimation of High-Frequency Volatility: An Autoregressive Conditional Duration Approach

Estimation of High-Frequency Volatility: An Autoregressive Conditional Duration Approach Estimation of High-Frequency Volatility: An Autoregressive Conditional Duration Approach Yiu-Kuen Tse School of Economics, Singapore Management University Thomas Tao Yang Department of Economics, Boston

More information

Assicurazioni Generali: An Option Pricing Case with NAGARCH

Assicurazioni Generali: An Option Pricing Case with NAGARCH Assicurazioni Generali: An Option Pricing Case with NAGARCH Assicurazioni Generali: Business Snapshot Find our latest analyses and trade ideas on bsic.it Assicurazioni Generali SpA is an Italy-based insurance

More information

On modelling of electricity spot price

On 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 information

Parametric Inference and Dynamic State Recovery from Option Panels. Torben G. Andersen

Parametric Inference and Dynamic State Recovery from Option Panels. Torben G. Andersen Parametric Inference and Dynamic State Recovery from Option Panels Torben G. Andersen Joint work with Nicola Fusari and Viktor Todorov The Third International Conference High-Frequency Data Analysis in

More information

Economics 201FS: Variance Measures and Jump Testing

Economics 201FS: Variance Measures and Jump Testing 1/32 : Variance Measures and Jump Testing George Tauchen Duke University January 21 1. Introduction and Motivation 2/32 Stochastic volatility models account for most of the anomalies in financial price

More information

Testing for non-correlation between price and volatility jumps and ramifications

Testing for non-correlation between price and volatility jumps and ramifications Testing for non-correlation between price and volatility jumps and ramifications Claudia Klüppelberg Technische Universität München cklu@ma.tum.de www-m4.ma.tum.de Joint work with Jean Jacod, Gernot Müller,

More information

Applications of short-time asymptotics to the statistical estimation and option pricing of Lévy-driven models

Applications of short-time asymptotics to the statistical estimation and option pricing of Lévy-driven models Applications of short-time asymptotics to the statistical estimation and option pricing of Lévy-driven models José Enrique Figueroa-López 1 1 Department of Statistics Purdue University CIMAT and Universidad

More information

BROWNIAN MOTION Antonella Basso, Martina Nardon

BROWNIAN MOTION Antonella Basso, Martina Nardon BROWNIAN MOTION Antonella Basso, Martina Nardon basso@unive.it, mnardon@unive.it Department of Applied Mathematics University Ca Foscari Venice Brownian motion p. 1 Brownian motion Brownian motion plays

More information

An Introduction to Stochastic Calculus

An Introduction to Stochastic Calculus An Introduction to Stochastic Calculus Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 2-3 Haijun Li An Introduction to Stochastic Calculus Week 2-3 1 / 24 Outline

More information

Jumps in Equilibrium Prices. and Market Microstructure Noise

Jumps in Equilibrium Prices. and Market Microstructure Noise Jumps in Equilibrium Prices and Market Microstructure Noise Suzanne S. Lee and Per A. Mykland Abstract Asset prices we observe in the financial markets combine two unobservable components: equilibrium

More information

Optimal 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 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 information

Pricing and hedging with rough-heston models

Pricing and hedging with rough-heston models Pricing and hedging with rough-heston models Omar El Euch, Mathieu Rosenbaum Ecole Polytechnique 1 January 216 El Euch, Rosenbaum Pricing and hedging with rough-heston models 1 Table of contents Introduction

More information

Exact Sampling of Jump-Diffusion Processes

Exact Sampling of Jump-Diffusion Processes 1 Exact Sampling of Jump-Diffusion Processes and Dmitry Smelov Management Science & Engineering Stanford University Exact Sampling of Jump-Diffusion Processes 2 Jump-Diffusion Processes Ubiquitous in finance

More information

The Black-Scholes Model

The 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 information

Optimal Kernel Estimation of Spot Volatility

Optimal Kernel Estimation of Spot Volatility Optimal Kernel Estimation of Spot Volatility José E. Figueroa-López Department of Mathematics and Statistics Washington University in St. Louis figueroa@math.wustl.edu Joint work with Cheng Li from Purdue

More information

Basic Concepts in Mathematical Finance

Basic Concepts in Mathematical Finance Chapter 1 Basic Concepts in Mathematical Finance In this chapter, we give an overview of basic concepts in mathematical finance theory, and then explain those concepts in very simple cases, namely in the

More information

Lecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.

Lecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r. Lecture 7 Overture to continuous models Before rigorously deriving the acclaimed Black-Scholes pricing formula for the value of a European option, we developed a substantial body of material, in continuous

More information

Using Lévy Processes to Model Return Innovations

Using Lévy Processes to Model Return Innovations Using Lévy Processes to Model Return Innovations Liuren Wu Zicklin School of Business, Baruch College Option Pricing Liuren Wu (Baruch) Lévy Processes Option Pricing 1 / 32 Outline 1 Lévy processes 2 Lévy

More information

Option pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard

Option pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard Option pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard Indifference pricing and the minimal entropy martingale measure Fred Espen Benth Centre of Mathematics for Applications

More information

Control. Econometric Day Mgr. Jakub Petrásek 1. Supervisor: RSJ Invest a.s.,

Control. Econometric Day Mgr. Jakub Petrásek 1. Supervisor: RSJ Invest a.s., and and Econometric Day 2009 Petrásek 1 2 1 Department of Probability and Mathematical Statistics, Charles University, RSJ Invest a.s., email:petrasek@karlin.mff.cuni.cz 2 Department of Probability and

More information

Option Pricing Modeling Overview

Option Pricing Modeling Overview Option Pricing Modeling Overview Liuren Wu Zicklin School of Business, Baruch College Options Markets Liuren Wu (Baruch) Stochastic time changes Options Markets 1 / 11 What is the purpose of building a

More information

Chapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets

Chapter 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 information

Volume and volatility in European electricity markets

Volume and volatility in European electricity markets Volume and volatility in European electricity markets Roberto Renò reno@unisi.it Dipartimento di Economia Politica, Università di Siena Commodities 2007 - Birkbeck, 17-19 January 2007 p. 1/29 Joint work

More information

Calculating Implied Volatility

Calculating Implied Volatility Statistical Laboratory University of Cambridge University of Cambridge Mathematics and Big Data Showcase 20 April 2016 How much is an option worth? A call option is the right, but not the obligation, to

More information

The Black-Scholes Model

The 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 information

Advanced 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 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 information

Financial Econometrics and Volatility Models Estimating Realized Variance

Financial Econometrics and Volatility Models Estimating Realized Variance Financial Econometrics and Volatility Models Estimating Realized Variance Eric Zivot June 2, 2010 Outline Volatility Signature Plots Realized Variance and Market Microstructure Noise Unbiased Estimation

More information

Rough Heston models: Pricing, hedging and microstructural foundations

Rough Heston models: Pricing, hedging and microstructural foundations Rough Heston models: Pricing, hedging and microstructural foundations Omar El Euch 1, Jim Gatheral 2 and Mathieu Rosenbaum 1 1 École Polytechnique, 2 City University of New York 7 November 2017 O. El Euch,

More information

Pricing of a European Call Option Under a Local Volatility Interbank Offered Rate Model

Pricing 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 information

Parametric Inference and Dynamic State Recovery from Option Panels. Nicola Fusari

Parametric Inference and Dynamic State Recovery from Option Panels. Nicola Fusari Parametric Inference and Dynamic State Recovery from Option Panels Nicola Fusari Joint work with Torben G. Andersen and Viktor Todorov July 2012 Motivation Under realistic assumptions derivatives are nonredundant

More information

Tangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.

Tangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford. Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey

More information

Large Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models

Large Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models Large Deviations and Stochastic Volatility with Jumps: TU Berlin with A. Jaquier and A. Mijatović (Imperial College London) SIAM conference on Financial Mathematics, Minneapolis, MN July 10, 2012 Implied

More information

Lecture 11: Stochastic Volatility Models Cont.

Lecture 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 information

Optimum thresholding using mean and conditional mean square error

Optimum thresholding using mean and conditional mean square error Optimum tresolding using mean and conditional mean square error José E. Figueroa-López and Cecilia Mancini August 7, 7 Abstract We consider a univariate semimartingale model for te logaritm of an asset

More information

Smile in the low moments

Smile 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 information

IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS

IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS DENIS BELOMESTNY AND MARKUS REISS 1. Introduction The aim of this report is to describe more precisely how the spectral calibration method

More information

Efficient multipowers

Efficient multipowers Efficient multipowers Kolokolov, Aleksey; Reno, Roberto 2016 Link to publication Citation for published version (APA): Kolokolov, A., & Reno, R. (2016). Efficient multipowers. (Working Papers in Statistics;

More information

Remarks on rough Bergomi: asymptotics and calibration

Remarks on rough Bergomi: asymptotics and calibration Department of Mathematics, Imperial College London Advances in Financial Mathematics, Paris, January 2017 Based on joint works with C Martini, A Muguruza, M Pakkanen and H Stone January 11, 2017 Implied

More information

Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias

Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias José E. Figueroa-López 1 1 Department of Statistics Purdue University Computational Finance Seminar Purdue University

More information

IEOR E4703: Monte-Carlo Simulation

IEOR E4703: Monte-Carlo Simulation IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

Quadratic hedging in affine stochastic volatility models

Quadratic hedging in affine stochastic volatility models Quadratic hedging in affine stochastic volatility models Jan Kallsen TU München Pittsburgh, February 20, 2006 (based on joint work with F. Hubalek, L. Krawczyk, A. Pauwels) 1 Hedging problem S t = S 0

More information

On Using Shadow Prices in Portfolio optimization with Transaction Costs

On Using Shadow Prices in Portfolio optimization with Transaction Costs On Using Shadow Prices in Portfolio optimization with Transaction Costs Johannes Muhle-Karbe Universität Wien Joint work with Jan Kallsen Universidad de Murcia 12.03.2010 Outline The Merton problem The

More information

Pricing Dynamic Guaranteed Funds Under a Double Exponential. Jump Diffusion Process. Chuang-Chang Chang, Ya-Hui Lien and Min-Hung Tsay

Pricing Dynamic Guaranteed Funds Under a Double Exponential. Jump Diffusion Process. Chuang-Chang Chang, Ya-Hui Lien and Min-Hung Tsay Pricing Dynamic Guaranteed Funds Under a Double Exponential Jump Diffusion Process Chuang-Chang Chang, Ya-Hui Lien and Min-Hung Tsay ABSTRACT This paper complements the extant literature to evaluate the

More information

Weierstrass Institute for Applied Analysis and Stochastics Maximum likelihood estimation for jump diffusions

Weierstrass Institute for Applied Analysis and Stochastics Maximum likelihood estimation for jump diffusions Weierstrass Institute for Applied Analysis and Stochastics Maximum likelihood estimation for jump diffusions Hilmar Mai Mohrenstrasse 39 1117 Berlin Germany Tel. +49 3 2372 www.wias-berlin.de Haindorf

More information

Hedging 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 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 information

Exploring 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 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 information

Université de Montréal. Rapport de recherche. Empirical Analysis of Jumps Contribution to Volatility Forecasting Using High Frequency Data

Université de Montréal. Rapport de recherche. Empirical Analysis of Jumps Contribution to Volatility Forecasting Using High Frequency Data Université de Montréal Rapport de recherche Empirical Analysis of Jumps Contribution to Volatility Forecasting Using High Frequency Data Rédigé par : Imhof, Adolfo Dirigé par : Kalnina, Ilze Département

More information

1 The continuous time limit

1 The continuous time limit Derivative Securities, Courant Institute, Fall 2008 http://www.math.nyu.edu/faculty/goodman/teaching/derivsec08/index.html Jonathan Goodman and Keith Lewis Supplementary notes and comments, Section 3 1

More information

Modeling the dependence between a Poisson process and a continuous semimartingale

Modeling the dependence between a Poisson process and a continuous semimartingale 1 / 28 Modeling the dependence between a Poisson process and a continuous semimartingale Application to electricity spot prices and wind production modeling Thomas Deschatre 1,2 1 CEREMADE, University

More information

Universität Regensburg Mathematik

Universität Regensburg Mathematik Universität Regensburg Mathematik Modeling financial markets with extreme risk Tobias Kusche Preprint Nr. 04/2008 Modeling financial markets with extreme risk Dr. Tobias Kusche 11. January 2008 1 Introduction

More information

Lecture 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 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 information

Leverage Effect, Volatility Feedback, and Self-Exciting Market Disruptions 11/4/ / 24

Leverage 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 information

Statistical methods for financial models driven by Lévy processes

Statistical methods for financial models driven by Lévy processes Statistical methods for financial models driven by Lévy processes José Enrique Figueroa-López Department of Statistics, Purdue University PASI Centro de Investigación en Matemátics (CIMAT) Guanajuato,

More information

Implied Lévy Volatility

Implied Lévy Volatility Joint work with José Manuel Corcuera, Peter Leoni and Wim Schoutens July 15, 2009 - Eurandom 1 2 The Black-Scholes model The Lévy models 3 4 5 6 7 Delta Hedging at versus at Implied Black-Scholes Volatility

More information

Option P&L Attribution and Pricing

Option P&L Attribution and Pricing Option P&L Attribution and Pricing Liuren Wu joint with Peter Carr Baruch College March 23, 2018 Stony Brook University Carr and Wu (NYU & Baruch) P&L Attribution and Option Pricing March 23, 2018 1 /

More information

Continuous Processes. Brownian motion Stochastic calculus Ito calculus

Continuous Processes. Brownian motion Stochastic calculus Ito calculus Continuous Processes Brownian motion Stochastic calculus Ito calculus Continuous Processes The binomial models are the building block for our realistic models. Three small-scale principles in continuous

More information

Asset Pricing Models with Underlying Time-varying Lévy Processes

Asset Pricing Models with Underlying Time-varying Lévy Processes Asset Pricing Models with Underlying Time-varying Lévy Processes Stochastics & Computational Finance 2015 Xuecan CUI Jang SCHILTZ University of Luxembourg July 9, 2015 Xuecan CUI, Jang SCHILTZ University

More information

Extended Libor Models and Their Calibration

Extended 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 information

Statistical Analysis of Short-time Option Prices Based on a Levy Model

Statistical Analysis of Short-time Option Prices Based on a Levy Model Washington University in St. Louis Washington University Open Scholarship Arts & Sciences Electronic Theses and Dissertations Arts & Sciences Spring 5-2018 Statistical Analysis of Short-time Option Prices

More information

Rough volatility models: When population processes become a new tool for trading and risk management

Rough volatility models: When population processes become a new tool for trading and risk management Rough volatility models: When population processes become a new tool for trading and risk management Omar El Euch and Mathieu Rosenbaum École Polytechnique 4 October 2017 Omar El Euch and Mathieu Rosenbaum

More information

Numerical valuation for option pricing under jump-diffusion models by finite differences

Numerical valuation for option pricing under jump-diffusion models by finite differences Numerical valuation for option pricing under jump-diffusion models by finite differences YongHoon Kwon Younhee Lee Department of Mathematics Pohang University of Science and Technology June 23, 2010 Table

More information

Stochastic Dynamical Systems and SDE s. An Informal Introduction

Stochastic Dynamical Systems and SDE s. An Informal Introduction Stochastic Dynamical Systems and SDE s An Informal Introduction Olav Kallenberg Graduate Student Seminar, April 18, 2012 1 / 33 2 / 33 Simple recursion: Deterministic system, discrete time x n+1 = f (x

More information

Numerical 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 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 information

Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing

Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014

More information

Unified Credit-Equity Modeling

Unified Credit-Equity Modeling Unified Credit-Equity Modeling Rafael Mendoza-Arriaga Based on joint research with: Vadim Linetsky and Peter Carr The University of Texas at Austin McCombs School of Business (IROM) Recent Advancements

More information

The Implied Volatility Index

The Implied Volatility Index The Implied Volatility Index Risk Management Institute National University of Singapore First version: October 6, 8, this version: October 8, 8 Introduction This document describes the formulation and

More information

Variance derivatives and estimating realised variance from high-frequency data. John Crosby

Variance derivatives and estimating realised variance from high-frequency data. John Crosby Variance derivatives and estimating realised variance from high-frequency data John Crosby UBS, London and Centre for Economic and Financial Studies, Department of Economics, Glasgow University Presentation

More information

GARCH Options in Incomplete Markets

GARCH Options in Incomplete Markets GARCH Options in Incomplete Markets Giovanni Barone-Adesi a, Robert Engle b and Loriano Mancini a a Institute of Finance, University of Lugano, Switzerland b Dept. of Finance, Leonard Stern School of Business,

More information

Simple Robust Hedging with Nearby Contracts

Simple 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 information

Credit Risk using Time Changed Brownian Motions

Credit Risk using Time Changed Brownian Motions Credit Risk using Time Changed Brownian Motions Tom Hurd Mathematics and Statistics McMaster University Joint work with Alexey Kuznetsov (New Brunswick) and Zhuowei Zhou (Mac) 2nd Princeton Credit Conference

More information

NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours

NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 Stochastic Financial Modelling Time allowed: 2 hours Candidates should attempt all questions. Marks for each question

More information

Hedging with Life and General Insurance Products

Hedging with Life and General Insurance Products Hedging with Life and General Insurance Products June 2016 2 Hedging with Life and General Insurance Products Jungmin Choi Department of Mathematics East Carolina University Abstract In this study, a hybrid

More information