Some Mathematical Models from Population Genetics 4: Spatial Models

Total Page:16

File Type:pdf, Size:1020Kb

Some Mathematical Models from Population Genetics 4: Spatial Models Some mathematical models from population genetics 4: Spatial models Alison Etheridge University of Oxford joint work with Jochen Blath (TU Berlin), Mark Meredith (Oxford) Nick Barton (Edinburgh), Frantz Depaulis (Paris VI) New York, Sept. 07 – p.1 Lessons learned so far The key message of our balancing selection example is that if we want to study the genealogy of a sample from a structured population, then fluctuations in background frequencies matter. A central question then is how should we model fluctuations of spatially distributed populations. New York, Sept. 07 – p.2 Large population, long timescales, measured in units of size N. Write Zn for the population size after n generations. 0 N If Φ (1) = 1 + a, then E[ZN ] = (1 + a) Z0, so for non-trivial limit a assume a = N . If Z0 converges, so does ZbNtc . N N≥1 N N≥1 Limit process: dXt = aXtdt + pγXtdBt: Reminder: Feller’s rescaling Galton Watson process, offspring generating function Φ(s). Assume Φ00(1) < . 1 New York, Sept. 07 – p.3 0 N If Φ (1) = 1 + a, then E[ZN ] = (1 + a) Z0, so for non-trivial limit a assume a = N . If Z0 converges, so does ZbNtc . N N≥1 N N≥1 Limit process: dXt = aXtdt + pγXtdBt: Reminder: Feller’s rescaling Galton Watson process, offspring generating function Φ(s). Assume Φ00(1) < . 1 Large population, long timescales, measured in units of size N. Write Zn for the population size after n generations. New York, Sept. 07 – p.3 If Z0 converges, so does ZbNtc . N N≥1 N N≥1 Limit process: dXt = aXtdt + pγXtdBt: Reminder: Feller’s rescaling Galton Watson process, offspring generating function Φ(s). Assume Φ00(1) < . 1 Large population, long timescales, measured in units of size N. Write Zn for the population size after n generations. 0 N If Φ (1) = 1 + a, then E[ZN ] = (1 + a) Z0, so for non-trivial limit a assume a = N . New York, Sept. 07 – p.3 Reminder: Feller’s rescaling Galton Watson process, offspring generating function Φ(s). Assume Φ00(1) < . 1 Large population, long timescales, measured in units of size N. Write Zn for the population size after n generations. 0 N If Φ (1) = 1 + a, then E[ZN ] = (1 + a) Z0, so for non-trivial limit a assume a = N . If Z0 converges, so does ZbNtc . N N≥1 N N≥1 Limit process: dXt = aXtdt + pγXtdBt: New York, Sept. 07 – p.3 Galton-Watson branching process branching Brownian motion/ branching random walk. Offspring born where parent died. 1 Feller rescaling: individual represented by atom of mass N , time in units of size N. Spatially distributed populations Populations dispersed in Rd or Zd, New York, Sept. 07 – p.4 1 Feller rescaling: individual represented by atom of mass N , time in units of size N. Spatially distributed populations Populations dispersed in Rd or Zd, Galton-Watson branching process branching Brownian motion/ branching random walk. Offspring born where parent died. New York, Sept. 07 – p.4 Spatially distributed populations Populations dispersed in Rd or Zd, Galton-Watson branching process branching Brownian motion/ branching random walk. Offspring born where parent died. 1 Feller rescaling: individual represented by atom of mass N , time in units of size N. New York, Sept. 07 – p.4 Super-random walk: dX (t) = m (X (t) X (t)) dt+aX (t)dt+ γX (t)dW (t); i Zd; i ij j − i i i i 2 Xj p where W (t); t 0 Zd is a collection of independent Brownian mo- f i ≥ gi2 tions and Xi(t) is the size of the population in deme i at time t. The limiting processes The Dawson-Watanabe superprocess: For positive, twice differentiable test functions φ, t t φ, Xt φ, X0 D∆φ, Xs ds aφ, Xs ds h i − h i − Z0 h i − Z0 h i is a martingale with quadratic variation t γφ2; X ds: 0 h si R New York, Sept. 07 – p.5 The limiting processes The Dawson-Watanabe superprocess: For positive, twice differentiable test functions φ, t t φ, Xt φ, X0 D∆φ, Xs ds aφ, Xs ds h i − h i − Z0 h i − Z0 h i is a martingale with quadratic variation t γφ2; X ds: 0 h si Super-random walk: R dX (t) = m (X (t) X (t)) dt+aX (t)dt+ γX (t)dW (t); i Zd; i ij j − i i i i 2 Xj p where W (t); t 0 Zd is a collection of independent Brownian mo- f i ≥ gi2 tions and Xi(t) is the size of the population in deme i at time t. New York, Sept. 07 – p.5 Clumping and extinction. E [ φ, X ] = eat T φ, X ; h ti h t 0i Take a = 0. t 2 var ( φ, Xt ) = γTt−s (Tsφ) ; X0 ds: h i Z h i 0 In one and two dimensions grows without bound. New York, Sept. 07 – p.6 As we saw in the first lecture, in Zd can specify the population size locally The Classical Stepping Stone Model. Populations should be regulated by local rules. Individuals living in locally crowded regions will have a lower reproduc- tive success than those living in sparsely populated regions. Controlling the population Exogenously specify total population size Fleming-Viot superprocess. New York, Sept. 07 – p.7 Populations should be regulated by local rules. Individuals living in locally crowded regions will have a lower reproduc- tive success than those living in sparsely populated regions. Controlling the population Exogenously specify total population size Fleming-Viot superprocess. As we saw in the first lecture, in Zd can specify the population size locally The Classical Stepping Stone Model. New York, Sept. 07 – p.7 Individuals living in locally crowded regions will have a lower reproduc- tive success than those living in sparsely populated regions. Controlling the population Exogenously specify total population size Fleming-Viot superprocess. As we saw in the first lecture, in Zd can specify the population size locally The Classical Stepping Stone Model. Populations should be regulated by local rules. New York, Sept. 07 – p.7 Controlling the population Exogenously specify total population size Fleming-Viot superprocess. As we saw in the first lecture, in Zd can specify the population size locally The Classical Stepping Stone Model. Populations should be regulated by local rules. Individuals living in locally crowded regions will have a lower reproduc- tive success than those living in sparsely populated regions. New York, Sept. 07 – p.7 For infinite initial measures, to prevent immediate catastrophe, 1 h(r)rd−1dr < . 0 1 R The stepping-stone version of the Bolker-Pacala model: In the super-random walk setting the corresponding model is dXt(i) = mij (Xt(j) Xt(i)) dt + α 0M λij Xt(j)1 Xt(i)dt − − Xj Xj @ A (i) + γXt(i)dBt : p Note that moment equations not closed. Locally regulated populations a(s; x) = α M h(x; y); X (dy) : − h s i For simplicity h(x; y) = h( x y ). k − k New York, Sept. 07 – p.8 The stepping-stone version of the Bolker-Pacala model: In the super-random walk setting the corresponding model is dXt(i) = mij (Xt(j) Xt(i)) dt + α 0M λij Xt(j)1 Xt(i)dt − − Xj Xj @ A (i) + γXt(i)dBt : p Note that moment equations not closed. Locally regulated populations a(s; x) = α M h(x; y); X (dy) : − h s i For simplicity h(x; y) = h( x y ). For infinite initial measures, to k − k prevent immediate catastrophe, 1 h(r)rd−1dr < . 0 1 R New York, Sept. 07 – p.8 Note that moment equations not closed. Locally regulated populations a(s; x) = α M h(x; y); X (dy) : − h s i For simplicity h(x; y) = h( x y ). For infinite initial measures, to k − k prevent immediate catastrophe, 1 h(r)rd−1dr < . 0 1 R The stepping-stone version of the Bolker-Pacala model: In the super-random walk setting the corresponding model is dXt(i) = mij (Xt(j) Xt(i)) dt + α 0M λij Xt(j)1 Xt(i)dt − − Xj Xj @ A (i) + γXt(i)dBt : p New York, Sept. 07 – p.8 Locally regulated populations a(s; x) = α M h(x; y); X (dy) : − h s i For simplicity h(x; y) = h( x y ). For infinite initial measures, to k − k prevent immediate catastrophe, 1 h(r)rd−1dr < . 0 1 R The stepping-stone version of the Bolker-Pacala model: In the super-random walk setting the corresponding model is dXt(i) = mij (Xt(j) Xt(i)) dt + α 0M λij Xt(j)1 Xt(i)dt − − Xj Xj @ A (i) + γXt(i)dBt : Note that moment equations not closed. p New York, Sept. 07 – p.8 Survival and Extinction Theorem For each fixed interaction kernel h and γ; K > 0 there exists α0 = α0(K; γ; h) such that for α > α0, the superprocess version of the Bolker-Pacala model with parameters (h; K/α; α; γ) started from any finite initial measure dies out in finite time. If h also satisfies h(r)rd−1dr < , then when started from any tempered initial 1 measureR (with p > d) the process with these parameters suffers local extinction. New York, Sept. 07 – p.9 Let α > 0 be fixed. If r2−δh(r) is unbounded for some δ > 0, then for each fixed • γ > 0, there is an M0 > 0 such that for M < M0 the superprocess version of the Bolker-Pacala model with parameters (h; M; α; γ) started from any finite initial measure dies out in finite time.
Recommended publications
  • Poisson Representations of Branching Markov and Measure-Valued
    The Annals of Probability 2011, Vol. 39, No. 3, 939–984 DOI: 10.1214/10-AOP574 c Institute of Mathematical Statistics, 2011 POISSON REPRESENTATIONS OF BRANCHING MARKOV AND MEASURE-VALUED BRANCHING PROCESSES By Thomas G. Kurtz1 and Eliane R. Rodrigues2 University of Wisconsin, Madison and UNAM Representations of branching Markov processes and their measure- valued limits in terms of countable systems of particles are con- structed for models with spatially varying birth and death rates. Each particle has a location and a “level,” but unlike earlier con- structions, the levels change with time. In fact, death of a particle occurs only when the level of the particle crosses a specified level r, or for the limiting models, hits infinity. For branching Markov pro- cesses, at each time t, conditioned on the state of the process, the levels are independent and uniformly distributed on [0,r]. For the limiting measure-valued process, at each time t, the joint distribu- tion of locations and levels is conditionally Poisson distributed with mean measure K(t) × Λ, where Λ denotes Lebesgue measure, and K is the desired measure-valued process. The representation simplifies or gives alternative proofs for a vari- ety of calculations and results including conditioning on extinction or nonextinction, Harris’s convergence theorem for supercritical branch- ing processes, and diffusion approximations for processes in random environments. 1. Introduction. Measure-valued processes arise naturally as infinite sys- tem limits of empirical measures of finite particle systems. A number of ap- proaches have been developed which preserve distinct particles in the limit and which give a representation of the measure-valued process as a transfor- mation of the limiting infinite particle system.
    [Show full text]
  • Superprocesses and Mckean-Vlasov Equations with Creation of Mass
    Sup erpro cesses and McKean-Vlasov equations with creation of mass L. Overb eck Department of Statistics, University of California, Berkeley, 367, Evans Hall Berkeley, CA 94720, y U.S.A. Abstract Weak solutions of McKean-Vlasov equations with creation of mass are given in terms of sup erpro cesses. The solutions can b e approxi- mated by a sequence of non-interacting sup erpro cesses or by the mean- eld of multityp e sup erpro cesses with mean- eld interaction. The lat- ter approximation is asso ciated with a propagation of chaos statement for weakly interacting multityp e sup erpro cesses. Running title: Sup erpro cesses and McKean-Vlasov equations . 1 Intro duction Sup erpro cesses are useful in solving nonlinear partial di erential equation of 1+ the typ e f = f , 2 0; 1], cf. [Dy]. Wenowchange the p oint of view and showhowtheyprovide sto chastic solutions of nonlinear partial di erential Supp orted byanFellowship of the Deutsche Forschungsgemeinschaft. y On leave from the Universitat Bonn, Institut fur Angewandte Mathematik, Wegelerstr. 6, 53115 Bonn, Germany. 1 equation of McKean-Vlasovtyp e, i.e. wewant to nd weak solutions of d d 2 X X @ @ @ + d x; + bx; : 1.1 = a x; t i t t t t t ij t @t @x @x @x i j i i=1 i;j =1 d Aweak solution = 2 C [0;T];MIR satis es s Z 2 t X X @ @ a f = f + f + d f + b f ds: s ij s t 0 i s s @x @x @x 0 i j i Equation 1.1 generalizes McKean-Vlasov equations of twodi erenttyp es.
    [Show full text]
  • Local Conditioning in Dawson–Watanabe Superprocesses
    The Annals of Probability 2013, Vol. 41, No. 1, 385–443 DOI: 10.1214/11-AOP702 c Institute of Mathematical Statistics, 2013 LOCAL CONDITIONING IN DAWSON–WATANABE SUPERPROCESSES By Olav Kallenberg Auburn University Consider a locally finite Dawson–Watanabe superprocess ξ =(ξt) in Rd with d ≥ 2. Our main results include some recursive formulas for the moment measures of ξ, with connections to the uniform Brown- ian tree, a Brownian snake representation of Palm measures, continu- ity properties of conditional moment densities, leading by duality to strongly continuous versions of the multivariate Palm distributions, and a local approximation of ξt by a stationary clusterη ˜ with nice continuity and scaling properties. This all leads up to an asymptotic description of the conditional distribution of ξt for a fixed t> 0, given d that ξt charges the ε-neighborhoods of some points x1,...,xn ∈ R . In the limit as ε → 0, the restrictions to those sets are conditionally in- dependent and given by the pseudo-random measures ξ˜ orη ˜, whereas the contribution to the exterior is given by the Palm distribution of ξt at x1,...,xn. Our proofs are based on the Cox cluster representa- tions of the historical process and involve some delicate estimates of moment densities. 1. Introduction. This paper may be regarded as a continuation of [19], where we considered some local properties of a Dawson–Watanabe super- process (henceforth referred to as a DW-process) at a fixed time t> 0. Recall that a DW-process ξ = (ξt) is a vaguely continuous, measure-valued diffu- d ξtf µvt sion process in R with Laplace functionals Eµe− = e− for suitable functions f 0, where v = (vt) is the unique solution to the evolution equa- 1 ≥ 2 tion v˙ = 2 ∆v v with initial condition v0 = f.
    [Show full text]
  • Lecture 19: Polymer Models: Persistence and Self-Avoidance
    Lecture 19: Polymer Models: Persistence and Self-Avoidance Scribe: Allison Ferguson (and Martin Z. Bazant) Martin Fisher School of Physics, Brandeis University April 14, 2005 Polymers are high molecular weight molecules formed by combining a large number of smaller molecules (monomers) in a regular pattern. Polymers in solution (i.e. uncrystallized) can be characterized as long flexible chains, and thus may be well described by random walk models of various complexity. 1 Simple Random Walk (IID Steps) Consider a Rayleigh random walk (i.e. a fixed step size a with a random angle θ). Recall (from Lecture 1) that the probability distribution function (PDF) of finding the walker at position R~ after N steps is given by −3R2 2 ~ e 2a N PN (R) ∼ d (1) (2πa2N/d) 2 Define RN = |R~| as the end-to-end distance of the polymer. Then from previous results we q √ 2 2 ¯ 2 know hRN i = Na and RN = hRN i = a N. We can define an additional quantity, the radius of gyration GN as the radius from the centre of mass (assuming the polymer has an approximately 1 2 spherical shape ). Hughes [1] gives an expression for this quantity in terms of hRN i as follows: 1 1 hG2 i = 1 − hR2 i (2) N 6 N 2 N As an example of the type of prediction which can be made using this type of simple model, consider the following experiment: Pull on both ends of a polymer and measure the force f required to stretch it out2. Note that the force will be given by f = −(∂F/∂R) where F = U − TS is the Helmholtz free energy (T is the temperature, U is the internal energy, and S is the entropy).
    [Show full text]
  • 8 Convergence of Loop-Erased Random Walk to SLE2 8.1 Comments on the Paper
    8 Convergence of loop-erased random walk to SLE2 8.1 Comments on the paper The proof that the loop-erased random walk converges to SLE2 is in the paper \Con- formal invariance of planar loop-erased random walks and uniform spanning trees" by Lawler, Schramm and Werner. (arxiv.org: math.PR/0112234). This section contains some comments to give you some hope of being able to actually make sense of this paper. In the first sections of the paper, δ is used to denote the lattice spacing. The main theorem is that as δ 0, the loop erased random walk converges to SLE2. HOWEVER, beginning in section !2.2 the lattice spacing is 1, and in section 3.2 the notation δ reappears but means something completely different. (δ also appears in the proof of lemma 2.1 in section 2.1, but its meaning here is restricted to this proof. The δ in this proof has nothing to do with any δ's outside the proof.) After the statement of the main theorem, there is no mention of the lattice spacing going to zero. This is hidden in conditions on the \inner radius." For a domain D containing 0 the inner radius is rad0(D) = inf z : z = D (364) fj j 2 g Now fix a domain D containing the origin. Let be the conformal map of D onto U. We want to show that the image under of a LERW on a lattice with spacing δ converges to SLE2 in the unit disc. In the paper they do this in a slightly different way.
    [Show full text]
  • Non-Local Branching Superprocesses and Some Related Models
    Published in: Acta Applicandae Mathematicae 74 (2002), 93–112. Non-local Branching Superprocesses and Some Related Models Donald A. Dawson1 School of Mathematics and Statistics, Carleton University, 1125 Colonel By Drive, Ottawa, Canada K1S 5B6 E-mail: [email protected] Luis G. Gorostiza2 Departamento de Matem´aticas, Centro de Investigaci´ony de Estudios Avanzados, A.P. 14-740, 07000 M´exicoD. F., M´exico E-mail: [email protected] Zenghu Li3 Department of Mathematics, Beijing Normal University, Beijing 100875, P.R. China E-mail: [email protected] Abstract A new formulation of non-local branching superprocesses is given from which we derive as special cases the rebirth, the multitype, the mass- structured, the multilevel and the age-reproduction-structured superpro- cesses and the superprocess-controlled immigration process. This unified treatment simplifies considerably the proof of existence of the old classes of superprocesses and also gives rise to some new ones. AMS Subject Classifications: 60G57, 60J80 Key words and phrases: superprocess, non-local branching, rebirth, mul- titype, mass-structured, multilevel, age-reproduction-structured, superprocess- controlled immigration. 1Supported by an NSERC Research Grant and a Max Planck Award. 2Supported by the CONACYT (Mexico, Grant No. 37130-E). 3Supported by the NNSF (China, Grant No. 10131040). 1 1 Introduction Measure-valued branching processes or superprocesses constitute a rich class of infinite dimensional processes currently under rapid development. Such processes arose in appli- cations as high density limits of branching particle systems; see e.g. Dawson (1992, 1993), Dynkin (1993, 1994), Watanabe (1968). The development of this subject has been stimu- lated from different subjects including branching processes, interacting particle systems, stochastic partial differential equations and non-linear partial differential equations.
    [Show full text]
  • A Representation for Functionals of Superprocesses Via Particle Picture Raisa E
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector stochastic processes and their applications ELSEVIER Stochastic Processes and their Applications 64 (1996) 173-186 A representation for functionals of superprocesses via particle picture Raisa E. Feldman *,‘, Srikanth K. Iyer ‘,* Department of Statistics and Applied Probability, University oj’ Cull~ornia at Santa Barbara, CA 93/06-31/O, USA, and Department of Industrial Enyineeriny and Management. Technion - Israel Institute of’ Technology, Israel Received October 1995; revised May 1996 Abstract A superprocess is a measure valued process arising as the limiting density of an infinite col- lection of particles undergoing branching and diffusion. It can also be defined as a measure valued Markov process with a specified semigroup. Using the latter definition and explicit mo- ment calculations, Dynkin (1988) built multiple integrals for the superprocess. We show that the multiple integrals of the superprocess defined by Dynkin arise as weak limits of linear additive functionals built on the particle system. Keywords: Superprocesses; Additive functionals; Particle system; Multiple integrals; Intersection local time AMS c.lassijication: primary 60555; 60H05; 60580; secondary 60F05; 6OG57; 60Fl7 1. Introduction We study a class of functionals of a superprocess that can be identified as the multiple Wiener-It6 integrals of this measure valued process. More precisely, our objective is to study these via the underlying branching particle system, thus providing means of thinking about the functionals in terms of more simple, intuitive and visualizable objects. This way of looking at multiple integrals of a stochastic process has its roots in the previous work of Adler and Epstein (=Feldman) (1987), Epstein (1989), Adler ( 1989) Feldman and Rachev (1993), Feldman and Krishnakumar ( 1994).
    [Show full text]
  • Random Walk in Random Scenery (RWRS)
    IMS Lecture Notes–Monograph Series Dynamics & Stochastics Vol. 48 (2006) 53–65 c Institute of Mathematical Statistics, 2006 DOI: 10.1214/074921706000000077 Random walk in random scenery: A survey of some recent results Frank den Hollander1,2,* and Jeffrey E. Steif 3,† Leiden University & EURANDOM and Chalmers University of Technology Abstract. In this paper we give a survey of some recent results for random walk in random scenery (RWRS). On Zd, d ≥ 1, we are given a random walk with i.i.d. increments and a random scenery with i.i.d. components. The walk and the scenery are assumed to be independent. RWRS is the random process where time is indexed by Z, and at each unit of time both the step taken by the walk and the scenery value at the site that is visited are registered. We collect various results that classify the ergodic behavior of RWRS in terms of the characteristics of the underlying random walk (and discuss extensions to stationary walk increments and stationary scenery components as well). We describe a number of results for scenery reconstruction and close by listing some open questions. 1. Introduction Random walk in random scenery is a family of stationary random processes ex- hibiting amazingly rich behavior. We will survey some of the results that have been obtained in recent years and list some open questions. Mike Keane has made funda- mental contributions to this topic. As close colleagues it has been a great pleasure to work with him. We begin by defining the object of our study. Fix an integer d ≥ 1.
    [Show full text]
  • Immigration Structures Associated with Dawson-Watanabe Superprocesses
    View metadata, citation and similar papers at core.ac.uk brought to you COREby provided by Elsevier - Publisher Connector stochastic processes and their applications ELSEVIER StochasticProcesses and their Applications 62 (1996) 73 86 Immigration structures associated with Dawson-Watanabe superprocesses Zeng-Hu Li Department of Mathematics, Beijing Normal University, Be!ring 100875, Peoples Republic of China Received March 1995: revised November 1995 Abstract The immigration structure associated with a measure-valued branching process may be described by a skew convolution semigroup. For the special type of measure-valued branching process, the Dawson Watanabe superprocess, we show that a skew convolution semigroup corresponds uniquely to an infinitely divisible probability measure on the space of entrance laws for the underlying process. An immigration process associated with a Borel right superpro- cess does not always have a right continuous realization, but it can always be obtained by transformation from a Borel right one in an enlarged state space. Keywords: Dawson-Watanabe superprocess; Skew convolution semigroup; Entrance law; Immigration process; Borel right process AMS classifications: Primary 60J80; Secondary 60G57, 60J50, 60J40 I. Introduction Let E be a Lusin topological space, i.e., a homeomorphism of a Borel subset of a compact metric space, with the Borel a-algebra .N(E). Let B(E) denote the set of all bounded ~(E)-measurable functions on E and B(E) + the subspace of B(E) compris- ing non-negative elements. Denote by M(E) the space of finite measures on (E, ~(E)) equipped with the topology of weak convergence. Forf e B(E) and # e M(E), write ~l(f) for ~Efd/~.
    [Show full text]
  • Particle Filtering Within Adaptive Metropolis Hastings Sampling
    Particle filtering within adaptive Metropolis-Hastings sampling Ralph S. Silva Paolo Giordani School of Economics Research Department University of New South Wales Sveriges Riksbank [email protected] [email protected] Robert Kohn∗ Michael K. Pitt School of Economics Economics Department University of New South Wales University of Warwick [email protected] [email protected] October 29, 2018 Abstract We show that it is feasible to carry out exact Bayesian inference for non-Gaussian state space models using an adaptive Metropolis-Hastings sampling scheme with the likelihood approximated by the particle filter. Furthermore, an adaptive independent Metropolis Hastings sampler based on a mixture of normals proposal is computation- ally much more efficient than an adaptive random walk Metropolis proposal because the cost of constructing a good adaptive proposal is negligible compared to the cost of approximating the likelihood. Independent Metropolis-Hastings proposals are also arXiv:0911.0230v1 [stat.CO] 2 Nov 2009 attractive because they are easy to run in parallel on multiple processors. We also show that when the particle filter is used, the marginal likelihood of any model is obtained in an efficient and unbiased manner, making model comparison straightforward. Keywords: Auxiliary variables; Bayesian inference; Bridge sampling; Marginal likeli- hood. ∗Corresponding author. 1 Introduction We show that it is feasible to carry out exact Bayesian inference on the parameters of a general state space model by using the particle filter to approximate the likelihood and adaptive Metropolis-Hastings sampling to generate unknown parameters. The state space model can be quite general, but we assume that the observation equation can be evaluated analytically and that it is possible to generate from the state transition equation.
    [Show full text]
  • Resource Theories of Multi-Time Processes: a Window Into Quantum Non-Markovianity
    Resource theories of multi-time processes: A window into quantum non-Markovianity Graeme D. Berk1, Andrew J. P. Garner2,3, Benjamin Yadin4, Kavan Modi1, and Felix A. Pollock1 1School of Physics and Astronomy, Monash University, Clayton, Victoria 3800, Australia 2Institute for Quantum Optics and Quantum Information, Austrian Academy of Sciences, Boltzmanngasse 3, A-1090 Vienna, Austria 3School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, 637371, Singapore 4School of Mathematical Sciences, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom. April 8, 2021 We investigate the conditions under which an uncontrolled background process may be harnessed by an agent to perform a task that would otherwise be impossible within their operational framework. This situation can be understood from the perspective of resource theory: rather than harnessing ‘useful’ quantum states to perform tasks, we propose a resource theory of quantum processes across multiple points in time. Uncontrolled background processes fulfil the role of resources, and a new set of objects called superprocesses, corresponding to operationally implementable control of the system undergoing the process, constitute the transformations between them. After formally introducing a framework for deriving resource theories of multi-time processes, we present a hierarchy of examples induced by restricting quantum or classical communication within the superprocess – corresponding to a client-server scenario. The resulting nine resource theories have different notions of quantum or classical memory as the determinant of their utility. Furthermore, one of these theories has a strict correspondence between non- useful processes and those that are Markovian and, therefore, could be said to be a true ‘quantum resource theory of non-Markovianity’.
    [Show full text]
  • A Course in Interacting Particle Systems
    A Course in Interacting Particle Systems J.M. Swart January 14, 2020 arXiv:1703.10007v2 [math.PR] 13 Jan 2020 2 Contents 1 Introduction 7 1.1 General set-up . .7 1.2 The voter model . .9 1.3 The contact process . 11 1.4 Ising and Potts models . 14 1.5 Phase transitions . 17 1.6 Variations on the voter model . 20 1.7 Further models . 22 2 Continuous-time Markov chains 27 2.1 Poisson point sets . 27 2.2 Transition probabilities and generators . 30 2.3 Poisson construction of Markov processes . 31 2.4 Examples of Poisson representations . 33 3 The mean-field limit 35 3.1 Processes on the complete graph . 35 3.2 The mean-field limit of the Ising model . 36 3.3 Analysis of the mean-field model . 38 3.4 Functions of Markov processes . 42 3.5 The mean-field contact process . 47 3.6 The mean-field voter model . 49 3.7 Exercises . 51 4 Construction and ergodicity 53 4.1 Introduction . 53 4.2 Feller processes . 54 4.3 Poisson construction . 63 4.4 Generator construction . 72 4.5 Ergodicity . 79 4.6 Application to the Ising model . 81 4.7 Further results . 85 5 Monotonicity 89 5.1 The stochastic order . 89 5.2 The upper and lower invariant laws . 94 5.3 The contact process . 97 5.4 Other examples . 100 3 4 CONTENTS 5.5 Exercises . 101 6 Duality 105 6.1 Introduction . 105 6.2 Additive systems duality . 106 6.3 Cancellative systems duality . 113 6.4 Other dualities .
    [Show full text]