Source author record

A. D. Barbour

A. D. Barbour appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2020arXiv

Central moment inequalities using Stein's method

We derive explicit central moment inequalities for random variables that admit a Stein coupling, such as exchangeable pairs, size--bias couplings or local dependence, among others. The bounds are in terms of moments (not necessarily central) of variables in the Stein coupling, which are typically local in some sense, and therefore easier to bound. In cases where the Stein couplings have the kind of behaviour leading to good normal approximation, the central moments are closely bounded by those of a normal. We show how the bounds can be used to produce concentration inequalities, and compare them to those existing in related settings. Finally, we illustrate the power of the theory by bounding the central moments of sums of neighbourhood statistics in sparse Erdős--Rényi random graphs.

preprint2020arXiv

Local limit theorems for occupancy models

We present a rather general method for proving local limit theorems, with a good rate of convergence, for sums of dependent random variables. The method is applicable when a Stein coupling can be exhibited. Our approach involves both Stein's method for distributional approximation and Stein's method for concentration. As applications, we prove local central limit theorems with rate of convergence for the number of germs with $d$ neighbours in a germ--grain model, and the number of degree-$d$ vertices in an Erdős--Rényi random graph. In both cases, the error rate is optimal, up to logarithmic factors.

preprint2016arXiv

General random walk in a random environment defined on Galton-Watson trees

We consider the motion of a particle on a Galton Watson tree, when the probabilities of jumping from a vertex to any one of its neighbours is determined by a random process. Given the tree, positive weights are assigned to the edges in such a way that, viewed along any line of descent, they evolve as a random process. In order to introduce our method for proving transience or recurrence, we first suppose that the weights are i.i.d., reproving a result of Lyons and Pemantle. We then extend the argument to allow a Markovian environment, and finally to a random walk on a Markovian environment that changes the environment. Our approach involves studying the typical behaviour of processes on fixed lines of descent, which we then show determines the behaviour of the process on the whole tree.

preprint2016arXiv

Multivariate approximation in total variation, I: equilibrium distributions of Markov jump processes

For integer valued random variables, the translated Poisson distributions form a flexible family for approximation in total variation, in much the same way that the normal family is used for approximation in Kolmogorov distance. Using the Stein--Chen method, approximation can often be achieved with error bounds of the same order as those for the CLT. In this paper, an analogous theory, again based on Stein's method, is developed in the multivariate context. The approximating family consists of the equilibrium distributions of a collection of Markov jump processes, whose analogues in one dimension are the immigration--death processes with Poisson distributions as equilibria. The method is illustrated by providing total variation error bounds for the approximation of the equilibrium distribution of one Markov jump process by that of another. In a companion paper, it is shown how to use the method for discrete normal approximation in ${\mathbb Z}^d$.

preprint2016arXiv

Multivariate approximation in total variation, II: discrete normal approximation

The paper applies the theory developed in Part I to the discrete normal approximation in total variation of random vectors in ${\mathbb Z}^d$. We illustrate the use of the method for sums of independent integer valued random vectors, and for random vectors exhibiting an exchangeable pair. We conclude with an application to random colourings of regular graphs.

preprint2016arXiv

On the emergence of random initial conditions in fluid limits

The paper presents a phenomenon occurring in population processes that start near zero and have large carrying capacity. By the classical result of Kurtz~(1970), such processes, normalized by the carrying capacity, converge on finite intervals to the solutions of ordinary differential equations, also known as the fluid limit. When the initial population is small relative to carrying capacity, this limit is trivial. Here we show that, viewed at suitably chosen times increasing to infinity, the process converges to the fluid limit, governed by the same dynamics, but with a random initial condition. This random initial condition is related to the martingale limit of an associated linear birth and death process.

preprint2015arXiv

Connecting deterministic and stochastic metapopulation models

In this paper, we study the relationship between certain stochastic and deterministic versions of Hanski's incidence function model and the spatially realistic Levins model. We show that the stochastic version can be well approximated in a certain sense by the deterministic version when the number of habitat patches is large, provided that the presence or absence of individuals in a given patch is influenced by a large number of other patches. Explicit bounds on the deviation between the stochastic and deterministic models are given.

preprint2015arXiv

Stein factors for negative binomial approximation in Wasserstein distance

The paper gives the bounds on the solutions to a Stein equation for the negative binomial distribution that are needed for approximation in terms of the Wasserstein metric. The proofs are probabilistic, and follow the approach introduced in Barbour and Xia (Bernoulli 12 (2006) 943-954). The bounds are used to quantify the accuracy of negative binomial approximation to parasite counts in hosts. Since the infectivity of a population can be expected to be proportional to its total parasite burden, the Wasserstein metric is the appropriate choice.

preprint2014arXiv

Escape from the boundary in Markov population processes

Density dependent Markov population processes in large populations of size $N$ were shown by Kurtz (1970, 1971) to be well approximated over finite time intervals by the solution of the differential equations that describe their average drift, and to exhibit stochastic fluctuations about this deterministic solution on the scale $\sqrt N$ that can be approximated by a diffusion process. Here, motivated by an example from evolutionary biology, we are concerned with describing how such a process leaves an absorbing boundary. Initially, one or more of the populations is of size much smaller than $N$, and the length of time taken until all populations have sizes comparable to $N$ then becomes infinite as $N \to \infty$. Under suitable assumptions, we show that in the early stages of development, up to the time when all populations have sizes at least $N^{1-α}$, for $1/3 < α< 1$, the process can be accurately approximated in total variation by a Markov branching process. Thereafter, the process is well approximated by the deterministic solution starting from the original initial point, but with a random time delay. Analogous behaviour is also established for a Markov process approaching an equilibrium on a boundary, where one or more of the populations become extinct.

preprint2014arXiv

Individual and patch behaviour in structured metapopulation models

Density dependent Markov population processes with countably many types can often be well approximated over finite time intervals by the solution of the differential equations that describe their average drift, provided that the total population size is large. They also exhibit diffusive stochastic fluctuations on a smaller scale about this deterministic path. Here, it is shown that the individuals in such processes experience an almost deterministic environment. Small groups of individuals behave almost independently of one another, evolving as Markov jump processes, whose transition rates are prescribed functions of time. In the context of metapopulation models, we show that `individuals' can represent either patches or the individuals that migrate among the patches; in host--parasite systems, they can represent both hosts and parasites.

preprint2013arXiv

Approximating the epidemic curve

Many models of epidemic spread have a common qualitative structure. The numbers of infected individuals during the initial stages of an epidemic can be well approximated by a branching process, after which the proportion of individuals that are susceptible follows a more or less deterministic course. In this paper, we show that both of these features are consequences of assuming a locally branching structure in the models, and that the deterministic course can itself be determined from the distribution of the limiting random variable associated with the backward, susceptibility branching process. Examples considered include a stochastic version of the Kermack & McKendrick model, the Reed-Frost model, and the Volz configuration model.

preprint2013arXiv

Asymptotic behaviour of gossip processes and small world networks

Both small world models of random networks with occasional long range connections and gossip processes with occasional long range transmission of information have similar characteristic behaviour. The long range elements appreciably reduce the effective distances, measured in space or in time, between pairs of typical points. In this paper, we show that their common behaviour can be interpreted as a product of the locally branching nature of the models. In particular, it is shown that both typical distances between points and the proportion of space that can be reached within a given distance or time can be approximated by formulae involving the limit random variable of the branching process.

preprint2012arXiv

Central limit approximations for Markov population processes with countably many types

When modelling metapopulation dynamics, the influence of a single patch on the metapopulation depends on the number of individuals in the patch. Since there is usually no obvious natural upper limit on the number of individuals in a patch, this leads to systems in which there are countably infinitely many possible types of entity. Analogous considerations apply in the transmission of parasitic diseases. In this paper, we prove central limit theorems for quite general systems of this kind, together with bounds on the rate of convergence in an appropriately chosen weighted $\ell_1$ norm.

preprint2011arXiv

A law of large numbers approximation for Markov population processes with countably many types

When modelling metapopulation dynamics, the influence of a single patch on the metapopulation depends on the number of individuals in the patch. Since the population size has no natural upper limit, this leads to systems in which there are countably infinitely many possible types of individual. Analogous considerations apply in the transmission of parasitic diseases. In this paper, we prove a law of large numbers for rather general systems of this kind, together with a rather sharp bound on the rate of convergence in an appropriately chosen weighted $\ell_1$ norm.

preprint2011arXiv

Total variation approximation for quasi-equilibrium distributions, II

Quasi-stationary distributions, as discussed by Darroch & Seneta (1965), have been used in biology to describe the steady state behaviour of population models which, while eventually certain to become extinct, nevertheless maintain an apparent stochastic equilibrium for long periods. These distributions have some drawbacks: they need not exist, nor be unique, and their calculation can present problems. In an earlier paper, we gave biologically plausible conditions under which the quasi-stationary distribution is unique, and can be closely approximated by distributions that are simple to compute. In this paper, we consider conditions under which the quasi-stationary distribution, if it exists, need not be unique, but an apparent stochastic equilibrium can nonetheless be identified and computed; we call such a distribution a quasi-equilibrium distribution.

preprint2010arXiv

Approximation by the Dickman distribution and quasi-logarithmic combinatorial structures

Quasi-logarithmic combinatorial structures are a class of decomposable combinatorial structures which extend the logarithmic class considered by Arratia, Barbour and Tavaré (2003). In order to obtain asymptotic approximations to their component spectrum, it is necessary first to establish an approximation to the sum of an associated sequence of independent random variables in terms of the Dickman distribution. This in turn requires an argument that refines the Mineka coupling by incorporating a blocking construction, leading to exponentially sharper coupling rates for the sums in question. Applications include distributional limit theorems for the size of the largest component and for the vector of counts of the small components in a quasi-logarithmic combinatorial structure.

preprint2010arXiv

Assessing molecular variability in cancer genomes

The dynamics of tumour evolution are not well understood. In this paper we provide a statistical framework for evaluating the molecular variation observed in different parts of a colorectal tumour. A multi-sample version of the Ewens Sampling Formula forms the basis for our modelling of the data, and we provide a simulation procedure for use in obtaining reference distributions for the statistics of interest. We also describe the large-sample asymptotics of the joint distributions of the variation observed in different parts of the tumour. While actual data should be evaluated with reference to the simulation procedure, the asymptotics serve to provide theoretical guidelines, for instance with reference to the choice of possible statistics.

preprint2010arXiv

Couplings for irregular combinatorial assemblies

When approximating the joint distribution of the component counts of a decomposable combinatorial structure that is `almost' in the logarithmic class, but nonetheless has irregular structure, it is useful to be able first to establish that the distribution of a certain sum of non-negative integer valued random variables is smooth. This distribution is not like the normal, and individual summands can contribute a non-trivial amount to the whole, so its smoothness is somewhat surprising. In this paper, we consider two coupling approaches to establishing the smoothness, and contrast the results that are obtained.

preprint2010arXiv

The shortest distance in random multi-type intersection graphs

Using an associated branching process as the basis of our approximation, we show that typical inter-point distances in a multitype random intersection graph have a defective distribution, which is well described by a mixture of translated and scaled Gumbel distributions, the missing mass corresponding to the event that the vertices are not in the same component of the graph.

preprint2010arXiv

Total variation approximation for quasi-equilibrium distributions

Quasi-stationary distributions, as discussed by Darroch & Seneta (1965), have been used in biology to describe the steady state behaviour of population models which, while eventually certain to become extinct, nevertheless maintain an apparent stochastic equilibrium for long periods. These distributions have some drawbacks: they need not exist, nor be unique, and their calculation can present problems. In this paper, we give biologically plausible conditions under which the quasi-stationary distribution is unique, and can be closely approximated by distributions that are simple to compute.