arxiv: v9 [math.pr] 27 Jan 2014

Size: px
Start display at page:

Download "arxiv: v9 [math.pr] 27 Jan 2014"

Transcription

1 SENSITIVITIES VIA ROUGH PATHS NICOLAS MARIE arxiv: v9 math.pr 27 Jan 214 Abstract. Motivated by a problematic coming from mathematical finance, this paper is devoted to existing and additional results of continuity and differentiability of the Itô map associated to rough differential equations. These regularity results together with Malliavin calculus are applied to sensitivities analysis for stochastic differential equations driven by multidimensional Gaussian processes with continuous paths, especially fractional Brownian motions. Precisely, in that framework, results on computation of greeks for Itô s stochastic differential equations are extended. An application in mathematical finance, and simulations, are provided. Contents 1. Introduction 1 2. Regularity of the Itô map : existing and additional results Differentiability of the Itô map with respect to x and V Differentiability of the Itô map with respect to the driving signal Application to Gaussian stochastic analysis Sensitivity analysis for Gaussian rough differential equations Application to mathematical finance and simulations Computation of sensitivities in a fractional stochastic volatility model Simulations 31 Appendix A. Fractional Brownian motion 34 References 35 Acknowledgements. Many thanks to Laure Coutin and Laurent Decreusefond for their advices. That work was ported by ANR Masterie. 1. Introduction Motivated by a problematic coming from mathematical finance, this paper is devoted to existing and additional results of continuity and differentiability of the Itô map associated to rough differential equations RDEs. These regularity results together with Malliavin calculus are applied to sensitivities analysis for stochastic differential equations SDEs driven by multidimensional Gaussian processes with continuous paths, especially fractional Brownian motions. In order to state the problematic mentioned above, let s remind some notions of mathematical finance. Consider a probability space Ω, A, P, a d-dimensional Brownian motion B and F := A t ; t, T the natural filtration associated to B d N and T >. Key words and phrases. Rough paths, Rough differential equations, Malliavin calculus, Sensitivities, Mathematical finance, Gaussian processes. 1

2 2 NICOLAS MARIE Consider the financial market consisting of d + 1 assets one risk-free asset and d risky assets over the filtered probability space Ω, A, F, P. At time t, T, the deterministic price of the risk-free asset is denoted by St, and prices of the d risky assets are given by the random vector S t := St 1,..., St d. In a first place, assume that the process S is the solution of a stochastic differential equation, taken in the sense of Itô : S t = x + t µ S u du + t σ S u db u ; x R d where, µ : R d R d and σ : R d M d R are globally Lipschitz continuous functions. Let P P be the risk-neutral probability measure of the market i.e. such that S := S/S is a F, P -martingale. Theorem 1.1. Consider an option of payoff h L 2 Ω, A T, P. Then, there exists an admissible strategy ϕ such that : S t, T, V t ϕ = E t S h A t P -a.s. where V ϕ is the associated wealth process. Refer to D. Lamberton and B. Lapeyre 12, Theorem for a proof in the case of the Black-Scholes model. With notations of Theorem 1.1, V T ϕ = E h. It is the option s price and, when h := F S T with a function F from R d into R +, it is useful to get the existence and an expression of the sensitivities of V T ϕ to perturbations of the initial condition and of the volatility function σ in particular : := x E F S x T and V := σ E F S σ T. In finance, these sensitivities are called greeks. For instance, involves in a procedure called -hedging that allows to compute the admissible strategy of Theorem 1.1 cf. 12, Subsection However, there are no reasons to use these quantities in finance only. Indeed, they could also be used in pharmacokinetic as mentioned at 18, Section 5. Greeks have been studied by many authors. In 7, E. Fournié et al. have established the existence of greeks and have provided expressions of them via Malliavin calculus by assuming that σ satisfied the uniform elliptic condition cf. Theorem 1.2. In 1, E. Gobet and R. Münos have extended these results by assuming that σ only satisfied the hypoelliptic condition. On computation of greeks in the Black-Scholes model, refer to P. Malliavin and A. Thalmaier 16, Chapter 2. On sensitivities in models with jumps, refer to N. Privault et al. 21 and 6. Finally, via the cubature formula for the Brownian motion, J. Teichmann has provided estimators of Malliavin s weights for computation of greeks cf. J. Teichmann 24. On the regularity of the solution map associated to SDEs taken in the sense of Itô, refer to H. Kunita 11. At the following theorem, δ is the divergence operator associated to the Brownian motion B cf. D. Nualart 2, Section 1.3. Theorem 1.2. When b and σ are differentiable, of bounded and Lipschitz continuous derivatives, and F L 2 R d ; R + : 1 If σ satisfies the uniform elliptic condition i.e. there exists ε > such that for every a, b R d, b T σ T aσab ε b 2, then exists and there exists an adapted d-dimensional process h such that : = E F S T δh.

3 SENSITIVITIES VIA ROUGH PATHS 3 2 Let σ : R d M d R be a function such that, for every ε belonging to a closed neighborhood of, σ + ε σ satisfies the uniform elliptic condition. Then V exists and there exists an anticipative d-dimensional process h V such that : V = E F S T δh V. Refer to 7, propositions 3.2 and 3.3 for a proof. Under technical assumptions stated at Subsection 2.3, the main purpose of this paper is to extend Theorem 1.2 for the following SDE, taken in the sense of rough paths introduced by T. Lyons in 14 : X t = x + t µ X s ds + t σ X s dw s ; x R d where, W is a centered d-dimensional Gaussian process with continuous paths of finite p-variation p 1, and µ and σ satisfy the following assumption. Assumption 1.3. µ and σ are p + 1 times differentiable, bounded and of bounded derivatives. In order to solve that problematic, subsections 2.1 and 2.2 are devoted to existing and additional results of continuity and differentiability of the Itô map associated to rough differential equations. In particular, the continuous differentiability of the Itô map with respect to the collection of vector fields is proved in addition to existing regularity results for the initial condition and for the driving signal, coming from P. Friz and N. Victoir 9, chapters 4 and 11. In order to apply the probabilistic integrability results coming from T. Cass, C. Litterer and T. Lyons 2, tailor-made upper-bounds are provided for each derivative. Subsection 2.3 reminds definitions and results related to the natural geometric rough paths associated to a Gaussian process having a covariance function satisfying the technical Assumption 2.11, called enhanced Gaussian process by P. Friz and N. Victoir. Results of subsections 2.1 and 2.2 are applied together with results coming from 2, in order to get the probabilistic integrability of the solution of a Gaussian RDE and their derivatives. The main problematic is solved at Section 3 by using results of Section 2 together with Malliavin calculus arguments. Simulations of and V are provided at Subsection 4.2. The fractional Brownian motion fbm introduced in 17 by B.B. Mandelbrot and W. Van Ness has been studied by many authors in order to generalize the Brownian motion classically used to modelize risky assets prices process. For instance, the regularity of the process paths and its memory are both continuously controlled via the Hurst parameter H of the fbm. However, the fbm is not a semi-martingale if H 1/2 cf. 2, Proposition In 22, L.C.G. Rogers has shown the existence of arbitrages for assets having a prices process modelized by a fbm. In order to bypass that difficulty, several approaches have been studied. For instance, in 3, P. Cheridito assumed that the prices process was modelized by a mixed fractional Brownian motion that is a semi-martingale that involves a fbm. At Subsection 4.1, prices of the risky assets are modelized by a fractional SDE, in which the volatility is modelized by another one. Results of Section 3 are applied in order to show the existence and provide an expression of the sensitivity of the option s price with respect to collection of vector fields of the volatility s equation. This paper uses many results on rough paths and rough differential equations coming from 9 and, T. Lyons and Z. Qian 15. It also involves Malliavin calculus results coming from 2.

4 4 NICOLAS MARIE Notations and, short definitions and results, used throughout the paper are stated below. However, in a sake of rigor, original well stated results of the literature are cited throughout the paper. Notations general : R e and R d e, d N are equipped with their euclidean norms, both denoted by.. The canonical basis of R d is denoted by e 1,..., e d. With respect to that basis, for k = 1,..., d ; the k-th component of any vector u R d is denoted by u k. The closed ball of R d for., of center a R d and of radius r >, is denoted by Ba, r. The usual matrix resp. operator norm on M e,d R resp. LR e ; R d is denoted by. M resp.. L. Consider s < t T. The set of dissections of s, t is denoted by D s,t. In particular, D T := D,T. T := {u, v R 2 : u < v T }. The space of continuous resp. continuously differentiable functions from s, t into R d is denoted by C s, t; R d resp. C 1 s, t; R d, and is equipped with the uniform norm. ;s,t. Differentiability means differentiability in the sense of Fréchet in general cf. H. Cartan 1, Chapter I.2. Consider two Banach spaces E and F. Let ϕ : E F be a derivable map at point x E, in the direction h E ; the associated directional derivative is denoted by : ϕx + εh ϕx D h ϕx := lim ε ε Notations rough paths : in F. Consider p 1 and α, 1. The space of continuous functions of finite p- variation resp. α-hölder continuous functions from s, t into R d is denoted by C p-var s, t; R d := y C s, t; R d : D 1 D={r k } D s,t k=1 y rk+1 y rk p < resp. C α-höl s, t; R d, that is a subset of C 1/α-var s, t; R d, and is equipped with the p-variation distance d p-var;s,t resp. the α-hölder distance d α-höl;s,t. Refer to 9, chapters 5 and 8 about these spaces. Consider N N and y :, T R d a continuous function of finite 1- variation. The step-n tensor algebra over R d is denoted by T N R d := R R d R d N, the step-n signature of y is denoted by. S N y := 1, dy r,..., <r 1< <r N <. dy r1 dy rn, and the step-n free nilpotent group over R d is denoted by G N R d := { S N y 1 ; y C 1-var, 1; R d }. Refer to 9, Chapter 7. For k =,..., N, the k + 1-th component of any X T N R d is denoted by X k.

5 SENSITIVITIES VIA ROUGH PATHS 5 The space of geometric p-rough paths is denoted by GΩ p,t R d := { S p y; y C 1-var, T ; R d }d p-var;t, and is equipped with the p-variation distance d p-var;t, or with the uniform distance d ;T, associated to the Carnot-Carathéodory s distance. Refer to 9, Chapter 9. The closed ball of GΩ p,t R d for d p-var;t, of center Y GΩ p,t R d and of radius r >, is denoted by B p,t Y, r. For every Y GΩ p,t R d, ω Y,p : s, t T Y p p-var;s,t is a control. Refer to 9, Chapter 1 about properties of controls. Consider q 1 such that 1/p + 1/q > 1, Y GΩ p,t R d and h GΩ q,t R e. The geometric p-rough path over Y 1, h 1 provided at 9, Theorem 9.26 is denoted by S p Y h. The translation of Y by h provided at 9, Theorem 9.34 is denoted by T h Y. Consider γ >. The space of collections of γ-lipschitz resp. locally γ- Lipschitz vector fields on R e is denoted by Lip γ R e ; R d resp. Lip γ loc Re ; R d cf. 9, Definition 1.2. Lip γ R e ; R d is equipped with the γ-lipschitz norm. lip γ such that, for every V Lip γ R e ; R d, } V lip γ := max { V, DV,..., D γ V, D γ V {γ}-höl. The closed ball of Lip γ R e ; R d for. lip γ, of center V Lip γ R e ; R d and of radius r >, is denoted by B Lip γ V, r. Consider ε >, I a compact interval of, T, a control ω : T R + and Y GΩ p,t R d. Let s put : M ε,i,ω := D = {r k } D I ω r k, r k+1 ε D 1 k=1 ω r k, r k+1, M ε,i,p Y := M ε,i,ωy,p and N ε,i,p Y := {n N : τ n I} where, τ := infi and for every n N, τ n+1 := inf { t I : Y p p-var;τ n,t ε and t > τ n } I. In the sequel, I :=, T. Consider γ > p and V Lip γ loc Re ; R d satisfying the p-non explosion condition i.e. V and D p V are respectively globally Lipschitz continuous and γ p-hölder continuous on R e. The unique solution of dx = V XdW with initial condition X G p R e or X R e, is denoted by π V, X ; W. By 9, Exercice 1.55, if V is a collection of affine vector fields and ω : T R + is a control satisfying W p-var;s,t ω 1/p s, t for every s, t T, there exists a constant C 1 >, not depending on X R e and W, such that : π V, x ; W ;T C x e C1M 1,I,ω. By 9, Theorem 1.36, if V Lip γ R e ; R d, there exists a constant C 2 >, not depending on X G p R e, V and W, such that for every s, t T, π V, X ; W p-var;s,t C 2 V lip γ 1 W p-var;s,t V p W p lip γ 1 p-var;s,t.

6 6 NICOLAS MARIE By 9, Theorem 1.47, if V Lip γ R e ; R d, there exists a constant C 3 >, not depending on V and W, such that for every s, t T, V WdW C 3 V lip γ 1 W p-var;s,t W p p-var;s,t. p-var;s,t Notations Gaussian stochastic analysis : For every t, T,, t is equipped with the Borel σ-algebra B t associated to the usual topology on, t. R d is equipped with the Borel σ-algebra associated to the usual euclidean topology on R d, and G p R d is equipped with the Borel σ-algebra associated to the Carnot-Carathéodory s topology on G p R d. These σ-algebras are both denoted by B. Let W be a d-dimensional centered Gaussian process with continuous paths. Its Cameron-Martin s space is denoted by H 1 := { h C, T ; R : Z W s.t. t, T, h t = EW t Z } with W := span {W t, t, T } L2 cf. 9, Subsection and Section 15.3, its reproducing kernel Hilbert space is denoted by H, and the Wiener integral with respect to W defined on H is denoted by W cf. 2, Section 1.1. The Malliavin derivative associated to W is denoted by D for the R d - valued resp. H-valued random variables, and its domain in L 2 Ω resp. L 2 Ω; H is denoted by D 1,2 resp. D 1,2 H cf. 2, Section 1.2. For R d -valued random variables, the divergence associated to D is denoted by δ, and its domain in L 2 Ω; H is denoted by domδ cf. 2, Section 1.3. The vector space of R d -valued resp. H-valued random variables locally derivable in the sense of Malliavin is denoted by D 1,2 loc 2, Subsection resp. D1,2 loc H cf. 2. Regularity of the Itô map : existing and additional results This section is devoted to the regularity of the Itô map associated to a RDE. On one hand, results of continuity and differentiability of the Itô map with respect to the initial condition and the driving signal coming from 9, Chapter 11 are stated. In addition, the continuous differentiability of the Itô map with respect to the collection of vector fields is proved. On the other hand, in order to apply the integrability results coming from 2, tailor-made upper-bounds are provided for each derivative. Let s first synthesize existing continuity results for the Itô map and the rough integral. Theorem 2.1. Consider R > : 1 The Itô map X, W, V π V, X ; W is uniformly continuous from 2 The map G p R e B p,t 1, R Lip γ R e ; R d into GΩ p,t R d. is uniformly continuous from J : W, V V WdW B p,t 1, R Lip γ 1 R d ; R d into GΩ p,t R d.

7 SENSITIVITIES VIA ROUGH PATHS 7 In both cases, the uniform continuity holds when B p,t 1, R and GΩ p,t R d are equipped with the uniform distance d ;T. Refer to 9, corollaries 1.39,4,48 for a proof. Remark. Consider x R e, W GΩ p,t R d and V := V 1,..., V d a collection of affine vector fields on R e. By 9, Theorem 1.53, π V, x ; W t belongs to the ball B; Rx, W of R e for every t, T, where Rx, W := C1 + x e C W p p-var;t and C > is a constant not depending on x and W. Moreover, for every x R e and every W GΩ p,t R d, x x and W p-var;t W p-var;t = R x, W Rx, W. So, if ˆV Lip γ R e ; R d is the collection of vector fields satisfying ˆV V on B; Rx, W, then π V,. π ˆV,. on the set B, x B p,t 1, W p-var;t. Therefore, by Theorem 2.1, the map π V,. is uniformly continuous from B, x B p,t 1, W p-var;t into C p-var, T ; R e. The uniform continuity holds when B p,t 1, W p-var;t and C p-var, T ; R e are equipped with the uniform distance d ;T. The following technical corollary of 9, Theorem 9.26 allows to apply integrability results of 2 to differential equations having a drift term. Corollary 2.2. Consider p > q 1 such that 1/p + 1/q > 1, Y GΩ p,t R d, h GΩ q,t R e and ε >. There exists a constant C >, depending only on p and q, such that : M ε,i,p S p Y h C h p q-var;t + M ε,i,py. Proof. On one hand, for every s, t T, ω Y,p s, t = Y p-var;s,t S p Y h p-var;s,t. On the other hand, since p/q 1 and, ω Y,p and ω h,q are two controls : is also a control. ω = Y p p-var + h p q-var = ω Y,p + ω p/q h,q Then, by 9, Proposition 7.52, there exists a constant C 1, depending only on p and q, such that for every s, t T, S p Y h p Cωs, t. p-var;s,t In conclusion, M ε,i,p Sp Y h C C D = {r k } D I ω r k, r k+1 ε D 1 k=1 h p q-var;t + M ε,i,py ωr k, r k+1 by er-additivity of the control ω p/q h,q.

8 8 NICOLAS MARIE 2.1. Differentiability of the Itô map with respect to x and V. In order to prove the continuous differentiability of the Itô map for RDEs with respect to the collection of vector fields, let s first prove it for ODEs. Proposition 2.3. Consider γ 1, x R e and a continuous function w :, T R d of finite 1-variation. Then, the map V π V, x ; w is continuously differentiable from Lip γ R e ; R d into C 1-var, T ; R e. Proof. In a first step, the derivability of the Itô map with respect the collection of vector fields is established at every points and in every directions of Lip γ R d ; R e. In a second step, via 9, Proposition B.5, the continuous differentiability of that partial Itô map is proved. Step 1. Consider V, Ṽ Lipγ R e ; R d, ε, 1, x V := π V, x ; w and y V,Ṽ the solution of the following ODE : 1 y V,Ṽ t = For every t, T, x V +εṽ t x V t ε t y V,Ṽ t = where, t P t ε := ε 1 Q t ε := R t ε := t t DV x V s, y V,Ṽ s dw s + t t t V xs V +εṽ V x V s ε Ṽ x V +εṽ s Ṽ xv s = P t ε + Q t ε + R t ε V x V +εṽ s Ṽ x V +εṽ s Ṽ xv s DV x V s, ε 1 x V +εṽ s Ṽ x V s dw s. DV x V s, y V,Ṽ s dw s dw s + V x V s DV x V s, x V s +εṽ x V s dw s, dw s and x V s y V,Ṽ s dw s. Firstly, since V is continuously differentiable on R e, by 9, Lemma 4.2 : P t ε ε 1 w 1-var;T t,t ηεε 1 w 1-var;T x V +εṽ x V ;T where, ηε when ε. By 9, Theorem 3.18 : V xt V +εṽ V x V t DV x V t, x V t +εṽ x V t 2 P ε ;T M 3 ε := 2ηεe 3M1M2 M 2 Ṽ w 1-var;T with M 1 := V lip γ + Ṽ lip γ V + εṽ lip 1 V lip 1 and M 2 := w 1-var;T. Secondly, since Ṽ is continuously differentiable and of bounded derivative on Re, it is a collection of Lipschitz continuous vector fields. Then, by 9, Theorem 3.18 again : 3 Qε ;T M 4 ε := 2εe 3M1M2 M 2 Ṽ 2 lip γ w 1-var;T.

9 SENSITIVITIES VIA ROUGH PATHS 9 Thirdly, 4 R t ε V lip γ t x V +εṽ s x V s ε y V,Ṽ s dw s. Therefore, inequalities 2, 3 and 4 imply that : xt V +εṽ x V t y V,Ṽ t ε M 3ε+M 4 ε+ V lip γ t x V +εṽ s x V s ε y V,Ṽ s In conclusion, by Gronwall s lemma : x V +εṽ x V y V,Ṽ M 3 ε + M 4 ε e V ε lipγ w 1-var;T ;T ε. Step 2. The solution of equation 1 satisfies : DṼ x V = π A, ;. JF V, Ṽ.,. π V, x ;.,.w dw s. where, A : R e LLR e R e ; R e and F V, Ṽ : Re R d LR e R d ; LR e R e are two collections of vector fields, respectively defined by : AaL, b := L.a + b and F V, Ṽ a, a b, b := DV a,. b ; Ṽ ab for every a, b R e, a, b R d and L LR e. Firstly, by the second point of Theorem 2.1, the map J is uniformly continuous on every bounded sets of C 1-var, T ; R e R d C 1-var, T ; R e R d. Secondly, the map V, Ṽ, a F V,Ṽ a is uniformly continuous on every bounded sets of Lip γ R e ; R d Lip γ R e ; R d R e R d by construction. Thirdly, maps π A, ;. and V π V, x ; w are respectively uniformly continuous on every bounded sets of by Theorem 2.1 and its remark. C 1-var, T ; LR e R e and Lip γ R e ; R d Therefore, by composition, the map V, Ṽ D Ṽ xv is uniformly continuous on every bounded sets of Lip γ R e ; R d Lip γ R e ; R d. In conclusion, by 9, Proposition B.5, the map V π V, x ; w is continuously differentiable from Lip γ R e ; R d into C 1-var, T ; R e. Theorem 2.4. Consider W GΩ p,t R d :

10 1 NICOLAS MARIE 1 Let V := V 1,..., V d be a collection of γ-lipschtiz vector fields on R e. The map x π V, x ; W is continuously differentiable from R e into C p-var, T ; R e. For every t, T, the Jacobian matrix of π V,.; W t at point x R e is denoted by J x,w t. Moreover, for every ε >, there exists a constant C 1 > only depending on p, γ, ε and V lip γ, such that for every x R e, J x,w. ;T C 1 e C1M ε,i,pw. 2 For every t, T, J x,w t is an invertible matrix. Moreover, for every ε >, there exists a constant C 2 > only depending on p, γ, ε and V lip γ, such that for every x R e, J x,w. 1 ;T C 2 e C2M ε,i,pw. 3 Consider x R e. The map V π V, x ; W is continuously differentiable from Lip γ R e ; R d into C p-var, T ; R e. Moreover, for every R > and V, Ṽ B Lipγ, R, there exists two constants η > and C 3 >, depending on R continuously and not on W, such that : V π V, x ; W.Ṽ ;T C 3 e C3M η,i,pw. Proof. Refer to 9, theorems for a proof of the Itô map s continuous differentiability with respect to the initial condition, and refer to 2, Corollary 3.4 for the upper-bound provided at the first point for J x,w. ;T ; x R e. Let I be the identity matrix of M e R. Proofs of points 1 and 2 are similar because if w :, T R d is a continuous function of finite 1-variation, then J x,w t J x,w t 1 = I t = I + DV π V, x ; W s, J x,w s dw s and t as mentioned at the proof of 9, Proposition DV π V, x ; W s, J x,w s 1 dw s The proof of the third point is detailed. In a first step, the continuous differentiability of the Itô map with respect to the collection of vector fields is established. In a second step, in order to apply integrability results coming from 2, a tailor-maid upper-bound for the derivative of the Itô map with respect to V is provided. Step 1. Since W GΩ p,t R d, there exists a sequence w n, n N of functions belonging to C 1-var, T ; R d and satisfying : 5 lim d p-var;t Sp w n,., W =. n Consider n N, W n := S p w n,., x R e, a := x,, X := 1, a,..., a p T p R e+1 p!

11 SENSITIVITIES VIA ROUGH PATHS 11 and V, Ṽ Lipγ R e ; R d. By Proposition 2.3, the map π., x ; w n is continuously differentiable from Lip γ R e ; R d into C 1-var, T ; R e. In particular, V π V, x ; w n.ṽ = ϕwn, V, Ṽ with where, ϕ., V, Ṽ := π A, ;. J.; F V, Ṽ π F V, X ;. A : R e LLR e R e ; R e, F V, Ṽ : Re R d LR e R d ; LR e R e and F V : R e LR d ; R e R d are three collections of vector fields, respectively defined by : AaL, b := L.a + b, F V, Ṽ a, a b, b := DV a,. b ; Ṽ ab and F V ab := V ab, b for every a, b R e, a, b R d and L LR e. Consider ε, 1. By Taylor s formula applied to π., x ; W n between V and V + εṽ, and 9, Definition 1.17 RDE s solutions : 6 π V +ε Ṽ, x ; W π V, x ; W = lim n uniformly. ε ϕw n, V + θṽ, Ṽ dθ Via Lebesgue s theorem and 9, Proposition B.1, let show that the derivative of π., x ; W at point V, in the direction Ṽ, exists in Cp-var, T ; R e equipped with the norm. p-var;t and coincides with ϕw, V, Ṽ. On one hand, the continuity results of Theorem 2.1 imply that : θ, 1, ϕw n, V + θṽ, Ṽ n in C p-var, T ; R e equipped with. ;T. ϕw, V + θṽ, Ṽ On the other hand, by applying successively 9, theorems 1.47 and 1.36, for every θ, 1 and every s, t T, ω 1/p 1 s, t; n; θ := F V +θ Ṽ,Ṽ π, X FV +θ Ṽ ; W n dπ, X FV +θ Ṽ ; W n p-var;s,t with and ω 1/p 2 s, t; n ω 1/p 2 s, t; n := ω 1/p 3 s, t; n ω 3 s, t; n ω p 3 s, t; n ω 3 s, t; n := η 1 W n p p-var;s,t where η 1 > is depending on V and Ṽ, but not on Wn and θ.

12 12 NICOLAS MARIE By 9, Exercice 1.55, there exists a constant C 4 >, not depending on W n and θ, such that : D 1 ϕw n, V + θṽ, Ṽ ;T C 4 exp C 4 ω 2 r k, r k+1 ; n = C 4 exp C 4 D = {r k } D I ω 2 rk, r k+1 ; n 1 D = {r k } D I ω 3 rk, r k+1 ; n 1 k=1 D 1 k=1 because ω 2.; n ω 3.; n when ω 2.; n 1. By er-additivity of the control ω 3.; n : ϕw n, V + θṽ, Ṽ ;T C 4 e η1c4 Wn p p-var;t. ω 3 r k, r k+1 ; n, In the right-hand side of that inequality, since η 1 and C 4 are not depending on W n and θ, and since W n p p-var;t < n N by 5 ; ϕw n, V + θṽ, Ṽ ;T < θ,1 n N in C p-var, T ; R e equipped with. ;T. Therefore, by Lebesgue s theorem and inequality 6 : π V +ε Ṽ, x ; W π V, x ; W = Since θ ϕw, V + θṽ, Ṽ is continuous from ε ϕw, V + θṽ, Ṽ dθ., 1 into C p-var, T ; R e equipped with. p-var;t by Theorem 2.1 ; by 9, Proposition B.1, the derivative of π., x ; W at point V, in the direction Ṽ, exists in Cp-var, T ; R e equipped with. p-var;t and coincides with ϕw, V, Ṽ. Finally, as at the second step of the proof of Proposition 2.3, via 9, Proposition B.5 and Lemma 4.2 ; the map π., x ; W is continuously differentiable from Lip γ R e ; R d into C p-var, T ; R e. Step 2. Consider R > and V, Ṽ B Lipγ, R. By applying successively 9, theorems 1.47 and 1.36, for every s, t T, ω 1/p 4 s, t := F V, Ṽ π F V, X ; W dπ FV, X ; W ω 1/p 5 s, t with ω 1/p 5 s, t := ω 1/p 6 s, t ω 6 s, t ω p 6 s, t and ω 6 s, t := η 2 W p p-var;s,t where η 2 > is depending on R continuously, but not on W. p-var;s,t

13 SENSITIVITIES VIA ROUGH PATHS 13 By 9, Exercice 1.55, there exists a constant C 5 >, not depending on R and W, such that : D 1 V π V, x ; W.Ṽ C 5 exp C 5 ω 5 r k, r k+1 ;T = C 5 exp C 5 D = {r k } D I ω 5 rk, r k+1 1 D = {r k } D I ω 6 rk, r k+1 1 k=1 D 1 k=1 ω 6 r k, r k+1, because However, Therefore, D = {r k } D I ω 6 rk, r k+1 1 ω 5 ω 6 when ω 5 1. D 1 k=1 V π V, x ; W.Ṽ ω 6 r k, r k+1 = η 2 M η 1 2,I,pW. C 3 e C3M η,i,pw ;T with C 3 := C 5 1 η 2 and η := η 1 2. Notations. In the sequel, matrices J x,w t and J x,w t 1 will be respectively denoted by J W t and J W t in a sake of simplicity. Moreover, for every s, t T, let s put : Then, J W s t := J W s J W t and J W t s := J W t J W s. J W s tj W t s = J W t sj W s t = I. At the following corollary, upper-bounds provided at the previous theorem are extended to RDEs having a drift term. Corollary 2.5. Consider m N, p > q 1 such that 1/p + 1/q > 1, h :, T R m a continuous function of finite q-variation, W GΩ p,t R d and W h := S p W h : 1 Let V := V 1,..., V d+m be a collection of γ-lipschitz vector fields on R e. For every ε >, there exists a constant C 1 > depending only on p, q, γ, ε and V lip γ, such that for every x R e, J Wh. ;T C 1 exp C 1 h p q-var;t + M ε,i,pw. 2 Consider x R e. For every R > and V, Ṽ B Lipγ, R, there exists two constants ε > and C 2 >, depending on R but not on h and W, such that : V π V, x ; W h.ṽ ;T C 2 exp C 2 h p q-var;t + M ε,i,pw. Proof. By Corollary 2.2, there exists a constant C 3 >, depending only on p and q, such that for every ε >, M ε,i,p W h C 3 h p q-var;t + M ε,i,pw. Therefore, by Theorem 2.4 :

14 14 NICOLAS MARIE 1 For a given collection of vector fields V Lip γ R e ; R d+m ; for every ε >, there exists a constant C 4 > depending only on p, γ, ε and V lip γ, such that for every x R e, J Wh. ;T C 4 e C4M ε,i,pw h C 1 exp C 1 h p q-var;t + M ε,i,pw with C 1 := C 4 1 C 3. 2 For a given initial condition x R e ; for every R > and V, Ṽ B Lip γ, R, there exists two constants ε > and C 5 >, depending on R but not on W h, such that : V π V, x ; W h.ṽ ;T C 5 e C5M ε,i,pw h C 2 exp C 2 h p q-var;t + M ε,i,pw with C 2 := C 5 1 C Differentiability of the Itô map with respect to the driving signal. Let s first remind the notion of differentiability introduced by P. Friz and N. Victoir on GΩ p,t R d. Definition 2.6. Consider a Banach space F, p > q 1 such that 1/p + 1/q > 1, and an open set U of GΩ p,t R d. The map ϕ : GΩ p,t R d F is continuously differentiable in the sense of Friz-Victoir on U if and only if, for every Y U, is continuously differentiable. h C q-var, T ; R d ϕt h Y F With notations of Definition 2.6, if ϕ is continuously differentiable from U into F in the sense of Friz-Victoir, in particular : Y U, ψ Y : h C q-var, T ; R d ψ Y h = ϕt h Y is derivable at every points and in every directions of C q-var, T ; R d. Notation. For every continuous functions h :, T R d of finite q-variation, D FV h ϕy := D h ψ Y ϕt εh Y ϕt Y = lim. ε ε In the sequel, D FV is called the Friz-Victoir directional derivative operator. Theorem 2.7. Consider a collection V := V 1,..., V d of γ-lipschitz vector fields on R e and x R e. The map W π V, x ; W is continuously differentiable from in the sense of Friz-Victoir. GΩ p,t R d into C p-var, T ; R e Moreover, for every W GΩ p,t R d and every continous function h :, T R d of finite q-variation, D F V h π V, x ; W =. J Ẉ sv π V, x ; W s dh s.

15 SENSITIVITIES VIA ROUGH PATHS 15 Duhamel s principle. Finally, consider W GΩ p,t R d and a control ω : T R + satisfying : Then, s, t T, W p-var;s,t ω 1/p s, t. 1 There exists a constant C 1 >, not depending on W and ω, such that for every continous function h :, T R d of finite q-variation, Dh FV π V, x ; W ;T C 1 exp C 1 h p q-var;t + M 1,I,ω. 2 There exists a constant C 2 >, not depending on W and ω, such that for every continous function h :, T R d of finite q-variation, Dh FV π V, x ; W p-var;t C 2 exp C 2 h pq-var;t + ω, T. Proof. Refer to 9, theorems and Exercice 11.9 for a proof of the first part. Consider a continuous function h :, T R d of finite q-variation, W h := S p W h, a := x,, and X := 1, a,..., a p T p R e+2. p! By 9, Theorem 11.3 : where, D FV h π V, x ; W = π A, ;. J., F π G, X ;.W h A : R e LLR e R e ; R e, F : R e R d R d LR e R d R d ; LR e R e and G : R e LR d R d ; R e R d R d are three collections of vector fields, respectively defined by : AaL, b := L.a + b, F a, a, a b, b, b := DV a,. b ; V ab and Gab, b := V ab, b, b for every a, b R e, a, b, a, b R d and L LR e. By applying successively 9, theorems 1.47 and 1.36, for every s, t T, ω 1/p 1 s, t := F π G, X ; W h dπ G, X ; W h ω 1/p 2 s, t with ω 1/p 2 s, t := ω 1/p 3 s, t ω 3 s, t ω p 3 s, t and, by 9, Proposition 7.52 : p-var;s,t ω 3 s, t := ε 1 h p q-var;s,t + ωs, t ε 2 W h p p-var;s,t where, ε 1, ε 2 1 are two constants not depending on W, ω and h.

16 16 NICOLAS MARIE On one hand, by 9, Exercice 1.55, there exists a constant C 3 >, not depending on W, ω and h, such that : D h FV π V, x ; W ;T C 3 exp C 3 with C 1 := C 3 ε 1, because = C 3 exp C 3 C 1 exp D = {r k } D I ω 2 rk, r k+1 1 D = {r k } D I ω 3 rk, r k+1 1 D 1 k=1 D 1 k=1 C 1 h p q-var;t + M 1,I,ω 7 ω 2 ω 3 when ω 2 1 and s, t T, ωs, t ω 3 s, t. ω 2 r k, r k+1 ω 3 r k, r k+1 On the other hand, by 9, Theorem 1.53, there exists a constant C 4 >, not depending on W, ω and h, such that for every s, t T satisfying ω 2 s, t 1, Dh FV π V, x ; W s,t C D FV h π V, x ; W s ω 1/p 2 s, te C4ω2s,t by 7. C D FV h π V, x ; W ;T ω 1/p 3 s, te C4ω3,T, Therefore, by er-additivity of the control ω 3, there exists a constant C 2 >, not depending on W, ω and h, such that : Dh FV π V, x ; W p-var;t C 2 exp C 2 h pq-var;t + ω, T. At the following corollary, the upper-bound provided at the previous theorem are extended to RDEs having a drift term. Corollary 2.8. Consider m N, p > q 1 such that 1/p + 1/q > 1, r 1, p such that 1/p+1/r > 1, g :, T R m a continuous function of finite r-variation, W GΩ p,t R d, W g := S p W g, V := V 1,..., V d+m a collection of γ-lipschitz vector fields on R e and x R e. Then, there exists a constant C >, not depending on g and W, such that for every continuous function h :, T R d+m of finite q-variation, Dh FV π V, x ; W g ;T C exp C h p q-var;t + g p r-var;t + M 1,I,pW. Proof. Let h :, T R d+m be a continuous function of finite q-variation. By Corollary 2.2, there exists a constant C 1 >, depending only on p and r, such that : M 1,I,p W g C 1 g p r-var;t + M 1,I,pW. Then, by Theorem 2.7, there exists a constant C 2 >, not depending on W g and h, such that : Dh FV π V, x ; W g ;T C 2 exp C 2 h p q-var;t + M 1,I,pW g C exp C h p q-var;t + g p r-var;t + M 1,I,pW with C := C 2 1 C 1.

17 SENSITIVITIES VIA ROUGH PATHS 17 Finally, let study the second order continuous differentiability of the Itô map with respect to the driving signal. In order to apply Fernique s theorem at Section 3, tailor-maid upper-bounds are provided for p 1, 2. Lemma 2.9. Consider p 1, 2, x R e and A := A 1,..., A d the collection of affine vector fields on R e defined by : L LR e, a, b R e, AaL, b := L.a + b. The map L, w π A, x ; L, w is continuously differentiable from C p-var, T ; LR e R e into C p-var, T ; R e. Moreover, consider a continuous map L, w :, T LR e R e of finite p- variation, and ω : T R + a control satisfying : s, t T, L, w p p-var;s,t ωs, t. Then, 1 There exists a constant C 1 >, not depending on w L := L, w and ω, such that for every continuous map h H := H, h :, T LR e R e of finite p-variation, w Lπ A, x ; w L.h H ;T C x exp C 1 h H p p-var;t + M 1,I,ω. 2 There exists a constant C 2 >, not depending on w L := L, w and ω, such that for every continuous map h H := H, h :, T LR e R e of finite q-variation, w Lπ A, x ; w L.h H p-var;t C x exp C 2 h H pp-var;t + ω, T. Proof. Since A is a collection of affine vector fields on R e, they are locally γ- Lipschitz on R e. Then, by 9, Theorem 1.3, π A, x ;. is derivable at every points and in every directions of C p-var, T ; LR e R e. Consider a continuous map h H := H, h :, T LR e R e of finite p-variation. Firstly, D h H π A, x ; w L = =.. dl s.d h H π A, x ; w L s +. dh s.π A, x ; w L s + h A 1 πa, x ; w L s ; D h H π A, x ; w L s dw L s, dh H s where, A 1 : R e R e LLR e R e LR e R e ; R e is the collection of affine vector fields defined by : A 1 a, bc, d := c 1.b + d 1.a + d 2 for every a, b R e and c, d LR e R e. By putting Qw L, h H := π A, x ; w L ; D h H π A, x ; w L : Qw L, h H = x, +. A 2 Qs w L, h H dw L s, dh H s where, A 2 : R e R e LLR e R e LR e R e ; R e R e is the collection of affine vector fields defined by : Aac A 2 a, bc, d := A 1 a, bc, d c = c a c 2.b d 1 +.a d 2

18 18 NICOLAS MARIE for every a, b R e and c, d LR e R e. Then, Qw L, h H = π A2, x, ; w L, h H. Secondly, there exists a constant C 3 1, depending only on p, such that for every s, t T, w L, h H p p-var;s,t ω 1 s, t := C 3 h H p p-var;s,t + ωs, t. By construction, ω 1 : T R + is a control. By 9, Exercice 1.55, there exists a constant C 4 >, not depending on w L, ω and h H, such that : D h H π A, x ; w L ;T with C 1 = C 3 C 4, because Qw L, h H ;T C x exp C 4 D = {r k } D I ω 1 rk, r k+1 1 C x expc 1 h H p p-var;t + D = {r k } D I ω r k, r k+1 1 D 1 k=1 ωr k, r k+1 s, t T, ωs, t ω 1 s, t. D 1 k=1 ω 1 r k, r k+1 Thirdly, by 9, Theorem 1.53, there exists a constant C 5 >, not depending on w L, ω and h H, such that for every s, t T, D h H π A, x ; w L s,t C x ω 1/p 1 s, te C5ω1,T. Therefore, by er-additivity of the control ω 1, there exists a constant C 2 >, not depending on w L, ω and h H, such that : D h H π A, x ; w L p-var;t C x exp C 2 h H pp-var;t + ω, T. Finally, by the remark following Theorem 2.1, since Q = π A2, x, ;. and A 2 is a collection of affine vector fields on R e, Q and then w L, h H D h H π A, x ; w L are uniformly continuous on bounded set of C p-var, T ; LR e R e C p-var, T ; LR e R e. In conclusion, by 9, Proposition B.5, π A, x ;. is continuously differentiable as stated. Proposition 2.1. Consider 2 > p q 1, x R e and V, Ṽ Lipγ+1 R e ; R d. The maps { D 1 C : p-var, T ; R d C p-var, T ; R e e w D 1 w = J ẉ, { D 2 C : p-var, T ; R d C p-var, T ; R e e w D 2 w = J. w and { C D 3 p-var, T ; R d C p-var, T ; R e : w D 3 w = V π V, x ; w.ṽ are continuously differentiable. Moreover, let w :, T R d be a continuous map of finite p-variation. there

19 SENSITIVITIES VIA ROUGH PATHS 19 exists three strictly positive constants C 1, C 2 and C 3, not depending on w, such that for every continuous function h :, T R d of finite q-variation, w D 1 w.h p-var;t C 1 exp C 1 h pq-var;t + w pp-var;t, w D 2 w.h p-var;t C 2 exp C 2 h pq-var;t + w pp-var;t and w D 3 w.h p-var;t C 3 exp C 3 h pq-var;t + w pp-var;t. Proof. Proofs are similar for D 1, D 2 and D 3. It is detailed for D 3. At the proof of Theorem 2.4, with the same notations, it has been established that : where, D 3 w = I J Kw I := π A, ;., J := J.; F V, Ṽ and K := π F V, x, ;.. When p 1, 2, since F V and F V, Ṽ are respectively two collections of γ-lipschitz vector fields on R e and R e R d by construction, by Theorem 2.7 : J : C p-var, T ; R e R d C p-var, T ; LR e R e and are continuously differentiable. K : C p-var, T ; R d C p-var, T ; R e R d Moreover, by Proposition 2.9, I is a continuously differentiable map from C p-var, T ; LR e R e into C p-var, T ; R e. Then, by composition, D 3 is a continuously differentiable map from C p-var, T ; R d into C p-var, T ; R e. Let h :, T R d be a continuous function of finite q-variation : Then, w D 3 w.h = DIJ Kw; DJKw, w Kw.h. 8 w D 3 w.h p-var;t DIJ Kw L;p DJKw L;p where,. L;p is the usual operator norm on or on w Kw.h p-var;t LC p-var, T ; LR e R e ; C p-var, T ; R e LC p-var, T ; R e R d ; C p-var, T ; LR e R e. Let s find upper-bounds for each term of the product of the right-hand side of inequality 8 : 1 Since Kw = π FV, x ; w, by the second point of Theorem 2.7 for the control ω w,p, there exists a constant C 4 >, not depending on w and h, such that : 9 w Kw.h p-var;t C 4 exp C 4 h pq-var;t + w pp-var;t.

20 2 NICOLAS MARIE 2 Let g :, T R e R d be a continuous function of finite p-variation such that g p-var;t 1. On one hand, by construction of the rough integral : where, J Kw = π LR e R eπ G V, Ṽ, x, ; Kw G V, Ṽ : Re R d 2 LR e R d ; R d R e LR e R e is the collection of γ-lipschitz vector fields as F V, Ṽ defined by : a, b, c R e R d, G V, Ṽ a, bc = c, F V,Ṽ a.c. On the other hand, by 9, Theorem 1.36, for every s, t T, Kw p-var;s,t ω 1/p 1 s, t := ω 1/p 2 s, t ω 2 s, t with ω 2 s, t := ε 1 w p p-var;s,t where, ε 1 1 is a constant not depending on w. However, for every ε 1, 1 ω 3.; ε ε g p p-var + ω 2 when ω 3.; ε 1 where, for every s, t T, and ω 1/p 3 s, t; ε := ω 1/p 4 s, t; ε ω 4 s, t; ε ω p 4 s, t; ε ω 4 s, t; ε := ε g p p-var;s,t + ω 1 s, t. Then, as at the proof of point 2 of Theorem 2.7, by using 1 ; there exists a constant C 5 >, not depending on w and g, such that : DJKw.g p-var;t C 5 expc 5 g p p-var;t + w p p-var;t. Therefore, 11 DJKw L;p C 5 expc w p p-var;t. 3 Let g H := H, g :, T LR e R e be the continuous map of finite p-variation such that g p-var;t 1. On one hand, it has been established at Lemma 2.9 that there exists a collection A 1 of affine vector fields on R e1 R e2 e 1 = e 2 = e such that : DIJ Kw.g H = π R e 2 π A1,, ; J Kw, g H. On the other hand, by applying successively 9, theorems 1.36 and 1.47, for every s, t T, J Kw p-var;s,t ω 1/p 5 s, t := ω 1/p 6 s, t ω 6 s, t ω p 6 s, t with ω 6 s, t := ε 2 w p p-var;s,t where, ε 2 1 is a constant not depending on w. Then, as at the proof of point 2 of Lemma 2.9 ; there exists a constant C 6 >, not depending on w and g H, such that : DIJ Kw.g H p-var;t C 6 expc 6 g H p p-var;t + w p p-var;t. Therefore, 12 DIJ Kw L;p C 6 expc w p p-var;t.

21 SENSITIVITIES VIA ROUGH PATHS 21 In conclusion, inequalities 9, 11, 12 and 8 imply together that there exists a constant C 3 >, not depending on w and h, such that : w D 3 w.h p-var;t C 3 exp C 3 h pq-var;t + w pp-var;t Application to Gaussian stochastic analysis. Consider a d-dimensional stochastic process W and the probability space Ω, A, P, where Ω is the vector space of continuous functions from, T into R d, A is the σ-algebra generated by cylinder sets of Ω, and P is the probability measure induced by the process W on Ω, A. In order to prove Corollary 2.17 that is crucial at Section 3, let s first state existing results on Gaussian rough paths established by P. Friz and N. Victoir in 8, and by T. Cass, C. Litterer and T. Lyons in 2. Consider the two following technical assumptions on the stochastic process W. Assumption W is a d-dimensionnel centered Gaussian process with continuous paths. Moreover, its covariance function c W is of finite 2D ρ-variation with ρ 1, 2 cf. 9, Definition 5.5. Assumption There exists p > q 1 such that : 1 p + 1 q > 1 and H1 C q-var, T ; R d. Example. By 9, Proposition 15.5, Proposition 15.7 and Exercice 2.2, a fractional Brownian motion of Hurst parameter H 1/4, 1/2 satisfies assumptions 2.11 and Theorem Consider a stochastic process W satisfying Assumption 2.11, and p > 2ρ. For almost every ω Ω, there exists a geometric p-rough path Wω over W ω satisfying : 1 There exists a deterministic constant C >, only depending on ρ, p and c W ρ-var;,t 2, such that : E e C W 2 p-var;t <. generalized Fernique s theorem. 2 Let W n, n N be a sequence of linear approximations, or of mollifier approximations, of the process W. Then, W is the limit in p-variation, in L r Ω for every r 1, of the sequence S 3 W n, n N universality. W is the enhanced Gaussian process associated to W. Refer to 9, Theorem for a proof. Proposition Consider a stochastic process W satisfying assumptions 2.11 and 2.12, W the enhanced Gaussian process associated to W, and the Cameron- Martin s space H 1 Ω of the process W. Then, Refer to 9, Lemma for a proof. ω Ω, h H 1, Wω + h = T h Wω. Proposition For every geometric p-rough path Y and every ε >, Refer to 2, Proposition 4.6 for a proof. M ε,i,p Y ε 2N ε,i,p Y + 1.

22 22 NICOLAS MARIE Theorem Consider a stochastic process W satisfying assumptions 2.11 and 2.12, and W the enhanced Gaussian process associated to W. Then, C, ε, r >, Ce CN ε,i,pw L r Ω. Refer to 2, Theorem 6.4 and Remark 6.5 for a proof. Corollary Consider x R e, V := V 1,..., V d+1 and Ṽ := Ṽ1,..., Ṽd+1 two collections of γ-lipschitz vector fields on R e, a stochastic process W satisfying assumptions 2.11 and 2.12, W the enhanced Gaussian process associated to W, W g := S p W g with g := Id,T, and a continuous function h :, T R d+1 of finite q-variation. Then, J Wg. ;T, V π V, x ; W g.ṽ ;T and Dh FV π V, x ; W g ;T belong to L r Ω for every r >. Proof. It is a straightforward consequence of corollaries 2.5 and 2.8, of Proposition 2.15 deterministic results, and of Theorem 2.16 probabilistic result. 3. Sensitivity analysis for Gaussian rough differential equations This section solves the problematic stated in introduction of the paper by using the deterministic results on RDEs of subsections 2.1 and 2.2, the probabilistic results on Gaussian RDEs of Subsection 2.3 and Malliavin calculus arguments. Assume that W, µ and σ defined in introduction satisfy the following condition. Assumption 3.1. The process W satisfies assumptions 2.11 and 2.12, and, T ; R d H 1. Moreover, there exists a constant C > such that : h C 1, T ; R d, h H 1 C ḣ ;T. C 1 Functions µ and σ satisfy Assumption 1.3 and, for every a R d, σa is an invertible matrix. Moreover, the function σ 1 : R d M d R is bounded. Example. A fractional Brownian motion B H of Hurst parameter H 1/4, 1 satisfies Assumption 3.1. Indeed, it as been stated at Subsection 2.3 that B H satisfies assumptions 2.11 and Moreover, by the first point of L. Decreusefond and S. Ustunel 4, Theorem 3.3 :, T ; R d H 1. C 1 Consider h C 1, T ; R d. By the second point of 4, Theorem 3.3 : h H 1 = J H ḣ H 1 = ḣ L 2,T T 1/2 ḣ ;T. Assume also that the function F : R d R satisfies one of the two following assumptions. Assumption 3.2. The function F is continuously differentiable from R d into R. Moreover, there exists two constants C > and N N such that, for every a R d, F a C1 + a N and DF a L C1 + a N. Assumption 3.3. There exists two constants C > and N N such that, for every a R d, F a C1 + a N.

23 SENSITIVITIES VIA ROUGH PATHS 23 The following results are solving, at least partially, the problematic stated in introduction of the paper. Notations : Under Assumption 3.1, the enhanced Gaussian process associated to W is denoted by W, W g := S p W g with g := Id,T, and V := V 1,..., V d+1 is the collection of vector fields defined by : V ab, c := µac + σab for every a, b R d and c R. Let S p Lip γ R d ; R d be the space of functions from R d into M d R, p+1 times differentiable, bounded and of bounded derivatives. For every x R d, EF X T is denoted by f T x, σ. Lemma 3.4. Let I = I 1,..., I d be the map from H into H 1 such that : I i h := E W i h i W i H 1 for every h H = H 1 H d and i = 1,..., d. Then, I is an isometry from H into H 1. Proof. On one hand, the linearity of I as a map from H into H 1 is a straightforward consequence of the linearity of W as a map from H into L 2 Ω. On the other hand, by construction of W and of the scalar products on H and H 1 : d Ih, Ig H 1 = E W i h i W i ; E W i g i W i H 1 i for every functions h, g H. = i=1 d E W i h i W i g i = h, g H i=1 The following corollary extends 9, Proposition 2.5 to Gaussian RDEs having a drift term. Lemma 3.5. For every x R d and almost every ω Ω, the map h π V, x ; W g ω+ h is continuously differentiable from H 1 into C p-var, T ; R d. In particular, for every t, T, π V, x ; W g t D 1,2 loc and for every h H1, Dπ V, x ; W g t, I 1 h H = Dh, FV π V, x ; W g t = t J Wg t sσ π V, x ; W g s dh s. Proof. By Proposition 2.14, for almost every ω Ω and every h H 1, Then, almost surely : W g ω + h = S p Wω + h g = S p T h Wω g = T h, S p Wω g = T h, W g ω. 13 π V, x ; W g. + h = π V, x ; T h, W g.

24 24 NICOLAS MARIE Moreover, Assumpution 2.12 and Corollary 2.8 imply that h π V, x ; T h, W g is continuously differentiable from H 1 C q-var, T ; R d into C p-var, T ; R d. Then, by equality 13, the map h π V, x ; W g. + h is also continuously differentiable from H 1 into C p-var, T ; R d, and for almost every ω Ω and every h H 1, with F ω := π V, x ; W g ω +.. D FV h, π V, x ; W g ω = D h F ω Moreover, by Duhamel s principle cf. Theorem 2.7, for every t, T and every h H 1, D FV h, π V, x ; W g t = = t t J Wg t sv π V, x ; W g s dh s, J Wg t sσ π V, x ; W g s dh s. In conclusion, by 2, Proposition and Lemma 4.1.2, π V, x ; W g t is continuously H 1 -differentiable and then locally derivable in the sense of Malliavin, with : Dπ V, x ; W g t, I 1 h H = D h F ω. Theorem 3.6. Under assumptions 3.1 and 3.2 : 1 The function f T., σ is differentiable from R d into R and, for every x, v R d, x f T x, σ.v = E DF X x T, I 1 h x,v H where h x,v := 1 T. σ 1 X x s J Wg s vds. 2 For every x R d, the function f T x,. is differentiable from S p into R and, for every σ, σ S p satisfying Assumption 3.1, σ f T x, σ. σ = E DF X σ T, I 1 h σ, σ H where h σ, σ := 1 T. σ 1 X σ s J Wg s T σ X σ T. σds. Proof. Proofs of points 1 and 2 are similar : 1 On one hand, for every ε, 1, η > and x, v R d, by Taylor s formula, and the first point of Corollary 2.5 ; there exists a constant C 1 > depending only on p, γ, η and V lip γ such that : F X x+εv T F X x T ε = 1 DF X x+θεv T C 1 v e C1M η,i,pw, DX x+θεv T 1.v dθ DF X x+θεv T L dθ. Moreover, since F satisfies Assumption 3.2, there exists two constants C 2 > and N N, depending only on F, such that for every θ, 1, DF X x+θεv L T C X x+θεv N T.

25 SENSITIVITIES VIA ROUGH PATHS 25 Then, by the triangle inequality together with 9, Theorem 1.36, there exists a constant C 3 >, not depending on x, v, θ, ε, V and W, such that : DF X x+θεv T L C x + v + C 3 V lip γ 1 W g p-var;t V p lip γ 1 W g p p-var;t N. Since W satisfies assumptions 2.11 and 2.12, by Corollary 2.17, the generalized Fernique s theorem cf. Theorem 2.13 and Cauchy-Schwarz s inequality : ε, 1 F X x+εv T F X x T ε is bounded by an integrable random variable not depending on ε. Therefore, by Lebesgue s theorem, f T., σ is differentiable on R d and 14 x, v R d, x f T x, σ.v = E DF X x T, x X x T.v. On the other hand, consider x, v R d. By construction, paths of the process h x,v are continuously differentiable from, T into R d and h x,v =. Then, since W satisfies Assumption 3.1, h x,v is a H 1 -valued random variable. By Duhamel s principle cf. Theorem 2.7 : D FV h x,v, Xx T = T = x X x T.v. JT Wg sσ Xs x dh x,v s Therefore, by equality 14, Lemma 3.5 and 2, Proposition Malliavin derivative s chain rule : x f T x, σ.v = EDF X x T.D FV h x,v, Xx T = EDF X x T. DX x T, I 1 h x,v H = E DF X x T, I 1 h x,v H. 2 Let x R d be fixed. On one hand, for every ε, 1 and σ, σ S p, by Taylor s formula : F X σ+ε σ T F XT σ 1 = ε DF X σ+θε σ T, DX σ+θε σ T. σ dθ 1 C X σ+θε σ N T DX σ+θε σ. σ ;T dθ. At 9, Theorem 1.36, the constant involving in the upper-bound doesn t depend on the signal and on the collection of vector fields. At the second point of Corollary 2.5, the two constants involving in the upper-bound continuously depend on the γ-lipschitz norm of the collection of vector fields. Then, there exists a constant C 4 >, depending on σ and σ but not on ε and W, such that for every θ, 1, DX σ+θε σ. σ ;T C 4 e C4M C 4,I,pW and X σ+θε σ T C 4 W g p-var;t W g p p-var;t. Since W satisfies assumptions 2.11 and 2.12, by Proposition 2.15, Theorem 2.16, the generalized Fernique s theorem cf. Theorem 2.13 and Cauchy- Schwarz s inequality : ε, 1 F X σ+ε σ T F X σ T ε

26 26 NICOLAS MARIE is bounded by an integrable random variable not depending on ε. Therefore, by Lebesgue s theorem, f T x,. is differentiable on S p and 15 σ, σ S p, σ f T x, σ. σ = E DF X σ T, σ X σ T. σ. Remarks : On the other hand, consider σ, σ S p satisfying Assumption 3.1. By construction, paths of the process h σ, σ are continuously differentiable from, T into R d and h σ, σ =. Then, h σ, σ is a H 1 -valued random variable. By Duhamel s principle cf. Theorem 2.7 : D FV h σ, σ, Xσ T = T = σ X σ T. σ. JT Wg sσ Xs σ dh σ, σ s Therefore, by equality 15, Lemma 3.5 and 2, Proposition Malliavin derivative s chain rule : σ f T x, σ. σ = EDF X σ T.D FV h σ, σ, Xσ T = EDF X σ T. DX σ T, I 1 h σ, σ H = E DF X σ T, I 1 h σ, σ H. 1 Note that the differentiability of partial maps of f T only requires assumptions 2.11, 2.12, 1.3 and Assumption 2.11 only guarantees the existence of the enhanced Gaussian process W at Theorem Therefore, when p 1, 2, Theorem 3.6 holds under the following assumption. Assumption 3.7. The process W is d-dimensional, Gaussian, centered and with continuous paths of finite p-variation p 1, 2. Moreover, T ; R d H 1 C p-var, T ; R d, C 1 and there exists a constant C > such that : h C 1, T ; R d, h H 1 C ḣ ;T. Functions µ and σ are p+2 times differentiable, bounded and of bounded derivatives and, for every a R d, σa is an invertible matrix. Moreover, the function σ 1 : R d M d R is bounded. Example. A fractional Brownian motion B H of Hurst parameter H 1/2, 1 satisfies Assumption 3.7. Indeed, by Kolmogorov s continuity criterion, there exists ε > and C T L 1 Ω; R + such that for every s, t T, 16 B H t B H s C T t s 1/p P-a.s. where p := 1/H ε 1, 2. Corollary 3.8. Under assumptions 3.7 and, 3.2 or 3.3 : 1 For every x, v R d, I 1 h x,v D 1,2 H and x f T x, σ.v = E F X x T δ I 1 h x,v. 2 For every x R d and every σ, σ S p+1 satisfying Assumption 3.7, I 1 h σ, σ D 1,2 H and σ f T x, σ. σ = E F X σ T δ I 1 h σ, σ.

Sensitivities via Rough Paths

Sensitivities via Rough Paths Sensitivities via Rough Paths Nicolas Marie To cite this version: Nicolas Marie. Sensitivities via Rough Paths. ESAIM: Probability and Statistics, EDP Sciences, 215, 19, pp.515-543. .

More information

Rough paths methods 4: Application to fbm

Rough paths methods 4: Application to fbm Rough paths methods 4: Application to fbm Samy Tindel Purdue University University of Aarhus 2016 Samy T. (Purdue) Rough Paths 4 Aarhus 2016 1 / 67 Outline 1 Main result 2 Construction of the Levy area:

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

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

MULTIDIMENSIONAL WICK-ITÔ FORMULA FOR GAUSSIAN PROCESSES

MULTIDIMENSIONAL WICK-ITÔ FORMULA FOR GAUSSIAN PROCESSES MULTIDIMENSIONAL WICK-ITÔ FORMULA FOR GAUSSIAN PROCESSES D. NUALART Department of Mathematics, University of Kansas Lawrence, KS 6645, USA E-mail: nualart@math.ku.edu S. ORTIZ-LATORRE Departament de Probabilitat,

More information

Malliavin Calculus: Analysis on Gaussian spaces

Malliavin Calculus: Analysis on Gaussian spaces Malliavin Calculus: Analysis on Gaussian spaces Josef Teichmann ETH Zürich Oxford 2011 Isonormal Gaussian process A Gaussian space is a (complete) probability space together with a Hilbert space of centered

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

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

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

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

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

From Fractional Brownian Motion to Multifractional Brownian Motion

From Fractional Brownian Motion to Multifractional Brownian Motion From Fractional Brownian Motion to Multifractional Brownian Motion Antoine Ayache USTL (Lille) Antoine.Ayache@math.univ-lille1.fr Cassino December 2010 A.Ayache (USTL) From FBM to MBM Cassino December

More information

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539 Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory

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

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

GENERALIZED COVARIATION FOR BANACH SPACE VALUED PROCESSES, ITÔ FORMULA AND APPLICATIONS

GENERALIZED COVARIATION FOR BANACH SPACE VALUED PROCESSES, ITÔ FORMULA AND APPLICATIONS Di Girolami, C. and Russo, F. Osaka J. Math. 51 (214), 729 783 GENERALIZED COVARIATION FOR BANACH SPACE VALUED PROCESSES, ITÔ FORMULA AND APPLICATIONS CRISTINA DI GIROLAMI and FRANCESCO RUSSO (Received

More information

Some SDEs with distributional drift Part I : General calculus. Flandoli, Franco; Russo, Francesco; Wolf, Jochen

Some SDEs with distributional drift Part I : General calculus. Flandoli, Franco; Russo, Francesco; Wolf, Jochen Title Author(s) Some SDEs with distributional drift Part I : General calculus Flandoli, Franco; Russo, Francesco; Wolf, Jochen Citation Osaka Journal of Mathematics. 4() P.493-P.54 Issue Date 3-6 Text

More information

CUBATURE ON WIENER SPACE: PATHWISE CONVERGENCE

CUBATURE ON WIENER SPACE: PATHWISE CONVERGENCE CUBATUR ON WINR SPAC: PATHWIS CONVRGNC CHRISTIAN BAYR AND PTR K. FRIZ Abstract. Cubature on Wiener space [Lyons, T.; Victoir, N.; Proc. R. Soc. Lond. A 8 January 2004 vol. 460 no. 2041 169-198] provides

More information

An Introduction to Malliavin Calculus. Denis Bell University of North Florida

An Introduction to Malliavin Calculus. Denis Bell University of North Florida An Introduction to Malliavin Calculus Denis Bell University of North Florida Motivation - the hypoellipticity problem Definition. A differential operator G is hypoelliptic if, whenever the equation Gu

More information

Evaluation of the HJM equation by cubature methods for SPDE

Evaluation of the HJM equation by cubature methods for SPDE Evaluation of the HJM equation by cubature methods for SPDEs TU Wien, Institute for mathematical methods in Economics Kyoto, September 2008 Motivation Arbitrage-free simulation of non-gaussian bond markets

More information

Itô s formula. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Itô s formula. Samy Tindel. Purdue University. Probability Theory 2 - MA 539 Itô s formula Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Itô s formula Probability Theory

More information

ROUGH PATH THEORY AND STOCHASTIC CALCULUS

ROUGH PATH THEORY AND STOCHASTIC CALCULUS Unspecified Journal Volume 00, Number 0, Pages 000 000 S????-????(XX)0000-0 ROUGH PATH THEORY AND STOCHASTIC CALCULUS YUZURU INAHAMA Abstract. T. Lyons rough path theory is something like a deterministic

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

COVARIANCE IDENTITIES AND MIXING OF RANDOM TRANSFORMATIONS ON THE WIENER SPACE

COVARIANCE IDENTITIES AND MIXING OF RANDOM TRANSFORMATIONS ON THE WIENER SPACE Communications on Stochastic Analysis Vol. 4, No. 3 (21) 299-39 Serials Publications www.serialspublications.com COVARIANCE IDENTITIES AND MIXING OF RANDOM TRANSFORMATIONS ON THE WIENER SPACE NICOLAS PRIVAULT

More information

Integral representations in models with long memory

Integral representations in models with long memory Integral representations in models with long memory Georgiy Shevchenko, Yuliya Mishura, Esko Valkeila, Lauri Viitasaari, Taras Shalaiko Taras Shevchenko National University of Kyiv 29 September 215, Ulm

More information

Malliavin calculus and central limit theorems

Malliavin calculus and central limit theorems Malliavin calculus and central limit theorems David Nualart Department of Mathematics Kansas University Seminar on Stochastic Processes 2017 University of Virginia March 8-11 2017 David Nualart (Kansas

More information

The Wiener Itô Chaos Expansion

The Wiener Itô Chaos Expansion 1 The Wiener Itô Chaos Expansion The celebrated Wiener Itô chaos expansion is fundamental in stochastic analysis. In particular, it plays a crucial role in the Malliavin calculus as it is presented in

More information

Some Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2)

Some Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2) Some Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2) Statistical analysis is based on probability theory. The fundamental object in probability theory is a probability space,

More information

White noise generalization of the Clark-Ocone formula under change of measure

White noise generalization of the Clark-Ocone formula under change of measure White noise generalization of the Clark-Ocone formula under change of measure Yeliz Yolcu Okur Supervisor: Prof. Bernt Øksendal Co-Advisor: Ass. Prof. Giulia Di Nunno Centre of Mathematics for Applications

More information

Topics in fractional Brownian motion

Topics in fractional Brownian motion Topics in fractional Brownian motion Esko Valkeila Spring School, Jena 25.3. 2011 We plan to discuss the following items during these lectures: Fractional Brownian motion and its properties. Topics in

More information

Laplace s Equation. Chapter Mean Value Formulas

Laplace s Equation. Chapter Mean Value Formulas Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic

More information

Applications of Ito s Formula

Applications of Ito s Formula CHAPTER 4 Applications of Ito s Formula In this chapter, we discuss several basic theorems in stochastic analysis. Their proofs are good examples of applications of Itô s formula. 1. Lévy s martingale

More information

High order Cubature on Wiener space applied to derivative pricing and sensitivity estimations

High order Cubature on Wiener space applied to derivative pricing and sensitivity estimations 1/25 High order applied to derivative pricing and sensitivity estimations University of Oxford 28 June 2010 2/25 Outline 1 2 3 4 3/25 Assume that a set of underlying processes satisfies the Stratonovich

More information

Wiener Measure and Brownian Motion

Wiener Measure and Brownian Motion Chapter 16 Wiener Measure and Brownian Motion Diffusion of particles is a product of their apparently random motion. The density u(t, x) of diffusing particles satisfies the diffusion equation (16.1) u

More information

i=1 α i. Given an m-times continuously

i=1 α i. Given an m-times continuously 1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable

More information

arxiv: v1 [math.pr] 11 Jan 2019

arxiv: v1 [math.pr] 11 Jan 2019 Concentration inequalities for Stochastic Differential Equations with additive fractional noise Maylis Varvenne January 4, 9 arxiv:9.35v math.pr] Jan 9 Abstract In this paper, we establish concentration

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

5 December 2016 MAA136 Researcher presentation. Anatoliy Malyarenko. Topics for Bachelor and Master Theses. Anatoliy Malyarenko

5 December 2016 MAA136 Researcher presentation. Anatoliy Malyarenko. Topics for Bachelor and Master Theses. Anatoliy Malyarenko 5 December 216 MAA136 Researcher presentation 1 schemes The main problem of financial engineering: calculate E [f(x t (x))], where {X t (x): t T } is the solution of the system of X t (x) = x + Ṽ (X s

More information

A Concise Course on Stochastic Partial Differential Equations

A Concise Course on Stochastic Partial Differential Equations A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original

More information

Discrete approximation of stochastic integrals with respect to fractional Brownian motion of Hurst index H > 1 2

Discrete approximation of stochastic integrals with respect to fractional Brownian motion of Hurst index H > 1 2 Discrete approximation of stochastic integrals with respect to fractional Brownian motion of urst index > 1 2 Francesca Biagini 1), Massimo Campanino 2), Serena Fuschini 2) 11th March 28 1) 2) Department

More information

Quasi-sure Stochastic Analysis through Aggregation

Quasi-sure Stochastic Analysis through Aggregation Quasi-sure Stochastic Analysis through Aggregation H. Mete Soner Nizar Touzi Jianfeng Zhang March 7, 21 Abstract This paper is on developing stochastic analysis simultaneously under a general family of

More information

Some Properties of NSFDEs

Some Properties of NSFDEs Chenggui Yuan (Swansea University) Some Properties of NSFDEs 1 / 41 Some Properties of NSFDEs Chenggui Yuan Swansea University Chenggui Yuan (Swansea University) Some Properties of NSFDEs 2 / 41 Outline

More information

µ X (A) = P ( X 1 (A) )

µ X (A) = P ( X 1 (A) ) 1 STOCHASTIC PROCESSES This appendix provides a very basic introduction to the language of probability theory and stochastic processes. We assume the reader is familiar with the general measure and integration

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

An introduction to Mathematical Theory of Control

An introduction to Mathematical Theory of Control An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018

More information

Stochastic calculus without probability: Pathwise integration and functional calculus for functionals of paths with arbitrary Hölder regularity

Stochastic calculus without probability: Pathwise integration and functional calculus for functionals of paths with arbitrary Hölder regularity Stochastic calculus without probability: Pathwise integration and functional calculus for functionals of paths with arbitrary Hölder regularity Rama Cont Joint work with: Anna ANANOVA (Imperial) Nicolas

More information

-variation of the divergence integral w.r.t. fbm with Hurst parameter H < 1 2

-variation of the divergence integral w.r.t. fbm with Hurst parameter H < 1 2 /4 On the -variation of the divergence integral w.r.t. fbm with urst parameter < 2 EL ASSAN ESSAKY joint work with : David Nualart Cadi Ayyad University Poly-disciplinary Faculty, Safi Colloque Franco-Maghrébin

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

Mean-field SDE driven by a fractional BM. A related stochastic control problem

Mean-field SDE driven by a fractional BM. A related stochastic control problem Mean-field SDE driven by a fractional BM. A related stochastic control problem Rainer Buckdahn, Université de Bretagne Occidentale, Brest Durham Symposium on Stochastic Analysis, July 1th to July 2th,

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

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

The Azéma-Yor Embedding in Non-Singular Diffusions

The Azéma-Yor Embedding in Non-Singular Diffusions Stochastic Process. Appl. Vol. 96, No. 2, 2001, 305-312 Research Report No. 406, 1999, Dept. Theoret. Statist. Aarhus The Azéma-Yor Embedding in Non-Singular Diffusions J. L. Pedersen and G. Peskir Let

More information

A NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES

A NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES A NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES STEFAN TAPPE Abstract. In a work of van Gaans (25a) stochastic integrals are regarded as L 2 -curves. In Filipović and Tappe (28) we have shown the connection

More information

Fast-slow systems with chaotic noise

Fast-slow systems with chaotic noise Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David

More information

Normal approximation of Poisson functionals in Kolmogorov distance

Normal approximation of Poisson functionals in Kolmogorov distance Normal approximation of Poisson functionals in Kolmogorov distance Matthias Schulte Abstract Peccati, Solè, Taqqu, and Utzet recently combined Stein s method and Malliavin calculus to obtain a bound for

More information

Precise Asymptotics for Infinite Dimensional Itô-Lyons Maps of Brownian Rough Paths

Precise Asymptotics for Infinite Dimensional Itô-Lyons Maps of Brownian Rough Paths Precise Asymptotics for Infinite Dimensional Itô-Lyons Maps of Brownian Rough Paths Jointwork with Yuzuru INAHAMA (Nagoya Univ.) Hiroshi KAWABI (Department of Mathematics, Faculty of Science Okayama University,

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

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

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

Jae Gil Choi and Young Seo Park

Jae Gil Choi and Young Seo Park Kangweon-Kyungki Math. Jour. 11 (23), No. 1, pp. 17 3 TRANSLATION THEOREM ON FUNCTION SPACE Jae Gil Choi and Young Seo Park Abstract. In this paper, we use a generalized Brownian motion process to define

More information

Backward Stochastic Differential Equations with Infinite Time Horizon

Backward Stochastic Differential Equations with Infinite Time Horizon Backward Stochastic Differential Equations with Infinite Time Horizon Holger Metzler PhD advisor: Prof. G. Tessitore Università di Milano-Bicocca Spring School Stochastic Control in Finance Roscoff, March

More information

Infinite-dimensional calculus under weak spatial regularity of the processes

Infinite-dimensional calculus under weak spatial regularity of the processes Infinite-dimensional calculus under weak spatial regularity of the processes Franco Flandoli1, Francesco Russo2, Giovanni Zanco3 November 8, 215 Abstract Two generalizations of Itô formula to infinite-dimensional

More information

THE MALLIAVIN CALCULUS FOR SDE WITH JUMPS AND THE PARTIALLY HYPOELLIPTIC PROBLEM

THE MALLIAVIN CALCULUS FOR SDE WITH JUMPS AND THE PARTIALLY HYPOELLIPTIC PROBLEM Takeuchi, A. Osaka J. Math. 39, 53 559 THE MALLIAVIN CALCULUS FOR SDE WITH JUMPS AND THE PARTIALLY HYPOELLIPTIC PROBLEM ATSUSHI TAKEUCHI Received October 11, 1. Introduction It has been studied by many

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

Regularity of the density for the stochastic heat equation

Regularity of the density for the stochastic heat equation Regularity of the density for the stochastic heat equation Carl Mueller 1 Department of Mathematics University of Rochester Rochester, NY 15627 USA email: cmlr@math.rochester.edu David Nualart 2 Department

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

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R

More information

Harmonic Functions and Brownian motion

Harmonic Functions and Brownian motion Harmonic Functions and Brownian motion Steven P. Lalley April 25, 211 1 Dynkin s Formula Denote by W t = (W 1 t, W 2 t,..., W d t ) a standard d dimensional Wiener process on (Ω, F, P ), and let F = (F

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

Econ Lecture 3. Outline. 1. Metric Spaces and Normed Spaces 2. Convergence of Sequences in Metric Spaces 3. Sequences in R and R n

Econ Lecture 3. Outline. 1. Metric Spaces and Normed Spaces 2. Convergence of Sequences in Metric Spaces 3. Sequences in R and R n Econ 204 2011 Lecture 3 Outline 1. Metric Spaces and Normed Spaces 2. Convergence of Sequences in Metric Spaces 3. Sequences in R and R n 1 Metric Spaces and Metrics Generalize distance and length notions

More information

Letting p shows that {B t } t 0. Definition 0.5. For λ R let δ λ : A (V ) A (V ) be defined by. 1 = g (symmetric), and. 3. g

Letting p shows that {B t } t 0. Definition 0.5. For λ R let δ λ : A (V ) A (V ) be defined by. 1 = g (symmetric), and. 3. g 4 Contents.1 Lie group p variation results Suppose G, d) is a group equipped with a left invariant metric, i.e. Let a := d e, a), then d ca, cb) = d a, b) for all a, b, c G. d a, b) = d e, a 1 b ) = a

More information

Densities for the Navier Stokes equations with noise

Densities for the Navier Stokes equations with noise Densities for the Navier Stokes equations with noise Marco Romito Università di Pisa Universitat de Barcelona March 25, 2015 Summary 1 Introduction & motivations 2 Malliavin calculus 3 Besov bounds 4 Other

More information

Recent results in game theoretic mathematical finance

Recent results in game theoretic mathematical finance Recent results in game theoretic mathematical finance Nicolas Perkowski Humboldt Universität zu Berlin May 31st, 2017 Thera Stochastics In Honor of Ioannis Karatzas s 65th Birthday Based on joint work

More information

ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS

ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS Bendikov, A. and Saloff-Coste, L. Osaka J. Math. 4 (5), 677 7 ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS ALEXANDER BENDIKOV and LAURENT SALOFF-COSTE (Received March 4, 4)

More information

ON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME

ON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME ON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME SAUL D. JACKA AND ALEKSANDAR MIJATOVIĆ Abstract. We develop a general approach to the Policy Improvement Algorithm (PIA) for stochastic control problems

More information

Rough Burgers-like equations with multiplicative noise

Rough Burgers-like equations with multiplicative noise Rough Burgers-like equations with multiplicative noise Martin Hairer Hendrik Weber Mathematics Institute University of Warwick Bielefeld, 3.11.21 Burgers-like equation Aim: Existence/Uniqueness for du

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

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

arxiv: v1 [math.pr] 24 Sep 2018

arxiv: v1 [math.pr] 24 Sep 2018 A short note on Anticipative portfolio optimization B. D Auria a,b,1,, J.-A. Salmerón a,1 a Dpto. Estadística, Universidad Carlos III de Madrid. Avda. de la Universidad 3, 8911, Leganés (Madrid Spain b

More information

Contents. 1 Preliminaries 3. Martingales

Contents. 1 Preliminaries 3. Martingales Table of Preface PART I THE FUNDAMENTAL PRINCIPLES page xv 1 Preliminaries 3 2 Martingales 9 2.1 Martingales and examples 9 2.2 Stopping times 12 2.3 The maximum inequality 13 2.4 Doob s inequality 14

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

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

Weak solutions of mean-field stochastic differential equations

Weak solutions of mean-field stochastic differential equations Weak solutions of mean-field stochastic differential equations Juan Li School of Mathematics and Statistics, Shandong University (Weihai), Weihai 26429, China. Email: juanli@sdu.edu.cn Based on joint works

More information

Quasi-invariant measures on the path space of a diffusion

Quasi-invariant measures on the path space of a diffusion Quasi-invariant measures on the path space of a diffusion Denis Bell 1 Department of Mathematics, University of North Florida, 4567 St. Johns Bluff Road South, Jacksonville, FL 32224, U. S. A. email: dbell@unf.edu,

More information

Elliptic Operators with Unbounded Coefficients

Elliptic Operators with Unbounded Coefficients Elliptic Operators with Unbounded Coefficients Federica Gregorio Universitá degli Studi di Salerno 8th June 2018 joint work with S.E. Boutiah, A. Rhandi, C. Tacelli Motivation Consider the Stochastic Differential

More information

Quasi-invariant Measures on Path Space. Denis Bell University of North Florida

Quasi-invariant Measures on Path Space. Denis Bell University of North Florida Quasi-invariant Measures on Path Space Denis Bell University of North Florida Transformation of measure under the flow of a vector field Let E be a vector space (or a manifold), equipped with a finite

More information

FRACTAL DIMENSIONS OF ROUGH DIFFERENTIAL EQUATIONS DRIVEN BY FRACTIONAL BROWNIAN MOTIONS

FRACTAL DIMENSIONS OF ROUGH DIFFERENTIAL EQUATIONS DRIVEN BY FRACTIONAL BROWNIAN MOTIONS FRACTAL DIMENSIONS OF ROUGH DIFFERENTIAL EQUATIONS DRIVEN BY FRACTIONAL BROWNIAN MOTIONS SHUWEN LOU AND CHENG OUYANG ABSTRACT. In this work we study fractal properties of a d-dimensional rough differential

More information

Smoothness of the distribution of the supremum of a multi-dimensional diffusion process

Smoothness of the distribution of the supremum of a multi-dimensional diffusion process Smoothness of the distribution of the supremum of a multi-dimensional diffusion process Masafumi Hayashi Arturo Kohatsu-Higa October 17, 211 Abstract In this article we deal with a multi-dimensional diffusion

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

arxiv: v2 [math.pr] 27 Oct 2015

arxiv: v2 [math.pr] 27 Oct 2015 A brief note on the Karhunen-Loève expansion Alen Alexanderian arxiv:1509.07526v2 [math.pr] 27 Oct 2015 October 28, 2015 Abstract We provide a detailed derivation of the Karhunen Loève expansion of a stochastic

More information

Homogenization for chaotic dynamical systems

Homogenization for chaotic dynamical systems Homogenization for chaotic dynamical systems David Kelly Ian Melbourne Department of Mathematics / Renci UNC Chapel Hill Mathematics Institute University of Warwick November 3, 2013 Duke/UNC Probability

More information

Solutions to the Exercises in Stochastic Analysis

Solutions to the Exercises in Stochastic Analysis Solutions to the Exercises in Stochastic Analysis Lecturer: Xue-Mei Li 1 Problem Sheet 1 In these solution I avoid using conditional expectations. But do try to give alternative proofs once we learnt conditional

More information

13 The martingale problem

13 The martingale problem 19-3-2012 Notations Ω complete metric space of all continuous functions from [0, + ) to R d endowed with the distance d(ω 1, ω 2 ) = k=1 ω 1 ω 2 C([0,k];H) 2 k (1 + ω 1 ω 2 C([0,k];H) ), ω 1, ω 2 Ω. F

More information

Stochastic Processes II/ Wahrscheinlichkeitstheorie III. Lecture Notes

Stochastic Processes II/ Wahrscheinlichkeitstheorie III. Lecture Notes BMS Basic Course Stochastic Processes II/ Wahrscheinlichkeitstheorie III Michael Scheutzow Lecture Notes Technische Universität Berlin Sommersemester 218 preliminary version October 12th 218 Contents

More information

arxiv: v1 [math.pr] 10 Jan 2019

arxiv: v1 [math.pr] 10 Jan 2019 Gaussian lower bounds for the density via Malliavin calculus Nguyen Tien Dung arxiv:191.3248v1 [math.pr] 1 Jan 219 January 1, 219 Abstract The problem of obtaining a lower bound for the density is always

More information

1 Math 241A-B Homework Problem List for F2015 and W2016

1 Math 241A-B Homework Problem List for F2015 and W2016 1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let

More information

Brownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion

Brownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion Brownian Motion An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Background We have already seen that the limiting behavior of a discrete random walk yields a derivation of

More information

Robust Mean-Variance Hedging via G-Expectation

Robust Mean-Variance Hedging via G-Expectation Robust Mean-Variance Hedging via G-Expectation Francesca Biagini, Jacopo Mancin Thilo Meyer Brandis January 3, 18 Abstract In this paper we study mean-variance hedging under the G-expectation framework.

More information

Tools of stochastic calculus

Tools of stochastic calculus slides for the course Interest rate theory, University of Ljubljana, 212-13/I, part III József Gáll University of Debrecen Nov. 212 Jan. 213, Ljubljana Itô integral, summary of main facts Notations, basic

More information

An introduction to rough paths

An introduction to rough paths An introduction to rough paths Antoine LEJAY INRIA, Nancy, France From works from T. Lyons, Z. Qian, P. Friz, N. Victoir, M. Gubinelli, D. Feyel, A. de la Pradelle, A.M. Davie,... SPDE semester Isaac Newton

More information