Multivariate Generalized Ornstein-Uhlenbeck Processes

Size: px
Start display at page:

Download "Multivariate Generalized Ornstein-Uhlenbeck Processes"

Transcription

1 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, 2013

2 The Ornstein-Uhlenbeck process : Origin The Ornstein-Uhlenbeck process: Origin Einstein (1905) models the movement of a free particle in fluid by Brownian Motion. Alexander Lindner, 2

3 The Ornstein-Uhlenbeck process : Origin The Ornstein-Uhlenbeck process: Origin Einstein (1905) models the movement of a free particle in fluid by Brownian Motion. Ornstein and Uhlenbeck (1930) add the concept of friction to Einsteins model. Alexander Lindner, 2

4 The Ornstein-Uhlenbeck process : Origin The Ornstein-Uhlenbeck process: Origin Einstein (1905) models the movement of a free particle in fluid by Brownian Motion. Ornstein and Uhlenbeck (1930) add the concept of friction to Einsteins model: A particle moving from left to right gets hit by more particles from the right than from the left side which results in a slowdown. The velocity v(t) of the particle is given by m dv(t) = λv(t)dt + db(t) i.e. v(t) = e λt/m v(0) + e λt/m (0,t] e λs/m db(s). Alexander Lindner, 2

5 The Ornstein-Uhlenbeck process : Origin The Ornstein-Uhlenbeck process: Origin Einstein (1905) models the movement of a free particle in fluid by Brownian Motion. Ornstein and Uhlenbeck (1930) add the concept of friction to Einsteins model: A particle moving from left to right gets hit by more particles from the right than from the left side which results in a slowdown. The velocity v(t) of the particle is given by i.e. m dv(t) = λv(t)dt + db(t) v(t) = e λt/m v(0) + e λt/m (0,t] e λs/m db(s). This solution is called an Ornstein-Uhlenbeck (OU) process. Setting λ = 0 yields the original formula by Einstein. Alexander Lindner, 2

6 The Ornstein-Uhlenbeck process : From OU to GOU processes OU processes as AR(1) time series For every h > 0 the Ornstein-Uhlenbeck process V t = e λt V 0 + e λt e λs db s, t 0, (0,t] fulfills the random recurrence equation V nh = e λh V (n 1)h + e λnh ((n 1)h,nh] e λs db s, n N. Hence it can be seen as a natural generalization in continuous time of the AR(1) time series X n = e λ X n 1 + Z n, n N, with i.i.d. noise (Z n ) n N such that L(Z 1 ) = L( (0,1] e λ(1 s) db s ). Alexander Lindner, 3

7 The Ornstein-Uhlenbeck process : From OU to GOU processes A more general AR(1) time series By embedding the more general random sequence Y n = A n Y n 1 + B n, n N, with (A n, B n ) n N i.i.d., A 1 > 0 a.s., into a continuous time setting in 1989 De Haan and Karandikar introduced the generalized Ornstein-Uhlenbeck process ) V t = e (V ξt 0 + e ξ s dη s, t 0. (0,t] driven by a bivariate Lévy process (ξ t, η t ) t 0 with starting random variable V 0. Alexander Lindner, 4

8 Definition: Lévy processes The Ornstein-Uhlenbeck process : From OU to GOU processes A Lévy process in R d on a probability space (Ω, F, P) is a stochastic process X = (X t ) t 0, X t : Ω R d satisfying the following properties: X 0 = 0 a.s. X has independent increments, i.e. for all 0 t 0 t 1... t n the random variables X t0, X t1 X t0,..., X tn X tn 1 are independent. X has stationary increments, i.e. for all s, t 0 it holds X s+t X s d = Xt. X has a.s. càdlàg paths, i.e. for P-a.e. ω Ω the path t X t (ω) is right-continuous in t 0 and has left limits in t > 0. Alexander Lindner, 5

9 Definition: Lévy processes The Ornstein-Uhlenbeck process : From OU to GOU processes A Lévy process in R d on a probability space (Ω, F, P) is a stochastic process X = (X t ) t 0, X t : Ω R d satisfying the following properties: X 0 = 0 a.s. X has independent increments, i.e. for all 0 t 0 t 1... t n the random variables X t0, X t1 X t0,..., X tn X tn 1 are independent. X has stationary increments, i.e. for all s, t 0 it holds X s+t X s d = Xt. X has a.s. càdlàg paths, i.e. for P-a.e. ω Ω the path t X t (ω) is right-continuous in t 0 and has left limits in t > 0. Elementary examples of Lévy processes include linear deterministic processes, Brownian motions as well as compound Poisson processes. Alexander Lindner, 5

10 The Ornstein-Uhlenbeck process : The GOU process Definition: Generalized OU processes The generalized Ornstein-Uhlenbeck (GOU) process (V t ) t 0 driven by the bivariate Lévy process (ξ t, η t ) t 0 is given by ) V t = e (V ξt 0 + e ξ s dη s, t 0, (0,t] where V 0 is a finite random variable, usually chosen independent of (ξ, η). In the case that (ξ t, η t ) = (λt, η t ) with a Lévy process (η t ) t 0 and a constant λ 0 the process (V t ) t 0 is called Lévy-driven Ornstein-Uhlenbeck process or Ornstein-Uhlenbeck type process. Obviously, if additionally (η t ) t 0 is a Brownian motion, we get the classical Ornstein-Uhlenbeck process. Alexander Lindner, 6

11 The Ornstein-Uhlenbeck process : The GOU process The generalized Ornstein-Uhlenbeck process: Applications Example 1: ξ t = t deterministic ) V t = e (V t 0 + e s dη s (0,t] Lévy driven Ornstein-Uhlenbeck process, classical for η t = B t applications in storage theory stochastic volatility model of Barndorff-Nielsen and Shephard (2001): η t subordinator, V t squared volatility, price process G t defined by dg t = (µ + bv t )dt + V t db t for constants µ, b. Alexander Lindner, 7

12 The Ornstein-Uhlenbeck process : The GOU process The generalized Ornstein-Uhlenbeck process: Applications Example 1: ξ t = t deterministic ) V t = e (V t 0 + e s dη s (0,t] Lévy driven Ornstein-Uhlenbeck process, classical for η t = B t applications in storage theory stochastic volatility model of Barndorff-Nielsen and Shephard (2001): η t subordinator, V t squared volatility, price process G t defined by dg t = (µ + bv t )dt + V t db t for constants µ, b. Example 2: η t = t deterministic. Applications for Asian options, or COGARCH(1,1) model of Klüppelberg, L., Maller (2004). Alexander Lindner, 7

13 The corresponding SDE The Ornstein-Uhlenbeck process : The corresponding SDE The generalized Ornstein-Uhlenbeck process driven by (ξ, η) ) V t = e (V ξt 0 + e ξ s dη s, t 0 is the unique solution of the SDE (0,t] dv t = V t du t + dl t, t 0 where (U t, L t ) t 0 is a bivariate Lévy process completely determined by (ξ, η). Alexander Lindner, 8

14 The corresponding SDE The Ornstein-Uhlenbeck process : The corresponding SDE The generalized Ornstein-Uhlenbeck process driven by (ξ, η) ) V t = e (V ξt 0 + e ξ s dη s, t 0 is the unique solution of the SDE (0,t] dv t = V t du t + dl t, t 0 where (U t, L t ) t 0 is a bivariate Lévy process completely determined by (ξ, η). In particular we have ξ t = log(e(u) t ) Definition: The ( Doléans-Dade Exponential E(U) t = exp U t 1 ) 2 tσ2 U ((1 + U s ) exp( U s )) is the unique solution of the SDE dz t = Z t du t, Z 0 = 1 a.s. s t Alexander Lindner, 8

15 The corresponding SDE The Ornstein-Uhlenbeck process : The corresponding SDE The generalized Ornstein-Uhlenbeck process driven by (ξ, η) ) V t = e (V ξt 0 + e ξ s dη s, t 0 is the unique solution of the SDE (0,t] dv t = V t du t + dl t, t 0 where (U t, L t ) t 0 is a bivariate Lévy process completely determined by (ξ, η). In particular we have ξ t = log(e(u) t ) and η t = L t + 0<s t U s L s 1 + U s tcov (B U1, B L1 ). Alexander Lindner, 8

16 Multivariate generalized Ornstein-Uhlenbeck processes : Multivariate Generalized Ornstein-Uhlenbeck Processes Alexander Lindner, 9

17 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Recall: An AR(1) time series The generalized Ornstein-Uhlenbeck process ) V t = e (V ξt 0 + e ξ s dη s, t 0. (0,t] driven by a bivariate Lévy process (ξ t, η t ) t 0 with starting random variable V 0 had been derived by embedding the AR(1) time series V n = A n V n 1 + B n, n N, with (A n, B n ) n N i.i.d., A 1 > 0 a.s., into a continuous time setting. Alexander Lindner, 10

18 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Constructing a multivariate GOU We aim to embed the random sequence V n = A n V n 1 + B n, n N, with (A n, B n ) n N i.i.d., (A n, B n ) R m m R m, A 1 a.s. non-singular, into a continuous time setting. Alexander Lindner, 11

19 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Constructing a multivariate GOU We aim to embed the random sequence V n = A n V n 1 + B n, n N, with (A n, B n ) n N i.i.d., (A n, B n ) R m m R m, A 1 a.s. non-singular, into a continuous time setting. More precisely, we want to find all stochastic processes (V t ) t 0 such that V t = A s,t V s + B s,t, 0 s t and is i.i.d. for each h > 0. (A (n 1)h,nh, B (n 1)h,nh ) n N Alexander Lindner, 11

20 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Constructing a multivariate GOU We aim to embed the random sequence V n = A n V n 1 + B n, n N, with (A n, B n ) n N i.i.d., (A n, B n ) R m m R m, A 1 a.s. non-singular, into a continuous time setting. More precisely, we want to find all stochastic processes (V t ) t 0 such that V t = A s,t V s + B s,t, 0 s t and (A (n 1)h,nh, B (n 1)h,nh ) n N is i.i.d. for each h > 0. Assuming slightly more leads to the following requirements: Alexander Lindner, 11

21 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Assumptions For each 0 s t let (A s,t, B s,t ) GL(R, m) R m s.t.: Assumption (a) For all 0 u s t almost surely A u,t = A s,t A u,s and B u,t = A s,t B u,s + B s,t. Alexander Lindner, 12

22 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Assumptions For each 0 s t let (A s,t, B s,t ) GL(R, m) R m s.t.: Assumption (a) For all 0 u s t almost surely A u,t = A s,t A u,s and B u,t = A s,t B u,s + B s,t. Assumption (b) The families of random matrices {(A s,t, B s,t ), a s t b} and {(A s,t, B s,t ), c s t d} are independent for 0 a b c d. Alexander Lindner, 12

23 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Assumptions For each 0 s t let (A s,t, B s,t ) GL(R, m) R m s.t.: Assumption (a) For all 0 u s t almost surely A u,t = A s,t A u,s and B u,t = A s,t B u,s + B s,t. Assumption (b) The families of random matrices {(A s,t, B s,t ), a s t b} and {(A s,t, B s,t ), c s t d} are independent for 0 a b c d. Assumption (c) For all 0 s t it holds (A s,t, B s,t ) d = (A 0,t s, B 0,t s ). Alexander Lindner, 12

24 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Assumptions For each 0 s t let (A s,t, B s,t ) GL(R, m) R m s.t.: Assumption (a) For all 0 u s t almost surely A u,t = A s,t A u,s and B u,t = A s,t B u,s + B s,t. Assumption (b) The families of random matrices {(A s,t, B s,t ), a s t b} and {(A s,t, B s,t ), c s t d} are independent for 0 a b c d. Assumption (c) For all 0 s t it holds Assumption (d) It holds (A s,t, B s,t ) d = (A 0,t s, B 0,t s ). P lim t 0 A 0,t = A 0,0 = I and P lim t 0 B 0,t = B 0,0 = 0, where I denotes the identity matrix and 0 the vector (or matrix) only having zero entries. Alexander Lindner, 12

25 Multivariate generalized Ornstein-Uhlenbeck processes : Construction The Process A t := A 0,t Lemma: Every stochastic process A t := A 0,t in the autoregressive model above which fulfills Assumptions (a) to (d) has a version which is a multiplicative right Lévy process in the general linear group GL(R, m) of order m. Alexander Lindner, 13

26 Multivariate generalized Ornstein-Uhlenbeck processes : Construction The Process A t := A 0,t Lemma: Every stochastic process A t := A 0,t in the autoregressive model above which fulfills Assumptions (a) to (d) has a version which is a multiplicative right Lévy process in the general linear group GL(R, m) of order m. That means, (A t ) t 0 is a stochastic process with values in GL(R, m) with the following properties: A 0 = I a.s. it has independent left increments, i.e. for all 0 t 1... t n, the random variables A 0, A t1 A 1 0,..., A t n A 1 t n 1 are independent. it has stationary left increments, i.e. for all s, t 0 it holds A s+t A 1 d s = A t. it has a.s. càdlàg paths, i.e. for P-a.e. ω Ω the path t A t (ω) is right-continuous in t 0 and has left limits in t > 0. Alexander Lindner, 13

27 Multivariate generalized Ornstein-Uhlenbeck processes : Construction The Multivariate Stochastic Exponential I Lemma: By an observation due to Skorokhod every right Lévy process in (GL(R, m), ) is the right stochastic exponential of a Lévy process in (R m m, +). Alexander Lindner, 14

28 Multivariate generalized Ornstein-Uhlenbeck processes : Construction The Multivariate Stochastic Exponential I Lemma: By an observation due to Skorokhod every right Lévy process in (GL(R, m), ) is the right stochastic exponential of a Lévy process in (R m m, +). Definition: Let (X t ) t 0 be a semimartingale in (R m m, +). Then its left stochastic exponential E (X ) t is defined as the unique R m m -valued, adapted, càdlàg solution of the SDE Z t = I + Z s dx s, t 0, (0,t] while the unique adapted, càdlàg solution of the SDE Z t = I + dx s Z s, t 0, (0,t] will be called right stochastic exponential and denoted by E (X ) t. Alexander Lindner, 14

29 Multivariate generalized Ornstein-Uhlenbeck processes : Construction The Multivariate Stochastic Exponential II Let (X t ) t 0 be a semimartingale in (R m m, +). Then we observe: A stochastic exponential of X is invertible for all t 0 if and only if det(i + X t ) 0 for all t 0. ( ) Alexander Lindner, 15

30 Multivariate generalized Ornstein-Uhlenbeck processes : Construction The Multivariate Stochastic Exponential II Let (X t ) t 0 be a semimartingale in (R m m, +). Then we observe: A stochastic exponential of X is invertible for all t 0 if and only if det(i + X t ) 0 for all t 0. ( ) Suppose (X t ) t 0 fulfills ( ). Then for (U t ) t 0 given by U t := X t +[X, X ] c t + it holds 0<s t ( (I + Xs ) 1 I + X s ), t 0 [ E (X ) t ] 1 = E (U) t, t 0. Alexander Lindner, 15

31 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Choice of A s,t We have Every stochastic process A t := A 0,t in the autoregressive model V t = A s,t V s + B s,t, 0 s t, which fulfills Assumptions (a) to (d) is a multiplicative right Lévy process in GL(R, m). Alexander Lindner, 16

32 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Choice of A s,t We have Every stochastic process A t := A 0,t in the autoregressive model V t = A s,t V s + B s,t, 0 s t, which fulfills Assumptions (a) to (d) is a multiplicative right Lévy process in GL(R, m). Every right Lévy process in (GL(R, m), ) is the right stochastic exponential of a Lévy process in (R m m, +), i.e. we have A t = E (U) t. Alexander Lindner, 16

33 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Choice of A s,t We have Every stochastic process A t := A 0,t in the autoregressive model V t = A s,t V s + B s,t, 0 s t, which fulfills Assumptions (a) to (d) is a multiplicative right Lévy process in GL(R, m). Every right Lévy process in (GL(R, m), ) is the right stochastic exponential of a Lévy process in (R m m, +), i.e. we have A t = E (U) t. There exists another Lévy process X in (R m m, +) such that A t = E (X ) 1, t 0, and the increments A s,t = A t A 1 s of A t take the form A s,t = E (X ) 1 t E (X )s, 0 s t. Alexander Lindner, 16

34 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Choice of (A s,t, B s,t ) 0 s t Theorem: (i) Suppose (X t, Y t ) t 0 to be a Lévy process in (R m m R m, +) such that E (X ) is non-singular. For 0 s t define ( ) As,t := E (X ) 1 t E (X )s B. s,t E (X ) 1 t (s,t] E (X )u dy u Then (A s,t, B s,t ) 0 s t satisfies Assumptions (a) to (d) above and for any starting random variable V 0 the process ( ) V t := E (X ) 1 t V 0 + E (X )s dy s (0,t] satisfies V t = A s,t V s + B s,t, 0 s t. Alexander Lindner, 17

35 Multivariate generalized Ornstein-Uhlenbeck processes : Construction Choice of (A s,t, B s,t ) 0 s t Theorem: (i) Suppose (X t, Y t ) t 0 to be a Lévy process in (R m m R m, +) such that E (X ) is non-singular. For 0 s t define ( ) As,t := E (X ) 1 t E (X )s B. s,t E (X ) 1 t (s,t] E (X )u dy u Then (A s,t, B s,t ) 0 s t satisfies Assumptions (a) to (d) above and for any starting random variable V 0 the process ( ) V t := E (X ) 1 t V 0 + E (X )s dy s (0,t] satisfies V t = A s,t V s + B s,t, 0 s t. (ii) All processes satisfying V t = A s,t V s + B s,t, 0 s t, with (A s,t, B s,t ) 0 s t satisfying Assumptions (a) to (d), can be obtained in this way. Alexander Lindner, 17

36 Multivariate generalized Ornstein-Uhlenbeck processes : Definition The Multivariate Generalized Ornstein-Uhlenbeck Process Definition: Let (X t, Y t ) t 0 be a Lévy process in (R m m R m, +) such that det(i + X t ) 0 for all t 0 and let V 0 be a random variable in R m. Then the process (V t ) t 0 in R m given by ( ) V t := E (X ) 1 t V 0 + E (X )s dy s (0,t] will be called multivariate generalized Ornstein-Uhlenbeck (MGOU) process driven by (X t, Y t ) t 0. Alexander Lindner, 18

37 Multivariate generalized Ornstein-Uhlenbeck processes : Definition The Multivariate Generalized Ornstein-Uhlenbeck Process Definition: Let (X t, Y t ) t 0 be a Lévy process in (R m m R m, +) such that det(i + X t ) 0 for all t 0 and let V 0 be a random variable in R m. Then the process (V t ) t 0 in R m given by ( ) V t := E (X ) 1 t V 0 + E (X )s dy s (0,t] will be called multivariate generalized Ornstein-Uhlenbeck (MGOU) process driven by (X t, Y t ) t 0. Remark: V 0 not a priori independent of (X, Y ). Alexander Lindner, 18

38 Multivariate generalized Ornstein-Uhlenbeck processes : Definition The Multivariate Generalized Ornstein-Uhlenbeck Process Definition: Let (X t, Y t ) t 0 be a Lévy process in (R m m R m, +) such that det(i + X t ) 0 for all t 0 and let V 0 be a random variable in R m. Then the process (V t ) t 0 in R m given by ( ) V t := E (X ) 1 t V 0 + E (X )s dy s (0,t] will be called multivariate generalized Ornstein-Uhlenbeck (MGOU) process driven by (X t, Y t ) t 0. Remark: V 0 not a priori independent of (X, Y ). E (X ) 1 t may take negative (definite) values. Alexander Lindner, 18

39 Multivariate generalized Ornstein-Uhlenbeck processes : The SDE The corresponding SDE Theorem: The MGOU process ( V t := E (X ) 1 t V 0 + (0,t] E (X )s dy s ) driven by the Lévy process (X t, Y t ) t 0 in (R m m R m, +) is the unique solution of the SDE dv t = du t V t + dl t, t 0, for the Lévy process (U t, L t ) t 0 in (R m m R m, +) given by ( ) ( Ut Xt + [X, X ] c t + ( 0<s t (I + Xs ) 1 ) I + X s := L t Y t + ( 0<s t (I + Xs ) 1 I ) Y s [X, Y ] c t ), for t 0. skip subsection Alexander Lindner, 19

40 Multivariate generalized Ornstein-Uhlenbeck processes : Stationarity Stationary Solutions - Part 1 Theorem: Suppose (V t ) t 0 is a MGOU process driven by the Lévy process (X t, Y t ) t 0 in (R m m R m, +). Let (U t, L t ) t 0 be the Lévy process defined as above. (i) Suppose lim t E (U)t = 0 in probability, then: A finite random variable V 0 can be chosen such that (V t ) t 0 is strictly stationary E (U)s dl s converges in distribution. (0,t] Alexander Lindner, 20

41 Multivariate generalized Ornstein-Uhlenbeck processes : Stationarity Stationary Solutions - Part 1 Theorem: Suppose (V t ) t 0 is a MGOU process driven by the Lévy process (X t, Y t ) t 0 in (R m m R m, +). Let (U t, L t ) t 0 be the Lévy process defined as above. (i) Suppose lim t E (U)t = 0 in probability, then: A finite random variable V 0 can be chosen such that (V t ) t 0 is strictly stationary E (U)s dl s converges in distribution. (0,t] In this case, the distribution of the strictly stationary process (V t ) t 0 is uniquely determined and is obtained by choosing V 0 independent of (X t, Y t ) t 0 as the distributional limit of (0,t] E (U)s dl s as t. Alexander Lindner, 20

42 Multivariate generalized Ornstein-Uhlenbeck processes : Stationarity Stationary Solutions - Part 1 Theorem: Suppose (V t ) t 0 is a MGOU process driven by the Lévy process (X t, Y t ) t 0 in (R m m R m, +). Let (U t, L t ) t 0 be the Lévy process defined as above. (ii) Suppose lim t E (X )t = 0 in probability, then: A finite random variable V 0 can be chosen such that (V t ) t 0 is strictly stationary E (X )s dy s converges in probability. (0,t] Alexander Lindner, 20

43 Multivariate generalized Ornstein-Uhlenbeck processes : Stationarity Stationary Solutions - Part 1 Theorem: Suppose (V t ) t 0 is a MGOU process driven by the Lévy process (X t, Y t ) t 0 in (R m m R m, +). Let (U t, L t ) t 0 be the Lévy process defined as above. (ii) Suppose lim t E (X )t = 0 in probability, then: A finite random variable V 0 can be chosen such that (V t ) t 0 is strictly stationary E (X )s dy s converges in probability. (0,t] In this case the strictly stationary solution is unique and given by V t = E (X ) 1 t E (X )s dy s a.s. for all t 0. (t, ) skip now Alexander Lindner, 20

44 Multivariate generalized Ornstein-Uhlenbeck processes : Stationarity MGOU processes on affine subspaces Definition: Suppose (X t, Y t ) t 0 is a Lévy process in (R m m R m, +) such that E (X ) is non-singular and define (A s,t, B s,t ) 0 s t by ( ) As,t := E (X ) 1 t E (X )s B. s,t E (X ) 1 t (s,t] E (X )u dy u Then an affine subspace H of R m is called invariant under the autoregressive model V t = A s,t V s + B s,t, 0 s t, if A s,t H + B s,t H almost surely, holds for all 0 s t. If R m is the only invariant affine subspace, the model is called irreducible. Alexander Lindner, 21

45 Multivariate generalized Ornstein-Uhlenbeck processes : Stationarity MGOU processes on affine subspaces Theorem: The autoregressive model V t = A s,t V s + B s,t, 0 s t, is irreducible if and only if there exists no pair (O, K) of an orthogonal transformation O R m m and a constant K = (k 1,..., k d ) T R d, 1 d m, such that a.s. ( ) ( ) OX t O 1 X 1 = t 0 X 1 and OY t = t K X 2 t X 3 t where X 1 t R d d, t 0. With (U t, L t ) t 0 as defined above this is equivalent to ( ) ( ) OU t O 1 U 1 = t 0 U 1 and OL t = t K a.s. with U 1 t R d d. U 2 t U 3 t Y 2 t L 2 t Alexander Lindner, 22

46 Multivariate generalized Ornstein-Uhlenbeck processes : Stationarity Stationary Solutions of MGOU processes - Part 2 Theorem: Suppose (V t ) t 0 is a MGOU process driven by the Lévy process (X t, Y t ) t 0 in R m m R m such that the corresponding autoregressive model V t = A s,t V s + B s,t, 0 s t, with (A s,t, B s,t ) 0 s t as defined before is irreducible. Let (U t, L t ) t 0 be defined as above. Then A finite random variable V 0, independent of (X t, Y t ) t 0, can be chosen such that (V t ) t 0 is strictly stationary lim t E (U)t = 0 in probability and (0,t] E (U)s dl s converges in distribution. Alexander Lindner, 23

47 Multivariate generalized Ornstein-Uhlenbeck processes : Stationarity Stationary Solutions of MGOU processes - Part 2 Theorem: Suppose (V t ) t 0 is a MGOU process driven by the Lévy process (X t, Y t ) t 0 in R m m R m such that the corresponding autoregressive model V t = A s,t V s + B s,t, 0 s t, with (A s,t, B s,t ) 0 s t as defined before is irreducible. Let (U t, L t ) t 0 be defined as above. Then A finite random variable V 0, independent of (X t, Y t ) t 0, can be chosen such that (V t ) t 0 is strictly stationary lim t E (U)t = 0 in probability and (0,t] E (U)s dl s converges in distribution. Alexander Lindner, 23

48 Multivariate generalized Ornstein-Uhlenbeck processes : Stationarity Stationary Solutions of MGOU processes - Part 2 Theorem: Suppose (V t ) t 0 is a MGOU process driven by the Lévy process (X t, Y t ) t 0 in R m m R m such that the corresponding autoregressive model V t = A s,t V s + B s,t, 0 s t, with (A s,t, B s,t ) 0 s t as defined before is irreducible. Let (U t, L t ) t 0 be defined as above. Then A finite random variable V 0, independent of (X t, Y t ) t 0, can be chosen such that (V t ) t 0 is strictly stationary lim t E (U)t = 0 in probability and (0,t] E (U)s dl s converges in distribution. A similar result for strictly noncausal strictly stationary solutions of MGOU processes can be obtained, too. skip now Alexander Lindner, 23

49 Extensions : Multivariate volatilities Extensions Behme (2012) obtains various further results, in particular the moment structure of multivariate generalized Ornstein Uhlenbeck processes. She further considers matrix valued positive semidefinite generalized Ornstein Uhlenbeck processes: Alexander Lindner, 24

50 Extensions : Multivariate volatilities Extensions Behme (2012) obtains various further results, in particular the moment structure of multivariate generalized Ornstein Uhlenbeck processes. She further considers matrix valued positive semidefinite generalized Ornstein Uhlenbeck processes: Often, in the one dimensional case volatilities are modeled as square-root process of a generalized Ornstein-Uhlenbeck process. Hence to construct a multivariate volatility model similarly, we have to ensure our processes to be positive semidefinite. One possibility hereby is to consider processes which fulfill V t = A s,t V s A T s,t + B s,t, 0 s t with A s,t in GL(R, m) and B s,t R m m positive semidefinite. Alexander Lindner, 24

51 Extensions : Multivariate volatilities Extensions Behme (2012) obtains various further results, in particular the moment structure of multivariate generalized Ornstein Uhlenbeck processes. She further considers matrix valued positive semidefinite generalized Ornstein Uhlenbeck processes: Often, in the one dimensional case volatilities are modeled as square-root process of a generalized Ornstein-Uhlenbeck process. Hence to construct a multivariate volatility model similarly, we have to ensure our processes to be positive semidefinite. One possibility hereby is to consider processes which fulfill V t = A s,t V s A T s,t + B s,t, 0 s t with A s,t in GL(R, m) and B s,t R m m positive semidefinite. This is equivalent to vec V t = (A s,t A s,t )vec V s + vec B s,t. Alexander Lindner, 24

52 Extensions : Multivariate volatilities A Construction Arguing as above we see that the only process which fulfills the above random recurrence equation is given by ( ) V t = E (X ) 1 t V 0 + E (X )s dy s ( E (X ) s ) T (0,t] for a Lévy process (X, Y ) R m m R m m ( E (X ) 1 t ) T, Alexander Lindner, 25

53 Extensions : Multivariate volatilities A Construction Arguing as above we see that the only process which fulfills the above random recurrence equation is given by ( ) V t = E (X ) 1 t V 0 + E (X )s dy s ( E (X ) s ) T (0,t] for a Lévy process (X, Y ) R m m R m m and that ( vec V t = E (X ) 1 t t E (X ) 1 t vec V 0 + = E (X) 1 t ( vec V 0 + t is a MGOU process driven by the Lévy process (X, Y) R m2 m 2 R m2 with and Y t = vec (Y t ). 0 0 ( E (X ) 1 t ) T, E (X )s E (X ) s dy s ) E (X)s dy s ), t 0 X t = I X t + X t I + [X I, I X ] t, t 0 Alexander Lindner, 25

54 Extensions : Multivariate volatilities A Condition for Positive Semidefiniteness The process V t = E (X ) 1 t ( V 0 + (0,t] E (X )s dy s ( E (X ) s ) T ) ( E (X ) 1 t ) T, is positive semidefinite for all t 0 and all positive semidefinite starting random variables V 0 if and only if Y is a matrix subordinator. Alexander Lindner, 26

55 : Thank you for your attention! Alexander Lindner, 27

56 : Main references: A. Behme and A. Lindner (2012) Multivariate Generalized Ornstein-Uhlenbeck Processes. Stoch. Proc. Appl A. Behme (2012) Moments of MGOU Processes and Positive Semidefinite Matrix Processes. JMVA 111. P. Bougerol and N. Picard (1992) Strict stationarity of generalized autoregressive processes. Ann. Probab. 20. L. de Haan and R.L. Karandikar (1989) Embedding a stochastic difference equation into a continuous-time process. Stoch. Proc. Appl. 32. Alexander Lindner, 28

GARCH processes continuous counterparts (Part 2)

GARCH processes continuous counterparts (Part 2) GARCH processes continuous counterparts (Part 2) Alexander Lindner Centre of Mathematical Sciences Technical University of Munich D 85747 Garching Germany lindner@ma.tum.de http://www-m1.ma.tum.de/m4/pers/lindner/

More information

Continuous Time Approximations to GARCH(1, 1)-Family Models and Their Limiting Properties

Continuous Time Approximations to GARCH(1, 1)-Family Models and Their Limiting Properties Communications for Statistical Applications and Methods 214, Vol. 21, No. 4, 327 334 DOI: http://dx.doi.org/1.5351/csam.214.21.4.327 Print ISSN 2287-7843 / Online ISSN 2383-4757 Continuous Time Approximations

More information

GARCH processes probabilistic properties (Part 1)

GARCH processes probabilistic properties (Part 1) GARCH processes probabilistic properties (Part 1) Alexander Lindner Centre of Mathematical Sciences Technical University of Munich D 85747 Garching Germany lindner@ma.tum.de http://www-m1.ma.tum.de/m4/pers/lindner/

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

Ornstein-Uhlenbeck Processes and Extensions

Ornstein-Uhlenbeck Processes and Extensions Ornstein-Uhlenbeck Processes and Extensions Ross A. Maller 1, Gernot Müller 2, and Alex Szimayer 3 1 The Australian National University Ross.Maller@anu.edu.au 2 Munich University of Technology mueller@ma.tum.de

More information

Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals

Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico

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

A Note on the Non-negativity of Continuous-time ARMA and GARCH Processes

A Note on the Non-negativity of Continuous-time ARMA and GARCH Processes A Note on the Non-negativity of Continuous-time ARMA and GARCH Processes HENGHSIU TSAI Institute of Statistical Science, Academia Sinica, Taipei, Taiwan 115, R.O.C. htsai@stat.sinica.edu.tw K. S. CHAN

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

ELEMENTS OF PROBABILITY THEORY

ELEMENTS OF PROBABILITY THEORY ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable

More information

Stationary particle Systems

Stationary particle Systems Stationary particle Systems Kaspar Stucki (joint work with Ilya Molchanov) University of Bern 5.9 2011 Introduction Poisson process Definition Let Λ be a Radon measure on R p. A Poisson process Π is a

More information

Generalized fractional Lévy processes 1

Generalized fractional Lévy processes 1 Generalized fractional Lévy processes 1 Claudia Klüppelberg Center for Mathematical Sciences and Institute for Advanced Study Technische Universität München Sandbjerg, January 2010 1 1 joint work with

More information

Ornstein-Uhlenbeck processes for geophysical data analysis

Ornstein-Uhlenbeck processes for geophysical data analysis Ornstein-Uhlenbeck processes for geophysical data analysis Semere Habtemicael Department of Mathematics North Dakota State University May 19, 2014 Outline 1 Introduction 2 Model 3 Characteristic function

More information

Exponential functionals of Lévy processes

Exponential functionals of Lévy processes Exponential functionals of Lévy processes Víctor Rivero Centro de Investigación en Matemáticas, México. 1/ 28 Outline of the talk Introduction Exponential functionals of spectrally positive Lévy processes

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

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

Lower Tail Probabilities and Related Problems

Lower Tail Probabilities and Related Problems Lower Tail Probabilities and Related Problems Qi-Man Shao National University of Singapore and University of Oregon qmshao@darkwing.uoregon.edu . Lower Tail Probabilities Let {X t, t T } be a real valued

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

Bernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010

Bernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010 1 Stochastic Calculus Notes Abril 13 th, 1 As we have seen in previous lessons, the stochastic integral with respect to the Brownian motion shows a behavior different from the classical Riemann-Stieltjes

More information

Strictly Stationary Solutions of Autoregressive Moving Average Equations

Strictly Stationary Solutions of Autoregressive Moving Average Equations Strictly Stationary Solutions of Autoregressive Moving Average Equations Peter J. Brockwell Alexander Lindner Abstract Necessary and sufficient conditions for the existence of a strictly stationary solution

More information

Claudia Klüppelberg Technische Universität München Zentrum Mathematik. Oslo, September 2012

Claudia Klüppelberg Technische Universität München Zentrum Mathematik. Oslo, September 2012 Claudia Klüppelberg Technische Universität München Zentrum Mathematik Oslo, September 2012 Collaborators Fred E. Benth, Oslo U Christine Bernhard, former Diploma student TUM Peter Brockwell Vicky Fasen,

More information

1 Introduction. 2 Diffusion equation and central limit theorem. The content of these notes is also covered by chapter 3 section B of [1].

1 Introduction. 2 Diffusion equation and central limit theorem. The content of these notes is also covered by chapter 3 section B of [1]. 1 Introduction The content of these notes is also covered by chapter 3 section B of [1]. Diffusion equation and central limit theorem Consider a sequence {ξ i } i=1 i.i.d. ξ i = d ξ with ξ : Ω { Dx, 0,

More information

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

Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility 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,

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

A Short Introduction to Diffusion Processes and Ito Calculus

A Short Introduction to Diffusion Processes and Ito Calculus A Short Introduction to Diffusion Processes and Ito Calculus Cédric Archambeau University College, London Center for Computational Statistics and Machine Learning c.archambeau@cs.ucl.ac.uk January 24,

More information

Filtrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition

Filtrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition Filtrations, Markov Processes and Martingales Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition David pplebaum Probability and Statistics Department,

More information

On Minimal Entropy Martingale Measures

On Minimal Entropy Martingale Measures On Minimal Entropy Martingale Measures Andrii Andrusiv Friedrich Schiller University of Jena, Germany Marie Curie ITN Outline 1. One Period Model Definitions and Notations Minimal Entropy Martingale Measure

More information

Extremal behavior of stochastic volatility models

Extremal behavior of stochastic volatility models Extremal behavior of stochastic volatility models Vicky Fasen Claudia Klüppelberg Alexander Lindner Abstract Empirical volatility changes in time and exhibits tails, which are heavier than normal. Moreover,

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

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

Convoluted Brownian motions: a class of remarkable Gaussian processes

Convoluted Brownian motions: a class of remarkable Gaussian processes Convoluted Brownian motions: a class of remarkable Gaussian processes Sylvie Roelly Random models with applications in the natural sciences Bogotá, December 11-15, 217 S. Roelly (Universität Potsdam) 1

More information

HSC Research Report. The Lamperti transformation for self-similar processes HSC/97/02. Krzysztof Burnecki* Makoto Maejima** Aleksander Weron*, ***

HSC Research Report. The Lamperti transformation for self-similar processes HSC/97/02. Krzysztof Burnecki* Makoto Maejima** Aleksander Weron*, *** HSC Research Report HSC/97/02 The Lamperti transformation for self-similar processes Krzysztof Burnecki* Makoto Maeima** Aleksander Weron*, *** * Hugo Steinhaus Center, Wrocław University of Technology,

More information

Stochastic Processes

Stochastic Processes Stochastic Processes Stochastic Process Non Formal Definition: Non formal: A stochastic process (random process) is the opposite of a deterministic process such as one defined by a differential equation.

More information

16. M. Maejima, Some sojourn time problems for strongly dependent Gaussian processes, Z. Wahrscheinlichkeitstheorie verw. Gebiete 57 (1981), 1 14.

16. M. Maejima, Some sojourn time problems for strongly dependent Gaussian processes, Z. Wahrscheinlichkeitstheorie verw. Gebiete 57 (1981), 1 14. Publications 1. M. Maejima, Some limit theorems for renewal processes with non-identically distributed random variables, Keio Engrg. Rep. 24 (1971), 67 83. 2. M. Maejima, A generalization of Blackwell

More information

Gaussian Processes. 1. Basic Notions

Gaussian Processes. 1. Basic Notions Gaussian Processes 1. Basic Notions Let T be a set, and X : {X } T a stochastic process, defined on a suitable probability space (Ω P), that is indexed by T. Definition 1.1. We say that X is a Gaussian

More information

Homework #6 : final examination Due on March 22nd : individual work

Homework #6 : final examination Due on March 22nd : individual work Université de ennes Année 28-29 Master 2ème Mathématiques Modèles stochastiques continus ou à sauts Homework #6 : final examination Due on March 22nd : individual work Exercise Warm-up : behaviour of characteristic

More information

Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process

Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process Mátyás Barczy University of Debrecen, Hungary The 3 rd Workshop on Branching Processes and Related

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

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

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

Systems Driven by Alpha-Stable Noises

Systems Driven by Alpha-Stable Noises Engineering Mechanics:A Force for the 21 st Century Proceedings of the 12 th Engineering Mechanics Conference La Jolla, California, May 17-20, 1998 H. Murakami and J. E. Luco (Editors) @ASCE, Reston, VA,

More information

for all f satisfying E[ f(x) ] <.

for all f satisfying E[ f(x) ] <. . Let (Ω, F, P ) be a probability space and D be a sub-σ-algebra of F. An (H, H)-valued random variable X is independent of D if and only if P ({X Γ} D) = P {X Γ}P (D) for all Γ H and D D. Prove that if

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

Exercises in stochastic analysis

Exercises in stochastic analysis Exercises in stochastic analysis Franco Flandoli, Mario Maurelli, Dario Trevisan The exercises with a P are those which have been done totally or partially) in the previous lectures; the exercises with

More information

Affine Processes are regular

Affine Processes are regular Affine Processes are regular Martin Keller-Ressel kemartin@math.ethz.ch ETH Zürich Based on joint work with Walter Schachermayer and Josef Teichmann (arxiv:0906.3392) Forschungsseminar Stochastische Analysis

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

On pathwise stochastic integration

On pathwise stochastic integration On pathwise stochastic integration Rafa l Marcin Lochowski Afican Institute for Mathematical Sciences, Warsaw School of Economics UWC seminar Rafa l Marcin Lochowski (AIMS, WSE) On pathwise stochastic

More information

When is a Moving Average a Semimartingale?

When is a Moving Average a Semimartingale? 29 Barrett Lectures Ph.D.-student under supervision of Jan Pedersen, Thiele Centre, University of Aarhus, Denmark. 29 Barrett Lectures at The University of Tennessee: Stochastic Analysis and its Applications

More information

Max stable Processes & Random Fields: Representations, Models, and Prediction

Max stable Processes & Random Fields: Representations, Models, and Prediction Max stable Processes & Random Fields: Representations, Models, and Prediction Stilian Stoev University of Michigan, Ann Arbor March 2, 2011 Based on joint works with Yizao Wang and Murad S. Taqqu. 1 Preliminaries

More information

Asymptotic Irrelevance of Initial Conditions for Skorohod Reflection Mapping on the Nonnegative Orthant

Asymptotic Irrelevance of Initial Conditions for Skorohod Reflection Mapping on the Nonnegative Orthant Published online ahead of print March 2, 212 MATHEMATICS OF OPERATIONS RESEARCH Articles in Advance, pp. 1 12 ISSN 364-765X print) ISSN 1526-5471 online) http://dx.doi.org/1.1287/moor.112.538 212 INFORMS

More information

GAUSSIAN PROCESSES; KOLMOGOROV-CHENTSOV THEOREM

GAUSSIAN PROCESSES; KOLMOGOROV-CHENTSOV THEOREM GAUSSIAN PROCESSES; KOLMOGOROV-CHENTSOV THEOREM STEVEN P. LALLEY 1. GAUSSIAN PROCESSES: DEFINITIONS AND EXAMPLES Definition 1.1. A standard (one-dimensional) Wiener process (also called Brownian motion)

More information

Properties of the Autocorrelation Function

Properties of the Autocorrelation Function Properties of the Autocorrelation Function I The autocorrelation function of a (real-valued) random process satisfies the following properties: 1. R X (t, t) 0 2. R X (t, u) =R X (u, t) (symmetry) 3. R

More information

Long-Range Dependence and Self-Similarity. c Vladas Pipiras and Murad S. Taqqu

Long-Range Dependence and Self-Similarity. c Vladas Pipiras and Murad S. Taqqu Long-Range Dependence and Self-Similarity c Vladas Pipiras and Murad S. Taqqu January 24, 2016 Contents Contents 2 Preface 8 List of abbreviations 10 Notation 11 1 A brief overview of times series and

More information

Lecture 1: Pragmatic Introduction to Stochastic Differential Equations

Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic

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

Jump Processes. Richard F. Bass

Jump Processes. Richard F. Bass Jump Processes Richard F. Bass ii c Copyright 214 Richard F. Bass Contents 1 Poisson processes 1 1.1 Definitions............................. 1 1.2 Stopping times.......................... 3 1.3 Markov

More information

Bernstein-gamma functions and exponential functionals of Lévy processes

Bernstein-gamma functions and exponential functionals of Lévy processes Bernstein-gamma functions and exponential functionals of Lévy processes M. Savov 1 joint work with P. Patie 2 FCPNLO 216, Bilbao November 216 1 Marie Sklodowska Curie Individual Fellowship at IMI, BAS,

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

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

Some Tools From Stochastic Analysis

Some Tools From Stochastic Analysis W H I T E Some Tools From Stochastic Analysis J. Potthoff Lehrstuhl für Mathematik V Universität Mannheim email: potthoff@math.uni-mannheim.de url: http://ls5.math.uni-mannheim.de To close the file, click

More information

1. Stochastic Processes and filtrations

1. Stochastic Processes and filtrations 1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S

More information

Max stable processes: representations, ergodic properties and some statistical applications

Max stable processes: representations, ergodic properties and some statistical applications Max stable processes: representations, ergodic properties and some statistical applications Stilian Stoev University of Michigan, Ann Arbor Oberwolfach March 21, 2008 The focus Let X = {X t } t R be a

More information

Functional Limit Theorems for A New Class of Non-Stationary Shot Noise Processes

Functional Limit Theorems for A New Class of Non-Stationary Shot Noise Processes Functional Limit Theorems for A New Class of Non-Stationary Shot Noise Processes GUODONG PANG AND YUHANG ZHOU Abstract. We study a class of non-stationary shot noise processes which have a general arrival

More information

Limit theorems for multipower variation in the presence of jumps

Limit theorems for multipower variation in the presence of jumps Limit theorems for multipower variation in the presence of jumps Ole E. Barndorff-Nielsen Department of Mathematical Sciences, University of Aarhus, Ny Munkegade, DK-8 Aarhus C, Denmark oebn@imf.au.dk

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

Non-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings. Krzysztof Burdzy University of Washington

Non-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings. Krzysztof Burdzy University of Washington Non-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings Krzysztof Burdzy University of Washington 1 Review See B and Kendall (2000) for more details. See also the unpublished

More information

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3 Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................

More information

THE LENT PARTICLE FORMULA

THE LENT PARTICLE FORMULA THE LENT PARTICLE FORMULA Nicolas BOULEAU, Laurent DENIS, Paris. Workshop on Stochastic Analysis and Finance, Hong-Kong, June-July 2009 This is part of a joint work with Laurent Denis, concerning the approach

More information

Brownian Motion and Poisson Process

Brownian Motion and Poisson Process and Poisson Process She: What is white noise? He: It is the best model of a totally unpredictable process. She: Are you implying, I am white noise? He: No, it does not exist. Dialogue of an unknown couple.

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

AN EXTREME-VALUE ANALYSIS OF THE LIL FOR BROWNIAN MOTION 1. INTRODUCTION

AN EXTREME-VALUE ANALYSIS OF THE LIL FOR BROWNIAN MOTION 1. INTRODUCTION AN EXTREME-VALUE ANALYSIS OF THE LIL FOR BROWNIAN MOTION DAVAR KHOSHNEVISAN, DAVID A. LEVIN, AND ZHAN SHI ABSTRACT. We present an extreme-value analysis of the classical law of the iterated logarithm LIL

More information

B8.3 Mathematical Models for Financial Derivatives. Hilary Term Solution Sheet 2

B8.3 Mathematical Models for Financial Derivatives. Hilary Term Solution Sheet 2 B8.3 Mathematical Models for Financial Derivatives Hilary Term 18 Solution Sheet In the following W t ) t denotes a standard Brownian motion and t > denotes time. A partition π of the interval, t is a

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

Exponential martingales: uniform integrability results and applications to point processes

Exponential martingales: uniform integrability results and applications to point processes Exponential martingales: uniform integrability results and applications to point processes Alexander Sokol Department of Mathematical Sciences, University of Copenhagen 26 September, 2012 1 / 39 Agenda

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

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

Bessel-like SPDEs. Lorenzo Zambotti, Sorbonne Université (joint work with Henri Elad-Altman) 15th May 2018, Luminy

Bessel-like SPDEs. Lorenzo Zambotti, Sorbonne Université (joint work with Henri Elad-Altman) 15th May 2018, Luminy Bessel-like SPDEs, Sorbonne Université (joint work with Henri Elad-Altman) Squared Bessel processes Let δ, y, and (B t ) t a BM. By Yamada-Watanabe s Theorem, there exists a unique (strong) solution (Y

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

Logarithmic scaling of planar random walk s local times

Logarithmic scaling of planar random walk s local times Logarithmic scaling of planar random walk s local times Péter Nándori * and Zeyu Shen ** * Department of Mathematics, University of Maryland ** Courant Institute, New York University October 9, 2015 Abstract

More information

Stochastic integration. P.J.C. Spreij

Stochastic integration. P.J.C. Spreij Stochastic integration P.J.C. Spreij this version: April 22, 29 Contents 1 Stochastic processes 1 1.1 General theory............................... 1 1.2 Stopping times...............................

More information

Introduction to numerical simulations for Stochastic ODEs

Introduction to numerical simulations for Stochastic ODEs Introduction to numerical simulations for Stochastic ODEs Xingye Kan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616 August 9, 2010 Outline 1 Preliminaries 2 Numerical

More information

Semi-Lévy driven continuous-time GARCH process

Semi-Lévy driven continuous-time GARCH process Semi-Lévy driven continuous-time GARCH process M. Mohammadi S. Rezakhah N. Modarresi arxiv:1812.10777v1 [math.pr] 27 Dec 2018 Abstract In this paper, we study the class of semi-lévy driven continuous-time

More information

Modelling energy forward prices Representation of ambit fields

Modelling energy forward prices Representation of ambit fields Modelling energy forward prices Representation of ambit fields Fred Espen Benth Department of Mathematics University of Oslo, Norway In collaboration with Heidar Eyjolfsson, Barbara Rüdiger and Andre Süss

More information

Multivariate Generalized Ornstein-Uhlenbeck Processes

Multivariate Generalized Ornstein-Uhlenbeck Processes Mulivariae Generalized Ornsein-Uhlenbeck Processes Ania Behme and Alexander Lindner Absrac De Haan and Karandikar [12] inroduced generalized Ornsein Uhlenbeck processes as one-dimensional processes (V

More information

What s more chaotic than chaos itself? Brownian Motion - before, after, and beyond.

What s more chaotic than chaos itself? Brownian Motion - before, after, and beyond. Include Only If Paper Has a Subtitle Department of Mathematics and Statistics What s more chaotic than chaos itself? Brownian Motion - before, after, and beyond. Math Graduate Seminar March 2, 2011 Outline

More information

Change of Measure formula and the Hellinger Distance of two Lévy Processes

Change of Measure formula and the Hellinger Distance of two Lévy Processes Change of Measure formula and the Hellinger Distance of two Lévy Processes Erika Hausenblas University of Salzburg, Austria Change of Measure formula and the Hellinger Distance of two Lévy Processes p.1

More information

STATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, = (2n 1)(2n 3) 3 1.

STATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, = (2n 1)(2n 3) 3 1. STATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, 26 Problem Normal Moments (A) Use the Itô formula and Brownian scaling to check that the even moments of the normal distribution

More information

Local vs. Nonlocal Diffusions A Tale of Two Laplacians

Local vs. Nonlocal Diffusions A Tale of Two Laplacians Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2

More information

Stochastic Calculus February 11, / 33

Stochastic Calculus February 11, / 33 Martingale Transform M n martingale with respect to F n, n =, 1, 2,... σ n F n (σ M) n = n 1 i= σ i(m i+1 M i ) is a Martingale E[(σ M) n F n 1 ] n 1 = E[ σ i (M i+1 M i ) F n 1 ] i= n 2 = σ i (M i+1 M

More information

Bernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012

Bernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012 1 Stochastic Calculus Notes March 9 th, 1 In 19, Bachelier proposed for the Paris stock exchange a model for the fluctuations affecting the price X(t) of an asset that was given by the Brownian motion.

More information

SAMPLE CHAPTERS UNESCO EOLSS STATIONARY PROCESSES. E ξ t (1) K.Grill Institut für Statistik und Wahrscheinlichkeitstheorie, TU Wien, Austria

SAMPLE CHAPTERS UNESCO EOLSS STATIONARY PROCESSES. E ξ t (1) K.Grill Institut für Statistik und Wahrscheinlichkeitstheorie, TU Wien, Austria STATIONARY PROCESSES K.Grill Institut für Statistik und Wahrscheinlichkeitstheorie, TU Wien, Austria Keywords: Stationary process, weak sense stationary process, strict sense stationary process, correlation

More information

Generalized fractional Ornstein-Uhlenbeck processes

Generalized fractional Ornstein-Uhlenbeck processes Generalized fractional Ornstein-Uhlenbeck processes Kotaro ENDO and Muneya MATSUI Department of Mathematics, Keio University Center for Mathematical Science, Munich University of Technology and Department

More information

p 1 ( Y p dp) 1/p ( X p dp) 1 1 p

p 1 ( Y p dp) 1/p ( X p dp) 1 1 p Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p

More information

Lecture 4: Stochastic Processes (I)

Lecture 4: Stochastic Processes (I) Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 215 Lecture 4: Stochastic Processes (I) Readings Recommended: Pavliotis [214], sections 1.1, 1.2 Grimmett and Stirzaker [21], section 9.4 (spectral

More information

Infinite Divisibility for Stochastic Processes and Time Change

Infinite Divisibility for Stochastic Processes and Time Change Infinite Divisibility for Stochastic Processes and Time Change Ole E. Barndorff-Nielsen 1, Makoto Maejima 2 and Ken-iti Sato 3 4 Running head: Infinite Divisibility for Stochastic Processes and Time Change

More information

Tyler Hofmeister. University of Calgary Mathematical and Computational Finance Laboratory

Tyler Hofmeister. University of Calgary Mathematical and Computational Finance Laboratory JUMP PROCESSES GENERALIZING STOCHASTIC INTEGRALS WITH JUMPS Tyler Hofmeister University of Calgary Mathematical and Computational Finance Laboratory Overview 1. General Method 2. Poisson Processes 3. Diffusion

More information

Lecture 21 Representations of Martingales

Lecture 21 Representations of Martingales Lecture 21: Representations of Martingales 1 of 11 Course: Theory of Probability II Term: Spring 215 Instructor: Gordan Zitkovic Lecture 21 Representations of Martingales Right-continuous inverses Let

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

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

Lecture 22 Girsanov s Theorem

Lecture 22 Girsanov s Theorem Lecture 22: Girsanov s Theorem of 8 Course: Theory of Probability II Term: Spring 25 Instructor: Gordan Zitkovic Lecture 22 Girsanov s Theorem An example Consider a finite Gaussian random walk X n = n

More information