Exact Simulation of Diffusions and Jump Diffusions
|
|
- Garry Mills
- 5 years ago
- Views:
Transcription
1 Exact Simulation of Diffusions and Jump Diffusions A work by: Prof. Gareth O. Roberts Dr. Alexandros Beskos Dr. Omiros Papaspiliopoulos Dr. Bruno Casella 28 th May, 2008
2 Content 1 Exact Algorithm Construction EA1 EA2 EA3 2 Jump Exact Algorithm 3 Simulated Example 4 Work in progress
3 Construction Exact Algorithm (EA) Proposed by Beskos, Papaspiliopoulos and Roberts (2006). Exact simulation of a class of Itô s diffusions. Exact in the sense that there is no discretisation error. The EA performs rejection sampling by proposing paths from processes that we can simulate and accepting them according to appropriate probability density ratios. The novelty lies in the fact that the paths proposed are unveiled only at finite (but random) time instances and the decision whether to accept the path or not can be easily taken.
4 Construction Consider X := {X t : 0 t T } a one-dimensional diffusion process solving the SDE: dx t = b(x t )dt + σ(x t )db t, X 0 = x R, t [0, T ] We will restrict our attention to SDEs of the type dx t = α(x t )dt + db t, X 0 = x R, t [0, T ] since the transformation X t η(x t ) can be applied, where η(x) = x z 1 σ(u) du with z being some element of the state space of X.
5 Construction Consider X := {X t : 0 t T } a one-dimensional diffusion process solving the SDE: dx t = b(x t )dt + σ(x t )db t, X 0 = x R, t [0, T ] We will restrict our attention to SDEs of the type dx t = α(x t )dt + db t, X 0 = x R, t [0, T ] since the transformation X t η(x t ) can be applied, where η(x) = x z 1 σ(u) du with z being some element of the state space of X.
6 Construction Q: the probability law of X ; W: the probability law of a Brownian motion starting at x R. The Girsanov Formula implies that { dq T dw = exp α(x t )dx t 1 T } α 2 (X t )dt Under the condition that α is everywhere differentiable and applying Itô s lemma to A(X ), for A(u) := u R, we have u 0 α(y)dy, { dq dw = exp A(X T ) A(x) 1 T ( α 2 (X t ) + α (X t ) ) } dt 2 0
7 Construction Q: the probability law of X ; W: the probability law of a Brownian motion starting at x R. The Girsanov Formula implies that { dq T dw = exp α(x t )dx t 1 T } α 2 (X t )dt Under the condition that α is everywhere differentiable and applying Itô s lemma to A(X ), for A(u) := u R, we have u 0 α(y)dy, { dq dw = exp A(X T ) A(x) 1 T ( α 2 (X t ) + α (X t ) ) } dt 2 0
8 Construction The main idea is to perform rejection sampling. Problem: A(X T ) will typically be unbounded. Solution: To propose candidate paths from a process identical to BM starting at x except for the distribution of its ending point: biased Brownian motion. This biased BM is defined as Ŵ = (W W T h), where h exp{a(u) (u x) 2 /2T }, u R, has to be integrable.
9 Construction Let Z be the probability measure induced by the biased BM. Then It can be shown that dq dz = dq dw dw dz Implying dw dz = f N(X T ; x, T ) exp{ A(X T )} h { dq T ( 1 dz exp 0 2 α2 (X t ) + 1 ) } 2 α (X t ) dt
10 Construction We assume that (α 2 + α )/2 is bounded below implying that dq is bounded. dz It is possible to obtain a non-negative function φ such that Analytically, φ is defined as: { dq T } dz exp φ(x t )dt 1 0 φ(u) = α2 (u) + α (u) k, 2 u R, for k inf u R (α2 +α )(u)/2.
11 Construction Problem: It is not possible to draw complete continuous paths of Z on [0, T ] and calculate the integral involved in the acceptance probability analytically. Solution: Theorem Let X be any continuous mapping from [0, T ] to R, and M(X ) an upper bound for the mapping t φ(x t ), t [0, T ]. If Φ is a homogeneous Poisson process of unit intensity on [0, T ] [0, M(X )] and N is the number of points of Φ found below the graph {(t, φ(x t )); t [0, T ]}, then { T } P(N = 0 X ) = exp φ(x t )dt. 0
12 EA1 - The case when φ is bounded Suppose that M is an upper bound for φ(x t ), t [0, T ]. 1. Produce a realisation {x 1, x 2,..., x τ }, of a Poisson process on [0, T ] [0, M], where x i = (x i,1, x i,2 ), 1 i τ;
13 EA1 - The case when φ is bounded 2. Simulate a skeleton of X Z at the time instances {x 1,1, x 2,1,..., x τ,1 };
14 EA1 - The case when φ is bounded 3. Evaluate N; 4. If N = 0, go to step 5, else go to step 1; 5. Output the currently constructed skeleton S(X ) of X ; X is accepted and X is rejected.
15 EA1 - The case when φ is bounded The EA returns a Skeleton of the process X defined over the time interval [0, T ] with starting point x. We represent it as: S 0 (x; 0, T ) := {(t 0, x t0 ), (t 1, x t1 ),..., (t M, x tm )} with t 0 = 0 and t M = T. We can simulate further information on the process X using the following representation, immediately derived from elementary Brownian motion constructions: X S 0 (x; 0, T ) M 1 i=0 BB(t i, x ti ; t i+1, x ti+1 ) where BB(s 1, a; s 2, b) is the measure of a Brownian bridge starting in a at time s 1 and ending in b at time s 2.
16 EA2 - The case when either lim sup u φ(u) < or lim sup u φ(u) < If lim sup u φ(u) < or lim sup φ(u) <, it is possible u to identify an upper bound M(m) for φ(x t ), with t [0, T ], after decomposing the proposed path X t at its min or max. If lim sup φ(u) < and m is the minimum of X t in [0, T ] u M(m) = sup{φ(u); u m} If lim sup φ(u) < and m is the maximum of X t in [0, T ] u M(m) = sup{φ(u); u m}
17 EA2 - The case when either lim sup u φ(u) < or lim sup u φ(u) < The EA2 algorithm 1 Initiate a path X t Z on [0, T ] by drawing X T h; 2 Simulate its minimum or maximum m and the moment t m when it is achieved ; 3 Find an upper bound M(m) for φ(x t ); 4 Produce a realization {x 1, x 2,..., x τ }, of a Poisson process on [0, T ] [0, M(m)]; 5 Simulate a skeleton of X (m, t m ) at the time instances {x 1,1, x 2,1,..., x τ,1 } ; 6 Evaluate N; 7 If N = 0, go to step 8, else go to step 1; 8 Output the currently constructed skeleton S(X ) of X. The algorithms to perform steps 2 and 5 are described in Beskos, Papaspiliopoulos and Roberts (2006). ref
18 EA3 - The case when φ is only bounded below So far, we have made the following hypotheses to perform the Exact Algorithm. C0 The Radon-Nikodym derivative of Q w.r.t. W exists and it is given by Girsanov s formula; C1 α C 1 ; C2 α 2 + α is bounded below; C3 The function h(u) = exp{ (u x) 2 /(2T ) + A(u)} is integrable. EA3 requires no further hypotheses on φ. In this case we need to know upper and lower bounds for X t in [0, T ] in order to construct an upper bound for φ(x t ). The bounds of X t are obtained using the idea of a layered Brownian bridge.
19 EA3 - The case when φ is only bounded below Define two increasing sequences of positive real numbers {a i } i 1 and {b i } i 1 with a 0 = b 0 = 0. Given X 0 = x and X T = y, set x = x y and ȳ = x y and define the following events: U i = { sup X t [ȳ + b i 1, ȳ + b i ] 0 t T { L i = inf 0 t T Xt [ x a i, x a i 1 ] } { } { D i = U i Li, i 1. } inf 0 t T Xt > x a i, sup X t < ȳ + b i }, 0 t T Define the random variable I = I (X ) such that {I = i} = D i. {I = i} { x a i < X t < ȳ + b i, t [0, T ]}.
20 EA3 - The case when φ is only bounded below Example: a i = b i and I (X ) = 4.
21 EA3 - The case when φ is only bounded below The EA3 algorithm 1 Initiate a path X t Z on [0, T ] by drawing X T h; 2 Simulate I (X ) ; 3 Find an upper bound M(I ) for φ(x t ); 4 Produce a realization {x 1, x 2,..., x τ }, of a Poisson process on [0, T ] [0, M(I )]; 5 Simulate a skeleton of X I at the time instances {x 1,1, x 2,1,..., x τ,1 } ; 6 Evaluate N; 7 If N = 0, go to step 8, else go to step 1; 8 Output the currently constructed skeleton S(X ) of X. The algorithms to perform steps 2 and 5 are described in Beskos, Papaspiliopoulos and Roberts (2007). ref
22 EA3 - The case when φ is only bounded below If we want to simulate a diffusion conditioned on the ending point X T = y, the paths are proposed from a BB(0, x; T, y). The following paper proposes some methods for inference in diffusions based on the Exact Algorithm: Beskos, A., Papaspiliopoulos, O., Roberts, G, and Fearnhead, P. (2006), Exact and computationally efficient likelihood-based inference for discretely observed diffusion processes (with discussion), J. R. Stat. Soc. B 68(3),
23 Jump Exact Algorithm (JEA) Proposed by Casella and Roberts (2008). Generalises the Exact Algorithm for the simulation of jump diffusions.
24 Consider X := {X t : 0 t T } a one-dimensional jump diffusion process solving the SDE: dx t = b(x t )dt + σ(x t )db t + dj t, X 0 = x J t is a jump process where the jump times follow a Poisson process of rate λ(t, X t ) and the jump sizes are given by a function g(z t, X t ), where z t has a distribution f Zt.
25 Between any two jumps, the process X behaves as a homogeneous diffusion process with drift b and coefficient σ. Jump times and size may depend on time and on the state of the process. The jump times are generated by a Poisson process on [0, T ] with intensity function λ(t, X t ). A random variable Z i f Z (z; t i ) is associated to each jump time t i. Z i and the state of the process X t determine the amplitude g(z ti, X t ) of the jump in t i. Once more, we apply the transformation X t η(x t ) to work with dx t = α(x t )dt + db t + dj t, X 0 = x
26 Between any two jumps, the process X behaves as a homogeneous diffusion process with drift b and coefficient σ. Jump times and size may depend on time and on the state of the process. The jump times are generated by a Poisson process on [0, T ] with intensity function λ(t, X t ). A random variable Z i f Z (z; t i ) is associated to each jump time t i. Z i and the state of the process X t determine the amplitude g(z ti, X t ) of the jump in t i. Once more, we apply the transformation X t η(x t ) to work with dx t = α(x t )dt + db t + dj t, X 0 = x
27 Constructing the JEA We make the hypothesis that λ(t, X t ) is bounded. To break down the state-dependency of the jump times we use the following result. Poisson Thinning To simulate a Poisson process of rate λ(t, X t ) λ on [0, T ] do: 1 Simulate a Poisson process of rate λ on [0, T ]. 2 Keep each point with probability R(t i, X t ) = λ(t i, X t i λ i i = 1,..., n, where n is the number of points simulated on step 1. ), for
28 1. Simulate a Poisson process of rate λ on [0, T ];
29 2. Apply the Exact Algorithm to [0, T ];
30 3. Perform Poisson thinning until one point t 1 is kept or no point is kept;
31 4. Simulate the jump at t 1 ;
32 5. Apply the Exact Algorithm to [t 1, T ];
33 6. Perform Poisson thinning until one point t 2 is kept or no point is kept;
34 7. Simulate the jump at t 2 ;
35 8. Apply the Exact Algorithm to [t 2, T ];
36 9. Perform Poisson thinning until one point t 3 is kept or no point is kept;
37 10. If no point is kept...
38 10. Output skeleton S(X ) of X.
39 dx t = tanh(x t )dt + db t + dj t, X 0 = 0 X 2 t λ(x t ) = 1 + Xt 2 < 1, g(z ti, X t ) = Z Z N(0, 9) and T = 20 The JEA is compared to the Euler approximation scheme by comparing the distribution of the ending point.
40 One realization of JEA
41 Comparison with Euler approximation CPU time JEA 15.94s Euler t = s t = s t = s t = s t = s t = m 6s
42 Work in progress Extend the Jump Exact Algorithm to simulate jump diffusions with unbounded jump rate. Extend the Jump Exact Algorithm to simulate jump diffusions conditioned on the ending point. Develop methods for inference in jump diffusion processes based on the JEA. Unknown parameters in the drift and in the diffusion coefficient; Unknown parameters in the jump rate; Unknown parameters in the distribution of the jump sizes.
43 Thank you! Contact:
44 References Beskos, A; Papaspiliopoulos, O. and Roberts, G. O. (2006), Retrospective exact simulation of diffusion sample paths with applications. Bernoulli, 12(6), Return Beskos, A., Papaspiliopoulos, O., Roberts, G, and Fearnhead, P. (2006), Exact and computationally efficient likelihood-based inference for discretely observed diffusion processes (with discussion), J. R. Stat. Soc. B 68(3), Beskos, A; Papaspiliopoulos, O. and Roberts, G. O. (2007), A factorization of diffusion measure and finite sample path constructions. To appear in Methodology and Computing in Applied Probability. Return Casella, B. and Roberts, G. O. (2008), Exact simulation of jump-diffusion processes with Monte Carlo Applications. In preparation.
Retrospective Exact Simulation of Diffusion Sample Paths with Applications
Retrospective Exact Simulation of Diffusion Sample Paths with Applications Alexandros Beskos, Omiros Papaspiliopoulos and Gareth O. Roberts July 5, 2006 Abstract The objective of this paper is to present
More informationOn the exact and ε-strong simulation of (jump) diffusions
Bernoulli 22(2), 2016, 794 856 DOI: 10.3150/14-BEJ676 On the exact and ε-strong simulation of (jump) diffusions MURRAY POLLOCK *, ADAM M. JOHANSEN ** and GARETH O. ROBERTS Department of Statistics, University
More information(Exact and efficient) Simulation of Conditioned Markov Processes
(Exact and efficient) Simulation of Conditioned Markov Processes Omiros Papaspiliopoulos Universitat Pompeu Fabra, Barcelona http://www.econ.upf.edu/ omiros GrStoc α, ΛEYKA A 2009 Acknowledgement: K.Kalogeropoulos
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 informationRetrospective exact simulation of diffusion sample paths with applications
Bernoulli 12(6), 2006, 1077 1098 Retrospective exact simulation of diffusion sample paths with applications ALEXANDROS BESKOS*, OMIROS PAPASPILIOPOULOS** and GARETH O. ROBERTS Department of Mathematics
More informationBridging the Gap between Center and Tail for Multiscale Processes
Bridging the Gap between Center and Tail for Multiscale Processes Matthew R. Morse Department of Mathematics and Statistics Boston University BU-Keio 2016, August 16 Matthew R. Morse (BU) Moderate Deviations
More informationNonparametric Drift Estimation for Stochastic Differential Equations
Nonparametric Drift Estimation for Stochastic Differential Equations Gareth Roberts 1 Department of Statistics University of Warwick Brazilian Bayesian meeting, March 2010 Joint work with O. Papaspiliopoulos,
More informationDUBLIN CITY UNIVERSITY
DUBLIN CITY UNIVERSITY SAMPLE EXAMINATIONS 2017/2018 MODULE: QUALIFICATIONS: Simulation for Finance MS455 B.Sc. Actuarial Mathematics ACM B.Sc. Financial Mathematics FIM YEAR OF STUDY: 4 EXAMINERS: Mr
More informationIntroduction to Random Diffusions
Introduction to Random Diffusions The main reason to study random diffusions is that this class of processes combines two key features of modern probability theory. On the one hand they are semi-martingales
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 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 informationFast-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 informationSimulation Method for Solving Stochastic Differential Equations with Constant Diffusion Coefficients
Journal of mathematics and computer Science 8 (2014) 28-32 Simulation Method for Solving Stochastic Differential Equations with Constant Diffusion Coefficients Behrouz Fathi Vajargah Department of statistics,
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
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 informationExact and computationally efficient likelihood-based estimation for discretely observed diffusion processes
J. R. Statist. Soc. B (26) 68, Part 3, pp. 333 382 Exact and computationally efficient likelihood-based estimation for discretely observed diffusion processes Alexandros Beskos, Omiros Papaspiliopoulos,
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 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 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 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 informationarxiv: v1 [math.pr] 1 Jul 2013
ESTIMATION OF FIRST PASSAGE TIME DENSITIES OF DIFFUSIONS PROCESSES THROUGH TIME-VARYING BOUNDARIES arxiv:307.0336v [math.pr] Jul 03 Imene Allab and Francois Watier Department of Mathematics, Université
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 informationON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER
ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose
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 informationBayesian inference for stochastic differential mixed effects models - initial steps
Bayesian inference for stochastic differential ixed effects odels - initial steps Gavin Whitaker 2nd May 2012 Supervisors: RJB and AG Outline Mixed Effects Stochastic Differential Equations (SDEs) Bayesian
More informationParticle filters for partially observed diffusions
J. R. Statist. Soc. B (28) 7, Part 4, pp. 755 777 Particle filters for partially observed diffusions Paul Fearnhead Lancaster University, UK and Omiros Papaspiliopoulos and Gareth O. Roberts University
More informationSmall-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias
Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias José E. Figueroa-López 1 1 Department of Statistics Purdue University Seoul National University & Ajou University
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 informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationSome 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 informationLeast Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises
Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,
More informationlim n C1/n n := ρ. [f(y) f(x)], y x =1 [f(x) f(y)] [g(x) g(y)]. (x,y) E A E(f, f),
1 Part I Exercise 1.1. Let C n denote the number of self-avoiding random walks starting at the origin in Z of length n. 1. Show that (Hint: Use C n+m C n C m.) lim n C1/n n = inf n C1/n n := ρ.. Show that
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 informationPath Decomposition of Markov Processes. Götz Kersting. University of Frankfurt/Main
Path Decomposition of Markov Processes Götz Kersting University of Frankfurt/Main joint work with Kaya Memisoglu, Jim Pitman 1 A Brownian path with positive drift 50 40 30 20 10 0 0 200 400 600 800 1000-10
More informationOptimal transportation and optimal control in a finite horizon framework
Optimal transportation and optimal control in a finite horizon framework Guillaume Carlier and Aimé Lachapelle Université Paris-Dauphine, CEREMADE July 2008 1 MOTIVATIONS - A commitment problem (1) Micro
More informationLecture 4: Introduction to stochastic processes and stochastic calculus
Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London
More informationMore Empirical Process Theory
More Empirical Process heory 4.384 ime Series Analysis, Fall 2008 Recitation by Paul Schrimpf Supplementary to lectures given by Anna Mikusheva October 24, 2008 Recitation 8 More Empirical Process heory
More informationP. Fearnhead, O. Papaspiliopoulos, G.O. Roberts and A.M. Stuart, Random weight particle filtering of continuous time stochastic processes.
[84] P. Fearnhead, O. Papaspiliopoulos, G.O. Roberts and A.M. Stuart, Random weight particle filtering of continuous time stochastic processes. Journal of the Royal Statistical Society B. 72(4) (21) 497
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 informationSimulation and Parametric Estimation of SDEs
Simulation and Parametric Estimation of SDEs Jhonny Gonzalez School of Mathematics The University of Manchester Magical books project August 23, 2012 Motivation The problem Simulation of SDEs SDEs driven
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 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 informationHigher order weak approximations of stochastic differential equations with and without jumps
Higher order weak approximations of stochastic differential equations with and without jumps Hideyuki TANAKA Graduate School of Science and Engineering, Ritsumeikan University Rough Path Analysis and Related
More informationWeek 1 Quantitative Analysis of Financial Markets Distributions A
Week 1 Quantitative Analysis of Financial Markets Distributions A Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October
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 informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More informationfor all f satisfying E[ f(x) ] <.
. Let (Ω, F, P ) be a probability space and D be a sub-σ-algebra of F. An (H, H)-valued random variable X is independent of D if and only if P ({X Γ} D) = P {X Γ}P (D) for all Γ H and D D. Prove that if
More 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 informationOn 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 informationPerturbative Approaches for Robust Intertemporal Optimal Portfolio Selection
Perturbative Approaches for Robust Intertemporal Optimal Portfolio Selection F. Trojani and P. Vanini ECAS Course, Lugano, October 7-13, 2001 1 Contents Introduction Merton s Model and Perturbative Solution
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 informationUnderstanding Regressions with Observations Collected at High Frequency over Long Span
Understanding Regressions with Observations Collected at High Frequency over Long Span Yoosoon Chang Department of Economics, Indiana University Joon Y. Park Department of Economics, Indiana University
More informationHomework If the inverse T 1 of a closed linear operator exists, show that T 1 is a closed linear operator.
Homework 3 1 If the inverse T 1 of a closed linear operator exists, show that T 1 is a closed linear operator Solution: Assuming that the inverse of T were defined, then we will have to have that D(T 1
More informationI. ANALYSIS; PROBABILITY
ma414l1.tex Lecture 1. 12.1.2012 I. NLYSIS; PROBBILITY 1. Lebesgue Measure and Integral We recall Lebesgue measure (M411 Probability and Measure) λ: defined on intervals (a, b] by λ((a, b]) := b a (so
More informationMath 4381 / 6378 Symmetry Analysis
Math 438 / 6378 Smmetr Analsis Elementar ODE Review First Order Equations Ordinar differential equations of the form = F(x, ( are called first order ordinar differential equations. There are a variet of
More informationThe Chaotic Character of the Stochastic Heat Equation
The Chaotic Character of the Stochastic Heat Equation March 11, 2011 Intermittency The Stochastic Heat Equation Blowup of the solution Intermittency-Example ξ j, j = 1, 2,, 10 i.i.d. random variables Taking
More informationTools 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 informationDiscretization of Stochastic Differential Systems With Singular Coefficients Part II
Discretization of Stochastic Differential Systems With Singular Coefficients Part II Denis Talay, INRIA Sophia Antipolis joint works with Mireille Bossy, Nicolas Champagnat, Sylvain Maire, Miguel Martinez,
More informationThe 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 informationAnswers and expectations
Answers and expectations For a function f(x) and distribution P(x), the expectation of f with respect to P is The expectation is the average of f, when x is drawn from the probability distribution P E
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
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 informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
More informationPartial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula
Partial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula Group 4: Bertan Yilmaz, Richard Oti-Aboagye and Di Liu May, 15 Chapter 1 Proving Dynkin s formula
More informationNumerical Methods with Lévy Processes
Numerical Methods with Lévy Processes 1 Objective: i) Find models of asset returns, etc ii) Get numbers out of them. Why? VaR and risk management Valuing and hedging derivatives Why not? Usual assumption:
More informationIntroduction to numerical simulations for Stochastic ODEs
Introduction to numerical simulations for Stochastic ODEs Xingye Kan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616 August 9, 2010 Outline 1 Preliminaries 2 Numerical
More informationLecture on Parameter Estimation for Stochastic Differential Equations. Erik Lindström
Lecture on Parameter Estimation for Stochastic Differential Equations Erik Lindström Recap We are interested in the parameters θ in the Stochastic Integral Equations X(t) = X(0) + t 0 µ θ (s, X(s))ds +
More informationLecture 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 information1 Brownian Local Time
1 Brownian Local Time We first begin by defining the space and variables for Brownian local time. Let W t be a standard 1-D Wiener process. We know that for the set, {t : W t = } P (µ{t : W t = } = ) =
More informationLimit theorems for multipower variation in the presence of jumps
Limit theorems for multipower variation in the presence of jumps Ole E. Barndorff-Nielsen Department of Mathematical Sciences, University of Aarhus, Ny Munkegade, DK-8 Aarhus C, Denmark oebn@imf.au.dk
More informationAffine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009
Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine
More informationSupplement A: Construction and properties of a reected diusion
Supplement A: Construction and properties of a reected diusion Jakub Chorowski 1 Construction Assumption 1. For given constants < d < D let the pair (σ, b) Θ, where Θ : Θ(d, D) {(σ, b) C 1 (, 1]) C 1 (,
More informationApproximation of BSDEs using least-squares regression and Malliavin weights
Approximation of BSDEs using least-squares regression and Malliavin weights Plamen Turkedjiev (turkedji@math.hu-berlin.de) 3rd July, 2012 Joint work with Prof. Emmanuel Gobet (E cole Polytechnique) Plamen
More informationMarkov chain Monte Carlo algorithms for SDE parameter estimation
Markov chain Monte Carlo algorithms for SDE parameter estimation Andrew Golightly and Darren J. Wilkinson Abstract This chapter considers stochastic differential equations for Systems Biology models derived
More informationReview For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: = 0 : homogeneous equation.
Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: y y x y2 = 0 : homogeneous equation. x2 v = y dy, y = vx, and x v + x dv dx = v + v2. dx =
More informationAsymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½
University of Pennsylvania ScholarlyCommons Statistics Papers Wharton Faculty Research 1998 Asymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½ Lawrence D. Brown University
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 informationOptimal Stopping Problems and American Options
Optimal Stopping Problems and American Options Nadia Uys A dissertation submitted to the Faculty of Science, University of the Witwatersrand, in fulfilment of the requirements for the degree of Master
More informationAn adaptive numerical scheme for fractional differential equations with explosions
An adaptive numerical scheme for fractional differential equations with explosions Johanna Garzón Departamento de Matemáticas, Universidad Nacional de Colombia Seminario de procesos estocásticos Jointly
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationApplications of Optimal Stopping and Stochastic Control
Applications of and Stochastic Control YRM Warwick 15 April, 2011 Applications of and Some problems Some technology Some problems The secretary problem Bayesian sequential hypothesis testing the multi-armed
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 informationAsymptotic Coupling of an SPDE, with Applications to Many-Server Queues
Asymptotic Coupling of an SPDE, with Applications to Many-Server Queues Mohammadreza Aghajani joint work with Kavita Ramanan Brown University March 2014 Mohammadreza Aghajanijoint work Asymptotic with
More informationSolving the Poisson Disorder Problem
Advances in Finance and Stochastics: Essays in Honour of Dieter Sondermann, Springer-Verlag, 22, (295-32) Research Report No. 49, 2, Dept. Theoret. Statist. Aarhus Solving the Poisson Disorder Problem
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 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 informationContinuous Time Finance
Continuous Time Finance Lisbon 2013 Tomas Björk Stockholm School of Economics Tomas Björk, 2013 Contents Stochastic Calculus (Ch 4-5). Black-Scholes (Ch 6-7. Completeness and hedging (Ch 8-9. The martingale
More informationStochastic Viral Dynamics with Beddington-DeAngelis Functional Response
Stochastic Viral Dynamics with Beddington-DeAngelis Functional Response Junyi Tu, Yuncheng You University of South Florida, USA you@mail.usf.edu IMA Workshop in Memory of George R. Sell June 016 Outline
More informationBrownian 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 informationSome Tools From Stochastic Analysis
W H I T E Some Tools From Stochastic Analysis J. Potthoff Lehrstuhl für Mathematik V Universität Mannheim email: potthoff@math.uni-mannheim.de url: http://ls5.math.uni-mannheim.de To close the file, click
More informationConvergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics
Convergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics Meng Xu Department of Mathematics University of Wyoming February 20, 2010 Outline 1 Nonlinear Filtering Stochastic Vortex
More informationTyler Hofmeister. University of Calgary Mathematical and Computational Finance Laboratory
JUMP PROCESSES GENERALIZING STOCHASTIC INTEGRALS WITH JUMPS Tyler Hofmeister University of Calgary Mathematical and Computational Finance Laboratory Overview 1. General Method 2. Poisson Processes 3. Diffusion
More informationAsymptotical distribution free test for parameter change in a diffusion model (joint work with Y. Nishiyama) Ilia Negri
Asymptotical distribution free test for parameter change in a diffusion model (joint work with Y. Nishiyama) Ilia Negri University of Bergamo (Italy) ilia.negri@unibg.it SAPS VIII, Le Mans 21-24 March,
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 informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationWeak convergence and large deviation theory
First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov
More informationON 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 informationElementary ODE Review
Elementary ODE Review First Order ODEs First Order Equations Ordinary differential equations of the fm y F(x, y) () are called first der dinary differential equations. There are a variety of techniques
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 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 information