Asymptotic Methods in Financial Mathematics

Size: px
Start display at page:

Download "Asymptotic Methods in Financial Mathematics"

Transcription

1 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 (Joint work with Cecilia Mancini) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 1 / 29

2 Outline 1 Motivational Example 2 The Main Estimators Multipower Variations and Truncated Realized Variations 3 Optimal Threshold Selection via Expected number of jump misclassifications via conditional Mean Square Error (cmse) 4 Conclusions J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 2 / 29

3 Motivational Example Outline 1 Motivational Example 2 The Main Estimators Multipower Variations and Truncated Realized Variations 3 Optimal Threshold Selection via Expected number of jump misclassifications via conditional Mean Square Error (cmse) 4 Conclusions J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 3 / 29

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

5 Motivational Example Model Simulation 1 Fix the time horizon T and the number of sampling points n. Set the time mesh and the sampling times h = h n = T n, t i = t i,n = ih n, i = 0,..., n. 2 Generate i.i.d. i.i.d. Z i N(0, 1), γ i N(µ jmp, σjmp 2 ), and I i.i.d. i Bernoulli(λh) 3 Iteratively generate the process: X ti = X ti 1 + ah + σ hz i + I i γ i J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 5 / 29

6 The Main Estimators Outline 1 Motivational Example 2 The Main Estimators Multipower Variations and Truncated Realized Variations 3 Optimal Threshold Selection via Expected number of jump misclassifications via conditional Mean Square Error (cmse) 4 Conclusions J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 6 / 29

7 The Main Estimators 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 QV (X) n n σ 2 + j:τ j T ζ2 j J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 7 / 29

8 The Main Estimators Multipower Variations and Truncated Realized Variations Two main classes of estimators Precursor. Realized Quadratic Variation: QV (X) n := 1 T QV (X) n n 1 ( Xti+1 X ) 2 t i i=0 n σ 2 + j:τ j T ζ2 j 1 Multipower Realized Variations (Barndorff-Nielsen and Shephard): n 1 BPV (X) n := Xti+1 X ti Xti+2 X ti+1, i=0 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 7 / 29

9 The Main Estimators Multipower Variations and Truncated Realized Variations Two main classes of estimators Precursor. Realized Quadratic Variation: QV (X) n := 1 T QV (X) n n 1 ( Xti+1 X ) 2 t i i=0 n σ 2 + j:τ j T ζ2 j 1 Multipower Realized Variations (Barndorff-Nielsen and Shephard): n k MPV [X] (r 1,...,r k ) n := Xti+1 X ti r 1... X ti+k X ti+k 1 r k. i=0 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 7 / 29

10 The Main Estimators Multipower Variations and Truncated Realized Variations Two main classes of estimators Precursor. Realized Quadratic Variation: QV (X) n := 1 T QV (X) n n 1 ( Xti+1 X ) 2 t i i=0 n σ 2 + j:τ j T ζ2 j 1 Multipower Realized Variations (Barndorff-Nielsen and Shephard): n k MPV [X] (r 1,...,r k ) n := Xti+1 X ti r 1... X ti+k X ti+k 1 r k. 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 ε }, (ε [0, )). J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 7 / 29

11 The Main Estimators 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) Merton Model with σ = 0.4, σ jmp = 3(h), µ jmp = 0, λ = 200; (right) TRV performance wrt the truncation level J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 8 / 29

12 The Main Estimators 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) Merton Model 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 WUSTL Statistics Seminar 9 / 29

13 The Main Estimators 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. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 10 / 29

14 The Main Estimators 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. Power Threshold (Mancini (2003)) ε Prw α,ω := α h ω, for α > 0 and ω (0, 1/2). J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 10 / 29

15 The Main Estimators 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. Power Threshold (Mancini (2003)) ε Prw α,ω := α h ω, for α > 0 and ω (0, 1/2). Bonferroni Threshold (Bollerslev et al. (2007) and Gegler & Stadtmüller (2010)) ( B σ,c BF := σh1/2 Φ 1 1 C h ), for C > 0 and σ > 0. 2 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 10 / 29

16 Outline 1 Motivational Example 2 The Main Estimators Multipower Variations and Truncated Realized Variations 3 Optimal Threshold Selection via Expected number of jump misclassifications via conditional Mean Square Error (cmse) 4 Conclusions J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 11 / 29

17 Classical Approach 1 Fix a suitable and sensible metric for estimation error; typically, the MSE: MSE(ε) = E [ ( ) ] 1 2 T TRV n(ε) σ 2 2 Show the existence of the optimal threshold ε n minimizing the error function; 3 Analyze 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 calibration of ε by estimating those parameters (if possible). J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 12 / 29

18 via Expected number of jump misclassifications Expected number of jump misclassifications 1 Notation: N t = # of jumps by time t n i X = X t i X ti 1 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 WUSTL Statistics Seminar 13 / 29

19 via Expected number of jump misclassifications Expected number of jump misclassifications 1 Notation: N t = # of jumps by time t n i X = X t i X ti 1 n i N = N t i N ti 1 = # of jumps during (t i 1, t i ] 2 Estimation Error: (F-L & Nisen, 2013) n ) Loss n (ε) := E (1 [ ni X >ε, ni N=0] + 1 [ ni X ε, ni N 0]. 3 Motivation i=1 The underlying principle is that the estimation error is directly link to the ability of the threshold to detect jumps. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 13 / 29

20 via Expected number of jump misclassifications Existence and Infill Asymptotic Characterization Theorem (FL & Nisen (SPA, 2013)) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 14 / 29

21 via Expected number of jump misclassifications Existence and Infill Asymptotic Characterization Theorem (FL & Nisen (SPA, 2013)) 1 For n large enough, the loss function Loss n (ε) is convex and, moreover, possesses a unique global minimum ε n. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 14 / 29

22 via Expected number of jump misclassifications Existence and Infill Asymptotic Characterization Theorem (FL & Nisen (SPA, 2013)) 1 For n large enough, the loss function Loss n (ε) is convex and, moreover, possesses a unique global minimum ε n. 2 As n, the optimal threshold sequence (ε n) n is such that ( ) ( ) 1 log 2πσλC(fζ ) σh 1/2 ε n = 3σ 2 n h n log + h.o.t., h n 3 log(1/hn ) 1 δ where C(f ζ ) = lim δ 0 2δ δ f ζ(x)dx > 0 and f ζ is the density of the jumps ζ i (h.o.t. refers to high order terms). J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 14 / 29

23 via Expected number of jump misclassifications Remarks 1 The threshold sequences ε 1 n := ( 1 3σ 2 h n log h n ), ε 2 n := ε 1 n ( ) log 2πσλC(fζ ) 3 log(1/hn ) σh 1/2 n are the first and second-order approximations for B n, and the biases of their corresponding TRV estimators attain the optimal rate of O(h n ) as n., J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 15 / 29

24 via Expected number of jump misclassifications Remarks 1 The threshold sequences ε 1 n := ( 1 3σ 2 h n log h n ), ε 2 n := ε 1 n ( ) log 2πσλC(fζ ) 3 log(1/hn ) σh 1/2 n are the first and second-order approximations for B n, and the biases of their corresponding TRV estimators attain the optimal rate of O(h n ) as n. 3 They both provide blueprints" for devising good threshold sequences!, J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 15 / 29

25 via Expected number of jump misclassifications 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 WUSTL Statistics Seminar 16 / 29

26 via Expected number of jump misclassifications 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 WUSTL Statistics Seminar 16 / 29

27 via Expected number of jump misclassifications 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 WUSTL Statistics Seminar 16 / 29

28 via Expected number of jump misclassifications 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 WUSTL Statistics Seminar 16 / 29

29 via Expected number of jump misclassifications 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 WUSTL Statistics Seminar 17 / 29

30 via Expected number of jump misclassifications 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 WUSTL Statistics Seminar 18 / 29

31 via conditional Mean Square Error (cmse) Setting 1 Framework: t X t = X 0 + σ s dw s + J t 0 where J is a pure-jump process. 2 Key Assumptions: σ t > 0 for all t, and σ and J are independent of W. 3 Notation: m i = n i J = J t i J ti 1, σ 2 i = ti t i 1 σ 2 sds 4 Key Observation: Conditioning on σ and J, the increments n i X are i.i.d. normal with mean m i and variance σ 2 i J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 19 / 29

32 via conditional Mean Square Error (cmse) Key Relationships Let b i (ε) := E [ ( i X) 2 1 { i X ε} σ, J ] J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 20 / 29

33 via conditional Mean Square Error (cmse) Key Relationships Let b i (ε) := E [ ( i X) 2 1 { i X ε} σ, J ] 1 As it turns out, b i (ε) = (e (ε m i ) 2 2σ i 2 (ε + m i ) + e (ε+m i ) 2 2σ i 2 σi (ε m i )) 2π + m2 i + σi 2 π ( ε mi 2σi 0 e t 2 ε+m i 2σi ) dt + e t 2 dt 0 J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 20 / 29

34 via conditional Mean Square Error (cmse) Key Relationships Let b i (ε) := E [ ( i X) 2 1 { i X ε} σ, J ] 1 As it turns out, b i (ε) = (e (ε m i ) 2 2σ i 2 (ε + m i ) + e (ε+m i ) 2 2σ i 2 σi (ε m i )) 2π + m2 i + σi 2 ( ε mi 2σi e t 2 ε+m i 2σi ) dt + e t 2 dt π Furthermore, db i (ε) dε = ε 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 WUSTL Statistics Seminar 20 / 29

35 via conditional Mean Square Error (cmse) Conditional Mean Square Error (MSE c ) Theorem (F-L & Mancini (2017)) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 21 / 29

36 via conditional Mean Square Error (cmse) Conditional Mean Square Error (MSE c ) Theorem (F-L & Mancini (2017)) 1 Let IV := T 0 σ2 sds be the integrated variance of X and let MSE c (ε) := E [(TRV n (ε) IV ) 2 ] σ, J J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 21 / 29

37 via conditional Mean Square Error (cmse) Conditional Mean Square Error (MSE c ) Theorem (F-L & Mancini (2017)) 1 Let IV := T 0 σ2 sds be the integrated variance of X and let MSE c (ε) := E [(TRV n (ε) IV ) 2 ] σ, J 2 Then, MSE c (ε) is differential in (0, ) and, furthermore, d MSE c (ε) = ε 2 G(ε), dε where G(ε) := ( i a i(ε) ε ) j i b j(ε) 2IV. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 21 / 29

38 via conditional Mean Square Error (cmse) Conditional Mean Square Error (MSE c ) Theorem (F-L & Mancini (2017)) 1 Let IV := T 0 σ2 sds be the integrated variance of X and let MSE c (ε) := E [(TRV n (ε) IV ) 2 ] σ, J 2 Then, MSE c (ε) is differential in (0, ) and, furthermore, d MSE c (ε) = ε 2 G(ε), dε where G(ε) := ( i a i(ε) ε ) j i b j(ε) 2IV. 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 WUSTL Statistics Seminar 21 / 29

39 via conditional Mean Square Error (cmse) Asymptotics: FA process with constant variance Theorem (F-L & Mancini (2017)) J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 22 / 29

40 via conditional Mean Square Error (cmse) 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 WUSTL Statistics Seminar 22 / 29

41 via conditional Mean Square Error (cmse) 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 + 2 Then, as n, the optimal threshold ε n is such that ( ) 1 2σ 2 h n log ε n j=1 ζ j h n J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 22 / 29

42 via conditional Mean Square Error (cmse) 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 WUSTL Statistics Seminar 23 / 29

43 via conditional Mean Square Error (cmse) 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 WUSTL Statistics Seminar 23 / 29

44 via conditional Mean Square Error (cmse) 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 WUSTL Statistics Seminar 23 / 29

45 via conditional Mean Square Error (cmse) 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 WUSTL Statistics Seminar 23 / 29

46 via conditional Mean Square Error (cmse) 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 WUSTL Statistics Seminar 24 / 29

47 via conditional Mean Square Error (cmse) 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. Loss is the number of jump misclassifications and N is the number of iterations. bar is used to denote average. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 25 / 29

48 Conclusions Outline 1 Motivational Example 2 The Main Estimators Multipower Variations and Truncated Realized Variations 3 Optimal Threshold Selection via Expected number of jump misclassifications via conditional Mean Square Error (cmse) 4 Conclusions J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 26 / 29

49 Conclusions Ongoing and Future Research 1 In principle, we can apply the proposed method for varying volatility t σ t by localization; i.e., applying to periods where σ is approximately constant. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 27 / 29

50 Conclusions Ongoing and Future Research 1 In principle, we can apply the proposed method for varying volatility t σ t by localization; i.e., applying 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 WUSTL Statistics Seminar 27 / 29

51 Conclusions Ongoing and Future Research 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 WUSTL Statistics Seminar 28 / 29

52 Conclusions Ongoing and Future Research 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 which is a solution of the equation ε 2 + 2(n 1)E[b 1 (ε)] 2nh n σ 2 = 0. J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 28 / 29

53 Conclusions Ongoing and Future Research 3 Furthermore, J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 29 / 29

54 Conclusions Ongoing and Future Research 3 Furthermore, As expected, in the finite jump activity case, it turns out ( ) 1 ε n 2σ 2 h n log h n J.E. Figueroa-López (WUSTL) Asymptotics in Financial Mathematics WUSTL Statistics Seminar 29 / 29

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

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

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 Probability and Mathematical Finance Seminar Department of Mathematical

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

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

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

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

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

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

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

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

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. 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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Two hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER

Two hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS Answer any FOUR of the SIX questions.

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

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

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

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

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

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

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

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

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

Stochastic calculus Introduction I. Stochastic Finance. C. Azizieh VUB 1/91. C. Azizieh VUB Stochastic Finance

Stochastic calculus Introduction I. Stochastic Finance. C. Azizieh VUB 1/91. C. Azizieh VUB Stochastic Finance Stochastic Finance C. Azizieh VUB C. Azizieh VUB Stochastic Finance 1/91 Agenda of the course Stochastic calculus : introduction Black-Scholes model Interest rates models C. Azizieh VUB Stochastic Finance

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

Regression estimation in continuous time with a view towards pricing Bermudan options

Regression estimation in continuous time with a view towards pricing Bermudan options with a view towards pricing Bermudan options Tagung des SFB 649 Ökonomisches Risiko in Motzen 04.-06.06.2009 Financial engineering in times of financial crisis Derivate... süßes Gift für die Spekulanten

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

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

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

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

Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models

Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models David Prager 1 1 Associate Professor of Mathematics Anderson University (SC) Based on joint work with Professor Qing Zhang,

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

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

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

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

Local Volatility Dynamic Models

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

Homework Assignments

Homework Assignments Homework Assignments Week 1 (p. 57) #4.1, 4., 4.3 Week (pp 58 6) #4.5, 4.6, 4.8(a), 4.13, 4.0, 4.6(b), 4.8, 4.31, 4.34 Week 3 (pp 15 19) #1.9, 1.1, 1.13, 1.15, 1.18 (pp 9 31) #.,.6,.9 Week 4 (pp 36 37)

More information

induced by the Solvency II project

induced by the Solvency II project Asset Les normes allocation IFRS : new en constraints assurance induced by the Solvency II project 36 th International ASTIN Colloquium Zürich September 005 Frédéric PLANCHET Pierre THÉROND ISFA Université

More information

Optimal Dividend Policy of A Large Insurance Company with Solvency Constraints. Zongxia Liang

Optimal Dividend Policy of A Large Insurance Company with Solvency Constraints. Zongxia Liang Optimal Dividend Policy of A Large Insurance Company with Solvency Constraints Zongxia Liang Department of Mathematical Sciences Tsinghua University, Beijing 100084, China zliang@math.tsinghua.edu.cn Joint

More information

Using MCMC and particle filters to forecast stochastic volatility and jumps in financial time series

Using MCMC and particle filters to forecast stochastic volatility and jumps in financial time series Using MCMC and particle filters to forecast stochastic volatility and jumps in financial time series Ing. Milan Fičura DYME (Dynamical Methods in Economics) University of Economics, Prague 15.6.2016 Outline

More information

Realized Measures. Eduardo Rossi University of Pavia. November Rossi Realized Measures University of Pavia / 64

Realized Measures. Eduardo Rossi University of Pavia. November Rossi Realized Measures University of Pavia / 64 Realized Measures Eduardo Rossi University of Pavia November 2012 Rossi Realized Measures University of Pavia - 2012 1 / 64 Outline 1 Introduction 2 RV Asymptotics of RV Jumps and Bipower Variation 3 Realized

More information

Financial Econometrics Jeffrey R. Russell. Midterm 2014 Suggested Solutions. TA: B. B. Deng

Financial Econometrics Jeffrey R. Russell. Midterm 2014 Suggested Solutions. TA: B. B. Deng Financial Econometrics Jeffrey R. Russell Midterm 2014 Suggested Solutions TA: B. B. Deng Unless otherwise stated, e t is iid N(0,s 2 ) 1. (12 points) Consider the three series y1, y2, y3, and y4. Match

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

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

Risk Measurement in Credit Portfolio Models

Risk Measurement in Credit Portfolio Models 9 th DGVFM Scientific Day 30 April 2010 1 Risk Measurement in Credit Portfolio Models 9 th DGVFM Scientific Day 30 April 2010 9 th DGVFM Scientific Day 30 April 2010 2 Quantitative Risk Management Profit

More information

Neil Shephard Oxford-Man Institute of Quantitative Finance, University of Oxford

Neil Shephard Oxford-Man Institute of Quantitative Finance, University of Oxford Measuring the impact of jumps on multivariate price processes using multipower variation Neil Shephard Oxford-Man Institute of Quantitative Finance, University of Oxford 1 1 Introduction Review the econometrics

More information

Asymptotic methods in risk management. Advances in Financial Mathematics

Asymptotic methods in risk management. Advances in Financial Mathematics Asymptotic methods in risk management Peter Tankov Based on joint work with A. Gulisashvili Advances in Financial Mathematics Paris, January 7 10, 2014 Peter Tankov (Université Paris Diderot) Asymptotic

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

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

Sparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models

Sparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models Sparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models Norbert Hilber Seminar of Applied Mathematics ETH Zürich Workshop on Financial Modeling with Jump Processes p. 1/18 Outline

More information

2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises

2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises 96 ChapterVI. Variance Reduction Methods stochastic volatility ISExSoren5.9 Example.5 (compound poisson processes) Let X(t) = Y + + Y N(t) where {N(t)},Y, Y,... are independent, {N(t)} is Poisson(λ) with

More information

Introduction to Sequential Monte Carlo Methods

Introduction to Sequential Monte Carlo Methods Introduction to Sequential Monte Carlo Methods Arnaud Doucet NCSU, October 2008 Arnaud Doucet () Introduction to SMC NCSU, October 2008 1 / 36 Preliminary Remarks Sequential Monte Carlo (SMC) are a set

More information

Risk Neutral Valuation

Risk Neutral Valuation copyright 2012 Christian Fries 1 / 51 Risk Neutral Valuation Christian Fries Version 2.2 http://www.christian-fries.de/finmath April 19-20, 2012 copyright 2012 Christian Fries 2 / 51 Outline Notation Differential

More information

Chapter 8: Sampling distributions of estimators Sections

Chapter 8: Sampling distributions of estimators Sections Chapter 8 continued Chapter 8: Sampling distributions of estimators Sections 8.1 Sampling distribution of a statistic 8.2 The Chi-square distributions 8.3 Joint Distribution of the sample mean and sample

More information

Value at Risk and Self Similarity

Value at Risk and Self Similarity Value at Risk and Self Similarity by Olaf Menkens School of Mathematical Sciences Dublin City University (DCU) St. Andrews, March 17 th, 2009 Value at Risk and Self Similarity 1 1 Introduction The concept

More information

Multi-period mean variance asset allocation: Is it bad to win the lottery?

Multi-period mean variance asset allocation: Is it bad to win the lottery? Multi-period mean variance asset allocation: Is it bad to win the lottery? Peter Forsyth 1 D.M. Dang 1 1 Cheriton School of Computer Science University of Waterloo Guangzhou, July 28, 2014 1 / 29 The Basic

More information

Conditional Density Method in the Computation of the Delta with Application to Power Market

Conditional Density Method in the Computation of the Delta with Application to Power Market Conditional Density Method in the Computation of the Delta with Application to Power Market Asma Khedher Centre of Mathematics for Applications Department of Mathematics University of Oslo A joint work

More information

Polynomial processes in stochastic portofolio theory

Polynomial processes in stochastic portofolio theory Polynomial processes in stochastic portofolio theory Christa Cuchiero University of Vienna 9 th Bachelier World Congress July 15, 2016 Christa Cuchiero (University of Vienna) Polynomial processes in SPT

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

Strategies for Improving the Efficiency of Monte-Carlo Methods

Strategies for Improving the Efficiency of Monte-Carlo Methods Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful

More information

Final exam solutions

Final exam solutions EE365 Stochastic Control / MS&E251 Stochastic Decision Models Profs. S. Lall, S. Boyd June 5 6 or June 6 7, 2013 Final exam solutions This is a 24 hour take-home final. Please turn it in to one of the

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

Arbitrage of the first kind and filtration enlargements in semimartingale financial models. Beatrice Acciaio

Arbitrage of the first kind and filtration enlargements in semimartingale financial models. Beatrice Acciaio Arbitrage of the first kind and filtration enlargements in semimartingale financial models Beatrice Acciaio the London School of Economics and Political Science (based on a joint work with C. Fontana and

More information

Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty

Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty George Photiou Lincoln College University of Oxford A dissertation submitted in partial fulfilment for

More information

Importance sampling and Monte Carlo-based calibration for time-changed Lévy processes

Importance sampling and Monte Carlo-based calibration for time-changed Lévy processes Importance sampling and Monte Carlo-based calibration for time-changed Lévy processes Stefan Kassberger Thomas Liebmann BFS 2010 1 Motivation 2 Time-changed Lévy-models and Esscher transforms 3 Applications

More information

Recent Advances in Fractional Stochastic Volatility Models

Recent Advances in Fractional Stochastic Volatility Models Recent Advances in Fractional Stochastic Volatility Models Alexandra Chronopoulou Industrial & Enterprise Systems Engineering University of Illinois at Urbana-Champaign IPAM National Meeting of Women in

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

Posterior Inference. , where should we start? Consider the following computational procedure: 1. draw samples. 2. convert. 3. compute properties

Posterior Inference. , where should we start? Consider the following computational procedure: 1. draw samples. 2. convert. 3. compute properties Posterior Inference Example. Consider a binomial model where we have a posterior distribution for the probability term, θ. Suppose we want to make inferences about the log-odds γ = log ( θ 1 θ), where

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

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

STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL

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

Simulating Stochastic Differential Equations

Simulating Stochastic Differential Equations IEOR E4603: Monte-Carlo Simulation c 2017 by Martin Haugh Columbia University Simulating Stochastic Differential Equations In these lecture notes we discuss the simulation of stochastic differential equations

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

Convergence Analysis of Monte Carlo Calibration of Financial Market Models

Convergence Analysis of Monte Carlo Calibration of Financial Market Models Analysis of Monte Carlo Calibration of Financial Market Models Christoph Käbe Universität Trier Workshop on PDE Constrained Optimization of Certain and Uncertain Processes June 03, 2009 Monte Carlo Calibration

More information

Summary Sampling Techniques

Summary Sampling Techniques Summary Sampling Techniques MS&E 348 Prof. Gerd Infanger 2005/2006 Using Monte Carlo sampling for solving the problem Monte Carlo sampling works very well for estimating multiple integrals or multiple

More information

Structural Models of Credit Risk and Some Applications

Structural Models of Credit Risk and Some Applications Structural Models of Credit Risk and Some Applications Albert Cohen Actuarial Science Program Department of Mathematics Department of Statistics and Probability albert@math.msu.edu August 29, 2018 Outline

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

UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions.

UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. Random Variables 2 A random variable X is a numerical (integer, real, complex, vector etc.) summary of the outcome of the random experiment.

More information

A Multivariate Analysis of Intercompany Loss Triangles

A Multivariate Analysis of Intercompany Loss Triangles A Multivariate Analysis of Intercompany Loss Triangles Peng Shi School of Business University of Wisconsin-Madison ASTIN Colloquium May 21-24, 2013 Peng Shi (Wisconsin School of Business) Intercompany

More information

Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs.

Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs. Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs Andrea Cosso LPMA, Université Paris Diderot joint work with Francesco Russo ENSTA,

More information

STAT/MATH 395 PROBABILITY II

STAT/MATH 395 PROBABILITY II STAT/MATH 395 PROBABILITY II Distribution of Random Samples & Limit Theorems Néhémy Lim University of Washington Winter 2017 Outline Distribution of i.i.d. Samples Convergence of random variables The Laws

More information

Problem Set 3. Thomas Philippon. April 19, Human Wealth, Financial Wealth and Consumption

Problem Set 3. Thomas Philippon. April 19, Human Wealth, Financial Wealth and Consumption Problem Set 3 Thomas Philippon April 19, 2002 1 Human Wealth, Financial Wealth and Consumption The goal of the question is to derive the formulas on p13 of Topic 2. This is a partial equilibrium analysis

More information

Unit 5: Sampling Distributions of Statistics

Unit 5: Sampling Distributions of Statistics Unit 5: Sampling Distributions of Statistics Statistics 571: Statistical Methods Ramón V. León 6/12/2004 Unit 5 - Stat 571 - Ramon V. Leon 1 Definitions and Key Concepts A sample statistic used to estimate

More information

Unit 5: Sampling Distributions of Statistics

Unit 5: Sampling Distributions of Statistics Unit 5: Sampling Distributions of Statistics Statistics 571: Statistical Methods Ramón V. León 6/12/2004 Unit 5 - Stat 571 - Ramon V. Leon 1 Definitions and Key Concepts A sample statistic used to estimate

More information

Volatility Measurement

Volatility Measurement Volatility Measurement Eduardo Rossi University of Pavia December 2013 Rossi Volatility Measurement Financial Econometrics - 2012 1 / 53 Outline 1 Volatility definitions Continuous-Time No-Arbitrage Price

More information