Stochastic modeling of the transmission of respiratory syncytial virus (RSV) in the region of Valencia,Spain

In this paper, we study the dynamics of the transmission of respiratory syncytial virus (RSV) in the population using stochastic models. The stochastic models are developed introducing stochastic perturbations on the demographic parameter as well as on the transmission rate of the RSV. Numerical simulations of the deterministic and stochastic models are performed in order to understand the effect of fluctuating birth rate and transmission rate of the RSV on the population dynamics. The numerical solutions of stochastic models are calculated using Euler-Maruyama and Milstein schemes, and confidence intervals for stochastic solutions are given using Monte-Carlo method. Analysis of the numerical results reveals that perturbations on the transmission rate are more decisive in the dynamics of RSV than perturbations on demographic parameters. In addition, the stochastic models show the advantage of reproducing more effectively the noisy RSV hospitalization data. It is concluded that these stochastic models are a viable option to provide a realistic modeling of the RSV dynamics on the population.     

Numerical integration of the master equation in some models of stochastic epidemiology

The processes by which disease spreads in a population of individuals are inherently stochastic. The master equation has proven to be a useful tool for modeling such processes. Unfortunately, solving the master equation analytically is possible only in limited cases (e.g., when the model is linear), and thus numerical procedures or approximation methods must be employed. Available approximation methods, such as the system size expansion method of van Kampen, may fail to provide reliable solutions, whereas current numerical approaches can induce appreciable computational cost. In this paper, we propose a new numerical technique for solving the master equation. Our method is based on a more informative stochastic process than the population process commonly used in the literature. By exploiting the structure of the master equation governing this process, we develop a novel technique for calculating the exact solution of the master equation--up to a desired precision--in certain models of stochastic epidemiology. We demonstrate the potential of our method by solving the master equation associated with the stochastic SIR epidemic model. MATLAB software that implements the methods discussed in this paper is freely available as Supporting Information S1.     

Numerical solution of structured population models

Numerical methods are presented for a general mass-structured population model with demographic rates that depend on individual mass, time, and total population mass. Several types of numerical methods are described, Eulerian methods, implicit methods, and the method of characteristics. These methods are compared for a sample problem with an exact solution. The preferred numerical technique combines the method of characteristics with an adaptive grid. Numerical solution of model equations developed for mosquitofish illustrates this method and demonstrates how seasons can play a dominant role in shaping population development.     

Equilibrium and extinction in stochastic population dynamics

Stochastic models of interacting biological populations, with birth and death rates depending on the population size are studied in the quasi-stationary state. Confidence regions in the state space are constructed by a new method for the numerical, solution of the ray equations. The concept of extinction time, which is closely related to the concept of stability for stochastic systems, is discussed. Results of numerical calculations for two-dimensional stochastic population models are presented.     

Persistence of structured populations in random environments

Environmental fluctuations often have different impacts on individuals that differ in size, age, or spatial location. To understand how population structure, environmental fluctuations, and density-dependent interactions influence population dynamics, we provide a general theory for persistence for density-dependent matrix models in random environments. For populations with compensating density dependence, exhibiting “bounded” dynamics, and living in a stationary environment, we show that persistence is determined by the stochastic growth rate (alternatively, dominant Lyapunov exponent) when the population is rare. If this stochastic growth rate is negative, then the total population abundance goes to zero with probability one. If this stochastic growth rate is positive, there is a unique positive stationary distribution. Provided there are initially some individuals in the population, the population converges in distribution to this stationary distribution and the empirical measures almost surely converge to the distribution of the stationary distribution. For models with overcompensating density-dependence, weaker results are proven. Methods to estimate stochastic growth rates are presented. To illustrate the utility of these results, applications to unstructured, spatially structured, and stage-structured population models are given. For instance, we show that diffusively coupled sink populations can persist provided that within patch fitness is sufficiently variable in time but not strongly correlated across space.     

Using statistics to design and estimate vital rates in matrix population models for a perennial herb

Matrix population models are widely used to assess population status and to inform management decisions. Despite existing theories for building such models, model construction is often partially based on expert opinion. So far, model structure has received relatively little attention, although it may affect estimates of population dynamics. Here, we assessed the consequences of two published matrix structures (a 4 × 4 matrix based on expert opinion and a 10 × 10 matrix based on statistical modeling) for estimates of vital rates and stochastic population dynamics of the long-lived herb Astragalus scaphoides. We explored the ways in which choice of model structure alters the accuracy (i.e., mean) and precision (i.e., variance) of predicted population dynamics. We found that model structure had a negligible effect on the accuracy and precision of vital rates and stochastic stage distribution. However, the 10 × 10 matrix produced lower estimates of stochastic population growth rates than the 4 × 4 matrix, and more accurately predicted the observed trends in population abundance for three out of four study populations. Moreover, estimates of realized variation in population growth rate due to fluctuations in population stage structure over time were occasionally sensitive to matrix structure, suggesting differential roles of transient dynamics. Our study indicates that statistical modeling for choosing categories in matrix models might be preferable over expert opinion to accurately predict population trends and can provide a more objective way for model construction when the biological knowledge of the species is limited.     

A discrete,stochastic model of colonial phytoplankton population size structure: Development and application to in situ imaging-in-flow cytometer observations of Dinobryon

The scientific community lacks models for the dynamic changes in population size structure that occur in colonial phytoplankton. This is surprising, as size is a key trait affecting many aspects of phytoplankton ecology, and colonial forms are very common. We aim to fill this gap with a new discrete, stochastic model of dynamic changes in phytoplankton colonies' population size structure. We use the colonial phytoplankton Dinobryon as a proof-of-concept organism. The model includes four stochastic functions—division, stomatocyst production, colony breakage, and colony loss—to determine Dinobryon population size structure and populations counts. Although the functions presented here are tailored to Dinobryon, the model is readily adaptable to represent other colonial taxa. We demonstrate how fitting our model to in situ observations of colony population size structure can provide a powerful approach to explore colony size dynamics. Here, we have (1) collected high-frequency in situ observations of Dinobryon in Lac (Lake) Montjoie (Quebec, Canada) in 2013 with a moored Imaging FlowCytobot (IFCB) and (2) fit the model to those observations with a genetic algorithm solver that extracts parameter estimates for each of the four stochastic functions. As an example of the power of this model-data integration, we also highlight ecological insights into Dinobryon colony size and stomatocyst production. The Dinobryon population was enriched in larger, flagellate-rich colonies near bloom initiation and shifted to smaller and emptier colonies toward bloom decline.     

The Modeling of Global Epidemics: Stochastic Dynamics and Predictability

The global spread of emergent diseases is inevitably entangled with the structure of the population flows among different geographical regions. The airline transportation network in particular shrinks the geographical space by reducing travel time between the world's most populated areas and defines the main channels along which emergent diseases will spread. In this paper, we investigate the role of the large-scale properties of the airline transportation network in determining the global propagation pattern of emerging diseases. We put forward a stochastic computational framework for the modeling of the global spreading of infectious diseases that takes advantage of the complete International Air Transport Association 2002 database complemented with census population data. The model is analyzed by using for the first time an information theory approach that allows the quantitative characterization of the heterogeneity level and the predictability of the spreading pattern in presence of stochastic fluctuations. In particular we are able to assess the reliability of numerical forecast with respect to the intrinsic stochastic nature of the disease transmission and travel flows. The epidemic pattern predictability is quantitatively determined and traced back to the occurrence of epidemic pathways defining a backbone of dominant connections for the disease spreading. The presented results provide a general computational framework for the analysis of containment policies and risk forecast of global epidemic outbreaks. On leave from CEA-Centre d'Etudes de Bruyères-Le-Chatel, France.     

On the comparative static properties of the expected population density in the presence of stochastic fluctuations

We consider the impact of increased stochastic fluctuations on the expected density of an unstructured population evolving according to a regular diffusion process subject to a concave expected growth rate. By relying on the flow nature of the solutions of stochastic differential equations and Girsanovs theorem, we demonstrate that typically increased volatility decreases the expected future population density. As a consequence, we are able to characterize the sensitivity of the expected population density with respect to changes in the diffusion coefficient measuring the size of the stochastic fluctuations. We provide both qualitative and quantitative information about the consequences of a mis-specified volatility structure and, especially, of a deterministic approximation to stochastic population growth. We also consider the effect of uncertainty in the initial density and demonstrate that the sign of the relationship between the expected population density and initial uncertainty is unambiguosly negative. Received: 15 February 1999 / Revised version: 29 September 1999 / Published online: 5 May 2000     

Linking demographic responses and life history tactics from longitudinal data in mammals

In stochastic environments, a change in a demographic parameter can influence the population growth rate directly or via a resulting impact on age structure. Stochastic elasticity of the long‐run stochastic growth rate λ_s to a demographic parameter offers a suitable way to measure the overall demographic response because it includes both the direct effect of changing the demographic parameter and its indirect effect through changes in the age structure. From 25 mammalian populations with contrasting life histories, we investigated how pace of life and population growth rate influence the demographic responses (measured as the relative contributions of the direct and indirect components of stochastic elasticity on λ_s). We found that in short‐lived species, the change in population structure resulting from an increase in yearling survival leads to an additional increase in λ_s, whereas in long‐lived species, the same change in population structure leads to a decrease. Short‐lived species thus display a boom‐bust life history strategy contrary to long‐lived species, for which the long lifespan dampens the demographic consequences of changing age structure. Irrespective of the species’ life history strategy, the change in population age structure resulting from an increase in adult survival leads to an additional increase in λ_s due to an increase of the proportion of mature individuals in the population. On the contrary, a change in population age structure resulting from an increase of reproductive performance leads to a decrease in λ_s that is due to the increase of the proportion of immature individuals in the population. Our comparative analysis of stochastic elasticity patterns in mammals shows the existence of different demographic responses to changes in age structure between short‐ and long‐lived species, which improves our understanding of population dynamics in variable environments in relation to the species‐specific pace of life.     

Spread rate for a nonlinear stochastic invasion

Despite the recognized importance of stochastic factors, models for ecological invasions are almost exclusively formulated using deterministic equations [29]. Stochastic factors relevant to invasions can be either extrinsic (quantities such as temperature or habitat quality which vary randomly in time and space and are external to the population itself) or intrinsic (arising from a finite population of individuals each reproducing, dying, and interacting with other individuals in a probabilistic manner). It has been long conjectured [27] that intrinsic stochastic factors associated with interacting individuals can slow the spread of a population or disease, even in a uniform environment. While this conjecture has been borne out by numerical simulations, we are not aware of a thorough analytical investigation. In this paper we analyze the effect of intrinsic stochastic factors when individuals interact locally over small neighborhoods. We formulate a set of equations describing the dynamics of spatial moments of the population. Although the full equations cannot be expressed in closed form, a mixture of a moment closure and comparison methods can be used to derive upper and lower bounds for the expected density of individuals. Analysis of the upper solution gives a bound on the rate of spread of the stochastic invasion process which lies strictly below the rate of spread for the deterministic model. The slow spread is most evident when invaders occur in widely spaced high density foci. In this case spatial correlations between individuals mean that density dependent effects are significant even when expected population densities are low. Finally, we propose a heuristic formula for estimating the true rate of spread for the full nonlinear stochastic process based on a scaling argument for moments. Received: 19 October 1998 / Revised version: 1 September 1999 / Published online: 4 October 2000     

A Stochastic Tick-Borne Disease Model: Exploring the Probability of Pathogen Persistence

We formulate and analyse a stochastic epidemic model for the transmission dynamics of a tick-borne disease in a single population using a continuous-time Markov chain approach. The stochastic model is based on an existing deterministic metapopulation tick-borne disease model. We compare the disease dynamics of the deterministic and stochastic models in order to determine the effect of randomness in tick-borne disease dynamics. The probability of disease extinction and that of a major outbreak are computed and approximated using the multitype Galton–Watson branching process and numerical simulations, respectively. Analytical and numerical results show some significant differences in model predictions between the stochastic and deterministic models. In particular, we find that a disease outbreak is more likely if the disease is introduced by infected deer as opposed to infected ticks. These insights demonstrate the importance of host movement in the expansion of tick-borne diseases into new geographic areas.     

Does environmental stochasticity matter? Analysis of red deer life-histories on Rum

Summary Most life-history theory assumes that short-term variation in an organism's environment does not affect the survivorships and fecundities of the organisms. This assumption is rarely met. Here we investigate the population and evolutionary biology of red deer,Cervus elephas, to see if relaxation of this assumption is likely to make significant differences to the predicted evolutionary biology of this species. To do this we used 21 years of data from a population of deer on Rum, Western Isles, Scotland. Population growth rates in a stochastic environment were estimated using Tuljapurkar's small noise approximation, confirmed by bootstrap simulation. Numerical differentiation was used to see if the selection pressures (i.e. sensitivities of population growth rate to changes in the vital rates) differ between the stochastic and deterministic cases. The data also allow the costs of reproduction to be estimated. These costs, incorporated as trade-offs into the sensitivity analysis, allow investigation of evolutionary benefits of different life-history tactics. Environmentally induced stochastic variation in the red deer vital rates causes a slight reduction ( 1%) in the predicted population growth rate and has little impact on the estimated selection pressures on the deer's life-history. We thus conclude that, even though density-independent stochastic effects on the population are marked, the deer's fitness is not markedly affected by these and they are adapted to the average conditions they experience. However, the selected life-history is sensitive to the trade-offs between current fecundity, survivorship and future fecundity and it is likely that the environmental variance will affect these trade-offs and, thus, affect the life-history favoured by selection. We also show that the current average life-history is non-optimal and suggest this is a result of selection pressures exerted by culling and predation, now much reduced. As the use of stochastic or deterministic methods provide similar estimates in this case, the use of the latter is justified. Thus,r (the annual per capita rate of population growth) is an appropriate measure of fitness in a population with stochastic numerical fluctuations. In a population of constant size lifetime reproductive success is the obvious measure of fitness to use.     

Stochastic epidemics: major outbreaks and the duration of the endemic period

A study is made of a two-dimensional stochastic system that models the spread of an infectious disease in a population. An asymptotic expression is derived for the probability that a major outbreak of the disease will occur in case the number of infectives is small. For the case that a major outbreak has occurred, an asymptotic approximation is derived for the expected time that the disease is in the population. The analytical expressions are obtained by asymptotically solving Dirichlet problems based on the Fokker-Planck equation for the stochastic system. Results of numerical calculations for the analytical expressions are compared with simulation results.     

A unifying framework for chaos and stochastic stability in discrete population models

 In this paper we propose a general framework for discrete time one-dimensional Markov population models which is based on two fundamental premises in population dynamics. We show that this framework incorporates both earlier population models, like the Ricker and Hassell models, and experimental observations concerning the structure of density dependence. The two fundamental premises of population dynamics are sufficient to guarantee that the model will exhibit chaotic behaviour for high values of the natural growth and the density-dependent feedback, and this observation is independent of the particular structure of the model. We also study these models when the environment of the population varies stochastically and address the question under what conditions we can find an invariant probability distribution for the population under consideration. The sufficient conditions for this stochastic stability that we derive are of some interest, since studying certain statistical characteristics of these stochastic population processes may only be possible if the process converges to such an invariant distribution. Received 15 May 1995; received in revised form 17 April 1996     

The rate of multi-step evolution in Moran and Wright-Fisher populations

Several groups have recently modeled evolutionary transitions from an ancestral allele to a beneficial allele separated by one or more intervening mutants. The beneficial allele can become fixed if a succession of intermediate mutants are fixed or alternatively if successive mutants arise while the previous intermediate mutant is still segregating. This latter process has been termed stochastic tunneling. Previous work has focused on the Moran model of population genetics. I use elementary methods of analyzing stochastic processes to derive the probability of tunneling in the limit of large population size for both Moran and Wright-Fisher populations. I also show how to efficiently obtain numerical results for finite populations. These results show that the probability of stochastic tunneling is twice as large under the Wright-Fisher model as it is under the Moran model.     

Study of an infant brain subjected to periodic motion via a custom experimental apparatus design and finite element modelling

This paper presents a rig that was specifically designed to simulate the shaking of mechanical models of biological systems, especially those related to shaken baby syndrome (SBS). The scope of this paper includes the testing of an anthropomorphic model that simulates an infant head and provides validation data for complex finite element (FE) modelling using three numerical methods (Lagrangian, Arbitrary-Lagrangian–Eulerian (ALE) and Eulerian method) for fluid structure coupling.The experiments for this study aim to provide an understanding of the influence of two factors on intracranial brain movement of the infant head during violent shaking: (1) the specific paediatric head structure: the anterior fontanelle and (2) the brain–skull interface.The results show that the Eulerian analysis method has significant advantages for the FSI modelling of brain–CSF–skull interactions over the more commonly used methods, e.g. the Lagrangian method. To the knowledge of the authors, this methodology has not been discussed in previous publication.The results indicate that the biomechanical investigation of SBS can provide more accurate results only if the skull with paediatric features and the brain–skull interface are correctly represented, which were overlooked in previous SBS studies.     

Food-signal theory: population regulation and the functional response

A fully stochastic food-signal model, a function of a patchy-preyencounter sequence, and a prey-processing function is described.The model shows how prey density and its second-order statisticalproperties can sequester prey from predators, questioning theuse of only numerical abundance of predator and prey organismsas a measure of preya — predator interactions. The modelhighlights the notion that patch structure can be generatedby relative velocity of predator and prey as well as by theirspatial distribution. The model extends ideas that include the‘biological pump’ and the downwelling of carbonfrom the upper ocean, the functional response, optimal-foragingtheory, and the connections between population dynamics andvariability in the physical environment.     

Perspectives in stochastic models of human reproduction: A review and analysis

The purpose of this paper is to review, to analyze, and to take steps toward synthesizing, two research areas in stochastic models of population dynamics. The first of these areas consists of stochastic models of human reproduction associated with the late Mindel C. Sheps and the second area is a class of stochastic models of population growth called generalized age-dependent branching processes, a class of models that shows promise of throwing more light on the classical mathematical demography of Lotka. The substantive material of the paper is arranged in eight sections ranging in content from a comparison of classical mathematical demography with generalized age-dependent branching processes, to suggestions for restructuring models of the Sheps school in quest of greater realism. The paper ends with a section on numerical examples illustrating applications in family planning evaluation and an appendix suggesting ways in which algebraic concepts are useful in short cutting computations.