LogFeller et Ray Knight
|
|
- Mercy Blake
- 5 years ago
- Views:
Transcription
1 LogFeller et Ray Knight Etienne Pardoux joint work with V. Le and A. Wakolbinger Etienne Pardoux (Marseille) MANEGE, 18/1/1 1 / 16
2 Feller s branching diffusion with logistic growth We consider the diffusion dz t = Z t (θ γz t )dt + 2 Z t dw t, t, (1) with Z = x, θ, γ >. Z t is not a branching process : the quadratic term introduces interactions between the branches. However there is still an interpretation in terms of the evolution of a population of the solution of the Fellog process Z t, see A. Lambert (25). Consider dh s = H =, ( ) θ 2 γl s(h s ) ds + db s dl s(), s, (2) where B is a standard Brownian motion. In this SDE, the term dl s ()/2 takes care of the reflection of H at the origin. Etienne Pardoux (Marseille) MANEGE, 18/1/1 2 / 16
3 Feller s branching diffusion with logistic growth We consider the diffusion dz t = Z t (θ γz t )dt + 2 Z t dw t, t, (1) with Z = x, θ, γ >. Z t is not a branching process : the quadratic term introduces interactions between the branches. However there is still an interpretation in terms of the evolution of a population of the solution of the Fellog process Z t, see A. Lambert (25). Consider dh s = H =, ( ) θ 2 γl s(h s ) ds + db s dl s(), s, (2) where B is a standard Brownian motion. In this SDE, the term dl s ()/2 takes care of the reflection of H at the origin. Etienne Pardoux (Marseille) MANEGE, 18/1/1 2 / 16
4 Feller s branching diffusion with logistic growth We consider the diffusion dz t = Z t (θ γz t )dt + 2 Z t dw t, t, (1) with Z = x, θ, γ >. Z t is not a branching process : the quadratic term introduces interactions between the branches. However there is still an interpretation in terms of the evolution of a population of the solution of the Fellog process Z t, see A. Lambert (25). Consider dh s = H =, ( ) θ 2 γl s(h s ) ds + db s dl s(), s, (2) where B is a standard Brownian motion. In this SDE, the term dl s ()/2 takes care of the reflection of H at the origin. Etienne Pardoux (Marseille) MANEGE, 18/1/1 2 / 16
5 Results Using Girsanov s theorem, one can show Proposition The SDE (2) has a unique weak solution. H being the solution of (2), L s (t) denoting its local time, let S x := inf{s > : L s () > x}. Our main result is the Theorem L Sx (t), t, solves the Fellog SDE (1). Etienne Pardoux (Marseille) MANEGE, 18/1/1 3 / 16
6 Results Using Girsanov s theorem, one can show Proposition The SDE (2) has a unique weak solution. H being the solution of (2), L s (t) denoting its local time, let S x := inf{s > : L s () > x}. Our main result is the Theorem L Sx (t), t, solves the Fellog SDE (1). Etienne Pardoux (Marseille) MANEGE, 18/1/1 3 / 16
7 Results Using Girsanov s theorem, one can show Proposition The SDE (2) has a unique weak solution. H being the solution of (2), L s (t) denoting its local time, let S x := inf{s > : L s () > x}. Our main result is the Theorem L Sx (t), t, solves the Fellog SDE (1). Etienne Pardoux (Marseille) MANEGE, 18/1/1 3 / 16
8 Discrete approximation Our proof of the above extended Ray Knight theorem is based on an approximation by finite population. For N N, let Z N,x be the total mass of a population of individuals, each of which has mass 1/N. The initial mass is Z N,x = Nx /N, and Z N,x t follows a Markovian jump dynamics : from its current state k/n, { (k + 1)/N at rate 2kN + kθ Z N,x jumps to (k 1)/N at rate 2kN + k(k 1)γ/N. For γ =, this is a GW process in cont. time : each individual independently spawns a child at rate 2N + θ, and dies (childless) at rate 2N. For γ, the quadratic death rate destroys independence. Viewing the individuals alive at time t as being arranged from left to right, and by decreeing that each of the pairwise fights (which happens at rate 2γ) is won by the individual to the left, we arrive at the additional death rate 2γL i (t)/n for individual i, where L i (t) denotes the number of indiv. currently living to the left of i at time t. Etienne Pardoux (Marseille) MANEGE, 18/1/1 4 / 16
9 Discrete approximation Our proof of the above extended Ray Knight theorem is based on an approximation by finite population. For N N, let Z N,x be the total mass of a population of individuals, each of which has mass 1/N. The initial mass is Z N,x = Nx /N, and Z N,x t follows a Markovian jump dynamics : from its current state k/n, { (k + 1)/N at rate 2kN + kθ Z N,x jumps to (k 1)/N at rate 2kN + k(k 1)γ/N. For γ =, this is a GW process in cont. time : each individual independently spawns a child at rate 2N + θ, and dies (childless) at rate 2N. For γ, the quadratic death rate destroys independence. Viewing the individuals alive at time t as being arranged from left to right, and by decreeing that each of the pairwise fights (which happens at rate 2γ) is won by the individual to the left, we arrive at the additional death rate 2γL i (t)/n for individual i, where L i (t) denotes the number of indiv. currently living to the left of i at time t. Etienne Pardoux (Marseille) MANEGE, 18/1/1 4 / 16
10 The just described reproduction dynamics gives rise to a forest F N of trees of descent. At any branch point, we imagine the new branch being placed to the right of the mother branch. Because of the asymmetric killing, the trees further to the right have a tendency to stay smaller : they are under attack by the trees to their left. For a given realization of F N, we read off a continuous and piecewise linear R + -valued path H N (called the exploration path associated with F N ) in the following way : Starting (, ), H N goes upwards at speed 2N until the top of the first mother branch (which is the leftmost leaf of the tree) is hit. There H N turns and goes downwards, again at speed 2N, until arriving at the closest branch point (which is the last birth time of a child of the first ancestor before his death). From there one goes upwards into the (yet unexplored) next branch, and proceeds in a similar fashion until being back at height, which means that the exploration of the leftmost tree is completed. Then explore the next tree, etc., until the exploration of the forest F N is completed. Etienne Pardoux (Marseille) MANEGE, 18/1/1 5 / 16
11 The just described reproduction dynamics gives rise to a forest F N of trees of descent. At any branch point, we imagine the new branch being placed to the right of the mother branch. Because of the asymmetric killing, the trees further to the right have a tendency to stay smaller : they are under attack by the trees to their left. For a given realization of F N, we read off a continuous and piecewise linear R + -valued path H N (called the exploration path associated with F N ) in the following way : Starting (, ), H N goes upwards at speed 2N until the top of the first mother branch (which is the leftmost leaf of the tree) is hit. There H N turns and goes downwards, again at speed 2N, until arriving at the closest branch point (which is the last birth time of a child of the first ancestor before his death). From there one goes upwards into the (yet unexplored) next branch, and proceeds in a similar fashion until being back at height, which means that the exploration of the leftmost tree is completed. Then explore the next tree, etc., until the exploration of the forest F N is completed. Etienne Pardoux (Marseille) MANEGE, 18/1/1 5 / 16
12 The just described reproduction dynamics gives rise to a forest F N of trees of descent. At any branch point, we imagine the new branch being placed to the right of the mother branch. Because of the asymmetric killing, the trees further to the right have a tendency to stay smaller : they are under attack by the trees to their left. For a given realization of F N, we read off a continuous and piecewise linear R + -valued path H N (called the exploration path associated with F N ) in the following way : Starting (, ), H N goes upwards at speed 2N until the top of the first mother branch (which is the leftmost leaf of the tree) is hit. There H N turns and goes downwards, again at speed 2N, until arriving at the closest branch point (which is the last birth time of a child of the first ancestor before his death). From there one goes upwards into the (yet unexplored) next branch, and proceeds in a similar fashion until being back at height, which means that the exploration of the leftmost tree is completed. Then explore the next tree, etc., until the exploration of the forest F N is completed. Etienne Pardoux (Marseille) MANEGE, 18/1/1 5 / 16
13 M 3 M 2 m 2 M 1 * M 4 m 1 * * m 3 Etienne Pardoux (Marseille) MANEGE, 18/1/1 6 / 16
14 The exploration process H N in case θ = γ = We have the following SDE ((H N, V N ) takes values in IR + { 1, 1}) dh N s ds = 2NVs N, H N =, dvs N = 21 {V N s = 1} dpn s 21 {V N r =1} dpn r + 2NdL N s (), V N = 1, where {P N s, s } is a Poisson processes with intensity 4N 2, i. e. M N s = P N s 4N 2 s is a martingale, and L N s (t) := the local time accumulated by H N at level t up to time s 1 := lim ε ε s 1 {t H N u <t+ε}du. Etienne Pardoux (Marseille) MANEGE, 18/1/1 7 / 16
15 The exploration process H N in case θ = γ = We have the following SDE ((H N, V N ) takes values in IR + { 1, 1}) dh N s ds = 2NVs N, H N =, dvs N = 21 {V N s = 1} dpn s 21 {V N r =1} dpn r + 2NdL N s (), V N = 1, where {P N s, s } is a Poisson processes with intensity 4N 2, i. e. M N s = P N s 4N 2 s is a martingale, and L N s (t) := the local time accumulated by H N at level t up to time s 1 := lim ε ε s 1 {t H N u <t+ε}du. Etienne Pardoux (Marseille) MANEGE, 18/1/1 7 / 16
16 If we let S N x It is easily seen that := inf{s > : L N s () [Nx]/N}, Zt N,x = L N S (t), t. x N Etienne Pardoux (Marseille) MANEGE, 18/1/1 8 / 16
17 If we let S N x It is easily seen that := inf{s > : L N s () [Nx]/N}, Zt N,x = L N S (t), t. x N Etienne Pardoux (Marseille) MANEGE, 18/1/1 8 / 16
18 Taking the limit As N, V N oscillates faster and faster. In order to study the limit of H N, we consider for f C 2 (IR), the perturbed test function, see e. g. Ethier, Kurtz (1986) f N (h, v) = f (h) + v 4N f (h). Implementing this with f (h) = h, we get where M 1,N s = 1 2N s H N s + V s N 4N = M1,N s Ms 2,N LN s (), 1 {V N r = 1} dmn r and Ms 2,N = 1 s 1 2N {V N r =1} dmn r are two mutually orthogonal martingales. Etienne Pardoux (Marseille) MANEGE, 18/1/1 9 / 16
19 Taking the limit As N, V N oscillates faster and faster. In order to study the limit of H N, we consider for f C 2 (IR), the perturbed test function, see e. g. Ethier, Kurtz (1986) f N (h, v) = f (h) + v 4N f (h). Implementing this with f (h) = h, we get where M 1,N s = 1 2N s H N s + V s N 4N = M1,N s Ms 2,N LN s (), 1 {V N r = 1} dmn r and Ms 2,N = 1 s 1 2N {V N r =1} dmn r are two mutually orthogonal martingales. Etienne Pardoux (Marseille) MANEGE, 18/1/1 9 / 16
20 Taking the limit as N We obtain the following joint convergence Theorem For any x >, as n, ( {H N s, Ms 1,N, Ms 2,N, s }, {L N s (t), s, t }, Sx N ) ( {Hs, Bs 1, Bs 2 ), s }, {L s (t), s, t }, S x, where B 1 and B 2 are two mutually independent B. M.s, if B = ( 2) 1 (B 1 B 2 ), H is B reflected above, L its local time, and S x = inf{s > ; L s () > x}. Etienne Pardoux (Marseille) MANEGE, 18/1/1 1 / 16
21 A Girsanov transformation 1 Let s Xs N,1 := X N,2 s := s X N := X N,1 + X N,2, θ 2N 1 {V N r = 1}dMN r, γl N r (Hr N ) 1 N {V N r =1} dmn r, and Y N := E(X N ) denote the Doléans exponential of X N, i. e. Y N solves the SDE s Ys N = 1 + Yr dx N r N, s. We prove Proposition Y N is a martingale ( EY N s = 1, s > ). Etienne Pardoux (Marseille) MANEGE, 18/1/1 11 / 16
22 A Girsanov transformation 1 Let s Xs N,1 := X N,2 s := s X N := X N,1 + X N,2, θ 2N 1 {V N r = 1}dMN r, γl N r (Hr N ) 1 N {V N r =1} dmn r, and Y N := E(X N ) denote the Doléans exponential of X N, i. e. Y N solves the SDE s Ys N = 1 + Yr dx N r N, s. We prove Proposition Y N is a martingale ( EY N s = 1, s > ). Etienne Pardoux (Marseille) MANEGE, 18/1/1 11 / 16
23 A Girsanov transformation 1 Let s Xs N,1 := X N,2 s := s X N := X N,1 + X N,2, θ 2N 1 {V N r = 1}dMN r, γl N r (Hr N ) 1 N {V N r =1} dmn r, and Y N := E(X N ) denote the Doléans exponential of X N, i. e. Y N solves the SDE s Ys N = 1 + Yr dx N r N, s. We prove Proposition Y N is a martingale ( EY N s = 1, s > ). Etienne Pardoux (Marseille) MANEGE, 18/1/1 11 / 16
24 A Girsanov transformation 2 Define new probability measures P N s. t. for all s >, d P N Fs dp Fs = Y N s, s Consider the 2-variate point process ( r (Qr 1,N, Qr 2,N ) = 1 {V N u = 1} dpn u, Under P N, r 1 {V N u =1} dpn u Q 1,N r has intensity (4N 2 + 2θN)1 {V N r = 1} dr Q 2,N r has intensity 4[N 2 + γnl N r (H N r )]1 {V N r =1} dr. ), r, This means that under P N, H N has the requested approximate Log Feller dynamics. Etienne Pardoux (Marseille) MANEGE, 18/1/1 12 / 16
25 A Girsanov transformation 2 Define new probability measures P N s. t. for all s >, d P N Fs dp Fs = Y N s, s Consider the 2-variate point process ( r (Qr 1,N, Qr 2,N ) = 1 {V N u = 1} dpn u, Under P N, r 1 {V N u =1} dpn u Q 1,N r has intensity (4N 2 + 2θN)1 {V N r = 1} dr Q 2,N r has intensity 4[N 2 + γnl N r (H N r )]1 {V N r =1} dr. ), r, This means that under P N, H N has the requested approximate Log Feller dynamics. Etienne Pardoux (Marseille) MANEGE, 18/1/1 12 / 16
26 A Girsanov transformation 2 Define new probability measures P N s. t. for all s >, d P N Fs dp Fs = Y N s, s Consider the 2-variate point process ( r (Qr 1,N, Qr 2,N ) = 1 {V N u = 1} dpn u, Under P N, r 1 {V N u =1} dpn u Q 1,N r has intensity (4N 2 + 2θN)1 {V N r = 1} dr Q 2,N r has intensity 4[N 2 + γnl N r (H N r )]1 {V N r =1} dr. ), r, This means that under P N, H N has the requested approximate Log Feller dynamics. Etienne Pardoux (Marseille) MANEGE, 18/1/1 12 / 16
27 A Girsanov transformation 2 We have that as N, X N X = θ B 1 + 2γ L r (H r )dbr 2 2 ( s [ ] ) θ Y N 2 Y = exp X s 4 + γ2 L 2 r (H r ) dr. Define the probability P s. t. for all s >, Under P, d P Fs dp Fs = Y s, s. H s = θ s 2 s γ L r (H r )dr + B r L s(), s. Etienne Pardoux (Marseille) MANEGE, 18/1/1 13 / 16
28 A Girsanov transformation 2 We have that as N, X N X = θ B 1 + 2γ L r (H r )dbr 2 2 ( s [ ] ) θ Y N 2 Y = exp X s 4 + γ2 L 2 r (H r ) dr. Define the probability P s. t. for all s >, Under P, d P Fs dp Fs = Y s, s. H s = θ s 2 s γ L r (H r )dr + B r L s(), s. Etienne Pardoux (Marseille) MANEGE, 18/1/1 13 / 16
29 A Girsanov transformation 2 We have that as N, X N X = θ B 1 + 2γ L r (H r )dbr 2 2 ( s [ ] ) θ Y N 2 Y = exp X s 4 + γ2 L 2 r (H r ) dr. Define the probability P s. t. for all s >, Under P, d P Fs dp Fs = Y s, s. H s = θ s 2 s γ L r (H r )dr + B r L s(), s. Etienne Pardoux (Marseille) MANEGE, 18/1/1 13 / 16
30 Conclusion Our main theorem follows taking the limit as N in the identity Z N,x t = L N S N x (t), Lemma under the measures P N and P, thanks to the following elementary Let (ξ N, η N ), (ξ, η) be random pairs defined on a probability space (Ω, F, P), with η N, η nonnegative scalar random variables, and ξ N, ξ taking values in some complete separable metric space X. Assume that E[η N ] = E[η] = 1. Write ( ξ N, η N ) for the random pair (ξ N, η N ) defined under the probability measure P N which has density η N with respect to P, and ( η, ξ) for the random pair (η, ξ) defined under the probability measure P which has density η with respect to P. Then ( ξ N, η N ) converges in distribution to ( η, ξ), provided that (ξ N, η N ) converges in distribution to (ξ, η). Etienne Pardoux (Marseille) MANEGE, 18/1/1 14 / 16
31 Conclusion Our main theorem follows taking the limit as N in the identity Z N,x t = L N S N x (t), Lemma under the measures P N and P, thanks to the following elementary Let (ξ N, η N ), (ξ, η) be random pairs defined on a probability space (Ω, F, P), with η N, η nonnegative scalar random variables, and ξ N, ξ taking values in some complete separable metric space X. Assume that E[η N ] = E[η] = 1. Write ( ξ N, η N ) for the random pair (ξ N, η N ) defined under the probability measure P N which has density η N with respect to P, and ( η, ξ) for the random pair (η, ξ) defined under the probability measure P which has density η with respect to P. Then ( ξ N, η N ) converges in distribution to ( η, ξ), provided that (ξ N, η N ) converges in distribution to (ξ, η). Etienne Pardoux (Marseille) MANEGE, 18/1/1 14 / 16
32 Appendix : Girsanov s theorem for Poisson processes 1 Let {(Q (1) s,..., Q s (d) ), s } be a d-variate point process adapted to, s } be the predictable some filtration F, and let {λ (i) s (P, F) intensity of Q (i), 1 i d. In other words, M (i) r := Q (i) s s λ(i) r dr is a (P, F) martingale, 1 i d. Assume that none of the Q (i), Q (j), i j, jump simultaneously (so that the M (i) s are mutually orthogonal). Let {µ r (i), r }, 1 i d, be nonnegative F-predictable processes such that for all s and all 1 i d s µ (i) r λ (i) r dr < P -a.s. Etienne Pardoux (Marseille) MANEGE, 18/1/1 15 / 16
33 Appendix : Girsanov s theorem for Poisson processes 1 Let {(Q (1) s,..., Q s (d) ), s } be a d-variate point process adapted to, s } be the predictable some filtration F, and let {λ (i) s (P, F) intensity of Q (i), 1 i d. In other words, M (i) r := Q (i) s s λ(i) r dr is a (P, F) martingale, 1 i d. Assume that none of the Q (i), Q (j), i j, jump simultaneously (so that the M (i) s are mutually orthogonal). Let {µ r (i), r }, 1 i d, be nonnegative F-predictable processes such that for all s and all 1 i d s µ (i) r λ (i) r dr < P -a.s. Etienne Pardoux (Marseille) MANEGE, 18/1/1 15 / 16
34 Appendix : Girsanov s theorem for Poisson processes 1 Let {(Q (1) s,..., Q s (d) ), s } be a d-variate point process adapted to, s } be the predictable some filtration F, and let {λ (i) s (P, F) intensity of Q (i), 1 i d. In other words, M (i) r := Q (i) s s λ(i) r dr is a (P, F) martingale, 1 i d. Assume that none of the Q (i), Q (j), i j, jump simultaneously (so that the M (i) s are mutually orthogonal). Let {µ r (i), r }, 1 i d, be nonnegative F-predictable processes such that for all s and all 1 i d s µ (i) r λ (i) r dr < P -a.s. Etienne Pardoux (Marseille) MANEGE, 18/1/1 15 / 16
35 Appendix : Girsanov s theorem for Poisson processes 2 Theorem {T i k, k = 1, 2...} denoting the jump times of Q(i), let for s ( Y s (i) ) := µ (i) exp T i k 1:Tk i s k { s (1 µ (i) r )λ (i) r } dr and Y s = d j=1 Y (j) s. If E[Y s ] = 1, s, then, for each 1 i d, the process Q (i) has the (i) ( P, F)-intensity λ r = µ (i) r λ (i) r, r, where the probability measure P is defined by d P Fs = Y s, s. dp Fs Note that Y (i) = E(X (i) ), with X (i) s := s (µ (i) r 1)dM r (i), 1 i d, s. Etienne Pardoux (Marseille) MANEGE, 18/1/1 16 / 16
Branching 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 informationDynamics of the evolving Bolthausen-Sznitman coalescent. by Jason Schweinsberg University of California at San Diego.
Dynamics of the evolving Bolthausen-Sznitman coalescent by Jason Schweinsberg University of California at San Diego Outline of Talk 1. The Moran model and Kingman s coalescent 2. The evolving Kingman 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 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 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 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 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 informationErdős-Renyi random graphs basics
Erdős-Renyi random graphs basics Nathanaël Berestycki U.B.C. - class on percolation We take n vertices and a number p = p(n) with < p < 1. Let G(n, p(n)) be the graph such that there is an edge between
More informationWXML Final Report: Chinese Restaurant Process
WXML Final Report: Chinese Restaurant Process Dr. Noah Forman, Gerandy Brito, Alex Forney, Yiruey Chou, Chengning Li Spring 2017 1 Introduction The Chinese Restaurant Process (CRP) generates random partitions
More informationApplications of Ito s Formula
CHAPTER 4 Applications of Ito s Formula In this chapter, we discuss several basic theorems in stochastic analysis. Their proofs are good examples of applications of Itô s formula. 1. Lévy s martingale
More 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 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 informationExact Simulation of Diffusions and Jump Diffusions
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 Content 1 Exact Algorithm Construction
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 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 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 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 informationLook-down model and Λ-Wright-Fisher SDE
Look-down model and Λ-Wright-Fisher SDE B. Bah É. Pardoux CIRM April 17, 2012 B. Bah, É. Pardoux ( CIRM ) Marseille, 1 Febrier 2012 April 17, 2012 1 / 15 Introduction We consider a new version of the look-down
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 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 informationThe multidimensional Ito Integral and the multidimensional Ito Formula. Eric Mu ller June 1, 2015 Seminar on Stochastic Geometry and its applications
The multidimensional Ito Integral and the multidimensional Ito Formula Eric Mu ller June 1, 215 Seminar on Stochastic Geometry and its applications page 2 Seminar on Stochastic Geometry and its applications
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 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 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 informationGENERALIZED RAY KNIGHT THEORY AND LIMIT THEOREMS FOR SELF-INTERACTING RANDOM WALKS ON Z 1. Hungarian Academy of Sciences
The Annals of Probability 1996, Vol. 24, No. 3, 1324 1367 GENERALIZED RAY KNIGHT THEORY AND LIMIT THEOREMS FOR SELF-INTERACTING RANDOM WALKS ON Z 1 By Bálint Tóth Hungarian Academy of Sciences We consider
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 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 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 informationAsymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process
Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process Mátyás Barczy University of Debrecen, Hungary The 3 rd Workshop on Branching Processes and Related
More informationPoint Process Control
Point Process Control The following note is based on Chapters I, II and VII in Brémaud s book Point Processes and Queues (1981). 1 Basic Definitions Consider some probability space (Ω, F, P). A real-valued
More informationAppendix B for The Evolution of Strategic Sophistication (Intended for Online Publication)
Appendix B for The Evolution of Strategic Sophistication (Intended for Online Publication) Nikolaus Robalino and Arthur Robson Appendix B: Proof of Theorem 2 This appendix contains the proof of Theorem
More informationThe Contour Process of Crump-Mode-Jagers Branching Processes
The Contour Process of Crump-Mode-Jagers Branching Processes Emmanuel Schertzer (LPMA Paris 6), with Florian Simatos (ISAE Toulouse) June 24, 2015 Crump-Mode-Jagers trees Crump Mode Jagers (CMJ) branching
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 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 informationMultivariate Generalized Ornstein-Uhlenbeck Processes
Multivariate Generalized Ornstein-Uhlenbeck Processes Anita Behme TU München Alexander Lindner TU Braunschweig 7th International Conference on Lévy Processes: Theory and Applications Wroclaw, July 15 19,
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 informationStochastic Processes III/ Wahrscheinlichkeitstheorie IV. Lecture Notes
BMS Advanced Course Stochastic Processes III/ Wahrscheinlichkeitstheorie IV Michael Scheutzow Lecture Notes Technische Universität Berlin Wintersemester 218/19 preliminary version November 28th 218 Contents
More informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
More 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 informationA. Bovier () Branching Brownian motion: extremal process and ergodic theorems
Branching Brownian motion: extremal process and ergodic theorems Anton Bovier with Louis-Pierre Arguin and Nicola Kistler RCS&SM, Venezia, 06.05.2013 Plan 1 BBM 2 Maximum of BBM 3 The Lalley-Sellke conjecture
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 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 Pedestrian s Guide to Local Time
The Pedestrian s Guide to Local Time Tomas Björk, Department of Finance, Stockholm School of Economics, Box 651, SE-113 83 Stockholm, SWEDEN tomas.bjork@hhs.se November 19, 213 Preliminary version Comments
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 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 informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
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 informationUPPER DEVIATIONS FOR SPLIT TIMES OF BRANCHING PROCESSES
Applied Probability Trust 7 May 22 UPPER DEVIATIONS FOR SPLIT TIMES OF BRANCHING PROCESSES HAMED AMINI, AND MARC LELARGE, ENS-INRIA Abstract Upper deviation results are obtained for the split time of a
More informationChapter 2 Event-Triggered Sampling
Chapter Event-Triggered Sampling In this chapter, some general ideas and basic results on event-triggered sampling are introduced. The process considered is described by a first-order stochastic differential
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 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 informatione - c o m p a n i o n
OPERATIONS RESEARCH http://dx.doi.org/1.1287/opre.111.13ec e - c o m p a n i o n ONLY AVAILABLE IN ELECTRONIC FORM 212 INFORMS Electronic Companion A Diffusion Regime with Nondegenerate Slowdown by Rami
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 informationStochastic Calculus. Kevin Sinclair. August 2, 2016
Stochastic Calculus Kevin Sinclair August, 16 1 Background Suppose we have a Brownian motion W. This is a process, and the value of W at a particular time T (which we write W T ) is a normally distributed
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 informationLarge Deviations for Perturbed Reflected Diffusion Processes
Large Deviations for Perturbed Reflected Diffusion Processes Lijun Bo & Tusheng Zhang First version: 31 January 28 Research Report No. 4, 28, Probability and Statistics Group School of Mathematics, The
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 informationanalysis for transport equations and applications
Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg
More informationCIMPA SCHOOL, 2007 Jump Processes and Applications to Finance Monique Jeanblanc
CIMPA SCHOOL, 27 Jump Processes and Applications to Finance Monique Jeanblanc 1 Jump Processes I. Poisson Processes II. Lévy Processes III. Jump-Diffusion Processes IV. Point Processes 2 I. Poisson Processes
More informationThe tree-valued Fleming-Viot process with mutation and selection
The tree-valued Fleming-Viot process with mutation and selection Peter Pfaffelhuber University of Freiburg Joint work with Andrej Depperschmidt and Andreas Greven Population genetic models Populations
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 informationThe Cameron-Martin-Girsanov (CMG) Theorem
The Cameron-Martin-Girsanov (CMG) Theorem There are many versions of the CMG Theorem. In some sense, there are many CMG Theorems. The first version appeared in ] in 944. Here we present a standard version,
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 STATE ESTIMATION FROM DOUBLY STOCHASTIC POINT PROCESS OBSERVATION. 0. Introduction
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLVI, Number 1, March 1 A NOTE ON STATE ESTIMATION FROM DOUBLY STOCHASTIC POINT PROCESS OBSERVATION. Introduction In this note we study a state estimation
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 informationRandom measures, intersections, and applications
Random measures, intersections, and applications Ville Suomala joint work with Pablo Shmerkin University of Oulu, Finland Workshop on fractals, The Hebrew University of Jerusalem June 12th 2014 Motivation
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 informationBranching Brownian motion seen from the tip
Branching Brownian motion seen from the tip J. Berestycki 1 1 Laboratoire de Probabilité et Modèles Aléatoires, UPMC, Paris 09/02/2011 Joint work with Elie Aidekon, Eric Brunet and Zhan Shi J Berestycki
More informationAN EXTREME-VALUE ANALYSIS OF THE LIL FOR BROWNIAN MOTION 1. INTRODUCTION
AN EXTREME-VALUE ANALYSIS OF THE LIL FOR BROWNIAN MOTION DAVAR KHOSHNEVISAN, DAVID A. LEVIN, AND ZHAN SHI ABSTRACT. We present an extreme-value analysis of the classical law of the iterated logarithm LIL
More 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 informationBrownian Motion on Manifold
Brownian Motion on Manifold QI FENG Purdue University feng71@purdue.edu August 31, 2014 QI FENG (Purdue University) Brownian Motion on Manifold August 31, 2014 1 / 26 Overview 1 Extrinsic construction
More informationRichard F. Bass Krzysztof Burdzy University of Washington
ON DOMAIN MONOTONICITY OF THE NEUMANN HEAT KERNEL Richard F. Bass Krzysztof Burdzy University of Washington Abstract. Some examples are given of convex domains for which domain monotonicity of the Neumann
More informationELEMENTS OF STOCHASTIC CALCULUS VIA REGULARISATION. A la mémoire de Paul-André Meyer
ELEMENTS OF STOCHASTIC CALCULUS VIA REGULARISATION A la mémoire de Paul-André Meyer Francesco Russo (1 and Pierre Vallois (2 (1 Université Paris 13 Institut Galilée, Mathématiques 99 avenue J.B. Clément
More informationBrownian Motion and lorentzian manifolds
The case of Jürgen Angst Institut de Recherche Mathématique Avancée Université Louis Pasteur, Strasbourg École d été de Probabilités de Saint-Flour juillet 2008, Saint-Flour 1 Construction of the diffusion,
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 informationSolutions to the Exercises in Stochastic Analysis
Solutions to the Exercises in Stochastic Analysis Lecturer: Xue-Mei Li 1 Problem Sheet 1 In these solution I avoid using conditional expectations. But do try to give alternative proofs once we learnt conditional
More informationSpatial diffusions with singular drifts: The construction of Super Brownian motion with immigration at unoccupied sites
Spatial diffusions with singular drifts: The construction of Super Brownian motion with immigration at unoccupied sites by Saifuddin Syed B.Math Hons., University of Waterloo, 24 AN ESSAY SUBMITTED IN
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 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 informationFunctions of Several Variables: Chain Rules
Functions of Several Variables: Chain Rules Calculus III Josh Engwer TTU 29 September 2014 Josh Engwer (TTU) Functions of Several Variables: Chain Rules 29 September 2014 1 / 30 PART I PART I: MULTIVARIABLE
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 informationOther properties of M M 1
Other properties of M M 1 Přemysl Bejda premyslbejda@gmail.com 2012 Contents 1 Reflected Lévy Process 2 Time dependent properties of M M 1 3 Waiting times and queue disciplines in M M 1 Contents 1 Reflected
More informationAsymptotic Properties of an Approximate Maximum Likelihood Estimator for Stochastic PDEs
Asymptotic Properties of an Approximate Maximum Likelihood Estimator for Stochastic PDEs M. Huebner S. Lototsky B.L. Rozovskii In: Yu. M. Kabanov, B. L. Rozovskii, and A. N. Shiryaev editors, Statistics
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 informationBrownian survival and Lifshitz tail in perturbed lattice disorder
Brownian survival and Lifshitz tail in perturbed lattice disorder Ryoki Fukushima Kyoto niversity Random Processes and Systems February 16, 2009 6 B T 1. Model ) ({B t t 0, P x : standard Brownian motion
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 informationSolving Stochastic Partial Differential Equations as Stochastic Differential Equations in Infinite Dimensions - a Review
Solving Stochastic Partial Differential Equations as Stochastic Differential Equations in Infinite Dimensions - a Review L. Gawarecki Kettering University NSF/CBMS Conference Analysis of Stochastic Partial
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 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 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 informationStochastic Demography, Coalescents, and Effective Population Size
Demography Stochastic Demography, Coalescents, and Effective Population Size Steve Krone University of Idaho Department of Mathematics & IBEST Demographic effects bottlenecks, expansion, fluctuating population
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 informationThe Moran Process as a Markov Chain on Leaf-labeled Trees
The Moran Process as a Markov Chain on Leaf-labeled Trees David J. Aldous University of California Department of Statistics 367 Evans Hall # 3860 Berkeley CA 94720-3860 aldous@stat.berkeley.edu http://www.stat.berkeley.edu/users/aldous
More informationA 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 informationTHE COALESCENT Lectures given at CIMPA school St Louis, Sénégal, April Étienne Pardoux
THE COALESCENT Lectures given at CIMPA school St Louis, Sénégal, April 2010 Étienne Pardoux 2 Contents 1 Kingman s coalescent 5 1.1 The Wright Fisher model.................... 5 1.2 Cannings model.........................
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 informationAn essay on the general theory of stochastic processes
Probability Surveys Vol. 3 (26) 345 412 ISSN: 1549-5787 DOI: 1.1214/1549578614 An essay on the general theory of stochastic processes Ashkan Nikeghbali ETHZ Departement Mathematik, Rämistrasse 11, HG G16
More informationLet (Ω, F) be a measureable space. A filtration in discrete time is a sequence of. F s F t
2.2 Filtrations Let (Ω, F) be a measureable space. A filtration in discrete time is a sequence of σ algebras {F t } such that F t F and F t F t+1 for all t = 0, 1,.... In continuous time, the second condition
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 informationExistence and Comparisons for BSDEs in general spaces
Existence and Comparisons for BSDEs in general spaces Samuel N. Cohen and Robert J. Elliott University of Adelaide and University of Calgary BFS 2010 S.N. Cohen, R.J. Elliott (Adelaide, Calgary) BSDEs
More information