Optimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error
|
|
- Amos Bradford
- 5 years ago
- Views:
Transcription
1 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 Indiana University, Bloomington, IN April 2, 2017 (Joint work with Cecilia Mancini from University of Florence) J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 1 / 35
2 Outline 1 Motivational Example 2 The Problem Multipower Variations and Truncated Realized Variations The problem of Parameter Tuning 3 Optimal Threshold Selection via Expected number of jump misclassifications via conditional Mean Square Error (cmse) 4 Ongoing and Future Research J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 2 / 35
3 Motivational Example Outline 1 Motivational Example 2 The Problem Multipower Variations and Truncated Realized Variations The problem of Parameter Tuning 3 Optimal Threshold Selection via Expected number of jump misclassifications via conditional Mean Square Error (cmse) 4 Ongoing and Future Research J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 3 / 35
4 Motivational Example Merton s 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 J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 4 / 35
5 Motivational Example Merton s 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 (σ = 0.4, λ = 200, σ jmp = 0.02, µ jmp = 0, t i t i 1 = 5 min, T = t n = 1/12) J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 4 / 35
6 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 as 2 Generate Z i i.i.d. N(0, 1), h = h n = T n, t i = t i,n = ih n, i = 0,..., n. γ i i.i.d. N(µ jmp, σ 2 jmp ), 3 Iteratively generate the process starting from X 0 : I i i.i.d. Bernoulli(λh) X ti = X ti 1 + ah + σ hz i + I i γ i J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 5 / 35
7 The Problem Outline 1 Motivational Example 2 The Problem Multipower Variations and Truncated Realized Variations The problem of Parameter Tuning 3 Optimal Threshold Selection via Expected number of jump misclassifications via conditional Mean Square Error (cmse) 4 Ongoing and Future Research J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 6 / 35
8 The Problem Multipower Variations and Truncated Realized Variations Two main classes of estimators Precursor. Realized Quadratic Variation: n 1 ( RQV [X] n := Xti+1 X ) 2 t i i=0 RQV [X] n n T σ 2 + j:τ j T ζ2 j J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 7 / 35
9 The Problem Multipower Variations and Truncated Realized Variations Two main classes of estimators Precursor. Realized Quadratic Variation: n 1 ( RQV [X] n := Xti+1 X ) 2 t i RQV [X] n i=0 n T σ 2 + j:τ j T ζ2 j 1 Bipower 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) Optimum Thresholding AMS Sectional Meeting 7 / 35
10 The Problem Multipower Variations and Truncated Realized Variations Two main classes of estimators Precursor. Realized Quadratic Variation: n 1 ( RQV [X] n := Xti+1 X ) 2 t i RQV [X] n i=0 n T σ 2 + j:τ j T ζ2 j 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 ε }, (ε [0, )). J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 7 / 35
11 The Problem The problem of Parameter Tuning Calibration or Tuning of the Estimator Problem: How do you choose the threshold parameter ε? J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 8 / 35
12 The Problem The problem of Parameter Tuning 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 performance wrt the truncation level J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 8 / 35
13 The Problem The problem of Parameter Tuning 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) Optimum Thresholding AMS Sectional Meeting 9 / 35
14 The Problem The problem of Parameter Tuning Truncation Selection 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) Optimum Thresholding AMS Sectional Meeting 10 / 35
15 Outline 1 Motivational Example 2 The Problem Multipower Variations and Truncated Realized Variations The problem of Parameter Tuning 3 Optimal Threshold Selection via Expected number of jump misclassifications via conditional Mean Square Error (cmse) 4 Ongoing and Future Research J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 11 / 35
16 Classical Approach 1 Fix a suitable and sensible metric of the 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 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) Optimum Thresholding AMS Sectional Meeting 12 / 35
17 via Expected number of jump misclassifications Via Expected number of jump misclassifications J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 13 / 35
18 via Expected number of jump misclassifications Via Expected number of jump misclassifications 1 Estimation Error: (F-L & Nisen, SPA 2013) n ) Loss n (ε) := E (1 [ ni X >ε, ni N=0] + 1 [ ni X ε, ni N 0]. 2 Notation: i=1 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) Optimum Thresholding AMS Sectional Meeting 13 / 35
19 via Expected number of jump misclassifications Via Expected number of jump misclassifications 1 Estimation Error: (F-L & Nisen, SPA 2013) n ) Loss n (ε) := E (1 [ ni X >ε, ni N=0] + 1 [ ni X ε, ni N 0]. 2 Notation: 3 Idea i=1 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 ] Intuitively, the estimation error should be directly linked to the ability of the threshold to detect jumps. J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 13 / 35
20 via Expected number of jump misclassifications Existence and Infill Asymptotic Characterization Theorem (FL & Nisen, SPA 2013) a J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 14 / 35
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 quasi-convex and, moreover, possesses a unique global minimum ε n. a J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 14 / 35
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 quasi-convex and, moreover, possesses 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., 1 δ where C(f ζ ) = lim δ 0 2δ δ f ζ(x)dx > 0 and f ζ is the density of the jumps ζ i. a a (throughout, h.o.t. refers to higher order terms ) h n J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 14 / 35
23 via Expected number of jump misclassifications Remarks 1 The threshold sequence ε 1 n := ( 1 3σ 2 h n log is the first-order approximations for ε n. h n ), J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 15 / 35
24 via Expected number of jump misclassifications Remarks 1 The threshold sequence ε 1 n := ( 1 3σ 2 h n log is the first-order approximations for ε n. 2 Why h log(1/h)? h n ), J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 15 / 35
25 via Expected number of jump misclassifications Remarks 1 The threshold sequence ε 1 n := ( 1 3σ 2 h n log is the first-order approximations for ε n. 2 Why h log(1/h)? This is proportional to modulus of continuity of the B.M.: lim sup h 0 1 2h log(1/h) h n ), sup W t W s = 1. s,t [0,1]: t s <h J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 15 / 35
26 via Expected number of jump misclassifications Remarks 1 The threshold sequence ε 1 n := ( 1 3σ 2 h n log is the first-order approximations for ε n. 2 Why h log(1/h)? This is proportional to modulus of continuity of the B.M.: lim sup h 0 1 2h log(1/h) h n ), sup W t W s = 1. s,t [0,1]: t s <h 3 Practically, ε 1 n provides us with a blueprint" for devising good threshold sequences! J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 15 / 35
27 via Expected number of jump misclassifications A Feasible Implementation based on ε 1 n (i) Get a rough" estimate of σ 2 via, e.g., the realized QV: ˆσ n,0 2 := 1 n X ti X ti 1 2 T i=1 J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 16 / 35
28 via Expected number of jump misclassifications A Feasible Implementation based on ε 1 n (i) Get a rough" estimate of σ 2 via, e.g., the 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) Optimum Thresholding AMS Sectional Meeting 16 / 35
29 via Expected number of jump misclassifications A Feasible Implementation based on ε 1 n (i) Get a rough" estimate of σ 2 via, e.g., the 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) Optimum Thresholding AMS Sectional Meeting 16 / 35
30 via Expected number of jump misclassifications A Feasible Implementation based on ε 1 n (i) Get a rough" estimate of σ 2 via, e.g., the 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): ˆσ 2 n,0 ˆε n,0 ˆσ2 n,1 ˆε n,1 ˆσ2 n,2 ˆσ2 n, ] J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 16 / 35
31 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) Optimum Thresholding AMS Sectional Meeting 17 / 35
32 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) Optimum Thresholding AMS Sectional Meeting 18 / 35
33 via conditional Mean Square Error (cmse) Via conditional Mean Square Error (cmse) J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 19 / 35
34 via conditional Mean Square Error (cmse) 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) Optimum Thresholding AMS Sectional Meeting 19 / 35
35 via conditional Mean Square Error (cmse) 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 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) Optimum Thresholding AMS Sectional Meeting 19 / 35
36 via conditional Mean Square Error (cmse) Key Relationships 1 Clearly, under the previous non-leverage assumption, i X σ, J independent N ( m i, σi 2 ), where m i := n i J = J t i J ti 1 and σ 2 i := t i t i 1 σ 2 s ds. J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 20 / 35
37 via conditional Mean Square Error (cmse) Key Relationships 1 Clearly, under the previous non-leverage assumption, i X σ, J independent N ( m i, σi 2 ), where m i := n i J = J t i J ti 1 and σi 2 := t i t i 1 σs 2 ds. 2 Then, b i (ε) := E [ ( i X) 2 ] 1 { i X ε} σ, J J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 20 / 35
38 via conditional Mean Square Error (cmse) Key Relationships 1 Clearly, under the previous non-leverage assumption, i X σ, J independent N ( m i, σi 2 ), where m i := n i J = J t i J ti 1 and σi 2 := t i t i 1 σs 2 ds. 2 Then, b i (ε) := E [ ( i X) 2 ] 1 { i X ε} σ, J is such that b i (ε) = σ ( i 2π 2σ i 2 (ε + m i ) + e (ε+m i ) 2 ) 2σ i 2 (ε m i ) e (ε mi ) 2 ε+m i + (mi 2 + σi 2 σ i ) ε m i σ i e x 2 /2 dx 2π J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 20 / 35
39 via conditional Mean Square Error (cmse) Key Relationships 1 Clearly, under the previous non-leverage assumption, i X σ, J independent N ( m i, σi 2 ), where m i := n i J = J t i J ti 1 and σi 2 := t i t i 1 σs 2 ds. 2 Then, b i (ε) := E [ ( i X) 2 ] 1 { i X ε} σ, J is such that b i (ε) = σ ( i 2π 2σ i 2 (ε + m i ) + e (ε+m i ) 2 ) 2σ i 2 (ε m i ) e (ε mi ) 2 ε+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) Optimum Thresholding AMS Sectional Meeting 20 / 35
40 via conditional Mean Square Error (cmse) Key Relationships 1 Clearly, under the previous non-leverage assumption, i X σ, J independent N ( m i, σi 2 ), where m i := n i J = J t i J ti 1 and σi 2 := t i t i 1 σs 2 ds. 2 Then, b i (ε) := E [ ( i X) 2 ] 1 { i X ε} σ, J is such that b i (ε) = σ ( i 2π 2σ i 2 (ε + m i ) + e (ε+m i ) 2 ) 2σ i 2 (ε m i ) e (ε mi ) 2 ε+m i + (mi 2 + σi 2 σ i ) ε m i σ i e x 2 /2 dx 2π (ε m i ) 2 db i (ε) = ε 2 a i (ε), with a i (ε) := e 2σ i 2 + e (ε+m i ) 2σ i 2 dε σ i 2π J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 20 / 35 2
41 via conditional Mean Square Error (cmse) Conditional Mean Square Error (MSE c ) Theorem (F-L & Mancini (2017)) J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 21 / 35
42 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) Optimum Thresholding AMS Sectional Meeting 21 / 35
43 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) Optimum Thresholding AMS Sectional Meeting 21 / 35
44 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) Optimum Thresholding AMS Sectional Meeting 21 / 35
45 via conditional Mean Square Error (cmse) Asymptotics: FA jumps with constant variance Theorem (F-L & Mancini (2017)) J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 22 / 35
46 via conditional Mean Square Error (cmse) Asymptotics: FA jumps 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) Optimum Thresholding AMS Sectional Meeting 22 / 35
47 via conditional Mean Square Error (cmse) Asymptotics: FA jumps 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) Optimum Thresholding AMS Sectional Meeting 22 / 35
48 via conditional Mean Square Error (cmse) A Feasible Implementation of ε n (i) Get a rough" estimate of σ 2 via the realized QV: σ n,0 2 := 1 n X ti X ti 1 2 T i=1 J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 23 / 35
49 via conditional Mean Square Error (cmse) A Feasible Implementation of ε n (i) Get a rough" estimate of σ 2 via the 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 := (2 σ 2 n,0 h n log(1/h n ) ) 1/2 J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 23 / 35
50 via conditional Mean Square Error (cmse) A Feasible Implementation of ε n (i) Get a rough" estimate of σ 2 via the 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 := (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) Optimum Thresholding AMS Sectional Meeting 23 / 35
51 via conditional Mean Square Error (cmse) A Feasible Implementation of ε n (i) Get a rough" estimate of σ 2 via the 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 := (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) Optimum Thresholding AMS Sectional Meeting 23 / 35
52 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) Optimum Thresholding AMS Sectional Meeting 24 / 35
53 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) Optimum Thresholding AMS Sectional Meeting 25 / 35
54 via conditional Mean Square Error (cmse) Heuristics of Proof I: FA jumps with constant variance Recall that the optimum threshold ε n G(ε) := i ( a i (ε) ε j i is a solution of the equation: ) b j (ε) 2IV } {{ } g i (ε) = 0 J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 26 / 35
55 via conditional Mean Square Error (cmse) Heuristics of Proof I: FA jumps with constant variance Recall that the optimum threshold ε n G(ε) := i ( a i (ε) ε j i is a solution of the equation: ) b j (ε) 2IV } {{ } g i (ε) = 0 It can be shown that, as n, σ i 2 2σ i εe ε 2 2σ i 2 + h.o.t., if m 2π i = 0, b i (ε) = σ i m i 2π ε2 e ( m i ε) 2 2σ i 2 + h.o.t., if m i 0. J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 26 / 35
56 via conditional Mean Square Error (cmse) Heuristics of Proof II: FA jumps with constant variance These in turn allows to show that g i (ε) := ε n b j (ε) 2 σj 2 j i j=1 = ε [ ] b j (ε) σj = ε 2 4ε = ε 2 4ε j i:m j =0 j i:m j =0 σ j e ε 2 2σ j 2 2π n σ j e ε 2π j=1 2 2σ 2 j + h.o.t. + h.o.t. j i:m j 0 b j (ε) 2 j:m j 0 σ 2 j J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 27 / 35
57 via conditional Mean Square Error (cmse) Heuristics of Proof III: FA jumps with constant variance Then, up to higher-order terms, ε n ε 4 is a solution of the equation: n σ j e ε 2 2σ j 2 = 0. 2π j=1 In the case of constant σ, we have: ε n 4n σ h n e ε n 2 2σ 2 hn + h.o.t. = 0, 2π which implies ε n 2σ 2 h n ln(1/h n ), n. J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 28 / 35
58 via conditional Mean Square Error (cmse) Generalization I What if t σ t is not constant? J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 29 / 35
59 via conditional Mean Square Error (cmse) Generalization I What if t σ t is not constant? By the mean value theorem, σj 2 ε n = ε n is such that := t j t j 1 σ 2 sds = σ 2 s j h n and, thus, J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 29 / 35
60 via conditional Mean Square Error (cmse) Generalization I What if t σ t is not constant? By the mean value theorem, σj 2 ε n = ε n is such that := t j t j 1 σ 2 sds = σ 2 s j h n and, thus, where w n := ε 2 n/h n. 0 = ε n 4 n σ j e ε 2π j=1 = ε n 4 1 hn n j=1 = ε n 4 hn T 0 2 n 2σ j 2 σ sj 2π e + h.o.t. ε2 n 2σ 2 s j hn h n + h.o.t. σ s 2π e w n 2σ 2 s ds + h.o.t., J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 29 / 35
61 via conditional Mean Square Error (cmse) Generalization II It can be shown that w n = ε 2 n/h n. Then, by the Laplace method, J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 30 / 35
62 via conditional Mean Square Error (cmse) Generalization II It can be shown that w n = ε 2 n/h n. Then, by the Laplace method, T where σ s0 := max s [0,T ] σ s. 0 σ s e w n 2σs 2 ds σ5/2 s σ 0 2π s 0 1 e wn wn 2σ 2 s 0 J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 30 / 35
63 via conditional Mean Square Error (cmse) Generalization II It can be shown that w n = ε 2 n/h n. Then, by the Laplace method, T where σ s0 := max s [0,T ] σ s. 0 σ s e w n 2σs 2 ds σ5/2 s σ 0 2π s 0 1 e wn wn Then, for some constant K and σ := max s [0,T ] σ s, which again implies ε n 0 = ε n K 1 e ε 2 n 2 σ 2 hn ε n + h.o.t., 2 σ 2 h n ln(1/h n ), n. 2σ 2 s 0 J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 30 / 35
64 Ongoing and Future Research Outline 1 Motivational Example 2 The Problem Multipower Variations and Truncated Realized Variations The problem of Parameter Tuning 3 Optimal Threshold Selection via Expected number of jump misclassifications via conditional Mean Square Error (cmse) 4 Ongoing and Future Research J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 31 / 35
65 Ongoing and Future Research Ongoing and Future Research 1 In principle, we can apply the constant-volatility method for varying volatility t σ t by localization; i.e., applying it to periods where σ is approximately constant. J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 32 / 35
66 Ongoing and Future Research Ongoing and Future Research 1 In principle, we can apply the constant-volatility method for varying volatility t σ t by localization; i.e., applying it to periods where σ is approximately constant. 2 This could also help us to estimate σ = max s [0,T ] σ s and, then, apply the asymptotic equivalence ε n 2 σ 2 h n ln(1/h n ). J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 32 / 35
67 Ongoing and Future Research Ongoing and Future Research 3 As it turns, for a Lévy process J and constant σ, the expected [ mean square error, 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) Optimum Thresholding AMS Sectional Meeting 33 / 35
68 Ongoing and Future Research Ongoing and Future Research 3 As it turns, for a Lévy process J and constant σ, the expected [ mean square error, 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 the solution of the equation ε 2 + 2(n 1)E[b 1 (ε)] 2nh n σ 2 = 0. J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 33 / 35
69 Ongoing and Future Research Ongoing and Future Research 3 Furthermore, J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 34 / 35
70 Ongoing and Future Research 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) Optimum Thresholding AMS Sectional Meeting 34 / 35
71 Ongoing and Future Research 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. h n h n J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 34 / 35
72 Ongoing and Future Research 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. Does this phenomenon holds for the minimizer ε n of the cmse? h n h n J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 34 / 35
73 Ongoing and Future Research 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. 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) Optimum Thresholding AMS Sectional Meeting 34 / 35
74 Ongoing and Future Research Further Reading J.E. Figueroa-López & C. Mancini. Optimum thresholding using mean and conditional mean square error. Available at 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. Li, & J. Nisen. Optimal iterative threshold-kernel estimation of jump diffusion processes. In preparation, J.E. Figueroa-López (WUSTL) Optimum Thresholding AMS Sectional Meeting 35 / 35
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 informationShort-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 informationOptimally 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 informationOptimally 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 informationOptimally 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 informationOptimally 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 informationA 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 informationShort-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 informationVolatility. 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 informationShort-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 informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationShort-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 informationOptimal 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 informationI 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 informationEconomics 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 informationLecture 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 informationOptimal 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 informationNear-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 informationAsymptotic 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 informationBROWNIAN 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 informationWeierstrass 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 informationEstimation 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 informationMgr. 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 informationVolume 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 informationArbitrage 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 informationTesting 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 informationAn 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 informationIEOR 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 informationEconomics 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 informationEquity correlations implied by index options: estimation and model uncertainty analysis
1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to
More informationEconomics 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 informationRegression 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 informationParameters Estimation in Stochastic Process Model
Parameters Estimation in Stochastic Process Model A Quasi-Likelihood Approach Ziliang Li University of Maryland, College Park GEE RIT, Spring 28 Outline 1 Model Review The Big Model in Mind: Signal + Noise
More informationControl. 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 informationUltra 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 informationFinancial 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 informationOption 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 informationNumerical 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 informationSmall-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 informationMulti-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 informationJumps 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 informationLecture 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 informationHomework 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 informationSTOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL
STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL YOUNGGEUN YOO Abstract. Ito s lemma is often used in Ito calculus to find the differentials of a stochastic process that depends on time. This paper will introduce
More informationApplications 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 informationEfficient 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 informationUniversitä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 informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationDrunken Birds, Brownian Motion, and Other Random Fun
Drunken Birds, Brownian Motion, and Other Random Fun Michael Perlmutter Department of Mathematics Purdue University 1 M. Perlmutter(Purdue) Brownian Motion and Martingales Outline Review of Basic Probability
More informationExact 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 informationOptimal 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 informationModeling 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 informationA Continuity Correction under Jump-Diffusion Models with Applications in Finance
A Continuity Correction under Jump-Diffusion Models with Applications in Finance Cheng-Der Fuh 1, Sheng-Feng Luo 2 and Ju-Fang Yen 3 1 Institute of Statistical Science, Academia Sinica, and Graduate Institute
More informationLarge 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 informationConvergence 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 information1 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 informationPolynomial 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 informationStrategies 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 informationRealized 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 informationVolatility 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 informationStochastic 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 informationStochastic 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 informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
More informationSADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1. By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD
The Annals of Applied Probability 1999, Vol. 9, No. 2, 493 53 SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1 By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD The use of saddlepoint
More informationThere are no predictable jumps in arbitrage-free markets
There are no predictable jumps in arbitrage-free markets Markus Pelger October 21, 2016 Abstract We model asset prices in the most general sensible form as special semimartingales. This approach allows
More informationOptimum 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 informationLocal Volatility Dynamic Models
René Carmona Bendheim Center for Finance Department of Operations Research & Financial Engineering Princeton University Columbia November 9, 27 Contents Joint work with Sergey Nadtochyi Motivation 1 Understanding
More informationModeling 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 informationOptimal Hedging of Variance Derivatives. John Crosby. Centre for Economic and Financial Studies, Department of Economics, Glasgow University
Optimal Hedging of Variance Derivatives John Crosby Centre for Economic and Financial Studies, Department of Economics, Glasgow University Presentation at Baruch College, in New York, 16th November 2010
More informationIMPLEMENTING 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 informationSimulating 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 informationRecent 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 informationBasic Arbitrage Theory KTH Tomas Björk
Basic Arbitrage Theory KTH 2010 Tomas Björk Tomas Björk, 2010 Contents 1. Mathematics recap. (Ch 10-12) 2. Recap of the martingale approach. (Ch 10-12) 3. Change of numeraire. (Ch 26) Björk,T. Arbitrage
More informationBeyond the Black-Scholes-Merton model
Econophysics Lecture Leiden, November 5, 2009 Overview 1 Limitations of the Black-Scholes model 2 3 4 Limitations of the Black-Scholes model Black-Scholes model Good news: it is a nice, well-behaved model
More informationEE266 Homework 5 Solutions
EE, Spring 15-1 Professor S. Lall EE Homework 5 Solutions 1. A refined inventory model. In this problem we consider an inventory model that is more refined than the one you ve seen in the lectures. The
More informationGPD-POT and GEV block maxima
Chapter 3 GPD-POT and GEV block maxima This chapter is devoted to the relation between POT models and Block Maxima (BM). We only consider the classical frameworks where POT excesses are assumed to be GPD,
More informationOptimal stopping problems for a Brownian motion with a disorder on a finite interval
Optimal stopping problems for a Brownian motion with a disorder on a finite interval A. N. Shiryaev M. V. Zhitlukhin arxiv:1212.379v1 [math.st] 15 Dec 212 December 18, 212 Abstract We consider optimal
More informationVolatility Trading Strategies: Dynamic Hedging via A Simulation
Volatility Trading Strategies: Dynamic Hedging via A Simulation Approach Antai Collage of Economics and Management Shanghai Jiao Tong University Advisor: Professor Hai Lan June 6, 2017 Outline 1 The volatility
More informationPricing and hedging in incomplete markets
Pricing and hedging in incomplete markets Chapter 10 From Chapter 9: Pricing Rules: Market complete+nonarbitrage= Asset prices The idea is based on perfect hedge: H = V 0 + T 0 φ t ds t + T 0 φ 0 t ds
More informationLecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing
Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing We shall go over this note quickly due to time constraints. Key concept: Ito s lemma Stock Options: A contract giving
More informationTangent 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 informationAsymptotic results discrete time martingales and stochastic algorithms
Asymptotic results discrete time martingales and stochastic algorithms Bernard Bercu Bordeaux University, France IFCAM Summer School Bangalore, India, July 2015 Bernard Bercu Asymptotic results for discrete
More informationUsing 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 informationFunctional 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 informationStatistical 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 informationIntroduction 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 informationNeil 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 informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationValuation of Volatility Derivatives. Jim Gatheral Global Derivatives & Risk Management 2005 Paris May 24, 2005
Valuation of Volatility Derivatives Jim Gatheral Global Derivatives & Risk Management 005 Paris May 4, 005 he opinions expressed in this presentation are those of the author alone, and do not necessarily
More informationLecture Quantitative Finance Spring Term 2015
implied Lecture Quantitative Finance Spring Term 2015 : May 7, 2015 1 / 28 implied 1 implied 2 / 28 Motivation and setup implied the goal of this chapter is to treat the implied which requires an algorithm
More information2 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 informationOn 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 informationApproximations of Stochastic Programs. Scenario Tree Reduction and Construction
Approximations of Stochastic Programs. Scenario Tree Reduction and Construction W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany www.mathematik.hu-berlin.de/~romisch
More informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More informationLocal vs Non-local Forward Equations for Option Pricing
Local vs Non-local Forward Equations for Option Pricing Rama Cont Yu Gu Abstract When the underlying asset is a continuous martingale, call option prices solve the Dupire equation, a forward parabolic
More informationValue 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 informationFinal 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 informationQuadratic 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 informationAsymptotic 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 informationLecture 17: More on Markov Decision Processes. Reinforcement learning
Lecture 17: More on Markov Decision Processes. Reinforcement learning Learning a model: maximum likelihood Learning a value function directly Monte Carlo Temporal-difference (TD) learning COMP-424, Lecture
More information