Two viewpoints on measure valued processes
|
|
- Milton Watts
- 5 years ago
- Views:
Transcription
1 Two viewpoints on measure valued processes Olivier Hénard Université Paris-Est, Cermics
2 Contents 1 The classical framework : from no particle to one particle 2 The lookdown framework : many particles.
3 The basic ingredients Superprocesses form a specific class of measure valued processes, whose definition requires two basic ingredients : a Markov process valued in some space E, with law P x when started at x E. a branching mechanism ψ : ψ(x, λ) = α(x)λ + β(x)λ 2 + (0, ) [e λl 1 + λl]π(x, dl) for λ 0, α B, β B +, and π a kernel from E to (0, ) such that sup (l x R d (0, ) l2 )π(x, dl) <.
4 The inhomogeneous superprocesses Definition ((P, ψ)-superprocesses, Dynkin) There exists a unique M f (E) valued Markov process (Z = Z t, t 0) such that, for all f C + b (E), E µ ( e Z t (f ) ) = e µ(u t ), µ M F (E), t 0, where u t is the unique positive solution to the mild equation ( t ) ( u t (x) + E x ds ψ(x s, u t s (X s )) = E x f (Xt ) ). 0
5 The superprocess started from one individual There exists a sigma-finite measure N x called the canonical measure, such that, if i I δ Z i is a Poisson Point Measure with intensity N x, then Z under P δx equals Z i in distribution. i I Exponential formula for Poisson Point measure writes down : E δx ( e Z t (f ) ) = e N x (1 e Z t (f )) Recalling the definition of the superprocess : u (f ) t (x) = N x (1 e Z t (f ) ). Let H max = inf {t 0, Z t = 0} [0, ] denote the height of the process. Taking f = l1 allows us to see that : v t (x) := N x (H max > t) equals lim u (l1) t (x). l
6 We will assume almost sure extinction, meaning : P µ (H max < ) = 1 for every µ M f (E), where H max = inf {t 0, Z t = 0} is the height of the process. Questions Desintegration of the law of superprocess with respect to H max : that means, find a family of probability measure (N (h) x, 0 < h < ) such that : N x (A) = N x (H max dh)n (h) x (A). h (0, ) Weak convergence of this family (N (h) x, 0 < h < ).
7 Spinal decomposition under N (h) x (i) Sample (X s, 0 s < h) with law P (h) x : 0 t h, dp (h) x D t = hv h t (X t ) e t 0 ds λψ(x s,v h s (X s )). dp x Dt h v h (x) (ii) Conditionally on (X s, 0 s < h), define µ (h) s (dz ) = 2β(X s )N Xs (dz ) + (0, ) e l v h s(x s ) l π(x s, dl) P lδxs (dz ) (iii) Conditionally on (µ (h) s, 0 s < h), let i I δ (si,z i )(ds, dz ) denote a Poisson point measure with intensity 1 {0 s<h,hmax (Z )+s h} ds µ(h) s (dz ). Theorem (Williams decomposition) ( i I Z i (t s i ) +, 0 t < h) has law N (h) x.
8 Weak convergence of this spinal decomposition (i) Sample (X s, s 0) with law P ( ) x : t 0, dp ( ) x D t = φ 0(X t ) t dp x Dt φ 0 (x) e 0 ds (α(x s) λ 0 ) with λ 0 = sup {λ R, f, (L α)f = λf } the generalized eigenvalue and φ 0 the associated generalized eigenvector. (ii) Conditionally on (X s, s 0), define µ ( ) s (dz ) = 2β(X s )N Xs (dz ) + l π(x s, dl) P lδxs (dz ) (0, ) (iii) Conditionally on (µ ( ) s, s 0), let i I δ (si,z i )(ds, dz ) denote a Poisson point measure with intensity ds µ ( ) s (dz ). Theorem (N (h) x, 0 < h < ) is weakly converging as h, towards N ( ) x say. ( i I Z i (t s i ) +, t 0) has law N ( ) x.
9 Hint for the proof, h fixed First prove an absolute continuity property : For every G t = σ(z s, s t) measurable functional F and 0 t < h, ( ) N (h) Zt ( x (F (Z )) = N h v h t ) x e (Z t (v h t ) v h (x)) F (Z ). h v h (x) Using the definition of the superprocess and the Feynman Kac formula, we also have that, for every f and g positive measurable function on E, N x ( Zt (f ) e Z t (g) ) = E x ( e t 0 ds λψ(x s,u (g) t s (X s)) f (X t ) ). since lhs and rhs satisfy the same mild equation. Set f ( ) hv h t ( ) h v h (x) ev h(x) and g g + v h t : 1 the lhs is first formula for F (Z ) = e Z t (g). 2 the integrand on the rhs is the product This gives the Williams decomposition. dp (h) x D t e t 0 ds µ(h) s ((1 e Z t s (g) )1 Hmax h s) dp x Dt
10 Hint for the proof, h 4. For passing to the limit as h, we work directly on the Williams decomposition. no problem for the laws of the subtrees. weak convergence of the law of the spine amounts to check that : h v h t ( ) h v h (x) φ(x) φ( ) eλ 0t 1 in L 1 (P x ) as h. Interpretation We get a recursive decomposition of the superprocess as a trunk on which copies of the same superprocess are grafted.
11 From one particle to many So far, we have unearthed one particle, namely the most persistent. What about the second most persistent, and the third, and the fourth...? Idea : conditionally on the total mass process, in the homogeneous case, these particles form an exchangeable particle particle system, which, ponderated by the total mass process, allows to recover the whole superprocess.
12 Contents 1 The classical framework : from no particle to one particle 2 The lookdown framework : many particles.
13 Conversely, let us assume that we are given : some process Y (in the rôle of Z (1)). some increasing process U (in the rôle of the quadratic variation [Z (1)] of Z (1)). Question Question : How to define the genealogy of a population process with size Y and resampling events directed by U?
14 U being increasing : U t = Ut k + s t U s where U k is continuous and U s = U s U s. We define N ρ = δ (t,π) and, for each t > 0, the exchangeable partition Π of t 0, U t =0 N has a unique non trivial block with asymptotic frequency U(t)/Y t. N k is a Poisson point measure with intensity (du k (t)/y 2 (t)) µ k, and the Kingman measure µ k assigns mass one to partitions with a unique non trivial block consisting of two different integers. We let : N = N ρ + N k
15 Let R 0 be some random probability measure on E independent of N. Conditionally on (R 0, N), we can define a particle system X = (X t (n), t 0, n N) : the initial state (X 0 (n), n N) is an exchangeable sequence valued in E with de Finetti measure R 0. between the reproduction events, X t (n) mutates according to a Markov process in E with law P x when started at x independently for each n. each atom (t, Π) of N k + N ρ is a reproduction event : let j 1 < j 2 <... be the elements of the unique block of the partition Π which is not a singleton. Then put { Xt (j l ) = X t (j 1 ) for each integer l X t (n) = X t (n #{l > 1, j l n}) if n / {j l, l N}.
16 Theorem (Donnelly and Kurtz, 99) The sequence (X t (n), n N) is then exchangeable. We denote by R its de Finetti measure : and define 1 R t := lim n n n δ Xt (i). i=1 (Z t, t 0) := (Y t R t, t 0). Notice the following consequence of the de Finetti Theorem : Conditionally on R t, (X t (n), n N) is i.i.d. according to R t.
17 Two useful examples : taking Y to be a CB process and U t = [Y ] t yields the (homogeneous) superprocesses. taking Y = 1 and U a subordinator (with jumps not bigger than 1) yields a Generalized Fleming Viot process. We introduce the filtrations : F t = σ((y s, s t), (X s, s t)) knows everything about the genealogy. G t = σ(z s, s t) just knows the resulting type distribution of the population. D t will denote the filtration induced by the spatial motion with law P.
18 Definition of the additive h-transform Assume there is some deterministic m = m(t) such that (Y t /m(t), t 0) and (m(t)h t (X t (1)), t 0) are two martingales in their own filtrations. Then we have that : (Z t (h t ), t 0) is a G t - martingale We may now define the additive h-transform Z H of Z. ( ) A G t, E(Z H Zt (h A) = E t ) Z 0 (h 0 ) 1 A(Z )
19 A particle system for the h-transform Assume m = 1. We define the law of a new process by the following requirements : (Y h, U h, X h ) (i) Initial condition (if deterministic) is unchanged. (iv) (Y h, U h ) is distributed according to : A G t, P((Y h, U h ) A) = E ( ) Yt 1 Y A (Y, U). 0 (v) (Xt h(n), n N, t 0) is defined as (X t (n), n N, t 0), except that (Xt h (1), t 0) now follows : ( ) A D t, P(X h ht (X (1) A) = E t (1)) h 0 (X 0 (1)) 1 A(X (1)). Interpretation Factorization of the effect of the h-transform : first particle responsible for the change of the spatial motion/ total mass process for the change of the mass.
20 Noting that E(X t (1) G t ) = R t, we get by projection on the smaller filtration G t : Theorem We have that : (a) the limit of the empirical measure : R h t 1 := lim n n δ X h t (i) 1 i n exists a.s. (b) the process (Yt hrh t, t 0) is distributed as (Z H t, t 0).
21 The first particle is precursory Corollary Conditionally on {R h t = µ}, Xt h (1) is distributed according to : P(X h t (1) dx) = h t (x) µ(h t ) µ(dx), and (Xt h (n), n 2) is an independent exchangeable random sequence with de Finetti measure µ.
22 Application : the branching case. Question Is there an individual responsible for the additional mass? Assume Y is a is pure jump CB(ψ) with ψ(λ) = αλ + (0, ) [e λl 1 + λl]π(dl) and U = [Y ]. Jumps of Y form a PPM with intensity Y t π(dl) at time t. Y h is a CBI(ψ, ψ ψ (0)). Jumps of Y h form a PPM with intensity Yt h π(dl) + lπ(dl) at time t. Jump l of Y h at time t is associated to a reproduction of the first level particle with probability l/y h (t) = l/(y h (t ) + l). Reproduction rate of the first level particle : l/(yt h + l)(y t h π(dl) + lπ(dl)) = lπ(dl). Reproduction rate of the other particles : (1 l/(yt h + l))(y t h π(dl) + lπ(dl)) = Y h t π(dl). Corollary The additive h-transform of a homogeneous superprocess is just the usual process, plus some immigration along a spine moving as an h-transform.
23 Application : the GFV process The Generalized Fleming Viot process is the process Z (= R) obtained when setting : Y = 1, U is a subordinator with jumps no greater than 1, In that case, we denote by φ the Laplace exponent of the subordinator U : 1 φ(λ) = cλ + (1 e λx )π(dx) 0 where c 0 and 1 0 x π(dx) <. The genealogy of the lookdown particle system is by construction described by the Λ-coalescent of Pitman. The finite measure Λ may be recovered from φ : g(x)λ(dx) = cg(0) + g( x)π(dx). 0 Corollary The additive h-transform of a GFV is just the usual GFV process with the first level particle moving as an h-transform.
24 The GFV process in finite state space without mutation Let us assume furthermore that E = {1, 2,..., K }, and P x = δ x. Set ν(dx) = x 2 Λ(dx). Define the pushing rates. i(i 1) 1 r i = c + ν(dx) (1 (1 x) i ix(1 x) i 1) 2 0 i(i 1) 1 = c + ν(dx)p(bin(i, x) 2) 2 0 We have that, for 1 K K : ( ) K i=1 R t ({i}) E( K i=1 R 0({i})) er K t, t 0 is a non negative G-martingale. Thus we may define, for 1 K K, ( ) A G t, P(R H K i=1 A) := E 1 A (R) R t ({i}) E( K i=1 R 0({i})) er K t.
25 We define from X a new particle system X h as follows : (i) The random sequence X0 h is such that the finite sequence ) (X0 h(j), 1 j K is a uniform permutation of {1,..., K }, and, ) independently, the sequence (X0 h(j), j K + 1 is exchangeable with asymptotic frequencies R0 h given by : P(R0 h dx) := ( K i=1 R 0({i}) ) P(R 0 dx). E K i=1 R 0({i}) (ii) The reproduction events are driven { by the restriction of the Poisson Point } measure N = N k + N ρ to V := (s, Π), Π [K ] = {{1}, {2},..., {K }}, that is the atoms for which the reproductions events do not involve more than one of the first K levels. Interpretation Forget some reproduction events/ Introduce immigration.
26 We define the lowest level L(t) at which the first K types have appeared : L(t) = inf{i 1, {1,..., K } {X t (1),..., X t (i)}} The random variable L(t) is F t measurable, but not G t measurable. Lemma The law of the process X h is absolutely continuous with respect to the law of the process X on each F t and : ( ) A F t, P(X h 1 {L(t)=K } A) = E 1 A (X ) P(L(0) = K ) er K t.
27 Noting that E(L(t) = K G t ) = K! 1 i K R t ({i}), we get by projection on the smaller filtration G t : Theorem We have that : (a) the limit of the empirical measure : exists a.s. R h t 1 := lim n n δ X h t (i) 1 i n (b) the process (R h t, t 0) is distributed as (RH t, t 0).
28 Define the family of processes R ( t) by : A G t, P(R ( t) A) := P ( R A K R t ({i}) = 0 ). i=1 The process R ( t) thus corresponds to the process R conditioned on non extinction of each of the first K types before time t. We have the following extension of Kimura s Theorem for GFV processes : Theorem Let t 0 0 be fixed. Assume that lim l P K +1 (L(l)< ) P K (L(l)< ) = 0. The family of processes (R ( t), t t 0 ) weakly converges on G t0 as t towards the law of R h.
29 Let us assume that K = K = 2 : two type population conditioned on non absorbtion. We define, for f C 2 ([0, 1]), and x [0, 1] : and 1 L 0 f (x) := c(1 2x)f (x) + y(1 y)ν(dy)[f (x(1 y) + y) f (x)] 1 + y(1 y)ν(dy)[f (x(1 y)) f (x)], 0 L 1 f (x) := 1 2 cx(1 x)f (x) + x (1 y) 2 ν(dy)[f (x(1 y) + y)) f (x)] 1 + (1 x) (1 y) 2 ν(dy)[f (x(1 y)) f (x)]. 0 Proposition Let L h be the generator of R h {1}. We have that : L h = L 0 + L 1.
30 Assume K = 1, and K = 2 : two type population conditioned on absorbtion by 1. We define, for f C 2 ([0, 1]), and x [0, 1] : and 1 L 0 f (x) := c(1 x)f (x) + yν(dy)[f (x(1 y) + y) f (x)] 0 L 1 f (x) := 1 2 cx(1 x)f (x) + x 1 0 (1 y)ν(dy)[f (x(1 y) + y)) f (x)] 1 + (1 x) (1 y)ν(dy)[f (x(1 y)) f (x)]. 0 Proposition Let L h be the generator of R h {1}. We have that : L h = L 0 + L 1.
(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(This is a sample cover image for this issue. The actual cover is not yet available at this time.) This article appeared in a journal published by Elsevier. The attached copy is furnished to the author
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 informationContinuous-state branching processes, extremal processes and super-individuals
Continuous-state branching processes, extremal processes and super-individuals Clément Foucart Université Paris 13 with Chunhua Ma Nankai University Workshop Berlin-Paris Berlin 02/11/2016 Introduction
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 informationStochastic flows associated to coalescent processes
Stochastic flows associated to coalescent processes Jean Bertoin (1) and Jean-François Le Gall (2) (1) Laboratoire de Probabilités et Modèles Aléatoires and Institut universitaire de France, Université
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 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 informationSTOCHASTIC FLOWS ASSOCIATED TO COALESCENT PROCESSES III: LIMIT THEOREMS
Illinois Journal of Mathematics Volume 5, Number 1, Spring 26, Pages 147 181 S 19-282 STOCHASTIC FLOWS ASSOCIATED TO COALESCENT PROCESSES III: LIMIT THEOREMS JEAN BERTOIN AND JEAN-FRANÇOIS LE GALL Abstract.
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 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 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 informationYaglom-type limit theorems for branching Brownian motion with absorption. by Jason Schweinsberg University of California San Diego
Yaglom-type limit theorems for branching Brownian motion with absorption by Jason Schweinsberg University of California San Diego (with Julien Berestycki, Nathanaël Berestycki, Pascal Maillard) Outline
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 informationA Note on the Central Limit Theorem for a Class of Linear Systems 1
A Note on the Central Limit Theorem for a Class of Linear Systems 1 Contents Yukio Nagahata Department of Mathematics, Graduate School of Engineering Science Osaka University, Toyonaka 560-8531, Japan.
More 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 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 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 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 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 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 informationThe Λ-Fleming-Viot process and a connection with Wright-Fisher diffusion. Bob Griffiths University of Oxford
The Λ-Fleming-Viot process and a connection with Wright-Fisher diffusion Bob Griffiths University of Oxford A d-dimensional Λ-Fleming-Viot process {X(t)} t 0 representing frequencies of d types of individuals
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 informationEXCHANGEABLE COALESCENTS. Jean BERTOIN
EXCHANGEABLE COALESCENTS Jean BERTOIN Nachdiplom Lectures ETH Zürich, Autumn 2010 The purpose of this series of lectures is to introduce and develop some of the main aspects of a class of random processes
More informationExercises with solutions (Set D)
Exercises with solutions Set D. A fair die is rolled at the same time as a fair coin is tossed. Let A be the number on the upper surface of the die and let B describe the outcome of the coin toss, where
More informationON COMPOUND POISSON POPULATION MODELS
ON COMPOUND POISSON POPULATION MODELS Martin Möhle, University of Tübingen (joint work with Thierry Huillet, Université de Cergy-Pontoise) Workshop on Probability, Population Genetics and Evolution Centre
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 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 information18.175: Lecture 15 Characteristic functions and central limit theorem
18.175: Lecture 15 Characteristic functions and central limit theorem Scott Sheffield MIT Outline Characteristic functions Outline Characteristic functions Characteristic functions Let X be a random variable.
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 informationMulti-operator scaling random fields
Available online at www.sciencedirect.com Stochastic Processes and their Applications 121 (2011) 2642 2677 www.elsevier.com/locate/spa Multi-operator scaling random fields Hermine Biermé a, Céline Lacaux
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 informationAsymptotics for posterior hazards
Asymptotics for posterior hazards Igor Prünster University of Turin, Collegio Carlo Alberto and ICER Joint work with P. Di Biasi and G. Peccati Workshop on Limit Theorems and Applications Paris, 16th January
More informationStochastic Integration.
Chapter Stochastic Integration..1 Brownian Motion as a Martingale P is the Wiener measure on (Ω, B) where Ω = C, T B is the Borel σ-field on Ω. In addition we denote by B t the σ-field generated by x(s)
More informationarxiv: v3 [math.pr] 30 May 2013
FROM FLOWS OF Λ FLEMING-VIOT PROCESSES TO LOOKDOWN PROCESSES VIA FLOWS OF PARTITIONS Cyril LABBÉ Université Pierre et Marie Curie May 31, 2013 arxiv:1107.3419v3 [math.pr] 30 May 2013 Abstract The goal
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 information9.2 Branching random walk and branching Brownian motions
168 CHAPTER 9. SPATIALLY STRUCTURED MODELS 9.2 Branching random walk and branching Brownian motions Branching random walks and branching diffusions have a long history. A general theory of branching Markov
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 informationMandelbrot s cascade in a Random Environment
Mandelbrot s cascade in a Random Environment A joint work with Chunmao Huang (Ecole Polytechnique) and Xingang Liang (Beijing Business and Technology Univ) Université de Bretagne-Sud (Univ South Brittany)
More informationHeavy Tailed Time Series with Extremal Independence
Heavy Tailed Time Series with Extremal Independence Rafa l Kulik and Philippe Soulier Conference in honour of Prof. Herold Dehling Bochum January 16, 2015 Rafa l Kulik and Philippe Soulier Regular variation
More informationScale functions for spectrally negative Lévy processes and their appearance in economic models
Scale functions for spectrally negative Lévy processes and their appearance in economic models Andreas E. Kyprianou 1 Department of Mathematical Sciences, University of Bath 1 This is a review talk and
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 informationModelling populations under fluctuating selection
Modelling populations under fluctuating selection Alison Etheridge With Aleksander Klimek (Oxford) and Niloy Biswas (Harvard) The simplest imaginable model of inheritance A population of fixed size, N,
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
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 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 informationMath 2C03 - Differential Equations. Slides shown in class - Winter Laplace Transforms. March 4, 5, 9, 11, 12, 16,
Math 2C03 - Differential Equations Slides shown in class - Winter 2015 Laplace Transforms March 4, 5, 9, 11, 12, 16, 18... 2015 Laplace Transform used to solve linear ODEs and systems of linear ODEs with
More informationStrongly Consistent Multivariate Conditional Risk Measures
Strongly Consistent Multivariate Conditional Risk Measures annes offmann Thilo Meyer-Brandis regor Svindland January 11, 2016 Abstract We consider families of strongly consistent multivariate conditional
More informationStatistical test for some multistable processes
Statistical test for some multistable processes Ronan Le Guével Joint work in progress with A. Philippe Journées MAS 2014 1 Multistable processes First definition : Ferguson-Klass-LePage series Properties
More informationA concentration theorem for the equilibrium measure of Markov chains with nonnegative coarse Ricci curvature
A concentration theorem for the equilibrium measure of Markov chains with nonnegative coarse Ricci curvature arxiv:103.897v1 math.pr] 13 Mar 01 Laurent Veysseire Abstract In this article, we prove a concentration
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 informationOptional Stopping Theorem Let X be a martingale and T be a stopping time such
Plan Counting, Renewal, and Point Processes 0. Finish FDR Example 1. The Basic Renewal Process 2. The Poisson Process Revisited 3. Variants and Extensions 4. Point Processes Reading: G&S: 7.1 7.3, 7.10
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 informationRandom Process Lecture 1. Fundamentals of Probability
Random Process Lecture 1. Fundamentals of Probability Husheng Li Min Kao Department of Electrical Engineering and Computer Science University of Tennessee, Knoxville Spring, 2016 1/43 Outline 2/43 1 Syllabus
More informationFormulas for probability theory and linear models SF2941
Formulas for probability theory and linear models SF2941 These pages + Appendix 2 of Gut) are permitted as assistance at the exam. 11 maj 2008 Selected formulae of probability Bivariate probability Transforms
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 informationA modified lookdown construction for the Xi-Fleming-Viot process with mutation and populations with recurrent bottlenecks
Alea 6, 25 61 (2009) A modified lookdown construction for the Xi-Fleming-Viot process with mutation and populations with recurrent bottlenecks Matthias Birkner, Jochen Blath, Martin Möhle, Matthias Steinrücken
More informationMatemáticas II: Segundo del Grado en Ingeniería Aeroespacial
Matemáticas II: Segundo del Grado en Ingeniería Aeroespacial Sergio Blanes, Dolors Roselló Ecuaciones diferenciales y transformadas de Laplace con aplicaciones L.M. Sánchez, M.P. Legua Ref.: 211-798 Capítulo
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 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 informationOn the posterior structure of NRMI
On the posterior structure of NRMI Igor Prünster University of Turin, Collegio Carlo Alberto and ICER Joint work with L.F. James and A. Lijoi Isaac Newton Institute, BNR Programme, 8th August 2007 Outline
More informationLecture 22: Variance and Covariance
EE5110 : Probability Foundations for Electrical Engineers July-November 2015 Lecture 22: Variance and Covariance Lecturer: Dr. Krishna Jagannathan Scribes: R.Ravi Kiran In this lecture we will introduce
More informationHarnack Inequalities and Applications for Stochastic Equations
p. 1/32 Harnack Inequalities and Applications for Stochastic Equations PhD Thesis Defense Shun-Xiang Ouyang Under the Supervision of Prof. Michael Röckner & Prof. Feng-Yu Wang March 6, 29 p. 2/32 Outline
More informationA Dynamic Contagion Process with Applications to Finance & Insurance
A Dynamic Contagion Process with Applications to Finance & Insurance Angelos Dassios Department of Statistics London School of Economics Angelos Dassios, Hongbiao Zhao (LSE) A Dynamic Contagion Process
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 informationSolution of Additional Exercises for Chapter 4
1 1. (1) Try V (x) = 1 (x 1 + x ). Solution of Additional Exercises for Chapter 4 V (x) = x 1 ( x 1 + x ) x = x 1 x + x 1 x In the neighborhood of the origin, the term (x 1 + x ) dominates. Hence, the
More informationADVANCED PROBABILITY: SOLUTIONS TO SHEET 1
ADVANCED PROBABILITY: SOLUTIONS TO SHEET 1 Last compiled: November 6, 213 1. Conditional expectation Exercise 1.1. To start with, note that P(X Y = P( c R : X > c, Y c or X c, Y > c = P( c Q : X > c, Y
More informationRandom and Deterministic perturbations of dynamical systems. Leonid Koralov
Random and Deterministic perturbations of dynamical systems Leonid Koralov - M. Freidlin, L. Koralov Metastability for Nonlinear Random Perturbations of Dynamical Systems, Stochastic Processes and Applications
More informationLogFeller et Ray Knight
LogFeller et Ray Knight Etienne Pardoux joint work with V. Le and A. Wakolbinger Etienne Pardoux (Marseille) MANEGE, 18/1/1 1 / 16 Feller s branching diffusion with logistic growth We consider the diffusion
More 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 informationHigher-order ordinary differential equations
Higher-order ordinary differential equations 1 A linear ODE of general order n has the form a n (x) dn y dx n +a n 1(x) dn 1 y dx n 1 + +a 1(x) dy dx +a 0(x)y = f(x). If f(x) = 0 then the equation is called
More informationOn Ergodic Impulse Control with Constraint
On Ergodic Impulse Control with Constraint Maurice Robin Based on joint papers with J.L. Menaldi University Paris-Sanclay 9119 Saint-Aubin, France (e-mail: maurice.robin@polytechnique.edu) IMA, Minneapolis,
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationThe effects of a weak selection pressure in a spatially structured population
The effects of a weak selection pressure in a spatially structured population A.M. Etheridge, A. Véber and F. Yu CNRS - École Polytechnique Motivations Aim : Model the evolution of the genetic composition
More informationStochastic flows associated to coalescent processes II: Stochastic differential equations
Stochastic flows associated to coalescent processes II: Stochastic differential equations Dedicated to the memory of Paul-André MEYER Jean Bertoin (1) and Jean-François Le Gall (2) (1) Laboratoire de Probabilités
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 informationStochastic analysis of biochemical reaction networks with absolute concentration robustness
Stochastic analysis of biochemical reaction networks with absolute concentration robustness David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison Boston University
More informationWeak convergence and large deviation theory
First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov
More informationThird In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 3 December 2009 (1) [6] Given that 2 is an eigenvalue of the matrix
Third In-Class Exam Solutions Math 26, Professor David Levermore Thursday, December 2009 ) [6] Given that 2 is an eigenvalue of the matrix A 2, 0 find all the eigenvectors of A associated with 2. Solution.
More information4. Conditional risk measures and their robust representation
4. Conditional risk measures and their robust representation We consider a discrete-time information structure given by a filtration (F t ) t=0,...,t on our probability space (Ω, F, P ). The time horizon
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 informationExponential functionals of Lévy processes
Exponential functionals of Lévy processes Víctor Rivero Centro de Investigación en Matemáticas, México. 1/ 28 Outline of the talk Introduction Exponential functionals of spectrally positive Lévy processes
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 informationA representation for the semigroup of a two-level Fleming Viot process in terms of the Kingman nested coalescent
A representation for the semigroup of a two-level Fleming Viot process in terms of the Kingman nested coalescent Airam Blancas June 11, 017 Abstract Simple nested coalescent were introduced in [1] to model
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 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 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 informationMarkov Chain Approximation of Pure Jump Processes. Nikola Sandrić (joint work with Ante Mimica and René L. Schilling) University of Zagreb
Markov Chain Approximation of Pure Jump Processes Nikola Sandrić (joint work with Ante Mimica and René L. Schilling) University of Zagreb 8 th International Conference on Lévy processes Angers, July 25-29,
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 information8 Periodic Linear Di erential Equations - Floquet Theory
8 Periodic Linear Di erential Equations - Floquet Theory The general theory of time varying linear di erential equations _x(t) = A(t)x(t) is still amazingly incomplete. Only for certain classes of functions
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 informationHomework #6 : final examination Due on March 22nd : individual work
Université de ennes Année 28-29 Master 2ème Mathématiques Modèles stochastiques continus ou à sauts Homework #6 : final examination Due on March 22nd : individual work Exercise Warm-up : behaviour of characteristic
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 informationON THE TWO OLDEST FAMILIES FOR THE WRIGHT-FISHER PROCESS
ON THE TWO OLDEST FAMILIES FOR THE WRIGHT-FISHER PROCESS JEAN-FRANÇOIS DELMAS, JEAN-STÉPHANE DHERSIN, AND ARNO SIRI-JEGOUSSE Abstract. We extend some of the results of Pfaffelhuber and Walkobinger on the
More informationTranslation Invariant Experiments with Independent Increments
Translation Invariant Statistical Experiments with Independent Increments (joint work with Nino Kordzakhia and Alex Novikov Steklov Mathematical Institute St.Petersburg, June 10, 2013 Outline 1 Introduction
More informationEvolution in a spatial continuum
Evolution in a spatial continuum Drift, draft and structure Alison Etheridge University of Oxford Joint work with Nick Barton (Edinburgh) and Tom Kurtz (Wisconsin) New York, Sept. 2007 p.1 Kingman s Coalescent
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 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 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 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 informationConcentration inequalities and the entropy method
Concentration inequalities and the entropy method Gábor Lugosi ICREA and Pompeu Fabra University Barcelona what is concentration? We are interested in bounding random fluctuations of functions of many
More information