A Central Limit Theorem for Fleming-Viot Particle Systems Application to the Adaptive Multilevel Splitting Algorithm
|
|
- Rosemary Johns
- 5 years ago
- Views:
Transcription
1 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 1 3 Université Pierre et Marie Curie, Paris London 218
2 Rare event simulation Dynamical rare event problem s Y s denotes a Markov continuous (diffusion) process in F ; A and B denote two sets in F. Problem: Simulate a reactive trajectory, i.e. that starts from y close to A and reaches B before going back to A. Problem: Rare event simulation p 1.
3 Rare event simulation Use of a Reaction Coordinate Reaction coordinate: one-dimensional continuous function ξ : F R Notation: {ξ > t} :={y F, ξ(y) > z} S t := inf(s ξ(y s ) = t) S A := inf(s Y s A) For simplicity assume A {ξ < 1}, B {ξ > 1}, and ξ(y ) =. so that the considered rare event becomes {S 1 > S A }, p = P(S 1 > S A )
4 Rare event simulation Reaction Coordinate Define the score of each trajectory by score = τ = sup ξ(y s ) s S A Multilevel/Importance Splitting: clone the trajectories with highest score, kill the other ones [Kahn and Harris (1951)].
5 AMS algorithm Definition (AMS algorithm) At step j of the algorithm, we define a particle system (Y 1,j,, Y N,j ) [C(R +, F )] N defined by the following set of rules. Initialization j = : consider N i.i.d. particles with law (Y s ) s SA. At step j 1, kill the particle N j with minimal score τ j Splitting: the killed particle is given the state of particle M j, uniformly picked among the (N 1) remaining particles. The particle is then resampled starting from the entrance time in level τ j : S M j,j 1 τ j. Set j 1 j. And so on until all particles reach {ξ > 1}.
6
7
8
9
10
11
12
13
14 AMS yields Unbiased Estimator of Unormalized Quantity Denote by J t the first iteration of branchings/steps so that all particles reach {ξ > t}. Denote empirical measure of particles paths at iteration J t N η N t := 1 N j=1 δ Y n,jt s Unbiased Estimates: For any pathwise observable (test function) ϕ, we have γ t (ϕ) := E[ϕ(Y )1 St<S A ]= E[(1 1/N) Jt η N t (ϕ)]. Weak Law of Large Number? Central Limit Theorem?
15 Mapping to a Fleming-Viot particle system The level-indexed process Key Idea: Assuming Y is strong Markov, define a new Markov process called the level-indexed process in E := F { }, the state space enhanced by a cemetary, by setting (NB: ξ(y ) = ): { X t := Y St if S t < S A X t := else Can be extended to a time-homogenous process by setting X h := Y Sξ(x )+h if S ξ(x )+h < S A Can be extended to the pathwise case by setting E := C(R +, F ) { } X t := Y s St if S t < S A
16 Mapping to a Fleming-Viot particle system The level-indexed process is a killed Process Lemma The level-indexed process X = (X t ) t is a Markov process in F { }, where the cemetary / F is absorbing (a trap). Moreover, the killing time τ satisfies τ = inf{t, X t = } = sup ξ(y s ) = score. s S A Moreover the conditional distribution on state space F of interest becomes η t = L(X t τ > t) = L(Y S t < S A ), and the probability of the rare event p t = P(S t < S A ) = P(τ > t).
17 Mapping to a Fleming-Viot particle system Fleming-Viot Particle System Definition (Fleming-Viot particle system) The Fleming-Viot particle system (X 1 t,, X N t ) t [,T ] is the Markov process with state space F N defined by the rules Initialization: consider N i.i.d. particles X 1,..., X N i.i.d. η, Evolution and killing: each particle evolves independently according to the law of the underlying Markov process X until one of them hits, Splitting: the killed particle is taken from, and is given instantaneously the state of one of the (N 1) other particles (randomly uniformly chosen). and so on until final time T.
18 Mapping to a Fleming-Viot particle system Lemma The AMS algorithm mapped with the level-indexed process and stopped at iteration J t (the first iteration when all particles reach {ξ > t}) is actually a Fleming-Viot particle system. Proof Picture.
19 Mapping to a Fleming-Viot particle system Lemma The AMS algorithm mapped with the level-indexed process and stopped at iteration J t (the first iteration when all particles reach {ξ > t}) is actually a Fleming-Viot particle system. Proof Picture.
20 Mapping to a Fleming-Viot particle system Lemma The AMS algorithm mapped with the level-indexed process and stopped at iteration J t (the first iteration when all particles reach {ξ > t}) is actually a Fleming-Viot particle system. Proof Picture.
21 Mapping to a Fleming-Viot particle system Lemma The AMS algorithm mapped with the level-indexed process and stopped at iteration J t (the first iteration when all particles reach {ξ > t}) is actually a Fleming-Viot particle system. Proof Picture.
22 Mapping to a Fleming-Viot particle system Lemma The AMS algorithm mapped with the level-indexed process and stopped at iteration J t (the first iteration when all particles reach {ξ > t}) is actually a Fleming-Viot particle system. Proof Picture.
23 Mapping to a Fleming-Viot particle system Lemma The AMS algorithm mapped with the level-indexed process and stopped at iteration J t (the first iteration when all particles reach {ξ > t}) is actually a Fleming-Viot particle system. Proof Picture.
24 Mapping to a Fleming-Viot particle system Lemma The AMS algorithm mapped with the level-indexed process and stopped at iteration J t (the first iteration when all particles reach {ξ > t}) is actually a Fleming-Viot particle system. Proof Picture.
25 Mapping to a Fleming-Viot particle system Lemma The AMS algorithm mapped with the level-indexed process and stopped at iteration J t (the first iteration when all particles reach {ξ > t}) is actually a Fleming-Viot particle system. Proof Picture.
26 Mapping to a Fleming-Viot particle system Some Important Estimators Normalized distribution proba of rare event un-normalized distribution η t = L(X t X t ) p t = P(X t ) γ t := p t η t Empirical Empirical Empirical ηt N = 1 N N i=1 δ Xt i t = ( 1 1 ) NNt N, γt N = pt N ηt N. N T := J T /N is the average number of branchings per particle system.
27 Main Results Regularity Assumption For x F { }, t T, consider the sub-markovian semi-group Q t (ϕ)(x) := E[ϕ(X t )1 τ >t X = x]. Assumption (Non-synchronous jumps) For any initial condition X t = x F and any ϕ C b (F ): (i) the jump times of the càdlàg version of the martingale process t L t := Q T t (ϕ)(x t ) have an atomless distribution: P(L t L t X = x) = t. (ii) The killing time has also an atomless distribution.
28 Main Results Non-explosion Assumption The non-synchronous jumps Assumption is morally equivalent to: martingale jumps and branchings in the Fleming-Viot system are never simultaneous. In addition we ask: Assumption (Non-explosion) The Fleming-Viot system is non-explosive in the sense that the number of branching at any finite time is almost surely finite P(N T < + ) = 1.
29 Main Results Example with Hard Obstacle (The originality of our result!) t X t process in E, and let F be an open domain with boundary F = F \ F. Let τ be the hitting time of E \ F. Proposition (i) The hitting time τ on boundary has an atomless distribution. (ii) The process starting from boundary F exits in open E \ F immediately almost surely. Then Assumption no synchronous jumps is fullfilled. Proposition (Grigorescu and Kang, 212) Assume that F is open and bounded in R d with smooth boundary F, and that X is a diffusion with smooth and uniformly elliptic coefficients. Then Assumption non explosion holds true, as well as Assumption no synchronuous jumps.
30 Main Results Central Limit Theorems Theorem Under Assumptions non-synchronous jumps and non-explosion, for any ϕ C b (F ) one has the convergence in distribution ( ) N γt N (ϕ) γ D T (ϕ) N (, N σ2 T (ϕ)), where by definition σ 2 T (ϕ) =p2 T Var η T (ϕ) p 2 T ln(p T ) η T (ϕ) 2 2 T Var ηt (Q T t (ϕ))p t dp t.
31 Main Results Central Limit Theorems By Slutsky Lemma: Corollary Under Assumptions non-synchronous jumps and non-explosion, for any ϕ C b (F ), one has the convergence in distribution ( ) N ηt N (ϕ) η D T (ϕ) N (, N σ2 T (ϕ η T (ϕ))/pt 2 ). Besides, where ( ) N pt N p D T N (, N σ2 ), σ 2 = σ 2 T (1 F ) = p 2 T ln(p T ) 2 T Var ηt (Q T t (1 F ))p t dp t.
32 Main Results Remarks on Asymptotic Variances The variance is the limit of the CLT of a discrete time Fleming-Viot particle system (CLT in Del Moral s book). Bounds on the probability estimator p 2 T log(p T ) σ 2 2p T (1 p T ) + p 2 T log(p T ) Lower bound: (soft) killing time independent of the state. Upper bound: killing time mainly dependent on the initial condition. Dominant term 2p T (1 p T ) is twice the naive Monte Carlo variance.
33 Proof Stochastic calculus with jumps Recall that one can integrate with respect ot semi-martingales X = montonous processes + martingales as follows: Y t dx t Y t (X t+dt X t ) We then have the chain rule d(x t Y t ) = Y t dx t + X t dy t + d[x, Y ] t where t [X, X ] t is an increasing process, bilinear with respect to vector space structure on X called the quadratic variation. Broadly speaking [X, X ] t = P lim t i+1 t i (X ti+1 X ti ) 2 If X is monotonous, [X, X ] t is the sum of the squares of the jumps. i
34 Proof Stochastic calculus with jumps If t M t is a (local) martingale, then t M 2 t [M, M] t is again a local martingale. In the presence of jumps there ar plenty of quadratic variation -like increasing processes t i(m) t such that t M 2 t i(m) t is a local martingale. For instance there is a unique i(m) t =< M, M > t which is predictable. Example: let t M t { 1, 1} be Poisson-like Markov + Martingale random walk process jumping up or down with proba 1/2 at indep. expo. times. Then [M, M] t = jumps Var(jump) = 1 < M, M > t = dt
35 Proof CLT for martingales with jumps Theorem (Martinagle CLT (Ethier-Kurtz)) On a filtered probability space, let t mt N denote a sequence of càdlàg local martingales indexed by N 1. Assume moreover that (i) m N D µ, where µ is a given probability on R. N + (ii) Vanishing jumps: One has lim N + E[sup t [,T ] m N t mt N 2 ] =. (iii) For each N, there exists an increasing càdlàg process t it N such that ( ) t mt N m N 2 it N is a local martingale. (iv) Vanishing jump: The process t ] it N satisfies lim N + E [sup t [,T ] it N i Nt =.
36 Proof CLT for martingales with jumps Theorem (Martinagle CLT (Ethier-Kurtz)) For the increasing càdlàg process t it N (with vanishing jumps) such that ( ) t mt N m N 2 it N is a local martingale: (v)!! Main Assumption!!: There is a cont. and incr. det. function t i t s. t., t [, T ], it N P i t. N + Then (m N t ) t [,T ] converges in law (under the Skorokhod topology) to (M t ) t [,T ], where M µ and (M t M ) t [,T ] is a Gaussian martingale, independent of M, with independent increments and variance function i t (time changed Brownian motion).
37 Proof CLT for martingales with jumps In short, we need to construct martingales of order 1/ N from the particle system and ensure the convergence of a quadratic variation of those martingales of order 1/N.
38 Proof Overview of the proof Key object: the càdlàg martingale t γ N t (Q) := γ N t ( ) Q T t (ϕ). Initial condition treated separately (easy). We will handle the distribution of γt N (Q) in the limit N by using a Central Limit Theorem for continuous time martingales. Not straightforward: the convergence of the quadratic variation N[γ N (Q), γ N (Q)] t is difficult (lots of IPPs!!).
39 Proof Martingale decomposition [Villemonais 214] The key martingale decomposition is the following: γ N t (Q) = γ N (Q) + 1 N t p N u (dm u + dm u ). With M t := 1 N N n=1 Mn t and M t := 1 N N n=1 Mn t. M n t is the martingale contribution except for branching times of particle n. M n t is the martingale contribution at branchings only of particle n. No ambiguity, natural way to do this.
40 Proof Orthogonality The 2N martingales {M n t, M m t } 1 n,m N are mutually orthogonal. More specifically (i) [M, M] t is a local martingale, (ii) [M, M] t = 1 N N [M n, M n ] t, n=1 (iii) Moreover, if we note the intermediate quadratic variation (M, M) t = 1 N N [M n, M n ] t, n=1 then the process [M, M] t (M, M) t is also a local martingale.
41 Proof Ingredient (i): A key formula Lemma The quadratic variation of martingales associated with the particles dynamics outside branchings can be related to through the key formula γ N t (Q 2 ) = γ N t ([Q T t (ϕ)] 2 ) p N t d(m, M) t = dγ N t (Q 2 ) + Martingale
42 Proof Ingredient (ii): L 2 apriori estimates Proposition (Villemonais 214, CDGR 217) For any ϕ D, we have [( ) ] E γt N (ϕ) γ T (ϕ) 2 6 ϕ 2 N. Proof. γ N T (ϕ) γ T (ϕ) = 1 N T p N t dm t + 1 N T + γ N (Q T ϕ) γ (Q T ϕ), p N t dm t (i) Initial condition is OK by independence.
43 Proof L 2 estimates (ii) M-terms. Using Ito s isometry and d[m, M] t 4 ϕ 2 dn t, we obtain [ ( T ) 2 ] [ T ] E pt N ( ) dm t = E p N 2d[M, t M]t 4 ϕ 2 1 ( ) 1 1 2(j 1) N N 4 ϕ 2. (iii) M-terms. In the same way, applying Ito s isometry and the key formula pt N d(m, M) t = dγt N (Q 2 ) + Martingale, we get [ ( T ) 2 ] [ T ] E pt N ( ) dm t = E p N 2d[M, t M]t [ T E j=1 ] [ ] pt N d(m, M) t = E γt N (Q2 ) ϕ 2.
44 Proof INgredient (iii): Time uniform a priori estimate of p N t Lemma One has sup t [,T ] pt N P p t N. Proof. Independent of the context.t p t is continuous on [, T ] by construction, it is cleat that t p N t is decreasing for all N 2. The Lemma results of last Proposition and a from a probabilistic version of Second Dini (or Pólya) theorem: if a sequence of monotone functions converges pointwise on a compact interval and if the limit function is also continuous, then the convergence is uniform on that interval.
45 Proof Increasing Process in general CLT In order to use martingale CLT, we need a càdlàg increasing process it N such that (γ N (Q) t ) 2 it N is a martingale (quadratic variation - like). After tedious computations and trials and errors we chose a quadratic variation with only the branching jumps integrated. i N t = t n=1 k=1 ( p N u ) 2 d(m, M) u + 1 N t ( ) p N 2dRu u. t NB: with the rest term is O(1/N) with ( N + (1 ) R t = 1 2 N Var (n) η τ n,k Var η N u (Q)pN u dp N u (Q) Var η N (Q) τ n,k ) 1 t τn,k.
46 Proof Integration by parts formulas Let t zt N be any càdlàg semi-martingale, c > a deterministic constant, and assume for any branching time τ j, j 1: If zτ N j c/n,then t p N s dz N s = p N t z N t z N If z N c(1 1/N) j, τ j t z N s (p N s ) 1 dp N s = If z N τ j c(1 1/N) j /N, t t z N s dp N s z N s d ln p N s + O(1/N). + O(1/N). t z N s d ln p N s = z N t ln p N t t ln p N s dz N s + O(1/N).
47 Proof Integration by parts for i N t Using the intergation by parts formula abvoe and the key formula : Lemma The increasing process it N can be integrated by parts and be rewritten as it N = pt N γt N ( Q 2 ) γ N ( Q 2 ) [ + γt (Q)] N 2 ln pt N 2 t γ N u (Q 2 )dp N u + O ( 1 N ).
48 Proof Convergence of i N t and final proof of the CLT By the non synchronous jump Assumption all the vanishing jump assumptions in the martingale CLT are verified. The only remaining remaining part to prove is the following: Proposition For any t [, T ], one has i N t = p N t γ N t i N t ( Q 2 ) γ N P N i t(ϕ). where ( Q 2 ) +[ γ N t (Q)] 2 ln p N t 2 t γ N u (Q 2 )dp N u +O ( 1 N t i t (ϕ) = p t γ t (Q 2 ) γ (Q 2 ) + [γ t (Q)] 2 ( ln p t 2 γ u Q 2 ) dp u.
49 Proof Proof of Convergence of i N t After some calculations, all but one term can be treated using the L 2 a priori convergence estimate. The only remaing problem is the following: i N t i t (ϕ) =easy converging terms with L 2 -estimate and IPP+ 2 t (p N u p u ) d γ N u (Q 2 ) It is difficult to prove its convergence to because the L 2 estimate is only pointwise. Hence handling the integrator is cumbersome!!.
50 Proof Convergence of i N t This difficult term is then treated as follows. Use the key formula t t (pu N p u ) d γu N (Q 2 ) = (pu N p u )pu N d(m, M) u +O ( 1 ) N Since (M, M) is an increasing process, it comes t ( t (pu N p u )pu N d(m, M) u sup pu N p u u The key formula back again implies: [ t ] [ ] E pu N d(m, M) u = E γt N (Q 2 ) ϕ 2. The a priori uniform estimate implies convergence of sup u p N u p u. p N u d(m, M) u ).
51 Proof Arxiv Bernard Delyon, Frédéric Cérou, Arnaud Guyader, Mathias Rousset. A Central Limit Theorem for Fleming-Viot Particle Systems with Hard Killing. 217 Bernard Delyon, Frédéric Cérou, Arnaud Guyader, Mathias Rousset. Asymptotic normality of the AMS algorithm. 218
A Central Limit Theorem for Fleming-Viot Particle Systems 1
A Central Limit Theorem for Fleming-Viot Particle Systems 1 Frédéric Cérou 2 IRIA Rennes & IRMAR, France frederic.cerou@inria.fr Bernard Delyon Université Rennes 1 & IRMAR, France bernard.delyon@univ-rennes1.fr
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 informationImportance splitting for rare event simulation
Importance splitting for rare event simulation F. Cerou Frederic.Cerou@inria.fr Inria Rennes Bretagne Atlantique Simulation of hybrid dynamical systems and applications to molecular dynamics September
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 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 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 informationTwo viewpoints on measure valued processes
Two viewpoints on measure valued processes Olivier Hénard Université Paris-Est, Cermics Contents 1 The classical framework : from no particle to one particle 2 The lookdown framework : many particles.
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 informationStochastic 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 information6. Brownian Motion. Q(A) = P [ ω : x(, ω) A )
6. Brownian Motion. stochastic process can be thought of in one of many equivalent ways. We can begin with an underlying probability space (Ω, Σ, P) and a real valued stochastic process can be defined
More informationConvergence of Feller Processes
Chapter 15 Convergence of Feller Processes This chapter looks at the convergence of sequences of Feller processes to a iting process. Section 15.1 lays some ground work concerning weak convergence of processes
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 informationPathwise construction of tree-valued Fleming-Viot processes
Pathwise construction of tree-valued Fleming-Viot processes Stephan Gufler November 9, 2018 arxiv:1404.3682v4 [math.pr] 27 Dec 2017 Abstract In a random complete and separable metric space that we call
More informationBernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012
1 Stochastic Calculus Notes March 9 th, 1 In 19, Bachelier proposed for the Paris stock exchange a model for the fluctuations affecting the price X(t) of an asset that was given by the Brownian motion.
More 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 informationLECTURE 2: LOCAL TIME FOR BROWNIAN MOTION
LECTURE 2: LOCAL TIME FOR BROWNIAN MOTION We will define local time for one-dimensional Brownian motion, and deduce some of its properties. We will then use the generalized Ray-Knight theorem proved in
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 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 informationIntroduction to Rare Event Simulation
Introduction to Rare Event Simulation Brown University: Summer School on Rare Event Simulation Jose Blanchet Columbia University. Department of Statistics, Department of IEOR. Blanchet (Columbia) 1 / 31
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 informationErgodic Theorems. Samy Tindel. Purdue University. Probability Theory 2 - MA 539. Taken from Probability: Theory and examples by R.
Ergodic Theorems Samy Tindel Purdue University Probability Theory 2 - MA 539 Taken from Probability: Theory and examples by R. Durrett Samy T. Ergodic theorems Probability Theory 1 / 92 Outline 1 Definitions
More informationSTAT 331. Martingale Central Limit Theorem and Related Results
STAT 331 Martingale Central Limit Theorem and Related Results In this unit we discuss a version of the martingale central limit theorem, which states that under certain conditions, a sum of orthogonal
More informationStochastic integration. P.J.C. Spreij
Stochastic integration P.J.C. Spreij this version: April 22, 29 Contents 1 Stochastic processes 1 1.1 General theory............................... 1 1.2 Stopping times...............................
More informationFiltrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition
Filtrations, Markov Processes and Martingales Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition David pplebaum Probability and Statistics Department,
More 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 informationSome recent improvements to importance splitting
Some recent improvements to importance splitting F. Cerou IRISA / INRIA Campus de Beaulieu 35042 RENNES Cédex, France Frederic.Cerou@inria.fr F. Le Gland IRISA / INRIA Campus de Beaulieu 35042 RENNES Cédex,
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 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 informationAsymptotics for posterior hazards
Asymptotics for posterior hazards Pierpaolo De Blasi University of Turin 10th August 2007, BNR Workshop, Isaac Newton Intitute, Cambridge, UK Joint work with Giovanni Peccati (Université Paris VI) and
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 informationSTAT Sample Problem: General Asymptotic Results
STAT331 1-Sample Problem: General Asymptotic Results In this unit we will consider the 1-sample problem and prove the consistency and asymptotic normality of the Nelson-Aalen estimator of the cumulative
More informationLecture 2. We now introduce some fundamental tools in martingale theory, which are useful in controlling the fluctuation of martingales.
Lecture 2 1 Martingales We now introduce some fundamental tools in martingale theory, which are useful in controlling the fluctuation of martingales. 1.1 Doob s inequality We have the following maximal
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 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 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 informationMean field simulation for Monte Carlo integration. Part II : Feynman-Kac models. P. Del Moral
Mean field simulation for Monte Carlo integration Part II : Feynman-Kac models P. Del Moral INRIA Bordeaux & Inst. Maths. Bordeaux & CMAP Polytechnique Lectures, INLN CNRS & Nice Sophia Antipolis Univ.
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 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 informationLecture 17 Brownian motion as a Markov process
Lecture 17: Brownian motion as a Markov process 1 of 14 Course: Theory of Probability II Term: Spring 2015 Instructor: Gordan Zitkovic Lecture 17 Brownian motion as a Markov process Brownian motion is
More informationSTOCHASTIC CALCULUS JASON MILLER AND VITTORIA SILVESTRI
STOCHASTIC CALCULUS JASON MILLER AND VITTORIA SILVESTRI Contents Preface 1 1. Introduction 1 2. Preliminaries 4 3. Local martingales 1 4. The stochastic integral 16 5. Stochastic calculus 36 6. Applications
More informationLecture 5. If we interpret the index n 0 as time, then a Markov chain simply requires that the future depends only on the present and not on the past.
1 Markov chain: definition Lecture 5 Definition 1.1 Markov chain] A sequence of random variables (X n ) n 0 taking values in a measurable state space (S, S) is called a (discrete time) Markov chain, if
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 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 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 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 informationSome functional (Hölderian) limit theorems and their applications (II)
Some functional (Hölderian) limit theorems and their applications (II) Alfredas Račkauskas Vilnius University Outils Statistiques et Probabilistes pour la Finance Université de Rouen June 1 5, Rouen (Rouen
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 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 informationThe strictly 1/2-stable example
The strictly 1/2-stable example 1 Direct approach: building a Lévy pure jump process on R Bert Fristedt provided key mathematical facts for this example. A pure jump Lévy process X is a Lévy process such
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 informationViscosity Solutions of Path-dependent Integro-Differential Equations
Viscosity Solutions of Path-dependent Integro-Differential Equations Christian Keller University of Southern California Conference on Stochastic Asymptotics & Applications Joint with 6th Western Conference
More informationKrzysztof Burdzy Robert Ho lyst Peter March
A FLEMING-VIOT PARTICLE REPRESENTATION OF THE DIRICHLET LAPLACIAN Krzysztof Burdzy Robert Ho lyst Peter March Abstract: We consider a model with a large number N of particles which move according to independent
More informationModern Discrete Probability Branching processes
Modern Discrete Probability IV - Branching processes Review Sébastien Roch UW Madison Mathematics November 15, 2014 1 Basic definitions 2 3 4 Galton-Watson branching processes I Definition A Galton-Watson
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 informationItô s excursion theory and random trees
Itô s excursion theory and random trees Jean-François Le Gall January 3, 200 Abstract We explain how Itô s excursion theory can be used to understand the asymptotic behavior of large random trees. We provide
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 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 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 informationStochastic 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 informationSUMMARY OF RESULTS ON PATH SPACES AND CONVERGENCE IN DISTRIBUTION FOR STOCHASTIC PROCESSES
SUMMARY OF RESULTS ON PATH SPACES AND CONVERGENCE IN DISTRIBUTION FOR STOCHASTIC PROCESSES RUTH J. WILLIAMS October 2, 2017 Department of Mathematics, University of California, San Diego, 9500 Gilman Drive,
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 informationOn the submartingale / supermartingale property of diffusions in natural scale
On the submartingale / supermartingale property of diffusions in natural scale Alexander Gushchin Mikhail Urusov Mihail Zervos November 13, 214 Abstract Kotani 5 has characterised the martingale property
More informationMath 6810 (Probability) Fall Lecture notes
Math 6810 (Probability) Fall 2012 Lecture notes Pieter Allaart University of North Texas September 23, 2012 2 Text: Introduction to Stochastic Calculus with Applications, by Fima C. Klebaner (3rd edition),
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 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 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 informationPart III Stochastic Calculus and Applications
Part III Stochastic Calculus and Applications Based on lectures by R. Bauerschmidt Notes taken by Dexter Chua Lent 218 These notes are not endorsed by the lecturers, and I have modified them often significantly
More informationExistence, Uniqueness and Stability of Invariant Distributions in Continuous-Time Stochastic Models
Existence, Uniqueness and Stability of Invariant Distributions in Continuous-Time Stochastic Models Christian Bayer and Klaus Wälde Weierstrass Institute for Applied Analysis and Stochastics and University
More informationConvergence exponentielle uniforme vers la distribution quasi-stationnaire en dynamique des populations
Convergence exponentielle uniforme vers la distribution quasi-stationnaire en dynamique des populations Nicolas Champagnat, Denis Villemonais Colloque Franco-Maghrébin d Analyse Stochastique, Nice, 25/11/2015
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 informationLecture 19 L 2 -Stochastic integration
Lecture 19: L 2 -Stochastic integration 1 of 12 Course: Theory of Probability II Term: Spring 215 Instructor: Gordan Zitkovic Lecture 19 L 2 -Stochastic integration The stochastic integral for processes
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 informationLearning Session on Genealogies of Interacting Particle Systems
Learning Session on Genealogies of Interacting Particle Systems A.Depperschmidt, A.Greven University Erlangen-Nuremberg Singapore, 31 July - 4 August 2017 Tree-valued Markov processes Contents 1 Introduction
More informationLogarithmic scaling of planar random walk s local times
Logarithmic scaling of planar random walk s local times Péter Nándori * and Zeyu Shen ** * Department of Mathematics, University of Maryland ** Courant Institute, New York University October 9, 2015 Abstract
More 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 informationNon-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings. Krzysztof Burdzy University of Washington
Non-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings Krzysztof Burdzy University of Washington 1 Review See B and Kendall (2000) for more details. See also the unpublished
More 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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 22 12/09/2013. Skorokhod Mapping Theorem. Reflected Brownian Motion
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.7J Fall 213 Lecture 22 12/9/213 Skorokhod Mapping Theorem. Reflected Brownian Motion Content. 1. G/G/1 queueing system 2. One dimensional reflection mapping
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 informationn E(X t T n = lim X s Tn = X s
Stochastic Calculus Example sheet - Lent 15 Michael Tehranchi Problem 1. Let X be a local martingale. Prove that X is a uniformly integrable martingale if and only X is of class D. Solution 1. If If direction:
More informationAn Ergodic Theorem for Fleming-Viot Models in Random Environments 1
An Ergodic Theorem for Fleming-Viot Models in Random Environments 1 Arash Jamshidpey 2 216 Abstract The Fleming-Viot (FV) process is a measure-valued diffusion that models the evolution of type frequencies
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 informationSequential Monte Carlo Samplers for Applications in High Dimensions
Sequential Monte Carlo Samplers for Applications in High Dimensions Alexandros Beskos National University of Singapore KAUST, 26th February 2014 Joint work with: Dan Crisan, Ajay Jasra, Nik Kantas, Alex
More informationMath 735: Stochastic Analysis
First Prev Next Go To Go Back Full Screen Close Quit 1 Math 735: Stochastic Analysis 1. Introduction and review 2. Notions of convergence 3. Continuous time stochastic processes 4. Information and conditional
More informationA 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 informationA Note on Auxiliary Particle Filters
A Note on Auxiliary Particle Filters Adam M. Johansen a,, Arnaud Doucet b a Department of Mathematics, University of Bristol, UK b Departments of Statistics & Computer Science, University of British Columbia,
More informationOn Stopping Times and Impulse Control with Constraint
On Stopping Times and Impulse Control with Constraint Jose Luis Menaldi Based on joint papers with M. Robin (216, 217) Department of Mathematics Wayne State University Detroit, Michigan 4822, USA (e-mail:
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 informationBernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010
1 Stochastic Calculus Notes Abril 13 th, 1 As we have seen in previous lessons, the stochastic integral with respect to the Brownian motion shows a behavior different from the classical Riemann-Stieltjes
More 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 informationSequential Monte Carlo Methods for Bayesian Computation
Sequential Monte Carlo Methods for Bayesian Computation A. Doucet Kyoto Sept. 2012 A. Doucet (MLSS Sept. 2012) Sept. 2012 1 / 136 Motivating Example 1: Generic Bayesian Model Let X be a vector parameter
More informationItô 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 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 informationThe range of tree-indexed random walk
The range of tree-indexed random walk Jean-François Le Gall, Shen Lin Institut universitaire de France et Université Paris-Sud Orsay Erdös Centennial Conference July 2013 Jean-François Le Gall (Université
More informationA Backward Particle Interpretation of Feynman-Kac Formulae
A Backward Particle Interpretation of Feynman-Kac Formulae P. Del Moral Centre INRIA de Bordeaux - Sud Ouest Workshop on Filtering, Cambridge Univ., June 14-15th 2010 Preprints (with hyperlinks), joint
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 7 9/25/2013
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.070J Fall 013 Lecture 7 9/5/013 The Reflection Principle. The Distribution of the Maximum. Brownian motion with drift Content. 1. Quick intro to stopping times.
More informationDoléans measures. Appendix C. C.1 Introduction
Appendix C Doléans measures C.1 Introduction Once again all random processes will live on a fixed probability space (Ω, F, P equipped with a filtration {F t : 0 t 1}. We should probably assume the filtration
More informationConcentration inequalities for Feynman-Kac particle models. P. Del Moral. INRIA Bordeaux & IMB & CMAP X. Journées MAS 2012, SMAI Clermond-Ferrand
Concentration inequalities for Feynman-Kac particle models P. Del Moral INRIA Bordeaux & IMB & CMAP X Journées MAS 2012, SMAI Clermond-Ferrand Some hyper-refs Feynman-Kac formulae, Genealogical & Interacting
More informationWeakly interacting particle systems on graphs: from dense to sparse
Weakly interacting particle systems on graphs: from dense to sparse Ruoyu Wu University of Michigan (Based on joint works with Daniel Lacker and Kavita Ramanan) USC Math Finance Colloquium October 29,
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 information