An SIS epidemic model with individual variation.

We study an extension of the stochastic SIS (Susceptible-Infectious-Susceptible) model in continuous time that accounts for variation amongst individuals. By examining its limiting behaviour as the population size grows we are able to exhibit conditions for the infection to become endemic.

A Model for Cell Proliferation in a Developing Organism.

In mathematical biology, there is a great deal of interest in producing continuum models by scaling discrete agent-based models governed by local stochastic rules. We discuss a particular example of this approach: a model for the proliferation of neural crest cells that can help us understand the development of Hirschprung's disease, a potentially-fatal condition in which the enteric nervous system of a new-born child does not extend all the way through the intestine and colon. Our starting point is a discrete-state, continuous-time Markov chain model proposed by Hywood et al. (2013a) for the location of the neural crest cells that make up the enteric nervous system. Hywood et al. (2013a) scaled their model to derive an approximate second order partial differential equation describing how the limiting expected number of neural crest cells evolve in space and time. In contrast, we exploit the relationship between the above-mentioned Markov chain model and the well-known Yule-Furry process to derive the exact form of the scaled version of the process. Furthermore, we provide expressions for other features of the domain agent occupancy process, such as the variance of the marginal occupancy at a particular site, the distribution of the number of agents that are yet to reach a given site and a stochastic description of the process itself.

A model for a spatially structured metapopulation accounting for within patch dynamics.

We develop a stochastic metapopulation model that accounts for spatial structure as well as within patch dynamics. Using a deterministic approximation derived from a functional law of large numbers, we develop conditions for extinction and persistence of the metapopulation in terms of the birth, death and migration parameters. Interestingly, we observe the Allee effect in a metapopulation comprising two patches of greatly different sizes, despite there being decreasing patch specific per-capita birth rates. We show that the Allee effect is due to the way the migration rates depend on the population density of the patches.

Development and testing of a genetic marker-based pedigree reconstruction system 'PR-genie' incorporating size-class data.

For wildlife populations, it is often difficult to determine biological parameters that indicate breeding patterns and population mixing, but knowledge of these parameters is essential for effective management. A pedigree encodes the relationship between individuals and can provide insight into the dynamics of a population over its recent history. Here, we present a method for the reconstruction of pedigrees for wild populations of animals that live long enough to breed multiple times over their lifetime and that have complex or unknown generational structures. Reconstruction was based on microsatellite genotype data along with ancillary biological information: sex and observed body size class as an indicator of relative age of individuals within the population. Using body size-class data to infer relative age has not been considered previously in wildlife genealogy and provides a marked improvement in accuracy of pedigree reconstruction. Body size-class data are particularly useful for wild populations because it is much easier to collect noninvasively than absolute age data. This new pedigree reconstruction system, PR-genie, performs reconstruction using maximum likelihood with optimization driven by the cross-entropy method. We demonstrated pedigree reconstruction performance on simulated populations (comparing reconstructed pedigrees to known true pedigrees) over a wide range of population parameters and under assortative and intergenerational mating schema. Reconstruction accuracy increased with the presence of size-class data and as the amount and quality of genetic data increased. We provide recommendations as to the amount and quality of data necessary to provide insight into detailed familial relationships in a wildlife population using this pedigree reconstruction technique.

Chain binomial models and binomial autoregressive processes.

We establish a connection between a class of chain-binomial models of use in ecology and epidemiology and binomial autoregressive (AR) processes. New results are obtained for the latter, including expressions for the lag-conditional distribution and related quantities. We focus on two types of chain-binomial model, extinction-colonization and colonization-extinction models, and present two approaches to parameter estimation. The asymptotic distributions of the resulting estimators are studied, as well as their finite-sample performance, and we give an application to real data. A connection is made with standard AR models, which also has implications for parameter estimation.

Quantitative risk stratification in Markov chains with limiting conditional distributions.

BACKGROUND: Many clinical decisions require patient risk stratification. The authors introduce the concept of limiting conditional distributions, which describe the equilibrium proportion of surviving patients occupying each disease state in a Markov chain with death. Such distributions can quantitatively describe risk stratification. METHODS: The authors first establish conditions for the existence of a positive limiting conditional distribution in a general Markov chain and describe a framework for risk stratification using the limiting conditional distribution. They then apply their framework to a clinical example of a treatment indicated for high-risk patients, first to infer the risk of patients selected for treatment in clinical trials and then to predict the outcomes of expanding treatment to other populations of risk. RESULTS: For the general chain, a positive limiting conditional distribution exists only if patients in the earliest state have the lowest combined risk of progression or death. The authors show that in their general framework, outcomes and population risk are interchangeable. For the clinical example, they estimate that previous clinical trials have selected the upper quintile of patient risk for this treatment, but they also show that expanded treatment would weakly dominate this degree of targeted treatment, and universal treatment may be cost-effective. CONCLUSIONS: Limiting conditional distributions exist in most Markov models of progressive diseases and are well suited to represent risk stratification quantitatively. This framework can characterize patient risk in clinical trials and predict outcomes for other populations of risk.

Metapopulation persistence in a dynamic landscape: more habitat or better stewardship?

Habitat loss and fragmentation has created metapopulations where there were once continuous populations. Ecologists and conservation biologists have become interested in the optimal way to manage and conserve such metapopulations. Several authors have considered the effect of patch disturbance and recovery on metapopulation persistence, but almost all such studies assume that every patch is equally susceptible to disturbance. We investigated the influence of protecting patches from disturbance on metapopulation persistence, and used a stochastic metapopulation model to answer the question: How can we optimally trade off returns from protection of patches vs. creation of patches? We considered the problem of finding, under budgetary constraints, the optimal combination of increasing the number of patches in the metapopulation network vs. increasing the number of protected patches in the network. We discovered that the optimal trade-off is dependent upon all of the properties of the system: the species dynamics, the dynamics of the landscape, and the relative costs of each action. A stochastic model and accompanying methodology are provided allowing a manager to determine the optimal policy for small metapopulations. We also provide two approximations, including a rule of thumb, for determining the optimal policy for larger metapopulations. The method is illustrated with an example inspired by information for the greater bilby, Macrotis lagotis, inhabiting southwestern Queensland, Australia. We found that given realistic costs for each action, protection of patches should be prioritized over patch creation for improving the persistence of the greater bilby during the next 20 years.