Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models
|
|
- Cynthia Johnston
- 5 years ago
- Views:
Transcription
1 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 de Cargèse May 29th, 2014 (Joint work with Jeff Nisen)
2 Outline 1 The Statistical Problems and the Main Estimators 2 Optimally Thresholded Estimators for Finite-Jump Activity Models 3 Main Results 4 Extensions Additive Processes Stochastic Volatility Processes 5 Conclusions
3 Outline 1 The Statistical Problems and the Main Estimators 2 Optimally Thresholded Estimators for Finite-Jump Activity Models 3 Main Results 4 Extensions Additive Processes Stochastic Volatility Processes 5 Conclusions
4 Set-up 1 Itô Semimartingales: Continuous-time stochastic process t X t with dynamics where dx t = γ t dt + σ t dw t + dj t, t W t is a standard Brownian motion; t J t := N t j=1 ζ j is a piece-wise constant process of finite jump activity; t γ t and t σ t are adapted processes; 2 Finite-Jump Activity (FJA) Lévy Model: where N t X t = γt + σw t + ζ j, {N t} t 0 is a homogeneous Poisson process with jump intensity λ; {ζ j } j 0 are i.i.d. with density f ζ : R R +; j=1 the triplet ({W t}, {N t}, {ζ j }) are mutually independent.
5 Statistical Problems Given a discrete record of observations, X t0, X t1,..., X tn, 0 = t 0 < t 1 < < t n = T, the following problems are of interest in a high-frequency sampling setting (i.e., mesh(π) := max i {t i t i 1 } 0): 1 Estimating the integrated variance (or continuous quadratic variation): σ 2 T := T 0 σ 2 t dt. 2 Estimating the jump features of the process: Jump times τ 1 < τ 2 < < τ NT Jump sizes ζ 1 < ζ 2 < < ζ NT
6 Two main classes of estimators Precursor. Realized Quadratic Variation: n 1 ( QV (X) π := Xti+1 X ) 2 t i, (π : 0 = t0 < < t n = T ). i=0 Under very general conditions: RV (X) π mesh(π) 0 σ 2 T + N T j=1 ζ2 j. 1 Multipower Realized Variations (Barndorff-Nielsen and Shephard (2004)): n 1 BPV (X) T := Xti+1 X ti Xti+2 X ti+1, MPV (X) (r 1,...,r k ) T := i=0 n k X ti+1 X ti r 1... X ti+k X ti+k 1 r k. i=0 2 Threshold Realized Variations (Mancini (2003), Jacod(2007)): n 1 ( TRV (X)[B] π T := Xti+1 X ) 2 t i 1{ Xti+1 Xti B}, (B (0, )). i=0
7 Two main classes of estimators Precursor. Realized Quadratic Variation: n 1 ( QV (X) π := Xti+1 X ) 2 t i, (π : 0 = t0 < < t n = T ). i=0 Under very general conditions: RV (X) π mesh(π) 0 σ 2 T + N T j=1 ζ2 j. 1 Multipower Realized Variations (Barndorff-Nielsen and Shephard (2004)): n 1 BPV (X) T := Xti+1 X ti Xti+2 X ti+1, MPV (X) (r 1,...,r k ) T := i=0 n k X ti+1 X ti r 1... X ti+k X ti+k 1 r k. i=0 2 Threshold Realized Variations (Mancini (2003), Jacod(2007)): n 1 ( TRV (X)[B] π T := Xti+1 X ) 2 t i 1{ Xti+1 Xti B}, (B (0, )). i=0
8 Two main classes of estimators Precursor. Realized Quadratic Variation: n 1 ( QV (X) π := Xti+1 X ) 2 t i, (π : 0 = t0 < < t n = T ). i=0 Under very general conditions: RV (X) π mesh(π) 0 σ 2 T + N T j=1 ζ2 j. 1 Multipower Realized Variations (Barndorff-Nielsen and Shephard (2004)): n 1 BPV (X) T := Xti+1 X ti Xti+2 X ti+1, MPV (X) (r 1,...,r k ) T := i=0 n k X ti+1 X ti r 1... X ti+k X ti+k 1 r k. i=0 2 Threshold Realized Variations (Mancini (2003), Jacod(2007)): n 1 ( TRV (X)[B] π T := Xti+1 X ) 2 t i 1{ Xti+1 Xti B}, (B (0, )). i=0
9 Advantages and Drawbacks 1 Multipower Realized Variations (MPV) Pros: Easy to implement; Cons: Can exhibit high levels of bias and variability in the presence of jumps; e.g., for a FJA Lévy model and regular sampling (t i = ih n with h n = T /n), ] E [MPV (X) (r 1,...,r k ) T C r σ T 2 T C r h max i r i n, with Cr = k i=1 E N (0, 1) r i r r k = 2. 2 Threshold Realized Variations (TRV): Pros: Can be adapted for estimating other jump features: n 1 n 1 N[B] π T := 1 { >B} Xti+1, Ĵ[B] π ( ) 2 X T := Xti+1 X ti 1 { ti Xti+1 Xti i=0 i=0 } >B Cons: Performance strongly depends on a good" choice of the threshold B; e.g., for a FJA Lévy model and regular sampling (t i = ih n with h n = T /n), [ ( E TRV (X)[B n] π T ] σ T 2 Th n γ 2 λσ 2) ( ) Bn Bn 2T σφ +2 T λb3 nc(f ζ ). hn hn 3
10 Advantages and Drawbacks 1 Multipower Realized Variations (MPV) Pros: Easy to implement; Cons: Can exhibit high levels of bias and variability in the presence of jumps; e.g., for a FJA Lévy model and regular sampling (t i = ih n with h n = T /n), ] E [MPV (X) (r 1,...,r k ) T C r σ T 2 T C r h max i r i n, with Cr = k i=1 E N (0, 1) r i r r k = 2. 2 Threshold Realized Variations (TRV): Pros: Can be adapted for estimating other jump features: n 1 n 1 N[B] π T := 1 { >B} Xti+1, Ĵ[B] π ( ) 2 X T := Xti+1 X ti 1 { ti Xti+1 Xti i=0 i=0 } >B Cons: Performance strongly depends on a good" choice of the threshold B; e.g., for a FJA Lévy model and regular sampling (t i = ih n with h n = T /n), [ ( E TRV (X)[B n] π T ] σ T 2 Th n γ 2 λσ 2) ( ) Bn Bn 2T σφ +2 T λb3 nc(f ζ ). hn hn 3
11 Advantages and Drawbacks 1 Multipower Realized Variations (MPV) Pros: Easy to implement; Cons: Can exhibit high levels of bias and variability in the presence of jumps; e.g., for a FJA Lévy model and regular sampling (t i = ih n with h n = T /n), ] E [MPV (X) (r 1,...,r k ) T C r σ T 2 T C r h max i r i n, with Cr = k i=1 E N (0, 1) r i r r k = 2. 2 Threshold Realized Variations (TRV): Pros: Can be adapted for estimating other jump features: n 1 n 1 N[B] π T := 1 { >B} Xti+1, Ĵ[B] π ( ) 2 X T := Xti+1 X ti 1 { ti Xti+1 Xti i=0 i=0 } >B Cons: Performance strongly depends on a good" choice of the threshold B; e.g., for a FJA Lévy model and regular sampling (t i = ih n with h n = T /n), [ ( E TRV (X)[B n] π T ] σ T 2 Th n γ 2 λσ 2) ( ) Bn Bn 2T σφ +2 T λb3 nc(f ζ ). hn hn 3
12 Advantages and Drawbacks 1 Multipower Realized Variations (MPV) Pros: Easy to implement; Cons: Can exhibit high levels of bias and variability in the presence of jumps; e.g., for a FJA Lévy model and regular sampling (t i = ih n with h n = T /n), ] E [MPV (X) (r 1,...,r k ) T C r σ T 2 T C r h max i r i n, with Cr = k i=1 E N (0, 1) r i r r k = 2. 2 Threshold Realized Variations (TRV): Pros: Can be adapted for estimating other jump features: n 1 n 1 N[B] π T := 1 { >B} Xti+1, Ĵ[B] π ( ) 2 X T := Xti+1 X ti 1 { ti Xti+1 Xti i=0 i=0 } >B Cons: Performance strongly depends on a good" choice of the threshold B; e.g., for a FJA Lévy model and regular sampling (t i = ih n with h n = T /n), [ ( E TRV (X)[B n] π T ] σ T 2 Th n γ 2 λσ 2) ( ) Bn Bn 2T σφ +2 T λb3 nc(f ζ ). hn hn 3
13 Advantages and Drawbacks 1 Multipower Realized Variations (MPV) Pros: Easy to implement; Cons: Can exhibit high levels of bias and variability in the presence of jumps; e.g., for a FJA Lévy model and regular sampling (t i = ih n with h n = T /n), ] E [MPV (X) (r 1,...,r k ) T C r σ T 2 T C r h max i r i n, with Cr = k i=1 E N (0, 1) r i r r k = 2. 2 Threshold Realized Variations (TRV): Pros: Can be adapted for estimating other jump features: n 1 n 1 N[B] π T := 1 { >B} Xti+1, Ĵ[B] π ( ) 2 X T := Xti+1 X ti 1 { ti Xti+1 Xti i=0 i=0 } >B Cons: Performance strongly depends on a good" choice of the threshold B; e.g., for a FJA Lévy model and regular sampling (t i = ih n with h n = T /n), [ ( E TRV (X)[B n] π T ] σ T 2 Th n γ 2 λσ 2) ( ) Bn Bn 2T σφ +2 T λb3 nc(f ζ ). hn hn 3 Remark: Bias can attain the rate h n for some suitable B n s
14 Selection of threshold parameter Literature consists of mostly somewhat ad hoc" selection methods for B, which satisfy sufficient conditions for the consistency of the associated estimators. Power Threshold (Mancini (2003)). Pow(α, ω) := α mesh(π) ω, for α > 0 and ω (0, 1/2). Bonferroni Threshold (Gegler & Stadtmüller (2010) and Bollerslev et al. (2007)) ( BF n ( σ, C) := σmesh(π) 1/2 Φ 1 1 C mesh(π) ), for C > 0 and σ > 0. 2
15 Numerical illustration I Diffusion Volatility Parameter (DVP) Estimates RMPV(1, 1) RMPV(2 3, 2 3, 2 3) RMPV(1 2,, 1 2) RMPV(2 5,, 2 5) RMPV(1 3,, 1 3) Min RV(2) Med RV(2) TBPV(Pow(0.05)) TBPV(Pow(0.15)) TBPV(Pow(0.25)) TBPV(Pow(0.35)) TBPV(Pow(0.45)) TBPV(Pow(0.495)) TBPV(B opt) TBPV(BF(0.05)) TBPV(BH(0.05)) Merton Model: Diffusion Volatility Parameter Estimates Actual DVP Multi Power Variation Style Estimators Thresholded Multi Power Style Estimators Multiple Testing Style Estimators Figure: Volatility Estimation Boxplots: Based on 2,500 sample paths (T = 1 year, 5 min sampling frequency). Parameters: σ = 0.3, λ = 20, ζ = D N ( 0.1, ).
16 Numerical illustration II Diffusion Volatility Parameter Estimates TRV(GE) TRV(Hist) TRV(BF) TRV(Pow(0.35)) TRV(Pow(0.4)) TRV(Pow(0.45)) TRV(Pow(0.495)) MPV(1, 1) MPV(2 3,, 2 3) MPV(1 2,, 1 2) MPV(2 5,, 2 5) MPV(1 3,, 1 3) MinRV MedRV Merton Model: Diffusion Volatility Parameter Estimates Calibrated TRV Estimators Uncalibrated TRV Estimators Power Variation Estimators T = 6 months Freq. 5 min. σ = 0.35 λ = δ = Figure: Volatility Estimation Boxplots: Based on 1,000 sample paths. Parameters: σ = 0.35, λ = , ζ = D N (0, ), T = 6-months, sampling frequency = 5-min.
17 Outline 1 The Statistical Problems and the Main Estimators 2 Optimally Thresholded Estimators for Finite-Jump Activity Models 3 Main Results 4 Extensions Additive Processes Stochastic Volatility Processes 5 Conclusions
18 Optimal Threshold Realized Estimators 1 Question: Can the threshold parameter be chosen in a meaningfully objective manner? 2 Aims Develop a well-posed optimal selection criterion for the threshold B, that minimizes a suitable statistical loss function of estimation. Develop a feasible implementation method for B. 3 Assumptions Finite-Jump Activity Lévy Model: N t i.i.d. X t = γt + σw t + ζ j, {N t} t 0 Poisson(λ), ζ j f ζ, ; j=1 Regular sampling scheme with mesh h n := T n ; i.e., π : t i = it n. The jump density function f ζ takes the mixture form: f ζ (x) = pf +(x)1 {x 0} + qf ( x)1 {x<0} with p + q = 1, p, q 0, f ± : [0, ) R + C 1 b(0, )
19 Optimal Threshold Realized Estimators 1 Question: Can the threshold parameter be chosen in a meaningfully objective manner? 2 Aims Develop a well-posed optimal selection criterion for the threshold B, that minimizes a suitable statistical loss function of estimation. Develop a feasible implementation method for B. 3 Assumptions Finite-Jump Activity Lévy Model: N t i.i.d. X t = γt + σw t + ζ j, {N t} t 0 Poisson(λ), ζ j f ζ, ; j=1 Regular sampling scheme with mesh h n := T n ; i.e., π : t i = it n. The jump density function f ζ takes the mixture form: f ζ (x) = pf +(x)1 {x 0} + qf ( x)1 {x<0} with p + q = 1, p, q 0, f ± : [0, ) R + C 1 b(0, )
20 Optimal Threshold Realized Estimators 1 Question: Can the threshold parameter be chosen in a meaningfully objective manner? 2 Aims Develop a well-posed optimal selection criterion for the threshold B, that minimizes a suitable statistical loss function of estimation. Develop a feasible implementation method for B. 3 Assumptions Finite-Jump Activity Lévy Model: N t i.i.d. X t = γt + σw t + ζ j, {N t} t 0 Poisson(λ), ζ j f ζ, ; j=1 Regular sampling scheme with mesh h n := T n ; i.e., π : t i = it n. The jump density function f ζ takes the mixture form: f ζ (x) = pf +(x)1 {x 0} + qf ( x)1 {x<0} with p + q = 1, p, q 0, f ± : [0, ) R + C 1 b(0, )
21 Loss Functions 1 Natural Loss Function Loss (1) n (B) := E [ TRV (X)[B]T T σ 2 2] + E [ N[B]T N T 2 ]. 2 Alternative Loss Function Loss (2) n (B) := E n i=1 ( ) 1 [ n i X >B, n N=0] + 1 i [ n i X B, n N 0], i where, as usual, n i X := X t i X ti 1 and n i N := N t i N ti 1. 3 Interpretation Loss (1) n (B) favors sequences that minimizes the estimation errors of both the continuous and the jump component. Loss (2) n (B) favors sequences that minimizes the total number of miss-classifications: flag a jump within [t i 1, t i ], when there is no jump, and fail to identify a jump within [t i 1, t i ] when there is a jump. Loss (2) n (B) is much more tractable than Loss (1) n (B).
22 Outline 1 The Statistical Problems and the Main Estimators 2 Optimally Thresholded Estimators for Finite-Jump Activity Models 3 Main Results 4 Extensions Additive Processes Stochastic Volatility Processes 5 Conclusions
23 Asymptotic Comparison of Loss Functions Theorem (FL & Nisen (2013)) Given a threshold sequence (B n ) n satisfying B n 0 and B n n, there exists a positive sequence (ε n ) n 1, with lim n ε n = 0, such that Loss (2) n (B) + R n (B) Loss (1) n (B) (1 + ε n )Loss (2) n (B) + R n (B) + R n (B), where, as n, R n (B) T ( ( ) 2 λ2 2σ 2 h n + T 2 Bn φ B n hn σ 2λB n C(f ζ )), h n R n (B) 6 [ T σ 4 h n + 2B 6 nt 2 λ 2 C(f ζ ) 2], C(f ζ ) := pf + (0) + qf (0). In particular, lim inf n B>0 Loss(1) n (B)/ inf B>0 Loss(2) n (B) = 1. Note: Loss (2) n is more analytically tractable: B n := arg min B>0 Loss (2) n (B).
24 Optimal Threshold Parameter Results Theorem: (FL & Nisen (2013)) (1) There exists an N N such that for all n N, the loss function Loss (2) n (B) is quasi-convex and, hence, possesses a unique global minimum Bn. (2) Furthermore, as n, the optimal threshold sequence (Bn) n is such that ( ) ( ) 1 ln 2πσλC(fζ ) σh 1/2 ( ) Bn = 3σ 2 n h h n log n + o, h n 3 log(1/hn ) log(1/h n ) where C(f ζ ) = p + f + (0) + qf (0). Corollary: For the optimally thresholded TRV estimator, as n, Bias(TRV (X)[B n]) h n (γ 2 λσ 2 ), MSE(TRV (X)[B n]) 2σ 4 h n.
25 Optimal Threshold Parameter Results Theorem: (FL & Nisen (2013)) (1) There exists an N N such that for all n N, the loss function Loss (2) n (B) is quasi-convex and, hence, possesses a unique global minimum Bn. (2) Furthermore, as n, the optimal threshold sequence (Bn) n is such that ( ) ( ) 1 ln 2πσλC(fζ ) σh 1/2 ( ) Bn = 3σ 2 n h h n log n + o, h n 3 log(1/hn ) log(1/h n ) where C(f ζ ) = p + f + (0) + qf (0). Corollary: For the optimally thresholded TRV estimator, as n, Bias(TRV (X)[B n]) h n (γ 2 λσ 2 ), MSE(TRV (X)[B n]) 2σ 4 h n.
26 Optimal Threshold Parameter Results Theorem: (FL & Nisen (2013)) (1) There exists an N N such that for all n N, the loss function Loss (2) n (B) is quasi-convex and, hence, possesses a unique global minimum Bn. (2) Furthermore, as n, the optimal threshold sequence (Bn) n is such that ( ) ( ) 1 ln 2πσλC(fζ ) σh 1/2 ( ) Bn = 3σ 2 n h h n log n + o, h n 3 log(1/hn ) log(1/h n ) where C(f ζ ) = p + f + (0) + qf (0). Corollary: For the optimally thresholded TRV estimator, as n, Bias(TRV (X)[B n]) h n (γ 2 λσ 2 ), MSE(TRV (X)[B n]) 2σ 4 h n.
27 Remarks 1 The leading term of the optimal sequence is proportional to the Lévy modulus of the Brownian motion: lim sup h 0 1 2h ln(1/h) 2 The leading order sequence B,1 n := sup W t W s = 1, a.s. t s <h,s,t [0,1] ( 1 3σ 2 h n ln also attains the same rates as B n and provides a blueprint" for devising a good threshold sequence. 3 The second-order term accounts for the volatility σ, the intensity of jumps λ, and the mass concentration of f ζ near 0, as measured by C(f ζ ). h n ),
28 A Feasible Iterative Algorithm to Find B n 1 Key Issue: The optimal threshold B would allow us to find an optimal estimate ˆσ for σ 2 of the form but B depends on precisely σ 2. ˆσ 2 := 1 T TRV (X)[B (σ 2 )] n, 2 The previous issue suggests a fixed-point" type of implementation: Set σ n,0 2 := 1 n T i=1 X t i X ti 1 2 and B ( ) 1/2 n,0 := 3 σ n,0 2 h n ln(1/h n ) while σ 2 n,k 1 > σ2 n,k do σ n,k T TRV (X)[ B n,k ] and B n,k+1 (3 σ n,k+1 2 h n ln(1/h n ) end while { } Let kn := inf k 1 : σ n,k+1 2 = σ2 n,k and take σ n,k 2 as the final n estimate for σ and the corresponding B n,k as an estimate for n B n. 3 The previous algorithm generates a non-increasing sequence of estimators { σ 2 n,k } k and finish in finite time. ) 1/2
29 A Feasible Iterative Algorithm to Find B n 1 Key Issue: The optimal threshold B would allow us to find an optimal estimate ˆσ for σ 2 of the form but B depends on precisely σ 2. ˆσ 2 := 1 T TRV (X)[B (σ 2 )] n, 2 The previous issue suggests a fixed-point" type of implementation: Set σ n,0 2 := 1 n T i=1 X t i X ti 1 2 and B ( ) 1/2 n,0 := 3 σ n,0 2 h n ln(1/h n ) while σ 2 n,k 1 > σ2 n,k do σ n,k T TRV (X)[ B n,k ] and B n,k+1 (3 σ n,k+1 2 h n ln(1/h n ) end while { } Let kn := inf k 1 : σ n,k+1 2 = σ2 n,k and take σ n,k 2 as the final n estimate for σ and the corresponding B n,k as an estimate for n B n. 3 The previous algorithm generates a non-increasing sequence of estimators { σ 2 n,k } k and finish in finite time. ) 1/2
30 A Feasible Iterative Algorithm to Find B n 1 Key Issue: The optimal threshold B would allow us to find an optimal estimate ˆσ for σ 2 of the form but B depends on precisely σ 2. ˆσ 2 := 1 T TRV (X)[B (σ 2 )] n, 2 The previous issue suggests a fixed-point" type of implementation: Set σ n,0 2 := 1 n T i=1 X t i X ti 1 2 and B ( ) 1/2 n,0 := 3 σ n,0 2 h n ln(1/h n ) while σ 2 n,k 1 > σ2 n,k do σ n,k T TRV (X)[ B n,k ] and B n,k+1 (3 σ n,k+1 2 h n ln(1/h n ) end while { } Let kn := inf k 1 : σ n,k+1 2 = σ2 n,k and take σ n,k 2 as the final n estimate for σ and the corresponding B n,k as an estimate for n B n. 3 The previous algorithm generates a non-increasing sequence of estimators { σ 2 n,k } k and finish in finite time. ) 1/2
31 A numerical illustration Merton Model: 4-year / 1-day σ = 0.3 λ = 5 µ = 0, δ = 0.6 Method TRV S TRV Loss S Loss B n,k n Pow ω=0.495;α= BF Table: Finite-sample performance of the threshold realized variation (TRV) estimators i.i.d. based on K = 5, 000 sample paths for the Merton model ζ i N (µ, δ 2 ). Loss represents the total number of Jump Misclassification Errors, while TRV, Loss, S TRV, and S Loss denote the corresponding sample means and standard deviations, respectively.
32 A numerical illustration (S2) Kou Model: 1-week / 5-minute σ = 0.5 λ = 50 p = 0.45, α + = 0.05, α = 0.1 Method TRV S TRV Loss S Loss B n,k n Pow ω=0.495;α= BF Table: Finite-sample performance of the threshold realized variation (TRV) estimators based on K = 5, 000 sample paths for the Kou model: f Kou (x) = p α + e x/α+ 1 [x 0] + (1 p) α e x /α 1 [x<0]. Loss represents the total number of Jump Misclassification Errors, while TRV, Loss, S TRV, and S Loss denote the corresponding sample means and standard deviations, respectively.
33 A numerical illustration (S3) Kou Model: 1-year / 5-minute σ = 0.4 λ = 1000 p = 0.5, α + = α = 0.1 Method TRV S TRV Loss S Loss B n,k n Pow ω=0.495;α= BF Table: Finite-sample performance of the threshold realized variation (TRV) estimators based on K = 5, 000 sample paths for the Kou model: f ζ (x) = p α + e x/α+ 1 [x 0] + q α e x /α 1 [x<0]. Loss represents the total number of Jump Misclassification Errors, while TRV, Loss, S TRV, and S Loss denote the corresponding sample means and standard deviations, respectively.
34 Outline 1 The Statistical Problems and the Main Estimators 2 Optimally Thresholded Estimators for Finite-Jump Activity Models 3 Main Results 4 Extensions Additive Processes Stochastic Volatility Processes 5 Conclusions
35 Outline 1 The Statistical Problems and the Main Estimators 2 Optimally Thresholded Estimators for Finite-Jump Activity Models 3 Main Results 4 Extensions Additive Processes Stochastic Volatility Processes 5 Conclusions
36 Additive Processes 1 The model X s := s 0 γ(u)du + s 0 N s σ(u)dw u + ζ j =: Xs c + J s, where (N s ) s 0 Poiss ({λ(s)} s 0 ), independent of W, and deterministic smooth functions σ, λ : [0, ) R + and γ : [0, ) R with σ and λ bounded away from 0. 2 Optimal Threshold Problem Given a sampling scheme π : t 0 < < t n = T, determine the vector B π, = (B π, t 1 inf E B=(Bt1,...,B tn ) R m + n = i=1 j=1,..., B π, t n ) that minimizes the problem n i=1 ( 1 [ Xti X ti 1 >B ti,n ti N ti 1 =0] + 1 [ Xti X ti 1 B ti,n ti N ti 1 0] inf {P ( i X > B ti, i N = 0) + P ( i X B ti, i N 0)}, B ti ( i X := X ti X ti 1, i N := N ti N ti 1 ) )
37 Additive Processes 1 The model X s := s 0 γ(u)du + s 0 N s σ(u)dw u + ζ j =: Xs c + J s, where (N s ) s 0 Poiss ({λ(s)} s 0 ), independent of W, and deterministic smooth functions σ, λ : [0, ) R + and γ : [0, ) R with σ and λ bounded away from 0. 2 Optimal Threshold Problem Given a sampling scheme π : t 0 < < t n = T, determine the vector B π, = (B π, t 1 inf E B=(Bt1,...,B tn ) R m + n = i=1 j=1,..., B π, t n ) that minimizes the problem n i=1 ( 1 [ Xti X ti 1 >B ti,n ti N ti 1 =0] + 1 [ Xti X ti 1 B ti,n ti N ti 1 0] inf {P ( i X > B ti, i N = 0) + P ( i X B ti, i N 0)}, B ti ( i X := X ti X ti 1, i N := N ti N ti 1 ) )
38 Additive Processes 1 The model X s := s 0 γ(u)du + s 0 N s σ(u)dw u + ζ j =: Xs c + J s, where (N s ) s 0 Poiss ({λ(s)} s 0 ), independent of W, and deterministic smooth functions σ, λ : [0, ) R + and γ : [0, ) R with σ and λ bounded away from 0. 2 Optimal Threshold Problem Given a sampling scheme π : t 0 < < t n = T, determine the vector B π, = (B π, t 1 inf E B=(Bt1,...,B tn ) R m + n = i=1 j=1,..., B π, t n ) that minimizes the problem n i=1 ( 1 [ Xti X ti 1 >B ti,n ti N ti 1 =0] + 1 [ Xti X ti 1 B ti,n ti N ti 1 0] inf {P ( i X > B ti, i N = 0) + P ( i X B ti, i N 0)}, B ti ( i X := X ti X ti 1, i N := N ti N ti 1 ) )
39 Well-posedness and Asymptotic Characterization Theorem (FL & Nisen (2013)) For any fixed T > 0, there exists h 0 := h 0 (T ) > 0 such that, for all t [0, T ] and h (0, h 0 ], the functions L t,h (B) := P( X t+h X t > B, N t+h N t = 0)+P( X t+h X t B, N t+h N t 0), are quasi-convex and possess a unique global minimum, Bt,h, such that, as h 0 +, B t,h = 3σ 2 (t)h ln(1/h) log( 2πσ(t)λC(f ζ ))σ(t)h 1/2 3 log(1/h) + o( h/ ln(1/h)).
40 Spot Volatility Estimation via Kernel Methods Notation: h i = t i t i 1 (Mesh), K θ (t) = 1 θ K ( t θ ) (Kernel), θ = Bandwidth Algorithm: [Based on kernel type estimators by Kristensen(2010)] For each i {1, 2,..., n}, set σ 2 0(t i ) := l j= l 1 h i+j i+j X 2 K θ (t i t i+j ) and B,0 t i := [ 3 σ 2 0(t i )h i ln(1/h i ) ] 1/2 while there exists i {1, 2,..., m} such that σ k 1 2 (t i) > σ k 2(t i) do σ k+1 2 (t i) l 1 j= l h i+j i+j X 2 1 [ ] K i+j X B,k θ (t i t i+j ) and t i+j B,k+1 t i [ 3 σ k+1 2 (t i)h i ln(1/h i ) ] 1/2 end while Let k := inf { k 1 : σ k+1 2 (t i) = σ k 2(t i); for all i = 1, 2,..., n } and take σ k 2 (t,k i) as the final estimate for σ(t i ) and the corresponding B t i as an estimate for Bt i. The previous algorithm generates a non-increasing sequence of estimators { σ 2 k (t i)} k, for each i, and finish in finite time.
41 Spot Volatility Estimation via Kernel Methods Notation: h i = t i t i 1 (Mesh), K θ (t) = 1 θ K ( t θ ) (Kernel), θ = Bandwidth Algorithm: [Based on kernel type estimators by Kristensen(2010)] For each i {1, 2,..., n}, set σ 2 0(t i ) := l j= l 1 h i+j i+j X 2 K θ (t i t i+j ) and B,0 t i := [ 3 σ 2 0(t i )h i ln(1/h i ) ] 1/2 while there exists i {1, 2,..., m} such that σ k 1 2 (t i) > σ k 2(t i) do σ k+1 2 (t i) l 1 j= l h i+j i+j X 2 1 [ ] K i+j X B,k θ (t i t i+j ) and t i+j B,k+1 t i [ 3 σ k+1 2 (t i)h i ln(1/h i ) ] 1/2 end while Let k := inf { k 1 : σ k+1 2 (t i) = σ k 2(t i); for all i = 1, 2,..., n } and take σ k 2 (t,k i) as the final estimate for σ(t i ) and the corresponding B t i as an estimate for Bt i. The previous algorithm generates a non-increasing sequence of estimators { σ 2 k (t i)} k, for each i, and finish in finite time.
42 Spot Volatility Estimation via Kernel Methods Notation: h i = t i t i 1 (Mesh), K θ (t) = 1 θ K ( t θ ) (Kernel), θ = Bandwidth Algorithm: [Based on kernel type estimators by Kristensen(2010)] For each i {1, 2,..., n}, set σ 2 0(t i ) := l j= l 1 h i+j i+j X 2 K θ (t i t i+j ) and B,0 t i := [ 3 σ 2 0(t i )h i ln(1/h i ) ] 1/2 while there exists i {1, 2,..., m} such that σ k 1 2 (t i) > σ k 2(t i) do σ k+1 2 (t i) l 1 j= l h i+j i+j X 2 1 [ ] K i+j X B,k θ (t i t i+j ) and t i+j B,k+1 t i [ 3 σ k+1 2 (t i)h i ln(1/h i ) ] 1/2 end while Let k := inf { k 1 : σ k+1 2 (t i) = σ k 2(t i); for all i = 1, 2,..., n } and take σ k 2 (t,k i) as the final estimate for σ(t i ) and the corresponding B t i as an estimate for Bt i. The previous algorithm generates a non-increasing sequence of estimators { σ 2 k (t i)} k, for each i, and finish in finite time.
43 Spot Volatility Estimation via Kernel Methods Notation: h i = t i t i 1 (Mesh), K θ (t) = 1 θ K ( t θ ) (Kernel), θ = Bandwidth Algorithm: [Based on kernel type estimators by Kristensen(2010)] For each i {1, 2,..., n}, set σ 2 0(t i ) := l j= l 1 h i+j i+j X 2 K θ (t i t i+j ) and B,0 t i := [ 3 σ 2 0(t i )h i ln(1/h i ) ] 1/2 while there exists i {1, 2,..., m} such that σ k 1 2 (t i) > σ k 2(t i) do σ k+1 2 (t i) l 1 j= l h i+j i+j X 2 1 [ ] K i+j X B,k θ (t i t i+j ) and t i+j B,k+1 t i [ 3 σ k+1 2 (t i)h i ln(1/h i ) ] 1/2 end while Let k := inf { k 1 : σ k+1 2 (t i) = σ k 2(t i); for all i = 1, 2,..., n } and take σ k 2 (t,k i) as the final estimate for σ(t i ) and the corresponding B t i as an estimate for Bt i. The previous algorithm generates a non-increasing sequence of estimators { σ 2 k (t i)} k, for each i, and finish in finite time.
44 Illustration of Opt. Thresh. Spot Vol. Estimation Alg. (A) Initial Estimates (B) Intermediate Estimates Actual Spot Volatility Est. Spot Vol. (Uniform) Est. Spot Vol. (Quad) Sample Increments Actual Spot Volatility Est. Spot Vol. (Uniform) Est. Spot Vol. (Quad) Sample Increments Time Horizon Time Horizon Figure: Estimation of Spot Volatility using Adaptive Kernel Weighted Realized Volatility. (A) The initial estimates. (B) Intermediate estimates. Parameters: γ(t) = 0.1t, σ(t) = 4.5t sin(2πe t2 ) , λ(t) = 25(e 3t 1), ζ i i.i.d. = D N (µ = 0.025, δ = 0.025).
45 Illustration of Opt. Thresh. Spot Vol. Estimation Alg. (C) Terminal Estimates Actual Spot Volatility Est. Spot Vol. (Uniform) Est. Spot Vol. (Quad) Sample Increments Time Horizon Figure: Estimation of Spot Volatility using Adaptive Kernel Weighted Realized Volatility. (C) The terminal estimates. (D) Estimation variability, based on 100 generated sample paths, for the Quadratic Kernel based estimator. Parameters: γ(t) = 0.1t, σ(t) = 4.5t sin(2πe t2 ) , λ(t) = 25(e 3t 1), ζ i i.i.d. = D N (µ = 0.025, δ = 0.025).
46 Outline 1 The Statistical Problems and the Main Estimators 2 Optimally Thresholded Estimators for Finite-Jump Activity Models 3 Main Results 4 Extensions Additive Processes Stochastic Volatility Processes 5 Conclusions
47 Stochastic Volatility Models with FJA 1 Motivation In practice (e.g., finance), one usually encounters models like X s := s 0 γ u du + s 0 N s σ u dw u + ζ j =: Xs c + J s, where {σ t } t 0 is itself stochastic, independent of (W, N, ζ). 2 Prototypical Example: Mean-reverting square-root process (CIR Model): dσ 2 t = κ ( α σ 2 t 3 Pitfalls of the Estimation Methods for σ: j=1 ) dt + γσt dw (σ) t, (2κα γ 2 > 0, W (σ) W ) Spot volatility estimation has received limited" attention in the literature: e.g., Kristensen(2010), Reno & Mancini(2013) consider Kernel Estimators; Fan & Wang(2008), Alvarez et al.(2010), and others consider estimators for [X, X] t = t 0 σ2 s ds and a finite-difference approx. of its derivative; Available kernel estimation methods for {σt 2 } t 0 (say, Kristensen s method) are quite sensitive to the bandwidth. Few Bandwidth Selection Methods (cf. Kristensen(2010))
48 Stochastic Volatility Models with FJA 1 Motivation In practice (e.g., finance), one usually encounters models like X s := s 0 γ u du + s 0 N s σ u dw u + ζ j =: Xs c + J s, where {σ t } t 0 is itself stochastic, independent of (W, N, ζ). 2 Prototypical Example: Mean-reverting square-root process (CIR Model): dσ 2 t = κ ( α σ 2 t 3 Pitfalls of the Estimation Methods for σ: j=1 ) dt + γσt dw (σ) t, (2κα γ 2 > 0, W (σ) W ) Spot volatility estimation has received limited" attention in the literature: e.g., Kristensen(2010), Reno & Mancini(2013) consider Kernel Estimators; Fan & Wang(2008), Alvarez et al.(2010), and others consider estimators for [X, X] t = t 0 σ2 s ds and a finite-difference approx. of its derivative; Available kernel estimation methods for {σt 2 } t 0 (say, Kristensen s method) are quite sensitive to the bandwidth. Few Bandwidth Selection Methods (cf. Kristensen(2010))
49 Stochastic Volatility Models with FJA 1 Motivation In practice (e.g., finance), one usually encounters models like X s := s 0 γ u du + s 0 N s σ u dw u + ζ j =: Xs c + J s, where {σ t } t 0 is itself stochastic, independent of (W, N, ζ). 2 Prototypical Example: Mean-reverting square-root process (CIR Model): dσ 2 t = κ ( α σ 2 t 3 Pitfalls of the Estimation Methods for σ: j=1 ) dt + γσt dw (σ) t, (2κα γ 2 > 0, W (σ) W ) Spot volatility estimation has received limited" attention in the literature: e.g., Kristensen(2010), Reno & Mancini(2013) consider Kernel Estimators; Fan & Wang(2008), Alvarez et al.(2010), and others consider estimators for [X, X] t = t 0 σ2 s ds and a finite-difference approx. of its derivative; Available kernel estimation methods for {σt 2 } t 0 (say, Kristensen s method) are quite sensitive to the bandwidth. Few Bandwidth Selection Methods (cf. Kristensen(2010))
50 Stochastic Volatility Models with FJA 1 Motivation In practice (e.g., finance), one usually encounters models like X s := s 0 γ u du + s 0 N s σ u dw u + ζ j =: Xs c + J s, where {σ t } t 0 is itself stochastic, independent of (W, N, ζ). 2 Prototypical Example: Mean-reverting square-root process (CIR Model): dσ 2 t = κ ( α σ 2 t 3 Pitfalls of the Estimation Methods for σ: j=1 ) dt + γσt dw (σ) t, (2κα γ 2 > 0, W (σ) W ) Spot volatility estimation has received limited" attention in the literature: e.g., Kristensen(2010), Reno & Mancini(2013) consider Kernel Estimators; Fan & Wang(2008), Alvarez et al.(2010), and others consider estimators for [X, X] t = t 0 σ2 s ds and a finite-difference approx. of its derivative; Available kernel estimation methods for {σt 2 } t 0 (say, Kristensen s method) are quite sensitive to the bandwidth. Few Bandwidth Selection Methods (cf. Kristensen(2010))
51 Stochastic Volatility Models with FJA 1 Motivation In practice (e.g., finance), one usually encounters models like X s := s 0 γ u du + s 0 N s σ u dw u + ζ j =: Xs c + J s, where {σ t } t 0 is itself stochastic, independent of (W, N, ζ). 2 Prototypical Example: Mean-reverting square-root process (CIR Model): dσ 2 t = κ ( α σ 2 t 3 Pitfalls of the Estimation Methods for σ: j=1 ) dt + γσt dw (σ) t, (2κα γ 2 > 0, W (σ) W ) Spot volatility estimation has received limited" attention in the literature: e.g., Kristensen(2010), Reno & Mancini(2013) consider Kernel Estimators; Fan & Wang(2008), Alvarez et al.(2010), and others consider estimators for [X, X] t = t 0 σ2 s ds and a finite-difference approx. of its derivative; Available kernel estimation methods for {σt 2 } t 0 (say, Kristensen s method) are quite sensitive to the bandwidth. Few Bandwidth Selection Methods (cf. Kristensen(2010))
52 A Proposed Model Selection Method 1 Technical Assumption (Kristensen, 2010): σ 2 t+δ σ 2 t = L t (δ)δ υ + o P (δ υ ), a.s. (δ 0) ( ) where υ (0, 1] and δ L t (δ) is slowly varying (random) function at 0 and t L t (0) := lim δ 0 + L t (δ) is continuous. 2 Under ( ), Kristensen(2010) argues that ( ) bw loc opt,t = n 1 σ 4 2υ+1 t K 1 2 2υ+1 2 υl 2 t (0) ( ) 1 bw glb opt = n 1 T 2υ+1 0 σ4 t dt K 2 2υ+1 2 υ T 0 L2 t (0)dt ( ) 3 Pitfall: In general, it is hard to check ( ) with explicit (nonzero) L t (0) and υ. 4 A (heuristic) alternative method: Suppose that [ (σ 2 E t+δ σt 2 ) 2 ] Ft = L 2 t (0)δ 2υ + o P (δ υ ), a.s. (δ 0), ( ) for some positive adapted {L t (0)} t 0. Then, use ( ).
53 Computation for the CIR model 1 Recall σ 2 t+δ σ 2 t = t+δ t κ ( α σu 2 ) t+δ du + γ t σ u dw (σ) u, 2 The leading term of ( σ 2 t+δ σt 2 2 ) is γ 2 t+δ 2, σ t u dw (σ) u which is such that ( ) 2 t+δ ( ) t+δ E γ 2 σ u dw (σ) u t F t = γ 2 E σudu 2 t F t = γ 2 σt 2 δ + o(δ). 3 The previous heuristics suggest that ( σ 2 E t+δ σt 2 2 ) Ft = γ 2 σt 2 δ + o(δ). 4 Hence, ( ) holds with υ = 1 2 and L2 t (0) = γ2 σ 2 t. 5 This implies the following local optimal estimator: ( 2σ bw loc 2 opt,t = t K 2 2 ) 1/2 n 1/2. γ 2
54 Outlined of Estimation Method 1 Get a rough" estimate t σt 2 ; e.g., using Alvarez et al. (2010), [X, ˆσ 0(t 2 X] ti + δ i ) = n [X, X] ti = 1 δn δ n 2 Refine ˆσ 2 0 (t i) using thresholding, with j:t j (t i,t i + δ n] ˆσ 1(t 2 i ) = 1 ( n δ j X ) 2 1 [ ] n j:t i (t i,t i + j X ˆB (0) t i δ n] B (0) t i := [ 3 σ 2 0 (t i)h i ln(1/h i ) ] 1/2. ( n j X ) 2. 3 Estimate γ using the realized variation of {ˆσ 1 (t i )} i=1,...,n, since σ, σ t = γ Apply iterated kernel estimation with thresholding and (estimated) optimal local bandwidth selection.
55 Numerical Results: Uniform Kernel Model: Normal jump sizes and CIR stochastic volatility. Parameters: γ drift = 0.05, λ = 60, ζ i N (0, 5(0.3)), κ = 5, α = 0.04; γ volatiity = 0.5. Sampling scheme: T=1/12 (one month) and h n = 5 min Monte Carlo results of MSE = n i=1 (ˆσ 2 (t i ) σ 2 (t i ) ) 2 based on 500 runs Method MSE Alvarez et al. Method Alvarez et al. Method with thresholding Kernel Est. with loc. bw selection (not thresholding) Kernel Est. with thresholding and loc. bw selection Oracle Kernel Est. with thresholding and loc. bw selection
56 Outline 1 The Statistical Problems and the Main Estimators 2 Optimally Thresholded Estimators for Finite-Jump Activity Models 3 Main Results 4 Extensions Additive Processes Stochastic Volatility Processes 5 Conclusions
57 Conclusions 1 Introduce an objective threshold selection procedure based on statistical optimality reasoning via a well-posed optimization problem. 2 Characterize precisely the infill asymptotic behavior of the optimal threshold sequence. 3 Proposed an iterative algorithm to find the optimal threshold sequence. 4 Extend the approach to more general stochastic models, which allows time-varying volatility and jump intensity.
58 For Further Reading I Figueroa-López & Nisen. Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models Stochastic Processes and their Applications 123(7), , Available at figueroa. Figueroa-López & Nisen. Optimality properties of thresholded multi power variation estimators. In preparation, 2014.
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 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 Financial Statistics Stevanovich Center The University
More informationAsymptotic 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 informationOptimum 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 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 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 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 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 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 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 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 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 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 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 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 informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November
More 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 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 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 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 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 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 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 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 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 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 informationOptimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing
Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Haindorf, 7 Februar 2008 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar
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 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 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 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 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 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 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 informationRisk 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 informationEstimation 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 informationModeling and Pricing of Variance Swaps for Local Stochastic Volatilities with Delay and Jumps
Modeling and Pricing of Variance Swaps for Local Stochastic Volatilities with Delay and Jumps Anatoliy Swishchuk Department of Mathematics and Statistics University of Calgary Calgary, AB, Canada QMF 2009
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 informationOption Pricing Modeling Overview
Option Pricing Modeling Overview Liuren Wu Zicklin School of Business, Baruch College Options Markets Liuren Wu (Baruch) Stochastic time changes Options Markets 1 / 11 What is the purpose of building a
More 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 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 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 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 informationRough volatility models: When population processes become a new tool for trading and risk management
Rough volatility models: When population processes become a new tool for trading and risk management Omar El Euch and Mathieu Rosenbaum École Polytechnique 4 October 2017 Omar El Euch and Mathieu Rosenbaum
More 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 informationChapter 3: Black-Scholes Equation and Its Numerical Evaluation
Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random
More informationImportance 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 informationNumerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps
Numerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps, Senior Quantitative Analyst Motivation: Swing Options An electricity or gas SUPPLIER needs to be capable,
More informationOptimal robust bounds for variance options and asymptotically extreme models
Optimal robust bounds for variance options and asymptotically extreme models Alexander Cox 1 Jiajie Wang 2 1 University of Bath 2 Università di Roma La Sapienza Advances in Financial Mathematics, 9th January,
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 informationEconomics 2010c: Lecture 4 Precautionary Savings and Liquidity Constraints
Economics 2010c: Lecture 4 Precautionary Savings and Liquidity Constraints David Laibson 9/11/2014 Outline: 1. Precautionary savings motives 2. Liquidity constraints 3. Application: Numerical solution
More informationValuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model
Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model 1(23) Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility
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 informationConditional 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 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 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 information3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors
3.4 Copula approach for modeling default dependency Two aspects of modeling the default times of several obligors 1. Default dynamics of a single obligor. 2. Model the dependence structure of defaults
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 informationParametric Inference and Dynamic State Recovery from Option Panels. Torben G. Andersen
Parametric Inference and Dynamic State Recovery from Option Panels Torben G. Andersen Joint work with Nicola Fusari and Viktor Todorov The Third International Conference High-Frequency Data Analysis in
More 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 informationinduced 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 informationChapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets
Chapter 5: Jump Processes and Incomplete Markets Jumps as One Explanation of Incomplete Markets It is easy to argue that Brownian motion paths cannot model actual stock price movements properly in reality,
More informationInterest rate models in continuous time
slides for the course Interest rate theory, University of Ljubljana, 2012-13/I, part IV József Gáll University of Debrecen Nov. 2012 Jan. 2013, Ljubljana Continuous time markets General assumptions, notations
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Generating Random Variables and Stochastic Processes Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationNon-semimartingales in finance
Non-semimartingales in finance Pricing and Hedging Options with Quadratic Variation Tommi Sottinen University of Vaasa 1st Northern Triangular Seminar 9-11 March 2009, Helsinki University of Technology
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 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 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 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 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 informationAsset Pricing Models with Underlying Time-varying Lévy Processes
Asset Pricing Models with Underlying Time-varying Lévy Processes Stochastics & Computational Finance 2015 Xuecan CUI Jang SCHILTZ University of Luxembourg July 9, 2015 Xuecan CUI, Jang SCHILTZ University
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Other Miscellaneous Topics and Applications of Monte-Carlo Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationLikelihood Estimation of Jump-Diffusions
Likelihood Estimation of Jump-Diffusions Extensions from Diffusions to Jump-Diffusions, Implementation with Automatic Differentiation, and Applications Berent Ånund Strømnes Lunde DEPARTMENT OF MATHEMATICS
More informationGeneralized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models
Generalized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models Bilkan Erkmen (joint work with Michael Coulon) Workshop on Stochastic Games, Equilibrium, and Applications
More informationMSc Financial Engineering CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL. To be handed in by monday January 28, 2013
MSc Financial Engineering 2012-13 CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL To be handed in by monday January 28, 2013 Department EMS, Birkbeck Introduction The assignment consists of Reading
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 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 informationRisk Neutral Measures
CHPTER 4 Risk Neutral Measures Our aim in this section is to show how risk neutral measures can be used to price derivative securities. The key advantage is that under a risk neutral measure the discounted
More informationM5MF6. Advanced Methods in Derivatives Pricing
Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................
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 informationEffectiveness of CPPI Strategies under Discrete Time Trading
Effectiveness of CPPI Strategies under Discrete Time Trading S. Balder, M. Brandl 1, Antje Mahayni 2 1 Department of Banking and Finance, University of Bonn 2 Department of Accounting and Finance, Mercator
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 informationHedging Under Jump Diffusions with Transaction Costs. Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo
Hedging Under Jump Diffusions with Transaction Costs Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo Computational Finance Workshop, Shanghai, July 4, 2008 Overview Overview Single factor
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets Liuren Wu ( c ) The Black-Merton-Scholes Model colorhmoptions Markets 1 / 18 The Black-Merton-Scholes-Merton (BMS) model Black and Scholes (1973) and Merton
More 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 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 informationQI SHANG: General Equilibrium Analysis of Portfolio Benchmarking
General Equilibrium Analysis of Portfolio Benchmarking QI SHANG 23/10/2008 Introduction The Model Equilibrium Discussion of Results Conclusion Introduction This paper studies the equilibrium effect of
More informationHedging under Arbitrage
Hedging under Arbitrage Johannes Ruf Columbia University, Department of Statistics Modeling and Managing Financial Risks January 12, 2011 Motivation Given: a frictionless market of stocks with continuous
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationLecture 3: Review of mathematical finance and derivative pricing models
Lecture 3: Review of mathematical finance and derivative pricing models Xiaoguang Wang STAT 598W January 21th, 2014 (STAT 598W) Lecture 3 1 / 51 Outline 1 Some model independent definitions and principals
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 informationControl Improvement for Jump-Diffusion Processes with Applications to Finance
Control Improvement for Jump-Diffusion Processes with Applications to Finance Nicole Bäuerle joint work with Ulrich Rieder Toronto, June 2010 Outline Motivation: MDPs Controlled Jump-Diffusion Processes
More informationLIBOR models, multi-curve extensions, and the pricing of callable structured derivatives
Weierstrass Institute for Applied Analysis and Stochastics LIBOR models, multi-curve extensions, and the pricing of callable structured derivatives John Schoenmakers 9th Summer School in Mathematical Finance
More informationExponential utility maximization under partial information
Exponential utility maximization under partial information Marina Santacroce Politecnico di Torino Joint work with M. Mania AMaMeF 5-1 May, 28 Pitesti, May 1th, 28 Outline Expected utility maximization
More informationRisk 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 informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationLogarithmic derivatives of densities for jump processes
Logarithmic derivatives of densities for jump processes Atsushi AKEUCHI Osaka City University (JAPAN) June 3, 29 City University of Hong Kong Workshop on Stochastic Analysis and Finance (June 29 - July
More information