Path-properties of the tree-valued Fleming-Viot process

Similar documents
Equivalence between Semimartingales and Itô Processes

Last Time. Martingale inequalities Martingale convergence theorem Uniformly integrable martingales. Today s lecture: Sections 4.4.1, 5.

Model-independent bounds for Asian options

RMSC 4005 Stochastic Calculus for Finance and Risk. 1 Exercises. (c) Let X = {X n } n=0 be a {F n }-supermartingale. Show that.

Dr. Maddah ENMG 625 Financial Eng g II 10/16/06

Sidney I. Resnick. A Probability Path. Birkhauser Boston Basel Berlin

CONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES

Robust Hedging of Options on a Leveraged Exchange Traded Fund

Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs.

Model-independent bounds for Asian options

Non-semimartingales in finance

Homework Assignments

M5MF6. Advanced Methods in Derivatives Pricing

The stochastic calculus

Stochastic calculus Introduction I. Stochastic Finance. C. Azizieh VUB 1/91. C. Azizieh VUB Stochastic Finance

Exponential martingales and the UI martingale property

Large Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models

Brownian Motion. Richard Lockhart. Simon Fraser University. STAT 870 Summer 2011

BROWNIAN MOTION II. D.Majumdar

Martingale representation theorem

Change of Measure (Cameron-Martin-Girsanov Theorem)

Optimal robust bounds for variance options and asymptotically extreme models

Logarithmic derivatives of densities for jump processes

Drunken Birds, Brownian Motion, and Other Random Fun

6. Martingales. = Zn. Think of Z n+1 as being a gambler s earnings after n+1 games. If the game if fair, then E [ Z n+1 Z n

STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL

Rohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)

Martingales. by D. Cox December 2, 2009

Martingale Transport, Skorokhod Embedding and Peacocks

Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes

Stochastic Calculus, Application of Real Analysis in Finance

Optimal Order Placement

1.1 Basic Financial Derivatives: Forward Contracts and Options

The Double Skorohod Map and Real-Time Queues

Economathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t

Chapter 3: Black-Scholes Equation and Its Numerical Evaluation

S t d with probability (1 p), where

Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies

Order book resilience, price manipulations, and the positive portfolio problem

Stochastic Dynamical Systems and SDE s. An Informal Introduction

PAPER 27 STOCHASTIC CALCULUS AND APPLICATIONS

Applications of Lévy processes

Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model

L 2 -theoretical study of the relation between the LIBOR market model and the HJM model Takashi Yasuoka

Asymptotic methods in risk management. Advances in Financial Mathematics

DRAFT. 1 exercise in state (S, t), π(s, t) = 0 do not exercise in state (S, t) Review of the Risk Neutral Stock Dynamics

The discounted portfolio value of a selffinancing strategy in discrete time was given by. δ tj 1 (s tj s tj 1 ) (9.1) j=1

Stochastic modelling of electricity markets Pricing Forwards and Swaps

Stochastic Calculus - An Introduction

PAPER 211 ADVANCED FINANCIAL MODELS

Fundamentals of Stochastic Filtering

Structural Models of Credit Risk and Some Applications

Drawdowns Preceding Rallies in the Brownian Motion Model

Continuous Time Finance. Tomas Björk

Brownian Motion and Ito s Lemma

Barrier Options Pricing in Uncertain Financial Market

Lecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.

Chapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets

AMH4 - ADVANCED OPTION PRICING. Contents

Corrections to the Second Edition of Modeling and Analysis of Stochastic Systems

Averaging fast subsystems

Hedging under Arbitrage

Kim Weston (Carnegie Mellon University) Market Stability and Indifference Prices. 1st Eastern Conference on Mathematical Finance.

Sum-Product: Message Passing Belief Propagation

Are the Azéma-Yor processes truly remarkable?

KØBENHAVNS UNIVERSITET (Blok 2, 2011/2012) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours

Blackwell Optimality in Markov Decision Processes with Partial Observation

Saddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance

MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS.

Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models

A No-Arbitrage Theorem for Uncertain Stock Model

Minimal Variance Hedging in Large Financial Markets: random fields approach

The ruin probabilities of a multidimensional perturbed risk model

Lecture 3. Sergei Fedotov Introduction to Financial Mathematics. Sergei Fedotov (University of Manchester) / 6

Replication under Price Impact and Martingale Representation Property

1 Mathematics in a Pill 1.1 PROBABILITY SPACE AND RANDOM VARIABLES. A probability triple P consists of the following components:

Volatility. Roberto Renò. 2 March 2010 / Scuola Normale Superiore. Dipartimento di Economia Politica Università di Siena

Exam Quantitative Finance (35V5A1)

Replication and Absence of Arbitrage in Non-Semimartingale Models

Lecture 14: Examples of Martingales and Azuma s Inequality. Concentration

Multiple Defaults and Counterparty Risks by Density Approach

Value at Risk Ch.12. PAK Study Manual

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 11 10/9/2013. Martingales and stopping times II

Chapter 1. Bond Pricing (continued)

Optimal stopping problems for a Brownian motion with a disorder on a finite interval

A note on the existence of unique equivalent martingale measures in a Markovian setting

Conditional Density Method in the Computation of the Delta with Application to Power Market

No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate

A new approach to LIBOR modeling

Indifference fee rate 1

Theoretical Problems in Credit Portfolio Modeling 2

Are the Azéma-Yor processes truly remarkable?

THE MARTINGALE METHOD DEMYSTIFIED

Constructing Markov models for barrier options

Valuation of derivative assets Lecture 8

e-companion ONLY AVAILABLE IN ELECTRONIC FORM

1 The continuous time limit

IEOR E4703: Monte-Carlo Simulation

Derivative Securities Section 9 Fall 2004 Notes by Robert V. Kohn, Courant Institute of Mathematical Sciences.

Errata, Mahler Study Aids for Exam 3/M, Spring 2010 HCM, 1/26/13 Page 1

Transcription:

Path-properties of the tree-valued Fleming-Viot process Peter Pfaffelhuber Joint work with Andrej Depperschmidt and Andreas Greven Luminy, 492012

The Moran model time t As every population model, the Moran model (of size N) gives rise to a tree-valued process (X N t ) t 0 Without proof, we assume that (X N t ) t 0 N === (X t ) t 0 in an appropriate sense We call (X t ) t 0 the tree-valued Fleming-Viot dynamics

The Moran model time t s As every population model, the Moran model (of size N) gives rise to a tree-valued process (X N t ) t 0 Without proof, we assume that (X N t ) t 0 N === (X t ) t 0 in an appropriate sense We call (X t ) t 0 the tree-valued Fleming-Viot dynamics

The Moran model time t s As every population model, the Moran model (of size N) gives rise to a tree-valued process (X N t ) t 0 Without proof, we assume that (X N t ) t 0 N === (X t ) t 0 in an appropriate sense We call (X t ) t 0 the tree-valued Fleming-Viot dynamics

The tree-valued Fleming-Viot dynamics Theorem The process (X t ) t 0 exists as limit of tree-valued Moran models and is unique Its state space is the set of (equivalence classes of) metric measure spaces {(X, r, µ) : µ M 1 (X )}, equipped with the Gromov-Prohorov topology It can be described by a martingale problem Almost surely, (Xt ) t 0 has continuous sample paths (Xt ) t 0 is compact for all t > 0 For many functions Φ, the quadratic variation of (Φ(X t )) t 0 can be computed

The Kingman coalescent as the equilibrium of (X t ) t 0 Theorem Let X be the Kingman coalescent, a random tree with: Start with many lines If there are k lines left, wait S k exp ( k 2) and merge two randomly chosen lines Stop upon reaching one line Then, X t t === X Important property: subtree with n leaves tree started with n lines

Goal Lift properties of Kingman coalescent X to the paths of (X t ) t 0 (when started in equilibrium) Example: Let N t ε be the number of ancestors of the time-t population at time t ε It is well-known that almost surely εn ε 2 ε 0 0 Is it also true that almost surely sup εnε t 2 ε 0 0? t 0

A law of large numbers for the number of ancestors Theorem Let N ε be the the number of ancestors of the Kingman coalescent X at time ε Then, almost surely, εn ε 2 ε 0 0

A law of large numbers for the number of ancestors εn ε 2 ε 0 0 Proof: Recall S k exp ( k 2) is the time the coalescent has k lines The assertion is the same as n (S n+1 + S n+2 + )n 2 0 }{{} =:T n With E[T n ] = 2/n and E[(T n E[T n ]) 4 ] 1 n 6, we find P( T n n 2 > ε) n4 E[(T n E[T n ]) 4 ] ε 4 1 ε 4 n 2 The result follows from the Borel-Cantelli-Lemma

A law of large numbers for the number of ancestors Theorem Let N t ε be the the number of ancestors of the time-t population X t, at time t ε Then, almost surely, sup εnε t 2 ε 0 0 t 0

A law of large numbers for the number of ancestors sup εnε t 2 ε 0 0 t 0 Proof: Let T t n = S t n + S t n+1 + and S t n is the time the time-t tree spends with n lines It suffices to show sup 0 t 1 T t nn 2 n 0 Using moment calculations, ( P sup Tn k/n2 k=0,,n 2 ) n 2 > ε 1 ε 8 n 2

A law of large numbers for the number of ancestors It suffices to show sup 0 t 1 Tnn t 2 n 0 We claim that {T t n > s} T t δ n time t time t δ > s δ} ( P sup 0 t 1 ( P ) Tnn t > 2 + ε sup 0 t 1 T tn2 /n 2 n >nancestors of time t population > 2 + ε n ( t tn2 )) }{{ n 2 } 1/n 2 1 ε 8 n 2

Small family sizes Theorem Let N t (x, ε) be the number of ancestors of the time-t population with families of size at least x Then, almost surely sup εn (xε, ε) 2e 2x ε 0 0 0 x< Open problem: Is it also true, that sup sup t 0 0 x< εn t (xε, ε) 2e 2x ε 0 0?

A law of large numbers for the tree metric Theorem Let F t 1 (ε),, F t N ε (ε) be the family sizes of the ancestors 1,, N t ε in X t Then, almost surely, 1 ε N t ε (Fi t (ε)) 2 = lim i=1 N 1 εn 2 N u,v=1 leaves in X t 1 {r t(u,v)<ε} ε 0 1

A law of large numbers for the tree metric Lemma lim N λ N 2 N u,v=1 leaves in X 1 {rt(u,v)<1/λ} λ 1 Ψ λ (X ) λ 1 with 1 Ψ λ (X ) := (λ + 1) lim N N 2 N u,v=1 leaves in X e λrt(u,v), where r t (u, v) is the time to the most recent common ancestor of u and v in X

A law of large numbers for the tree metric Theorem For Ψ λ as above, almost surely, Ψ λ (X ) 1 λ 0 ( ) Proof: For random leaves U, V, U, V from X, E[Ψ λ (X )] = E[(λ + 1)e λr(u,v ) ] = 1, E[(Ψ λ (X ) 1) 2 ] = (λ + 1) 2( E[e λ(r(u,v )+r(u,v )) ] 1 ) = some calculations on tree with 4 leaves 2λ 2 λ 1 = (λ + 3)(2λ + 1)(2λ + 3) 2λ, E[ ( Ψ λ (X ) 1) 4 ] = λ 3 4λ 2

A law of large numbers for the tree metric Theorem For Ψ λ as above, in probability, Ψ λ (X t ) 1 λ 0 sup 0 t T ( ) ( ) fdd-convergence in ( ) Lemma tightness in C R [0, ) sup E[(Ψ λ (X t ) Ψ λ (X 0 )) 4 ] t 2 λ>0

A law of large numbers for the tree metric Proof of Lemma: Recall E[Ψ λ (X t )] = (λ + 1)E[e λrt(u,v ) ] By the dynamics of the Moran model, d dt E[Ψ λ(x t )] = λ + 1 dt ( dt E[e λ 0 e λrt(u,v ) ] }{{} resampling between U and V + (1 dt) E[e λ(rt(u,v )+dt) e λrt(u,v ) ] }{{} tree growth = (λ + 1)(1 E[(λ + 1)e λrt(u,v ) ]) = (λ + 1)(E[Ψ λ (X t )] 1) )

A law of large numbers for the tree metric Proof of Lemma: Recall Similarly, d dt (E[Ψ λ(x t )] 1) = (λ + 1)(E[Ψ λ (X t )] 1) E[Ψ λ (X t ) 1 F s ] = e (λ+1)(t s) (Ψ λ (X s ) 1) and ( e (λ+1)t( )) Ψ λ (X t ) 1 is a martingale t 0

A law of large numbers for the tree metric Proof of Lemma: E[(Ψ λ (X t ) Ψ λ (X 0 )) 2 ] = E[2Ψ λ (X 0 )(Ψ λ (X t ) Ψ λ (X 0 ))] [ = 2E Ψ λ (X 0 ) ( [ e (λ+1)t E e (λ+1)t( ) ] Ψ λ (X t ) 1 F 0 }{{} martingale ( )] Ψ λ (X 0 ) 1 [ ] = 2E Ψ λ (X 0 )(Ψ λ (X 0 ) 1)(1 e (λ+1)t ) λ 1 e (λ+1)t λ t

A law of large numbers for the tree metric Proof of Lemma: E[(Ψ λ (X t ) Ψ λ (X 0 )) 4 ] = E[2Ψ λ (X 0 )(Ψ λ (X t ) Ψ λ (X 0 ))] [ = 2E Ψ λ (X 0 ) ( [ e (λ+1)t E e (λ+1)t( ) ] Ψ λ (X t ) 1 F 0 }{{} martingale ( )] Ψ λ (X 0 ) 1 [ = 2E Ψ λ (X 0 )(Ψ λ (X 0 ) 1)(1 e (λ+1)t )] λ 1 e (λ+1)t 2λ t 2

A Brownian motion in the Fleming-Viot dynamics Theorem Let W λ = (W λ (t)) t 0 be given by W λ (t) := λ t 0 (Ψ λ (X s ) 1)ds Then, W λ λ === W, where W = (W t ) t 0 is a Brownian motion

A Brownian motion in the Fleming-Viot dynamics W λ = λ Proof (part): E[W λ (t) 2 ] = 2λ 2 t 0 0 = 2λ 2 t λ 0 λ t 0 ( Ψ λ (X s ) 1 }{{} mean = 0 variance 1/2λ s 0 s 0 s 0 )ds λ === W E [ E[Ψ λ (X s ) 1 F r ](Ψ λ (X r ) 1) ] dr ds e (λ+1)(s r) E[(Ψ λ (X r ) 1) 2 ]dr ds e (λ+1)r dr ds λ t

Summary and outlook More about formalising trees (and Gromov-Prohorov convergence) and construction of tree-valued processes (via well-posed martingale problem) can be said All result also hold in models with mutation and selection (individuals also carry types which are (dis)favored to get offspring) All Theorems affect properties near the tree top do similar properties hold for branching trees?