Approximation of BSDEs using least-squares regression and Malliavin weights
|
|
- Katherine Russell
- 5 years ago
- Views:
Transcription
1 Approximation of BSDEs using least-squares regression and Malliavin weights Plamen Turkedjiev 3rd July, 2012 Joint work with Prof. Emmanuel Gobet (E cole Polytechnique) Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
2 FBSDEs T > 0, W q-dimensional Brownian motion, (Ω, F, P) filtered probability space with usual conditions, but filtration may be larger than that generated by W, ξ L 2 (F T ), T T Y t = ξ + f(s, X s, Y s, Z s )ds Z s dw s (L T L t ) t t where X is d-dimensional, (t, x, y, z) f(t, x, y, z) is Borel measurable. Typically, X is a jump-diffusion driven by W and a Poisson random measure, L is a martingale orthogonal to W, and ξ = Φ(X T ) X is a diffusion driven by W, L 0, and ξ = Φ(X t1,..., X T ) or ξ = Φ(X T, T 0 X tdt) Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
3 Local Lipschitz condition and Quadratic BSDEs We consider time-local Lipschitz continuous driver: f(t, x, y, z) f(t, x, y, z ) L f ( x x + y y + z z (T t) (1 θ)/2 ) MOTIVATION: Assume X is a diffusion (L 0) and driver satisfies quadratic growth condition f(t, x, y, z) c(1 + y + z 2 ) f(t, x, y, z) f(t, x, y, z ) c(1 + z + z )( y y + z z ) and x Φ(x) is Hölder continuous and bounded. Then Z t L f (T t) (1 θ)/2 holds P dt-a.e. for constants L f and θ independent of t. Locally Lipschitz can replace quadratic in this special problem! Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
4 Driver with exploding bound and variance reduction For all t > 0 and x R d, there exists α (0, 1] and C f > 0 such that f(t, x, 0, 0) C f (T t) 1 α MOTIVATION: ξ = Φ(X T ) and x Φ(x) is α-hölder continuous and bounded. X is a diffusion process (L 0) so that v t (x) = E[Φ(X T ) X t = x] is smooth. f(t, x, y, z) uniformly Lipschitz continuous and unif. bounded at (y, z) = (0, 0). (v t (X t ), v t (X t )σ(t, X t )) solves BSDE with data (ξ, 0). Suppose we can solve this BSDE! v t (x) C(T t) α 1 standard from PDE theory. (Y t v t (X t ), Z t v t (X t )σ(t, X t )) solves a BSDE with data (0, f 0 ), where f 0 (t, x, y, z) = f(t, x, y + v t (x), z + v t (x)σ(t, x)). This BSDE may be better behaved for simulation purposes. Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
5 Key property: discretizability of FBSDEs Time-grid: π = (0 = t 0 <... < t N = T ). Paritcularly important grid: for β (0, 1], π β for which t β i := T T (1 1 N )1/β. Theorem If α = 1, let β = 1; else let β < α. Under the given assumptions, there exists a positive constant C, independent of N, such that max sup 0 i N 1 t β i t<tβ i+1 E Y t Y t β 2 + i N 1 i=0 t β i+1 t β i E Z t Z t β 2 dt CN 1 i We say that O(N 1/2 ) is the optimal rate of convergence for a discrete-time approximation of the BSDE. Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
6 Algorithm 1: Multistep dynamical programming Let i = t i+1 t i and W i := W ti+1 W ti. Recurssively build approximation of the solution, starting at i = N 1: i Z i = E i [ Wi (ξ + N 1 k=i+1 f(t k, X k, Y k+1, Z k ) k )], Y i = E i [ξ + N 1 k=i f(t k, X k, Y k+1, Z k ) k ], Y N = ξ. Consistency conditions for the time-grid: sup k<n k (T t k ) 1 θ 0 as N, lim sup N sup k<n 1 k k+1. Theorem For N sufficiently large, there exists a positive constant C independent of the time-grid such that max E Y i Y 0 k N 1 t β i N i=0 E Z i Z t β 2 i CN 1 i Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
7 Assumptions and properties Markov structure Let ξ = Φ(X N ) and X be a Markov chain.this ensures (Y i, Z i) = (y i(x i), z i(x i)) for measurable (unknown) functions y i and z i. Almost sure bounds Let x Φ(x) be bounded. This ensures that C y > 0 such that, k, Y k C y and Z k Cy P-almost surely. k Basis functions For each 0 l q and 0 k N 1, take a finite number of functions p l,k ( ) = (p i l,k) 1 i K such that p l,k : R d R is deterministic and E[ p l,k (X k ) 2 ] <. Form basis of finite dimensional subspaces of L 2(F tk ). Simulations Take M independent simulations of the Brownian increments W and the explanatory Markov chain X. Denote these simulations by (X m k ) 1 m M and ( W m k ) 1 m M respectively. Let p m l,k := p l,k (X m k ). Definition For R > 0, the truncated Brownian increment is defined by [ W i] R = R i W i R i. Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
8 Emprical regression algorithm Set y R,M N ( ) = Φ( ). Then, for i < N, compute coefficients α M l,i α M 0,k = arg min α = arg min α 1 M 1 M M m=1 M m=1 [ W ( N 1 l,i] R Φ(XN m ) + i N 1 Φ(XN m ) + k=i+1 k=i+1 The coefficients are not independent of one another! Set f i (y R,M k+1 (Xm k+1 ), zr,m k (Xk m )) k f i (y R,M k+1 (Xm k+1 ), zr,m k (Xk m )) k α p m 0,k 2. y R,M i (x) = C y α M 0,i p 0,i(x) C y, z R,M l,i (x) = Cy αl,i M p 0,i(x) Cy. i i ) α p m l,k 2 Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
9 Key ingredient: concentration of measure inequalities Needed, amongst other things, to deal with the lack of independence between regression coefficients. The following example comes from [Györfi et al. 2002, Theorem 11.2]. Benefit: the estimates are distribution-free. Theorem Let F {f : R d [ B, B]} and (Z i ) 1 i n be i.i.d. Then, for all ε > 0. Proposition P( f F : (E[ f(z) 2 ]) 1/2 2( 1 n E[N 2 ( If F is in a K-dimensional vector space, n f(z i ) 2 ) 1/2 > ε) i=1 ) 2 24 ε, F, Z 1:n)] exp ( nε2 288B 2 ( ( )) 2eB 2 3eB 2 K N 2 (ε, F, z 1:n ) 3 ε 2 log ε 2 Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
10 Error estimates Norm: For function Ψ, define Ψ 2 k,m := 1 M M m=1 Ψ(Xm k ) 2. Theorem For N sufficiently large, there exists a possitive constant C independent of the time-grid, M and the basis functions such that N 1 k=0 { E[ y k y R,M k C N 1 k=0 + C 2 k,m ] + E[ z k z R,M k 2 k,m ] } k { min α E yr k (X k) α p 0,k 2 k + N 1 k=0 { KN k M N 1 + CNKR 2 3 CN q l=1 } + KNR2 + CN θconv M k=0 ( exp CMR 2 k KN 1+θconv } min α E zr l,k (X k) α p l,k 2 k ) N 1 i=k ( CKR 2 N θconv i k ) CK +... Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
11 Complexity analysis Aim: reduce error to O(N 2θconv ). Assume y i C κ+1+η b, z i C κ+η b. Local polynomials on disjoint hypercubes, degree κ + 1 for Y and κ for Z. Bias approximation: O(N 2θconv ) if δ z = cn θconv κ+η. K = cn d θconv κ+η up to log terms. θconv 2+2θconv+2d = cn κ+η Large deviation terms: M = ckn 2+2θconv terms. θconv 3+2θconv+2d Computational work C = cmn = cn κ+η N 2θconv cc 1 2(1+ 3 2θconv + κ+η d ). ODP scheme with θ conv = 1/2: N 1 cc if κ + η > 1, MDP has better performance. 1 2(4+ 2d κ+1+η ). up to log up to log terms. Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
12 BSDE numerics: Implicit vs explicit q = 3; f(z) = 1.5 z (1) ; Φ(x) = (x (3) 100) + ; X (j) t = 100e (σwt)(j) with σ 1 1 ρ 2 σ 0 σ ρ σ = 0 σ 2 1 ρ 2 σ 2ρ, σ 2 σ 3 = σ 3 ρ 0.1 N = 16; Basis q i=0 g i(ln(x i )) for g i Hermite polynomials with i deg(g i) 3. Explicit solution: Y t = BlackScholesCall(t, X t ; σ 3, 100), Z t = (0, 0, BlackScholesHedge(t, X t ; σ 3, 100)). Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
13 BSDE numerics Explicit vs Implicit Multistep forward scheme Explicit Implicit log(error) log(work) Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
14 Representation theorem due to Ma/Zhang X is a diffusion: dx t = b(t, X t )dt + σ(t, X t )dw t. Also need gradient process and it s inverse: d X t = b x (t, X t ) X t dt + σ(t, X t ) X t dw t, d X 1 t = ( b x (t, X t ) σ x (t, X t ) 2 ) Xt 1 dt + σ x (t, X t ) Xt 1 dw t. Representation theorem due to Ma/Zhang for Z: T Z t = E t [ξht t + f(r, X r, Y r, Z r )Hrdr] t t where (r t)h t r = ( t r [σ 1 (s, X s ) X s X 1 t σ(t, X t )] dw s ). H t r are the Malliavin weights; the representation formula is derived by means of Malliavin s calculus, but remains true in the Lipschitz case, even though the BSDE is not Mallivin differentiable. Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
15 Algorithm 2: Malliavin weights Let (t k t i )Hk i = ( k 1 j=i σ 1 (t j, X j ) X j X 1 i σ(t i, X i )] W j ). Recursively build the approximation starting at i = N 1: Z i = E i [ξhn i + N 1 k=i+1 f(t k, X k, Y k+1, Z k )Hk i k], Y i = E i [ξ + N 1 k=i f(t k, X k, Y k+1, Z k ) k ], Y N = ξ. Constraint on the time-grid: lim sup N sup i<n i+1 i <. Recall the special time-grid π β : Theorem For sufficiently high N, there exists a positive constant C independent of the time-grid such that max E Y i Y 0 k N 1 t β i N i=0 E Z i Z t β 2 i CN 1 i Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
16 Projection estimates Approximate conditional expectation E i by projection on finite subspace of L 2 (F ti ): Ẑ l,i = arg min α pl,i (X i) E[ Φ(X N )Hl,N i + N 1 k=i+1 f(t k, X k, Ŷk+1, Ẑk)Hl,k i k α p l,i (X i ) 2 ], Ŷ i = arg inf α p0,i(x i) E[ Φ(X N ) + N 1 k=i f(t k, X k, Ŷk+1, Ẑk) k α p 0,i (X i ) 2 ], Ŷ N = Φ(X N ). Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
17 Projection estimates Theorem There exists positive constant C independent of the time-grid such that E Y i Ŷi 2 2E Y i P Y i Y i 2 N 1 + C k=i E Z i Ẑi 2 2E Z i P Z i Z i 2 N 1 + C k=i {E Y k+1 P Y k+1 Y k E Z k P Z k Z k 2 } k (T t k ) 1 θ {E Y k+1 P Y k+1 Y k E Z k P Z k Z k 2 } k (T t k ) 1 θ Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
18 Almost sure bounds There exist positive constants C y and C z independent of the time-grid such that, i, Y i C y and Z i C z P-almost surely T ti Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
19 Finally... Thank You For Your Attention! Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
20 References [GT11] E. Gobet, T. Approximation for discrete BSDE using least-squares regression. [BD07] [CD11] [DG06] C. Bender and R. Denk. A forward scheme for backward SDEs. Stochastic Processes and their Applications, 117(12): , D. Crisan, F. Delarue. Sharp gradient bounds for solutions of semi-linear PDEs F. Delarue, G. Guatteri. Weak existence and uniqueness for forward-backward SDEs. Stochastic processes and their applications, 116(12): , [GGG11] C. Geiss, S. Geiss, and E. Gobet. Generalized fractional smoothness and L p -variation of BSDEs with non-lipschitz terminal condition. To appear in Stochastic Processes and their Applications, Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
21 References [GLW05] E. Gobet, J.P. Lemor, and X. Warin. A regression-based Monte Carlo method to solve backward stochastic differential equations. Annals of Applied Probability, 15(3): , [LGW06] J.P. Lemor, E. Gobet, and X. Warin. Rate of convergence of an empirical regression method for solving generalized BSDEs. Bernoulli, 12(5): , [GM10] E. Gobet and A. Makhlouf. L 2 -time regularity of BSDEs with irregular terminal functions. Stochastic Processes and their Applications, 120: , [GKKW02] L. Gyorfi, M. Kohler, A. Krzyzak, and H. Walk. A distribution-free theory of nonparametric regression. Springer Series in Statistics, Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
22 References [MZ02] [Mo10] J. Ma, J. Zhang. Representation theorems for BSDEs. Annals of Applied Probability, 12(4): , T. Moseler. A Picard-type iteration for BSDEs: Convergence and importance sampling. PhD thesis. Plamen Turkedjiev (HU Berlin) Approximation of BSDEs 3rd July, / 22
Approximation of backward stochastic differential equations using Malliavin weights and least-squares regression
Approximation of bacward stochastic differential equations using Malliavin weights and least-squares regression Emmanuel Gobet, Plamen Turedjiev To cite this version: Emmanuel Gobet, Plamen Turedjiev.
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationBrownian 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 informationUniformly Uniformly-ergodic Markov chains and BSDEs
Uniformly Uniformly-ergodic Markov chains and BSDEs Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) Centre Henri Lebesgue,
More informationNonlinear 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 informationFunctional 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 informationLinear regression MDP scheme for discrete backward stochastic differential equations under general conditions
Linear regression MDP scheme for discrete bacward stochastic differential equations under general conditions Emmanuel Gobet, Plamen Turedjiev To cite this version: Emmanuel Gobet, Plamen Turedjiev. Linear
More informationStrong Solutions and a Bismut-Elworthy-Li Formula for Mean-F. Formula for Mean-Field SDE s with Irregular Drift
Strong Solutions and a Bismut-Elworthy-Li Formula for Mean-Field SDE s with Irregular Drift Thilo Meyer-Brandis University of Munich joint with M. Bauer, University of Munich Conference for the 1th Anniversary
More informationThe 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 informationProperties of an infinite dimensional EDS system : the Muller s ratchet
Properties of an infinite dimensional EDS system : the Muller s ratchet LATP June 5, 2011 A ratchet source : wikipedia Plan 1 Introduction : The model of Haigh 2 3 Hypothesis (Biological) : The population
More informationMathematical 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 informationMean-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 information1. 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 informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More informationWeak 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 informationOn continuous time contract theory
Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem
More informationWiener 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 informationMarkov Chain BSDEs and risk averse networks
Markov Chain BSDEs and risk averse networks Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) 2nd Young Researchers in BSDEs
More informationMonte-Carlo MMD-MA, Université Paris-Dauphine. Xiaolu Tan
Monte-Carlo MMD-MA, Université Paris-Dauphine Xiaolu Tan tan@ceremade.dauphine.fr Septembre 2015 Contents 1 Introduction 1 1.1 The principle.................................. 1 1.2 The error analysis
More informationBrownian 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 informationEstimates for the density of functionals of SDE s with irregular drift
Estimates for the density of functionals of SDE s with irregular drift Arturo KOHATSU-HIGA a, Azmi MAKHLOUF a, a Ritsumeikan University and Japan Science and Technology Agency, Japan Abstract We obtain
More informationp 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 informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationStochastic 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 informationBrownian Motion and Conditional Probability
Math 561: Theory of Probability (Spring 2018) Week 10 Brownian Motion and Conditional Probability 10.1 Standard Brownian Motion (SBM) Brownian motion is a stochastic process with both practical and theoretical
More informationDynamic 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 informationBackward 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 informationUNCERTAINTY 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 informationPROBABILITY: 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 informationInformation and Credit Risk
Information and Credit Risk M. L. Bedini Université de Bretagne Occidentale, Brest - Friedrich Schiller Universität, Jena Jena, March 2011 M. L. Bedini (Université de Bretagne Occidentale, Brest Information
More informationLarge 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 informationLAN 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 informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationLAN 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 informationOn a class of stochastic differential equations in a financial network model
1 On a class of stochastic differential equations in a financial network model Tomoyuki Ichiba Department of Statistics & Applied Probability, Center for Financial Mathematics and Actuarial Research, University
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationMalliavin 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 informationOn 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 informationBSDEs and PDEs with discontinuous coecients Applications to homogenization K. Bahlali, A. Elouain, E. Pardoux. Jena, March 2009
BSDEs and PDEs with discontinuous coecients Applications to homogenization K. Bahlali, A. Elouain, E. Pardoux. Jena, 16-20 March 2009 1 1) L p viscosity solution to 2nd order semilinear parabolic PDEs
More informationExercises. 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 informationP (A G) dp G P (A G)
First homework assignment. Due at 12:15 on 22 September 2016. Homework 1. We roll two dices. X is the result of one of them and Z the sum of the results. Find E [X Z. Homework 2. Let X be a r.v.. Assume
More informationWEAK VERSIONS OF STOCHASTIC ADAMS-BASHFORTH AND SEMI-IMPLICIT LEAPFROG SCHEMES FOR SDES. 1. Introduction
WEAK VERSIONS OF STOCHASTIC ADAMS-BASHFORTH AND SEMI-IMPLICIT LEAPFROG SCHEMES FOR SDES BRIAN D. EWALD 1 Abstract. We consider the weak analogues of certain strong stochastic numerical schemes considered
More informationMA8109 Stochastic Processes in Systems Theory Autumn 2013
Norwegian University of Science and Technology Department of Mathematical Sciences MA819 Stochastic Processes in Systems Theory Autumn 213 1 MA819 Exam 23, problem 3b This is a linear equation of the form
More informationEULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS
Qiao, H. Osaka J. Math. 51 (14), 47 66 EULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS HUIJIE QIAO (Received May 6, 11, revised May 1, 1) Abstract In this paper we show
More informationA Central Limit Theorem for Fleming-Viot Particle Systems Application to the Adaptive Multilevel Splitting Algorithm
A Central Limit Theorem for Fleming-Viot Particle Systems Application to the Algorithm F. Cérou 1,2 B. Delyon 2 A. Guyader 3 M. Rousset 1,2 1 Inria Rennes Bretagne Atlantique 2 IRMAR, Université de Rennes
More informationOn the Multi-Dimensional Controller and Stopper Games
On the Multi-Dimensional Controller and Stopper Games Joint work with Yu-Jui Huang University of Michigan, Ann Arbor June 7, 2012 Outline Introduction 1 Introduction 2 3 4 5 Consider a zero-sum controller-and-stopper
More informationPROBABILITY: 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 informationPreliminary Exam: Probability 9:00am 2:00pm, Friday, January 6, 2012
Preliminary Exam: Probability 9:00am 2:00pm, Friday, January 6, 202 The exam lasts from 9:00am until 2:00pm, with a walking break every hour. Your goal on this exam should be to demonstrate mastery of
More informationStochastic differential equation models in biology Susanne Ditlevsen
Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential
More informationTopics 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 informationNumerical scheme for quadratic BSDEs and Dynamic Risk Measures
Numerical scheme for quadratic BSDEs, Numerics and Finance, Oxford Man Institute Julien Azzaz I.S.F.A., Lyon, FRANCE 4 July 2012 Joint Work In Progress with Anis Matoussi, Le Mans, FRANCE Outline Framework
More informationA 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 informationOPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS
APPLICATIONES MATHEMATICAE 29,4 (22), pp. 387 398 Mariusz Michta (Zielona Góra) OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS Abstract. A martingale problem approach is used first to analyze
More informationConvergence at first and second order of some approximations of stochastic integrals
Convergence at first and second order of some approximations of stochastic integrals Bérard Bergery Blandine, Vallois Pierre IECN, Nancy-Université, CNRS, INRIA, Boulevard des Aiguillettes B.P. 239 F-5456
More informationAn 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 informationA MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT
A MODEL FOR HE LONG-ERM OPIMAL CAPACIY LEVEL OF AN INVESMEN PROJEC ARNE LØKKA AND MIHAIL ZERVOS Abstract. We consider an investment project that produces a single commodity. he project s operation yields
More informationExponential 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 informationA class of globally solvable systems of BSDE and applications
A class of globally solvable systems of BSDE and applications Gordan Žitković Department of Mathematics The University of Texas at Austin Thera Stochastics - Wednesday, May 31st, 2017 joint work with Hao
More informationON THE PATHWISE UNIQUENESS OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS
PORTUGALIAE MATHEMATICA Vol. 55 Fasc. 4 1998 ON THE PATHWISE UNIQUENESS OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS C. Sonoc Abstract: A sufficient condition for uniqueness of solutions of ordinary
More informationSimulation methods for stochastic models in chemistry
Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation
More informationI 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 informationDeep learning-based numerical methods for high-dimensional parabolic partial differential equations and backward stochastic differential equations
Deep learning-based numerical methods for high-dimensional parabolic partial differential equations and backward stochastic differential equations arxiv:1706.04702v1 [math.na] 15 Jun 2017 Weinan E 1, Jiequn
More informationHigh 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 informationBrownian Motion and Stochastic Calculus
ETHZ, Spring 17 D-MATH Prof Dr Martin Larsson Coordinator A Sepúlveda Brownian Motion and Stochastic Calculus Exercise sheet 6 Please hand in your solutions during exercise class or in your assistant s
More informationMultiple points of the Brownian sheet in critical dimensions
Multiple points of the Brownian sheet in critical dimensions Robert C. Dalang Ecole Polytechnique Fédérale de Lausanne Based on joint work with: Carl Mueller Multiple points of the Brownian sheet in critical
More informationNumerical treatment of stochastic delay differential equations
Numerical treatment of stochastic delay differential equations Evelyn Buckwar Department of Mathematics Heriot-Watt University, Edinburgh Evelyn Buckwar Roma, 16.10.2007 Overview An example: machine chatter
More informationIntroduction to self-similar growth-fragmentations
Introduction to self-similar growth-fragmentations Quan Shi CIMAT, 11-15 December, 2017 Quan Shi Growth-Fragmentations CIMAT, 11-15 December, 2017 1 / 34 Literature Jean Bertoin, Compensated fragmentation
More informationSTATISTICS 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 informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationIn terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3.
1. GAUSSIAN PROCESSES A Gaussian process on a set T is a collection of random variables X =(X t ) t T on a common probability space such that for any n 1 and any t 1,...,t n T, the vector (X(t 1 ),...,X(t
More informationPoisson random measure: motivation
: motivation The Lévy measure provides the expected number of jumps by time unit, i.e. in a time interval of the form: [t, t + 1], and of a certain size Example: ν([1, )) is the expected number of jumps
More informationExact Simulation of Multivariate Itô Diffusions
Exact Simulation of Multivariate Itô Diffusions Jose Blanchet Joint work with Fan Zhang Columbia and Stanford July 7, 2017 Jose Blanchet (Columbia/Stanford) Exact Simulation of Diffusions July 7, 2017
More informationGAUSSIAN 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 informationExercises 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 informationMultilevel Monte Carlo for Stochastic McKean-Vlasov Equations
Multilevel Monte Carlo for Stochastic McKean-Vlasov Equations Lukasz Szpruch School of Mathemtics University of Edinburgh joint work with Shuren Tan and Alvin Tse (Edinburgh) Lukasz Szpruch (University
More informationLecture 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 informationA 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 informationSome 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 informationProf. Erhan Bayraktar (University of Michigan)
September 17, 2012 KAP 414 2:15 PM- 3:15 PM Prof. (University of Michigan) Abstract: We consider a zero-sum stochastic differential controller-and-stopper game in which the state process is a controlled
More informationMathematical 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 October 22nd, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationHarnack Inequalities and Applications for Stochastic Equations
p. 1/32 Harnack Inequalities and Applications for Stochastic Equations PhD Thesis Defense Shun-Xiang Ouyang Under the Supervision of Prof. Michael Röckner & Prof. Feng-Yu Wang March 6, 29 p. 2/32 Outline
More informationTheoretical Tutorial Session 2
1 / 36 Theoretical Tutorial Session 2 Xiaoming Song Department of Mathematics Drexel University July 27, 216 Outline 2 / 36 Itô s formula Martingale representation theorem Stochastic differential equations
More informationGARCH 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 informationRandom Fields: Skorohod integral and Malliavin derivative
Dept. of Math. University of Oslo Pure Mathematics No. 36 ISSN 0806 2439 November 2004 Random Fields: Skorohod integral and Malliavin derivative Giulia Di Nunno 1 Oslo, 15th November 2004. Abstract We
More informationLecture 12. F o s, (1.1) F t := s>t
Lecture 12 1 Brownian motion: the Markov property Let C := C(0, ), R) be the space of continuous functions mapping from 0, ) to R, in which a Brownian motion (B t ) t 0 almost surely takes its value. Let
More informationDensities 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 informationOptimal stopping for non-linear expectations Part I
Stochastic Processes and their Applications 121 (2011) 185 211 www.elsevier.com/locate/spa Optimal stopping for non-linear expectations Part I Erhan Bayraktar, Song Yao Department of Mathematics, University
More informationStatistical test for some multistable processes
Statistical test for some multistable processes Ronan Le Guével Joint work in progress with A. Philippe Journées MAS 2014 1 Multistable processes First definition : Ferguson-Klass-LePage series Properties
More information1. Stochastic Process
HETERGENEITY IN QUANTITATIVE MACROECONOMICS @ TSE OCTOBER 17, 216 STOCHASTIC CALCULUS BASICS SANG YOON (TIM) LEE Very simple notes (need to add references). It is NOT meant to be a substitute for a real
More informationA Note on the Central Limit Theorem for a Class of Linear Systems 1
A Note on the Central Limit Theorem for a Class of Linear Systems 1 Contents Yukio Nagahata Department of Mathematics, Graduate School of Engineering Science Osaka University, Toyonaka 560-8531, Japan.
More informationSkorokhod Embeddings In Bounded Time
Skorokhod Embeddings In Bounded Time Stefan Ankirchner Institute for Applied Mathematics University of Bonn Endenicher Allee 6 D-535 Bonn ankirchner@hcm.uni-bonn.de Philipp Strack Bonn Graduate School
More informationMalliavin 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 informationMan Kyu Im*, Un Cig Ji **, and Jae Hee Kim ***
JOURNAL OF THE CHUNGCHEONG MATHEMATICAL SOCIETY Volume 19, No. 4, December 26 GIRSANOV THEOREM FOR GAUSSIAN PROCESS WITH INDEPENDENT INCREMENTS Man Kyu Im*, Un Cig Ji **, and Jae Hee Kim *** Abstract.
More informationThe Real Option Approach to Plant Valuation: from Spread Options to Optimal Switching
The Real Option Approach to Plant Valuation: from Spread Options to René 1 and Michael Ludkovski 2 1 Bendheim Center for Finance Department of Operations Research & Financial Engineering Princeton University
More informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationQuantization of stochastic processes with applications on Euler-Maruyama schemes
UPTEC F15 063 Examensarbete 30 hp Oktober 2015 Quantization of stochastic processes with applications on Euler-Maruyama schemes Viktor Edward Abstract Quantization of stochastic processes with applications
More informationLogFeller 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 informationSolvability of G-SDE with Integral-Lipschitz Coefficients
Solvability of G-SDE with Integral-Lipschitz Coefficients Yiqing LIN joint work with Xuepeng Bai IRMAR, Université de Rennes 1, FRANCE http://perso.univ-rennes1.fr/yiqing.lin/ ITN Marie Curie Workshop
More informationSimulation of conditional diffusions via forward-reverse stochastic representations
Weierstrass Institute for Applied Analysis and Stochastics Simulation of conditional diffusions via forward-reverse stochastic representations Christian Bayer and John Schoenmakers Numerical methods for
More informationOn 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