Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility

Size: px
Start display at page:

Download "Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility"

Transcription

1 Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes, and Malliavin Calculus: Recent Applications Barcelona GSE Summer Forum June 26, 214 Joint work with Sveinn Ólafsson from Purdue University

2 Outline 1 The Setup Tempered-Stable-Like Processes A stochastic volatility financial model with Lévy jumps 2 Problem Formulation Some relevant literature 3 High-order Expansions for (near) At-the-money Options Tempered-stable-like Lévy with and without Brownian component Stochastic volatility model with tempered-stable-like jumps 4 Conclusions

3 The Setup Tempered-Stable-Like Processes Lévy Process 1 Lévy process {X t } t X = Independent Increments: t < t 1 < < t n = X t1 X t,..., X tn X tn 1 are independent Stationary Increments s < t = X t X D s = X t s The trajectories of the process, t X t(ω), are right-continuous with left-limits 2 The distribution law of {X t } t is determined by the distribution of X 1 : If L(X 1 ) N (, 1), then X t = W t is a standard Brownian Motion; If L(X 1 ) Poisson(λ), then X t = N t is a Poisson process with intensity λ;

4 The Setup Tempered-Stable-Like Processes Tempered-stable-like processes 1 A Lévy process {X t } t whose distribution at t = 1 has a characteristic function of the form: E ( ( ) e iux ) ( 1 = exp ibu + e iux ) 1 iux1 { x 1} s(x)dx, R\{} ( where s(x) = C x x ) q(x) x α 1 for some constants C(1), C( 1) (, ), α (, 2), and a function q : R\{} [, ) such that q(x) x 1, sup q(x) <. x 2 b, α, and s are called the drift", the index of jump activity, and the Lévy density of the process.

5 The Setup Tempered-Stable-Like Processes Connection to Stable Processes 1 If q(x) 1, the resulting Lévy process is a Stable Lévy Process {Z t } t ; 2 For a suitable c R, the process Z t := Z t ct is self-similar: {h 1/α Zht } t D = { Zt } t (h > ). If α > 1, c = EZ 1 ; If α < 1, c = xs(x)dx is the drift of Z. x 1 3 Distributions are often too fat" for some applications (e.g., finance): E( Z t p ) =, for any p > α. 4 The function q(x) can mitigate the intensity of large jumps so that E( X t p ) < x p α 1 q(x)dx <. x 1 5 Notation: Throughout, Z is called the strictly stable processes ( associated ) with X. This is characterized by its Lévy density s(x) := C x α 1 and that its center" is. x x

6 The Setup Tempered-Stable-Like Processes Short and long time behavior of TSP 1 In short-time, {X t } t behaves like a stable process (cf. Rosenbaum and Tankov (211)): 1 < α < 2: {h 1/α X ht } D { Z t} t, (h ), for the strictly α-stable process Z associated with X; < α < 1: {h 1/α (X ht cht)} D { Z t} t, (h ), for the strictly α-stable process Z associated with X; 2 In long-time, {X t } t behaves like a Brownian Motion: {h 1/2 X ht } D {B t } t, (h ), where {B t } t is a suitable Brownian motion.

7 The Setup Tempered-Stable-Like Processes Short and long time behavior of TSP 1 In short-time, {X t } t behaves like a stable process (cf. Rosenbaum and Tankov (211)): 1 < α < 2: {h 1/α X ht } D { Z t} t, (h ), for the strictly α-stable process Z associated with X; < α < 1: {h 1/α (X ht cht)} D { Z t} t, (h ), for the strictly α-stable process Z associated with X; 2 In long-time, {X t } t behaves like a Brownian Motion: {h 1/2 X ht } D {B t } t, (h ), where {B t } t is a suitable Brownian motion.

8 The Setup A stochastic volatility financial model with Lévy jumps Why stochastic volatility? 1 The Black-Scholes model (S t = S e Xt with X t = σw t + bt) offers tractable solutions, but is not in line with many stylized features observed in asset and option prices. 2 Exponential Lévy models (S t = S e Xt with a Lévy process X t ) generalize the Black-Scholes framework by allowing jumps in stock prices while preserving the independence and stationarity of returns. 3 They allow to generate implied volatility smiles similar to the ones observed in practice. 4 Exponential Lévy models still have their shortcomings: Independence of log-returns is not consistent with some fine stylized features of return s time series (e.g., volatility clustering and leverage) Exponential Lévy models are not capable of capturing the evolution of the implied volatility surface in time

9 The Setup A stochastic volatility financial model with Lévy jumps Introducing a stochastic volatility component 1 Consider a log-return process L t := log (S t /S ) such that dl t = dv t }{{} Continuous Component + dx t }{{} Jump Component = µ(y t )dt + σ(y t )dw t + dx t ; 2 Y = {Y t } t is a latent (hidden) risky factor with dynamics dy t = α(y t )dt + γ(y t )dw t, Y = y ; 3 X = {X t } t is a pure-jump tempered-stable-like process as before; 4 (W, W ) are correlated B.M. s independent of X; 5 If µ(y) and σ(y) σ(y ), L is a Lévy process with triplet (σ(y ), b, ν).

10 Problem Formulation The General Problems 1 Let us first take the drifts b and µ such that S t := e Lt is a martingale; in particular, q must be such that x 1 ex q(x)x α 1 dx < and E (S t ) = E ( e ) Lt = 1; 2 Consider the functional: [ (e L Π t := E t 1 ) ] + = E [ max ( e Lt 1, )]. 3 By Dominated Convergence Theorem, Π t when t. 4 The General Problems: We want to characterize the rate of convergence of Π t as t ; Characterize the effect of the different model s parameters in the short-time asymptotic behavior of Π t.

11 Problem Formulation Motivation 1 In finance, S E [ ( e Lt 1 ) + ] represents the price of an At-the-money (ATM) European Call Option with time-to-expiration t written on a stock whose price process is modeled by t S t := S e Lt. 2 Our results shed light on the behavior of ATM option prices close to expiration under a stochastic volatility model with Lévy type jumps. 3 In general, the call option price with time-to-expiry t and strike K = S e κ is [ E (S t K ) +] = S E [ ( e Xt e κ) + ]. 4 In mathematics, ϕ K (S) = (S K ) + are natural building blocks of convex functions f : R + R + : f (S) = f () + f +()S + (S K ) + µ(dk ).

12 Problem Formulation Some relevant literature Some relevant literature Two distinct regimes: Not ATM and ATM. Not ATM (K S κ ) 1 Tankov (211): Leading order term for general Lévy process: ( E e X t e κ) + (e = (1 e κ ) + + t x e κ)+ s(x)dx + o(t). 2 F-L & Forde (212): High-order term for quite general Lévy process; ( E e X t e κ) + (e = (1 e κ ) + + t x e κ)+ s(x)dx + t 2 2 d 2(κ) + o(t 2 ), where d 2 (κ) has an explicit form in terms of s.

13 Problem Formulation Some relevant literature Some relevant literature. Cont... ATM (K = S κ = ) 1 Tankov (211), Roper (211): Leading order term for a general bounded variation Lévy process (α < 1): ( + { } E e X t (e 1) = t max x 1 )+ (1 s(x)dx, e x )+ s(x)dx + o(t). 2 Tankov (211): Leading term for a pure-jump Lévy process with stable-like small-jump behavior with α > 1: ( ) + E e X t 1 = t 1/α E ( Z + ) ( + o t 1/α), (t ), where Z is the strictly α-stable process associated with X; 3 If C(1) = C( 1) =: C, 1 d 1 := E( Z + 1 ) = 1 π Γ (1 1/α) (2CΓ( α) cos (πα/2) )1/α. 4 Tankov (211): If L = σw + X with non-zero σ, ( ) + E e σw t +X t 1 = t 1/2 σ + o(t 1/2 ), (t ). 2π

14 Tempered-stable-like Lévy with and without Brownian component Second-Order Expansions Important question: What is the accuracy of the previous asymptotics? F-L, Gong, and Houdré (212): For a tempered-stable-like Lévy process X with α (1, 2), Π t = E ( e Xt 1 ) + = d1 t 1 α + d2 t + o(t), (t ), where d 1 = E( Z + 1 ) and d 2 = C(1)P( Z 1 < ) C( 1)P( Z 1 ) (e x q(x) q(x) x) x α 1 dx (e x q(x) q(x) x) x α 1 dx.

15 Tempered-stable-like Lévy with and without Brownian component Second-Order Expansions Important question: What is the accuracy of the previous asymptotics? F-L, Gong, and Houdré (212): For a tempered-stable-like Lévy process X with α (1, 2), Π t = E ( e Xt 1 ) + = d1 t 1 α + d2 t + o(t), (t ), where d 1 = E( Z + 1 ) and d 2 = C(1)P( Z 1 < ) C( 1)P( Z 1 ) (e x q(x) q(x) x) x α 1 dx (e x q(x) q(x) x) x α 1 dx.

16 Tempered-stable-like Lévy with and without Brownian component Second-Order Expansions Important question: What is the accuracy of the previous asymptotics? F-L, Gong, and Houdré (212): For a tempered-stable-like Lévy process X with α (1, 2), Π t = E ( e Xt 1 ) + = d1 t 1 α + d2 t + o(t), (t ), where d 1 = E( Z + 1 ) and d 2 = C(1)P( Z 1 < ) C( 1)P( Z 1 ) (e x q(x) q(x) x) x α 1 dx (e x q(x) q(x) x) x α 1 dx.

17 Tempered-stable-like Lévy with and without Brownian component Second-Order Expansions Important question: What is the accuracy of the previous asymptotics? F-L, Gong, and Houdré (212): For a tempered-stable-like Lévy process X with α (1, 2), Π t = E ( e Xt 1 ) + = d1 t 1 α + d2 t + o(t), (t ), where d 1 = E( Z + 1 ) and d 2 = C(1)P( Z 1 < ) C( 1)P( Z 1 ) (e x q(x) q(x) x) x α 1 dx (e x q(x) q(x) x) x α 1 dx.

18 Tempered-stable-like Lévy with and without Brownian component Second-Order Expansions. Continue... F-L, Gong, and Houdré (212): For a tempered-stable-like process X with α (1, 2) and a nonzero independent Brownian component σw t, where Π t = E ( e Xt +σwt 1 ) + = d1 t d2 t 3 α 2 + o(t 3 α 2 ), (t ), d 1 := σe ( W + ) σ 1 = 2π C(1) + C( 1) d 2 := σ 1 α E ( W 1 1 α) 2α(α 1) ( = 21 α Γ 1 α ) (C(1) + C( 1)) σ 1 α. π 2 2α(α 1)

19 Tempered-stable-like Lévy with and without Brownian component Second-Order Expansions. Continue... F-L, Gong, and Houdré (212): For a tempered-stable-like process X with α (1, 2) and a nonzero independent Brownian component σw t, where Π t = E ( e Xt +σwt 1 ) + = d1 t d2 t 3 α 2 + o(t 3 α 2 ), (t ), d 1 := σe ( W + ) σ 1 = 2π C(1) + C( 1) d 2 := σ 1 α E ( W 1 1 α) 2α(α 1) ( = 21 α Γ 1 α ) (C(1) + C( 1)) σ 1 α. π 2 2α(α 1)

20 Tempered-stable-like Lévy with and without Brownian component CGMY or Kobol Model 1 Lévy density: s(x) = C x Y 1 ( e x G 1x< + e x M 1x> ), G, M, C >, α (, 2); 2 Pure-jump Lévy model (σ = ): d 1 = E( Z + 1 ) = 1 ( π Γ 1 1 ) ( ( ) ) 2CΓ( Y ) πy 1/Y Y cos 2 d 2 = CΓ( Y ) 2 3 Mixed Lévy model (σ ): d 1 = σ 2π, ( (M 1) Y M Y (G + 1) Y + G Y ), ( d 2 := 21 Y Γ 1 Y ) Cσ 1 Y π 2 Y (Y 1).

21 Tempered-stable-like Lévy with and without Brownian component CGMY Model s(x) = C x α 1 ( e x /G 1 x< + e x /M 1 x> ), G, M, C > ATM Call Option Prices Pure Jump CGMY Model (C=.5,G=2,M=3.6,Y=1.5) IFT based Method MC based Method 1st order Approx. 2nd order Approx. ATM Call Option Prices General CGMY Model (σ=.4,c=.5,g=2,m=3.6,y= IFT based Method MC based Method 1st order Approx. 2nd order Approx Time to maturity, T (in years) x Time to maturity, T (in years) x 1 3 Figure: Comparisons of ATM call option prices computed by two methods (Inverse Fourier Transform and Monte-Carlo method) with the first- and second-order approximations.

22 Tempered-stable-like Lévy with and without Brownian component Additional Technical Conditions ( The Lévy density s(x) := C x x ) q(x) x α 1 is such that (i) α (1, 2), (ii) C(±1) (, ); (iii) q(x) x 1; (v) 1 q(x) lim =: ±q ( ± ) with q ( ± ) (, ); x ± x (vi) q(x) e x, for all x > ; (vii) q(x) 1, for all x < ; (viii) lim sup x ln q(x) x < ; Natural question: How necessary are these conditions? (iv) sup q(x) < ; x

23 Tempered-stable-like Lévy with and without Brownian component Additional Technical Conditions ( The Lévy density s(x) := C x x ) q(x) x α 1 is such that (i) α (1, 2), (ii) C(±1) (, ); (iii) q(x) x 1; (v) 1 q(x) lim =: ±q ( ± ) with q ( ± ) (, ); x ± x (vi) q(x) e x, for all x > ; (vii) q(x) 1, for all x < ; (viii) lim sup x ln q(x) x < ; Natural question: How necessary are these conditions? (iv) sup q(x) < ; x

24 Tempered-stable-like Lévy with and without Brownian component Additional Technical Conditions ( The Lévy density s(x) := C x x ) q(x) x α 1 is such that (i) α (1, 2), (ii) C(±1) (, ); (iii) q(x) x 1; (v) 1 q(x) lim =: ±q ( ± ) with q ( ± ) (, ); x ± x (vi) q(x) e x, for all x > ; (vii) q(x) 1, for all x < ; (viii) lim sup x ln q(x) x < ; Natural question: How necessary are these conditions? (iv) sup q(x) < ; x

25 Tempered-stable-like Lévy with and without Brownian component Additional Technical Conditions ( The Lévy density s(x) := C x x ) q(x) x α 1 is such that (i) α (1, 2), (ii) C(±1) (, ); (iii) q(x) x 1; (v) 1 q(x) lim =: ±q ( ± ) with q ( ± ) (, ); x ± x (vi) q(x) e x, for all x > ; (vii) q(x) 1, for all x < ; (viii) lim sup x ln q(x) x < ; Natural question: How necessary are these conditions? (iv) sup q(x) < ; x

26 Tempered-stable-like Lévy with and without Brownian component Additional Technical Conditions ( The Lévy density s(x) := C x x ) q(x) x α 1 is such that (i) α (1, 2), (ii) C(±1) (, ); (iii) q(x) x 1; (v) 1 q(x) lim =: ±q ( ± ) with q ( ± ) (, ); x ± x (vi) q(x) e x, for all x > ; (vii) q(x) 1, for all x < ; (viii) lim sup x ln q(x) x < ; Natural question: How necessary are these conditions? (iv) sup q(x) < ; x

27 Stochastic volatility model with tempered-stable-like jumps Specific Problems 1 Investigate the necessity of the technical conditions above. 2 The most liquid options are near" ATM: Π t (κ) := E [ ( e Lt e κ) + ], with κ ; Two approaches: Approximating Π t(κ) with Π t(); Consider another asymptotic regime of the form Π t(κ t) with κ t ; 3 Extend the result to handle a stochastic volatility component: Π t (κ t ) := E [ ( e Lt e κt ) + ], where dl t = dv t + dx t = µ(y t )dt + σ(y t )dw t + dx t, with dy t = α(y t )dt + γ(y t )dw t, Y = y;

28 Stochastic volatility model with tempered-stable-like jumps Specific Problems 1 Investigate the necessity of the technical conditions above. 2 The most liquid options are near" ATM: Π t (κ) := E [ ( e Lt e κ) + ], with κ ; Two approaches: Approximating Π t(κ) with Π t(); Consider another asymptotic regime of the form Π t(κ t) with κ t ; 3 Extend the result to handle a stochastic volatility component: Π t (κ t ) := E [ ( e Lt e κt ) + ], where dl t = dv t + dx t = µ(y t )dt + σ(y t )dw t + dx t, with dy t = α(y t )dt + γ(y t )dw t, Y = y;

29 Stochastic volatility model with tempered-stable-like jumps Specific Problems 1 Investigate the necessity of the technical conditions above. 2 The most liquid options are near" ATM: Π t (κ) := E [ ( e Lt e κ) + ], with κ ; Two approaches: Approximating Π t(κ) with Π t(); Consider another asymptotic regime of the form Π t(κ t) with κ t ; 3 Extend the result to handle a stochastic volatility component: Π t (κ t ) := E [ ( e Lt e κt ) + ], where dl t = dv t + dx t = µ(y t )dt + σ(y t )dw t + dx t, with dy t = α(y t )dt + γ(y t )dw t, Y = y;

30 Stochastic volatility model with tempered-stable-like jumps Specific Problems 1 Investigate the necessity of the technical conditions above. 2 The most liquid options are near" ATM: Π t (κ) := E [ ( e Lt e κ) + ], with κ ; Two approaches: Approximating Π t(κ) with Π t(); Consider another asymptotic regime of the form Π t(κ t) with κ t ; 3 Extend the result to handle a stochastic volatility component: Π t (κ t ) := E [ ( e Lt e κt ) + ], where dl t = dv t + dx t = µ(y t )dt + σ(y t )dw t + dx t, with dy t = α(y t )dt + γ(y t )dw t, Y = y;

31 Stochastic volatility model with tempered-stable-like jumps Specific Problems 1 Investigate the necessity of the technical conditions above. 2 The most liquid options are near" ATM: Π t (κ) := E [ ( e Lt e κ) + ], with κ ; Two approaches: Approximating Π t(κ) with Π t(); Consider another asymptotic regime of the form Π t(κ t) with κ t ; 3 Extend the result to handle a stochastic volatility component: Π t (κ t ) := E [ ( e Lt e κt ) + ], where dl t = dv t + dx t = µ(y t )dt + σ(y t )dw t + dx t, with dy t = α(y t )dt + γ(y t )dw t, Y = y;

32 Stochastic volatility model with tempered-stable-like jumps Pure-jump Case Theorem (F-L & Ólafsson, 213) ( Suppose the Lévy density s(x) := C x x ) q(x) x α 1 is such that (i) α (1, 2), (ii) C(±1) (, ); (iii) q(x) x 1; (v) 1 e x s(x)dx < ; (vi) x 1 x α 1 q(x) dx <. (iv) sup q(x) < ; x

33 Stochastic volatility model with tempered-stable-like jumps Pure-jump Case Theorem (F-L & Ólafsson, 213) ( Suppose the Lévy density s(x) := C x x ) q(x) x α 1 is such that (i) α (1, 2), (ii) C(±1) (, ); (iii) q(x) x 1; (v) 1 e x s(x)dx < ; (vi) x 1 x α 1 q(x) dx <. (iv) sup q(x) < ; x

34 Stochastic volatility model with tempered-stable-like jumps Pure-jump Case Theorem (F-L & Ólafsson, 213) ( Suppose the Lévy density s(x) := C x x ) q(x) x α 1 is such that (i) α (1, 2), (ii) C(±1) (, ); (iii) q(x) x 1; (v) 1 e x s(x)dx < ; (vi) x 1 x α 1 q(x) dx <. (iv) sup q(x) < ; x

35 Stochastic volatility model with tempered-stable-like jumps Pure-jump Case. Continuation... Theorem (F-L & Ólafsson, 213) Then, for κ t := θt + o(t), as t, for some θ R, Π t = E ( e Xt e κt ) + = d1 t 1 α + d2 t + o(t), (t ), where d 1 = E( Z + 1 ) and d 2 = C(1)P( Z 1 < ) C( 1)P( Z 1 ) Remark: The condition (vi) (e x q(x) q(x) x) x α 1 dx x 1 (e x q(x) q(x) x) x α 1 dx θp( Z 1 ). x α 1 q(x) dx <, is a necessary conditions for such a second-order expansion to exist.

36 Stochastic volatility model with tempered-stable-like jumps Nonzero Diffusion Component Theorem (F-L & Ólafsson, 213) Let dv t = µ(y t )dt + σ(y t )dw t, dy t = α(y t )dt + γ(y t )dw t, Y = y ; independent of X such that σ(y ) > and σ( ) is Lipschitz continuous at y. Then, for κ t := θt 3 α 2 + o(t 3 α 2 ), as t, for some θ R, Π t = E ( e Xt +Vt e ) ) κt + = d1 t d2 t 3 α 2 + o (t 3 α 2, (t ), where d 1 := σ(y ) 2π, d 2 := θ 2 + C(1) + C( 1) 2α(α 1) σ(y ) 1 α E ( W 1 1 α).

37 Stochastic volatility model with tempered-stable-like jumps Nonzero Diffusion Component Theorem (F-L & Ólafsson, 213) Let dv t = µ(y t )dt + σ(y t )dw t, dy t = α(y t )dt + γ(y t )dw t, Y = y ; independent of X such that σ(y ) > and σ( ) is Lipschitz continuous at y. Then, for κ t := θt 3 α 2 + o(t 3 α 2 ), as t, for some θ R, Π t = E ( e Xt +Vt e ) ) κt + = d1 t d2 t 3 α 2 + o (t 3 α 2, (t ), where d 1 := σ(y ) 2π, d 2 := θ 2 + C(1) + C( 1) 2α(α 1) σ(y ) 1 α E ( W 1 1 α).

38 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step1. Change of probability measure P P, where dp Ft = e Xt dp Ft : E ( e Xt e ) κt + = e x P (X t > κ t + x)dx (1) = t 1/α e κt = t 1/α e κt t 1/α κ t e t t 1/α κ t e t 1/αu P (t 1/α X t > u)du 1/αu P (t 1/α X t > u)du + t 1/α e κt e t 1/αu P (t 1/α X t > u)du (2) κ t P ( Z1 ) + t 1/α P ( Z 1 > u)du = κ t P ( Z1 ) + t 1/α E ( Z1 ) + (1) u = t 1/α (κ t + x), (2) t 1/α X t D Z1,

39 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step1. Change of probability measure P P, where dp Ft = e Xt dp Ft : E ( e Xt e ) κt + = e x P (X t > κ t + x)dx (1) = t 1/α e κt = t 1/α e κt t 1/α κ t e t t 1/α κ t e t 1/αu P (t 1/α X t > u)du 1/αu P (t 1/α X t > u)du + t 1/α e κt e t 1/αu P (t 1/α X t > u)du (2) κ t P ( Z1 ) + t 1/α P ( Z 1 > u)du = κ t P ( Z1 ) + t 1/α E ( Z1 ) + (1) u = t 1/α (κ t + x), (2) t 1/α X t D Z1,

40 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step1. Change of probability measure P P, where dp Ft = e Xt dp Ft : E ( e Xt e ) κt + = e x P (X t > κ t + x)dx (1) = t 1/α e κt = t 1/α e κt t 1/α κ t e t t 1/α κ t e t 1/αu P (t 1/α X t > u)du 1/αu P (t 1/α X t > u)du + t 1/α e κt e t 1/αu P (t 1/α X t > u)du (2) κ t P ( Z1 ) + t 1/α P ( Z 1 > u)du = κ t P ( Z1 ) + t 1/α E ( Z1 ) + (1) u = t 1/α (κ t + x), (2) t 1/α X t D Z1,

41 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step1. Change of probability measure P P, where dp Ft = e Xt dp Ft : E ( e Xt e ) κt + = e x P (X t > κ t + x)dx (1) = t 1/α e κt = t 1/α e κt t 1/α κ t e t t 1/α κ t e t 1/αu P (t 1/α X t > u)du 1/αu P (t 1/α X t > u)du + t 1/α e κt e t 1/αu P (t 1/α X t > u)du (2) κ t P ( Z1 ) + t 1/α P ( Z 1 > u)du = κ t P ( Z1 ) + t 1/α E ( Z1 ) + (1) u = t 1/α (κ t + x), (2) t 1/α X t D Z1,

42 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step1. Change of probability measure P P, where dp Ft = e Xt dp Ft : E ( e Xt e ) κt + = e x P (X t > κ t + x)dx (1) = t 1/α e κt = t 1/α e κt t 1/α κ t e t t 1/α κ t e t 1/αu P (t 1/α X t > u)du 1/αu P (t 1/α X t > u)du + t 1/α e κt e t 1/αu P (t 1/α X t > u)du (2) κ t P ( Z1 ) + t 1/α P ( Z 1 > u)du = κ t P ( Z1 ) + t 1/α E ( Z1 ) + (1) u = t 1/α (κ t + x), (2) t 1/α X t D Z1,

43 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step1. Change of probability measure P P, where dp Ft = e Xt dp Ft : E ( e Xt e ) κt + = e x P (X t > κ t + x)dx (1) = t 1/α e κt = t 1/α e κt t 1/α κ t e t t 1/α κ t e t 1/αu P (t 1/α X t > u)du 1/αu P (t 1/α X t > u)du + t 1/α e κt e t 1/αu P (t 1/α X t > u)du (2) κ t P ( Z1 ) + t 1/α P ( Z 1 > u)du = κ t P ( Z1 ) + t 1/α E ( Z1 ) + (1) u = t 1/α (κ t + x), (2) t 1/α X t D Z1,

44 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step 2. To study the convergence in e t 1/αu P (t 1/α X t > u)du E ( Z1 ) +, consider another change P P such that, under P, X is a stable process: e t 1/αu P (t 1/α X t > u)du = where d P = e Ut dp Ft. Since Ẽ ( Z ) + 1 = Ft Ũ t = U t ẼU t, t 1/α 1 ( = 1 t 1 t e t 1/α u Ẽ ( e Ut 1 {t 1/α X t >u}) du P(t 1/α Z t u)du and letting e t 1/αu P (t 1/α X t > u)du Ẽ ( Z ) ) + 1 ( e v 1 ) P( Z + t (e v 1) P( Z + t + Ũt v)dv + Ũt v)dv + o(1)

45 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step 2. To study the convergence in e t 1/αu P (t 1/α X t > u)du E ( Z1 ) +, consider another change P P such that, under P, X is a stable process: e t 1/αu P (t 1/α X t > u)du = where d P = e Ut dp Ft. Since Ẽ ( Z ) + 1 = Ft Ũ t = U t ẼU t, t 1/α 1 ( = 1 t 1 t e t 1/α u Ẽ ( e Ut 1 {t 1/α X t >u}) du P(t 1/α Z t u)du and letting e t 1/αu P (t 1/α X t > u)du Ẽ ( Z ) ) + 1 ( e v 1 ) P( Z + t (e v 1) P( Z + t + Ũt v)dv + Ũt v)dv + o(1)

46 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step 2. To study the convergence in e t 1/αu P (t 1/α X t > u)du E ( Z1 ) +, consider another change P P such that, under P, X is a stable process: e t 1/αu P (t 1/α X t > u)du = where d P = e Ut dp Ft. Since Ẽ ( Z ) + 1 = Ft Ũ t = U t ẼU t, t 1/α 1 ( = 1 t 1 t e t 1/α u Ẽ ( e Ut 1 {t 1/α X t >u}) du P(t 1/α Z t u)du and letting e t 1/αu P (t 1/α X t > u)du Ẽ ( Z ) ) + 1 ( e v 1 ) P( Z + t (e v 1) P( Z + t + Ũt v)dv + Ũt v)dv + o(1)

47 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step 2. In terms of the compensated jump measure µ of Z (under P), t t Z t = x µ(ds, dx), Ũ t = ( ln q(x) x) µ(ds, dx), Thus, as t, R\{} t 1/α 1 ( = R\{} e t 1/αu P (t 1/α X t > u)du Ẽ ( Z ) ) + 1 ( e v 1 ) 1 t P( Z + t + Ũt v)dv (e v 1) 1 t P( Z + t + Ũt v)dv ( e v 1 ) ( ) x 1 {x + ln q(x) x v}c x α 1 dxdv x ( ) x (e v 1) 1 {x + ln q(x) x v}c x α 1 dxdv x

48 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step 3. The above approach requires the the conditions: 1 q(x) = O(x), lim sup x ln q(x) x <. To relax those, the idea is to approximate X by a process X (δ) such that X (δ) D X, as δ X X (δ) is of bounded variation It turns out that the condition x 1 x α 1 q(x) dx < is exactly what is needed to do that.

49 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step 3. The above approach requires the the conditions: 1 q(x) = O(x), lim sup x ln q(x) x <. To relax those, the idea is to approximate X by a process X (δ) such that X (δ) D X, as δ X X (δ) is of bounded variation It turns out that the condition x 1 x α 1 q(x) dx < is exactly what is needed to do that. E ) (δ) + X (e t e e κt E X t e E ( e Xt e κt ) + E X (δ) t ) (δ) + X (e t e e κt + E X t e X (δ) t.

50 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step 3. The above approach requires the the conditions: 1 q(x) = O(x), lim sup x ln q(x) x <. To relax those, the idea is to approximate X by a process X (δ) such that X (δ) D X, as δ X X (δ) is of bounded variation It turns out that the condition x 1 x α 1 q(x) dx < is exactly what is needed to do that. t 1/α 1 (t 1/α E ) ) (δ) + X (e t e κt d1 1 e t E X t e t 1/α 1 ( t 1/α E ( e Xt e κt ) + d1 ) X (δ) t ) ) t (t 1/α 1 1/α (δ) + X E (e t e κt d1 + 1 e t E X t e X (δ) t

51 Stochastic volatility model with tempered-stable-like jumps Outline of the proof (pure-jump case) Step 3. Since X X (δ) is of bounded variation, it turns out that 1 e lim sup t t E X t e X (δ) t K 1 q(x) x α dx x <δ for some K <. Making t, d (δ) 2 K 1 q(x) x α dx lim inf t x <δ lim sup t d (δ) 2 + K t 1/α 1 ( t 1/α E ( e Xt e κt ) + d1 ) ( t 1/α 1 t 1/α E ( e Xt e ) ) κt + d1 1 q(x) x α dx x <δ Finally, the results follows from making δ since d (δ) 2 d 2.

52 Conclusions Conclusions 1 Obtained the second-order short-time expansions for near" ATM European call option prices under stable-like small jumps and a possible nonzero independent diffusion component. 2 Characterized explicitly the effects of the different parameters into the behavior of ATM option prices and implied volatility near expiration.

53 Appendix Bibliography For Further Reading I Figueroa-López, J.E., Gong, R., & Houdré, C. High-order short-time expansions for ATM option prices under a tempered stable Lévy model. To appear in Mathematical Finance, Figueroa-López, J.E., & Ólafsson, S. Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility. Arxiv 214. Tankov, P. Pricing and hedging in exponential Lévy models: review of recent results, Paris-Princeton Lecture Notes in Mathematical Finance, Springer 21. Rosenbaum, R., and Tankov, P. Asymptotic results for time-changed Lévy processes sampled at hitting times. Stochastic processes and their applications, 121, 211.

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

Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias José E. Figueroa-López 1 1 Department of Statistics Purdue University Seoul National University & Ajou University

More information

HJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011

HJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011 Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance

More information

Numerical Methods with Lévy Processes

Numerical Methods with Lévy Processes Numerical Methods with Lévy Processes 1 Objective: i) Find models of asset returns, etc ii) Get numbers out of them. Why? VaR and risk management Valuing and hedging derivatives Why not? Usual assumption:

More information

Jump-type Levy Processes

Jump-type Levy Processes Jump-type Levy Processes Ernst Eberlein Handbook of Financial Time Series Outline Table of contents Probabilistic Structure of Levy Processes Levy process Levy-Ito decomposition Jump part Probabilistic

More information

Generalised Fractional-Black-Scholes Equation: pricing and hedging

Generalised Fractional-Black-Scholes Equation: pricing and hedging Generalised Fractional-Black-Scholes Equation: pricing and hedging Álvaro Cartea Birkbeck College, University of London April 2004 Outline Lévy processes Fractional calculus Fractional-Black-Scholes 1

More information

Stochastic Calculus for Finance II - some Solutions to Chapter VII

Stochastic Calculus for Finance II - some Solutions to Chapter VII Stochastic Calculus for Finance II - some Solutions to Chapter VII Matthias hul Last Update: June 9, 25 Exercise 7 Black-Scholes-Merton Equation for the up-and-out Call) i) We have ii) We first compute

More information

Stochastic Volatility and Correction to the Heat Equation

Stochastic Volatility and Correction to the Heat Equation Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Introduction to Malliavin calculus and its applications Lecture 3: Clark-Ocone formula David Nualart Department of Mathematics Kansas University University of Wyoming Summer School 214 David Nualart

More information

Thomas Knispel Leibniz Universität Hannover

Thomas Knispel Leibniz Universität Hannover Optimal long term investment under model ambiguity Optimal long term investment under model ambiguity homas Knispel Leibniz Universität Hannover knispel@stochastik.uni-hannover.de AnStAp0 Vienna, July

More information

Obstacle problems for nonlocal operators

Obstacle problems for nonlocal operators Obstacle problems for nonlocal operators Camelia Pop School of Mathematics, University of Minnesota Fractional PDEs: Theory, Algorithms and Applications ICERM June 19, 2018 Outline Motivation Optimal regularity

More information

Statistical methods for financial models driven by Lévy processes

Statistical methods for financial models driven by Lévy processes Statistical methods for financial models driven by Lévy processes Part IV: Some Non-parametric Methods for Lévy Models José Enrique Figueroa-López 1 1 Department of Statistics Purdue University PASI CIMAT,

More information

Variance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models

Variance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models Variance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models Jean-Pierre Fouque North Carolina State University SAMSI November 3, 5 1 References: Variance Reduction for Monte

More information

Other properties of M M 1

Other properties of M M 1 Other properties of M M 1 Přemysl Bejda premyslbejda@gmail.com 2012 Contents 1 Reflected Lévy Process 2 Time dependent properties of M M 1 3 Waiting times and queue disciplines in M M 1 Contents 1 Reflected

More information

Optimal portfolio strategies under partial information with expert opinions

Optimal portfolio strategies under partial information with expert opinions 1 / 35 Optimal portfolio strategies under partial information with expert opinions Ralf Wunderlich Brandenburg University of Technology Cottbus, Germany Joint work with Rüdiger Frey Research Seminar WU

More information

Optimal investment strategies for an index-linked insurance payment process with stochastic intensity

Optimal investment strategies for an index-linked insurance payment process with stochastic intensity for an index-linked insurance payment process with stochastic intensity Warsaw School of Economics Division of Probabilistic Methods Probability space ( Ω, F, P ) Filtration F = (F(t)) 0 t T satisfies

More information

A new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009

A new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009 A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance

More information

Poisson random measure: motivation

Poisson random measure: motivation : motivation The Lévy measure provides the expected number of jumps by time unit, i.e. in a time interval of the form: [t, t + 1], and of a certain size Example: ν([1, )) is the expected number of jumps

More information

On the quantiles of the Brownian motion and their hitting times.

On the quantiles of the Brownian motion and their hitting times. On the quantiles of the Brownian motion and their hitting times. Angelos Dassios London School of Economics May 23 Abstract The distribution of the α-quantile of a Brownian motion on an interval [, t]

More information

Lecture 4: Introduction to stochastic processes and stochastic calculus

Lecture 4: Introduction to stochastic processes and stochastic calculus Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London

More information

Risk Bounds for Lévy Processes in the PAC-Learning Framework

Risk Bounds for Lévy Processes in the PAC-Learning Framework Risk Bounds for Lévy Processes in the PAC-Learning Framework Chao Zhang School of Computer Engineering anyang Technological University Dacheng Tao School of Computer Engineering anyang Technological University

More information

Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises

Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,

More information

Information and Credit Risk

Information and Credit Risk Information and Credit Risk M. L. Bedini Université de Bretagne Occidentale, Brest - Friedrich Schiller Universität, Jena Jena, March 2011 M. L. Bedini (Université de Bretagne Occidentale, Brest Information

More information

Convergence of price and sensitivities in Carr s randomization approximation globally and near barrier

Convergence of price and sensitivities in Carr s randomization approximation globally and near barrier Convergence of price and sensitivities in Carr s randomization approximation globally and near barrier Sergei Levendorskĭi University of Leicester Toronto, June 23, 2010 Levendorskĭi () Convergence of

More information

Definition: Lévy Process. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 2: Lévy Processes. Theorem

Definition: Lévy Process. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 2: Lévy Processes. Theorem Definition: Lévy Process Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 2: Lévy Processes David Applebaum Probability and Statistics Department, University of Sheffield, UK July

More information

Some Aspects of Universal Portfolio

Some Aspects of Universal Portfolio 1 Some Aspects of Universal Portfolio Tomoyuki Ichiba (UC Santa Barbara) joint work with Marcel Brod (ETH Zurich) Conference on Stochastic Asymptotics & Applications Sixth Western Conference on Mathematical

More information

Estimation for the standard and geometric telegraph process. Stefano M. Iacus. (SAPS VI, Le Mans 21-March-2007)

Estimation for the standard and geometric telegraph process. Stefano M. Iacus. (SAPS VI, Le Mans 21-March-2007) Estimation for the standard and geometric telegraph process Stefano M. Iacus University of Milan(Italy) (SAPS VI, Le Mans 21-March-2007) 1 1. Telegraph process Consider a particle moving on the real line

More information

Lecture 12: Diffusion Processes and Stochastic Differential Equations

Lecture 12: Diffusion Processes and Stochastic Differential Equations Lecture 12: Diffusion Processes and Stochastic Differential Equations 1. Diffusion Processes 1.1 Definition of a diffusion process 1.2 Examples 2. Stochastic Differential Equations SDE) 2.1 Stochastic

More information

Multivariate Generalized Ornstein-Uhlenbeck Processes

Multivariate Generalized Ornstein-Uhlenbeck Processes Multivariate Generalized Ornstein-Uhlenbeck Processes Anita Behme TU München Alexander Lindner TU Braunschweig 7th International Conference on Lévy Processes: Theory and Applications Wroclaw, July 15 19,

More information

Stochastic optimal control with rough paths

Stochastic optimal control with rough paths Stochastic optimal control with rough paths Paul Gassiat TU Berlin Stochastic processes and their statistics in Finance, Okinawa, October 28, 2013 Joint work with Joscha Diehl and Peter Friz Introduction

More information

First passage time for Brownian motion and piecewise linear boundaries

First passage time for Brownian motion and piecewise linear boundaries To appear in Methodology and Computing in Applied Probability, (2017) 19: 237-253. doi 10.1007/s11009-015-9475-2 First passage time for Brownian motion and piecewise linear boundaries Zhiyong Jin 1 and

More information

Gillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde

Gillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation

More information

On a class of stochastic differential equations in a financial network model

On a class of stochastic differential equations in a financial network model 1 On a class of stochastic differential equations in a financial network model Tomoyuki Ichiba Department of Statistics & Applied Probability, Center for Financial Mathematics and Actuarial Research, University

More information

Exercises. T 2T. e ita φ(t)dt.

Exercises. T 2T. e ita φ(t)dt. Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.

More information

The concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.

The concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt. The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes

More information

Optimal Execution Tracking a Benchmark

Optimal Execution Tracking a Benchmark Optimal Execution Tracking a Benchmark René Carmona Bendheim Center for Finance Department of Operations Research & Financial Engineering Princeton University Princeton, June 20, 2013 Optimal Execution

More information

A Change of Variable Formula with Local Time-Space for Bounded Variation Lévy Processes with Application to Solving the American Put Option Problem 1

A Change of Variable Formula with Local Time-Space for Bounded Variation Lévy Processes with Application to Solving the American Put Option Problem 1 Chapter 3 A Change of Variable Formula with Local Time-Space for Bounded Variation Lévy Processes with Application to Solving the American Put Option Problem 1 Abstract We establish a change of variable

More information

Malliavin Calculus in Finance

Malliavin Calculus in Finance Malliavin Calculus in Finance Peter K. Friz 1 Greeks and the logarithmic derivative trick Model an underlying assent by a Markov process with values in R m with dynamics described by the SDE dx t = b(x

More information

Estimation of arrival and service rates for M/M/c queue system

Estimation of arrival and service rates for M/M/c queue system Estimation of arrival and service rates for M/M/c queue system Katarína Starinská starinskak@gmail.com Charles University Faculty of Mathematics and Physics Department of Probability and Mathematical Statistics

More information

LAN property for sde s with additive fractional noise and continuous time observation

LAN property for sde s with additive fractional noise and continuous time observation LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,

More information

UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE

UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (2010), 275 284 UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Iuliana Carmen Bărbăcioru Abstract.

More information

Stochastic differential equation models in biology Susanne Ditlevsen

Stochastic differential equation models in biology Susanne Ditlevsen Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential

More information

Multi-Factor Lévy Models I: Symmetric alpha-stable (SαS) Lévy Processes

Multi-Factor Lévy Models I: Symmetric alpha-stable (SαS) Lévy Processes Multi-Factor Lévy Models I: Symmetric alpha-stable (SαS) Lévy Processes Anatoliy Swishchuk Department of Mathematics and Statistics University of Calgary Calgary, Alberta, Canada Lunch at the Lab Talk

More information

Ernesto Mordecki 1. Lecture III. PASI - Guanajuato - June 2010

Ernesto Mordecki 1. Lecture III. PASI - Guanajuato - June 2010 Optimal stopping for Hunt and Lévy processes Ernesto Mordecki 1 Lecture III. PASI - Guanajuato - June 2010 1Joint work with Paavo Salminen (Åbo, Finland) 1 Plan of the talk 1. Motivation: from Finance

More information

A problem of portfolio/consumption choice in a. liquidity risk model with random trading times

A problem of portfolio/consumption choice in a. liquidity risk model with random trading times A problem of portfolio/consumption choice in a liquidity risk model with random trading times Huyên PHAM Special Semester on Stochastics with Emphasis on Finance, Kick-off workshop, Linz, September 8-12,

More information

A numerical method for solving uncertain differential equations

A numerical method for solving uncertain differential equations Journal of Intelligent & Fuzzy Systems 25 (213 825 832 DOI:1.3233/IFS-12688 IOS Press 825 A numerical method for solving uncertain differential equations Kai Yao a and Xiaowei Chen b, a Department of Mathematical

More information

Affine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009

Affine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009 Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine

More information

Scale functions for spectrally negative Lévy processes and their appearance in economic models

Scale functions for spectrally negative Lévy processes and their appearance in economic models Scale functions for spectrally negative Lévy processes and their appearance in economic models Andreas E. Kyprianou 1 Department of Mathematical Sciences, University of Bath 1 This is a review talk and

More information

Duration-Based Volatility Estimation

Duration-Based Volatility Estimation A Dual Approach to RV Torben G. Andersen, Northwestern University Dobrislav Dobrev, Federal Reserve Board of Governors Ernst Schaumburg, Northwestern Univeristy CHICAGO-ARGONNE INSTITUTE ON COMPUTATIONAL

More information

Question 1. The correct answers are: (a) (2) (b) (1) (c) (2) (d) (3) (e) (2) (f) (1) (g) (2) (h) (1)

Question 1. The correct answers are: (a) (2) (b) (1) (c) (2) (d) (3) (e) (2) (f) (1) (g) (2) (h) (1) Question 1 The correct answers are: a 2 b 1 c 2 d 3 e 2 f 1 g 2 h 1 Question 2 a Any probability measure Q equivalent to P on F 2 can be described by Q[{x 1, x 2 }] := q x1 q x1,x 2, 1 where q x1, q x1,x

More information

Small-time expansions for the transition distributions of Lévy processes

Small-time expansions for the transition distributions of Lévy processes Small-time expansions for the transition distributions of Lévy processes José E. Figueroa-López and Christian Houdré Department of Statistics Purdue University W. Lafayette, IN 4796, USA figueroa@stat.purdue.edu

More information

EQUITY MARKET STABILITY

EQUITY MARKET STABILITY EQUITY MARKET STABILITY Adrian Banner INTECH Investment Technologies LLC, Princeton (Joint work with E. Robert Fernholz, Ioannis Karatzas, Vassilios Papathanakos and Phillip Whitman.) Talk at WCMF6 conference,

More information

Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise

Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise Robert J. Elliott 1 Samuel N. Cohen 2 1 Department of Commerce, University of South Australia 2 Mathematical Insitute, University

More information

Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )

Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( ) Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical

More information

Regular Variation and Extreme Events for Stochastic Processes

Regular Variation and Extreme Events for Stochastic Processes 1 Regular Variation and Extreme Events for Stochastic Processes FILIP LINDSKOG Royal Institute of Technology, Stockholm 2005 based on joint work with Henrik Hult www.math.kth.se/ lindskog 2 Extremes for

More information

Higher order weak approximations of stochastic differential equations with and without jumps

Higher order weak approximations of stochastic differential equations with and without jumps Higher order weak approximations of stochastic differential equations with and without jumps Hideyuki TANAKA Graduate School of Science and Engineering, Ritsumeikan University Rough Path Analysis and Related

More information

Utility Maximization in Hidden Regime-Switching Markets with Default Risk

Utility Maximization in Hidden Regime-Switching Markets with Default Risk Utility Maximization in Hidden Regime-Switching Markets with Default Risk José E. Figueroa-López Department of Mathematics and Statistics Washington University in St. Louis figueroa-lopez@wustl.edu pages.wustl.edu/figueroa

More information

Random Times and Their Properties

Random Times and Their Properties Chapter 6 Random Times and Their Properties Section 6.1 recalls the definition of a filtration (a growing collection of σ-fields) and of stopping times (basically, measurable random times). Section 6.2

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 7 9/25/2013

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 7 9/25/2013 MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.070J Fall 013 Lecture 7 9/5/013 The Reflection Principle. The Distribution of the Maximum. Brownian motion with drift Content. 1. Quick intro to stopping times.

More information

Multilevel Monte Carlo for Lévy Driven SDEs

Multilevel Monte Carlo for Lévy Driven SDEs Multilevel Monte Carlo for Lévy Driven SDEs Felix Heidenreich TU Kaiserslautern AG Computational Stochastics August 2011 joint work with Steffen Dereich Philipps-Universität Marburg supported within DFG-SPP

More information

1 Brownian Local Time

1 Brownian Local Time 1 Brownian Local Time We first begin by defining the space and variables for Brownian local time. Let W t be a standard 1-D Wiener process. We know that for the set, {t : W t = } P (µ{t : W t = } = ) =

More information

in Bounded Domains Ariane Trescases CMLA, ENS Cachan

in Bounded Domains Ariane Trescases CMLA, ENS Cachan CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline

More information

Annealed Brownian motion in a heavy tailed Poissonian potential

Annealed Brownian motion in a heavy tailed Poissonian potential Annealed Brownian motion in a heavy tailed Poissonian potential Ryoki Fukushima Research Institute of Mathematical Sciences Stochastic Analysis and Applications, Okayama University, September 26, 2012

More information

Infinitely divisible distributions and the Lévy-Khintchine formula

Infinitely divisible distributions and the Lévy-Khintchine formula Infinitely divisible distributions and the Cornell University May 1, 2015 Some definitions Let X be a real-valued random variable with law µ X. Recall that X is said to be infinitely divisible if for every

More information

Stochastic Processes and Advanced Mathematical Finance

Stochastic Processes and Advanced Mathematical Finance Steven R. Dunbar Department of Mathematics 23 Avery Hall University of Nebraska-Lincoln Lincoln, NE 68588-13 http://www.math.unl.edu Voice: 42-472-3731 Fax: 42-472-8466 Stochastic Processes and Advanced

More information

Point Process Control

Point Process Control Point Process Control The following note is based on Chapters I, II and VII in Brémaud s book Point Processes and Queues (1981). 1 Basic Definitions Consider some probability space (Ω, F, P). A real-valued

More information

On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem

On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem Koichiro TAKAOKA Dept of Applied Physics, Tokyo Institute of Technology Abstract M Yor constructed a family

More information

Stochastic Calculus. Kevin Sinclair. August 2, 2016

Stochastic Calculus. Kevin Sinclair. August 2, 2016 Stochastic Calculus Kevin Sinclair August, 16 1 Background Suppose we have a Brownian motion W. This is a process, and the value of W at a particular time T (which we write W T ) is a normally distributed

More information

Outline. A Central Limit Theorem for Truncating Stochastic Algorithms

Outline. A Central Limit Theorem for Truncating Stochastic Algorithms Outline A Central Limit Theorem for Truncating Stochastic Algorithms Jérôme Lelong http://cermics.enpc.fr/ lelong Tuesday September 5, 6 1 3 4 Jérôme Lelong (CERMICS) Tuesday September 5, 6 1 / 3 Jérôme

More information

Sequential Monte Carlo Methods for Bayesian Computation

Sequential Monte Carlo Methods for Bayesian Computation Sequential Monte Carlo Methods for Bayesian Computation A. Doucet Kyoto Sept. 2012 A. Doucet (MLSS Sept. 2012) Sept. 2012 1 / 136 Motivating Example 1: Generic Bayesian Model Let X be a vector parameter

More information

Statistical test for some multistable processes

Statistical test for some multistable processes Statistical test for some multistable processes Ronan Le Guével Joint work in progress with A. Philippe Journées MAS 2014 1 Multistable processes First definition : Ferguson-Klass-LePage series Properties

More information

Exact Simulation of Multivariate Itô Diffusions

Exact Simulation of Multivariate Itô Diffusions Exact Simulation of Multivariate Itô Diffusions Jose Blanchet Joint work with Fan Zhang Columbia and Stanford July 7, 2017 Jose Blanchet (Columbia/Stanford) Exact Simulation of Diffusions July 7, 2017

More information

I forgot to mention last time: in the Ito formula for two standard processes, putting

I forgot to mention last time: in the Ito formula for two standard processes, putting I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy

More information

SMALL-TIME EXPANSIONS FOR LOCAL JUMP-DIFFUSIONS MODELS WITH INFINITE JUMP ACTIVITY. 1. Introduction

SMALL-TIME EXPANSIONS FOR LOCAL JUMP-DIFFUSIONS MODELS WITH INFINITE JUMP ACTIVITY. 1. Introduction SMALL-TIME EXPANSIONS FOR LOCAL JUMP-DIFFUSIONS MODELS WITH INFINITE JUMP ACTIVITY JOSÉ E. FIGUEROA-LÓPEZ AND CHENG OUYANG Abstract. We consider a Markov process {X (x) t } t with initial condition X (x)

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

An introduction to Lévy processes

An introduction to Lévy processes with financial modelling in mind Department of Statistics, University of Oxford 27 May 2008 1 Motivation 2 3 General modelling with Lévy processes Modelling financial price processes Quadratic variation

More information

Inference for Lévy-Driven Continuous-Time ARMA Processes

Inference for Lévy-Driven Continuous-Time ARMA Processes Inference for Lévy-Driven Continuous-Time ARMA Processes Peter J. Brockwell Richard A. Davis Yu Yang Colorado State University May 23, 2007 Outline Background Lévy-driven CARMA processes Second order properties

More information

LAN property for ergodic jump-diffusion processes with discrete observations

LAN property for ergodic jump-diffusion processes with discrete observations LAN property for ergodic jump-diffusion processes with discrete observations Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Arturo Kohatsu-Higa (Ritsumeikan University, Japan) &

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

Effective dynamics for the (overdamped) Langevin equation

Effective dynamics for the (overdamped) Langevin equation Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath

More information

FE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology

FE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology FE610 Stochastic Calculus for Financial Engineers Lecture 3. Calculaus in Deterministic and Stochastic Environments Steve Yang Stevens Institute of Technology 01/31/2012 Outline 1 Modeling Random Behavior

More information

CIMPA SCHOOL, 2007 Jump Processes and Applications to Finance Monique Jeanblanc

CIMPA SCHOOL, 2007 Jump Processes and Applications to Finance Monique Jeanblanc CIMPA SCHOOL, 27 Jump Processes and Applications to Finance Monique Jeanblanc 1 Jump Processes I. Poisson Processes II. Lévy Processes III. Jump-Diffusion Processes IV. Point Processes 2 I. Poisson Processes

More information

A Class of Fractional Stochastic Differential Equations

A Class of Fractional Stochastic Differential Equations Vietnam Journal of Mathematics 36:38) 71 79 Vietnam Journal of MATHEMATICS VAST 8 A Class of Fractional Stochastic Differential Equations Nguyen Tien Dung Department of Mathematics, Vietnam National University,

More information

Stochastic Processes

Stochastic Processes Stochastic Processes A very simple introduction Péter Medvegyev 2009, January Medvegyev (CEU) Stochastic Processes 2009, January 1 / 54 Summary from measure theory De nition (X, A) is a measurable space

More information

Continuous Time Finance

Continuous Time Finance Continuous Time Finance Lisbon 2013 Tomas Björk Stockholm School of Economics Tomas Björk, 2013 Contents Stochastic Calculus (Ch 4-5). Black-Scholes (Ch 6-7. Completeness and hedging (Ch 8-9. The martingale

More information

Pathwise volatility in a long-memory pricing model: estimation and asymptotic behavior

Pathwise volatility in a long-memory pricing model: estimation and asymptotic behavior Pathwise volatility in a long-memory pricing model: estimation and asymptotic behavior Ehsan Azmoodeh University of Vaasa Finland 7th General AMaMeF and Swissquote Conference September 7 1, 215 Outline

More information

Exact and high order discretization schemes. Wishart processes and their affine extensions

Exact and high order discretization schemes. Wishart processes and their affine extensions for Wishart processes and their affine extensions CERMICS, Ecole des Ponts Paris Tech - PRES Université Paris Est Modeling and Managing Financial Risks -2011 Plan 1 Motivation and notations 2 Splitting

More information

Squared Bessel Process with Delay

Squared Bessel Process with Delay Southern Illinois University Carbondale OpenSIUC Articles and Preprints Department of Mathematics 216 Squared Bessel Process with Delay Harry Randolph Hughes Southern Illinois University Carbondale, hrhughes@siu.edu

More information

QMC methods in quantitative finance. and perspectives

QMC methods in quantitative finance. and perspectives , tradition and perspectives Johannes Kepler University Linz (JKU) WU Research Seminar What is the content of the talk Valuation of financial derivatives in stochastic market models using (QMC-)simulation

More information

Simulation methods for stochastic models in chemistry

Simulation methods for stochastic models in chemistry Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation

More information

Hardy-Stein identity and Square functions

Hardy-Stein identity and Square functions Hardy-Stein identity and Square functions Daesung Kim (joint work with Rodrigo Bañuelos) Department of Mathematics Purdue University March 28, 217 Daesung Kim (Purdue) Hardy-Stein identity UIUC 217 1 /

More information

LogFeller et Ray Knight

LogFeller et Ray Knight LogFeller et Ray Knight Etienne Pardoux joint work with V. Le and A. Wakolbinger Etienne Pardoux (Marseille) MANEGE, 18/1/1 1 / 16 Feller s branching diffusion with logistic growth We consider the diffusion

More information

Nonlinear representation, backward SDEs, and application to the Principal-Agent problem

Nonlinear representation, backward SDEs, and application to the Principal-Agent problem Nonlinear representation, backward SDEs, and application to the Principal-Agent problem Ecole Polytechnique, France April 4, 218 Outline The Principal-Agent problem Formulation 1 The Principal-Agent problem

More information

Discretization of SDEs: Euler Methods and Beyond

Discretization of SDEs: Euler Methods and Beyond Discretization of SDEs: Euler Methods and Beyond 09-26-2006 / PRisMa 2006 Workshop Outline Introduction 1 Introduction Motivation Stochastic Differential Equations 2 The Time Discretization of SDEs Monte-Carlo

More information

An Uncertain Control Model with Application to. Production-Inventory System

An Uncertain Control Model with Application to. Production-Inventory System An Uncertain Control Model with Application to Production-Inventory System Kai Yao 1, Zhongfeng Qin 2 1 Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China 2 School of Economics

More information

University Of Calgary Department of Mathematics and Statistics

University Of Calgary Department of Mathematics and Statistics University Of Calgary Department of Mathematics and Statistics Hawkes Seminar May 23, 2018 1 / 46 Some Problems in Insurance and Reinsurance Mohamed Badaoui Department of Electrical Engineering National

More information

Stationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos

Stationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos N Cufaro Petroni: CYCLOTRONS 2007 Giardini Naxos, 1 5 October, 2007 1 Stationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos CYCLOTRONS 2007 Giardini Naxos, 1 5 October Nicola

More information

Contagious default: application of methods of Statistical Mechanics in Finance

Contagious default: application of methods of Statistical Mechanics in Finance Contagious default: application of methods of Statistical Mechanics in Finance Wolfgang J. Runggaldier University of Padova, Italy www.math.unipd.it/runggaldier based on joint work with : Paolo Dai Pra,

More information

Jump-diffusion models driven by Lévy processes

Jump-diffusion models driven by Lévy processes Jump-diffusion models driven by Lévy processes José E. Figueroa-López Purdue University Department of Statistics West Lafayette, IN 4797-266 figueroa@stat.purdue.edu Abstract: During the past and this

More information

Branching Processes II: Convergence of critical branching to Feller s CSB

Branching Processes II: Convergence of critical branching to Feller s CSB Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied

More information

Large Deviations Principles for McKean-Vlasov SDEs, Skeletons, Supports and the law of iterated logarithm

Large Deviations Principles for McKean-Vlasov SDEs, Skeletons, Supports and the law of iterated logarithm Large Deviations Principles for McKean-Vlasov SDEs, Skeletons, Supports and the law of iterated logarithm Gonçalo dos Reis University of Edinburgh (UK) & CMA/FCT/UNL (PT) jointly with: W. Salkeld, U. of

More information