Some properties of an estimator for the basic reproduction number of the general epidemic model.

In this paper some properties of a convenient estimator, derived from a martingale estimating function, for the basic reproduction number of the general epidemic model are given for both finite and large samples. These properties give some guidelines for using this convenient estimator. It is shown that it underestimates the parameter and that the bias tends to zero when the population size and the initial number of infectives are increased simultaneously. The bias cannot be removed for a fixed number of introductory infectives. However, the estimator is asymptotically unbiased, conditional on a major outbreak. A simulation study shows that the central limit theorem applies for moderate population sizes.     

A stochastic SIS epidemic with demography: initial stages and time to extinction

We study an open population stochastic epidemic model from the time of introduction of the disease, through a possible outbreak and to extinction. The model describes an SIS (susceptible–infective–susceptible) epidemic where all individuals, including infectious ones, reproduce at a given rate. An approximate expression for the outbreak probability is derived using a coupling argument. Further, we analyse the behaviour of the model close to quasi-stationarity, and the time to disease extinction, with the aid of a diffusion approximation. In this situation the number of susceptibles and infectives behaves as an Ornstein–Uhlenbeck process, centred around the stationary point, for an exponentially distributed time before going extinct.     

Analysis of a stochastic SIR epidemic on a random network incorporating household structure

This paper is concerned with a stochastic SIR (susceptible → infective → removed) model for the spread of an epidemic amongst a population of individuals, with a random network of social contacts, that is also partitioned into households. The behaviour of the model as the population size tends to infinity in an appropriate fashion is investigated. A threshold parameter which determines whether or not an epidemic with few initial infectives can become established and lead to a major outbreak is obtained, as are the probability that a major outbreak occurs and the expected proportion of the population that are ultimately infected by such an outbreak, together with methods for calculating these quantities. Monte Carlo simulations demonstrate that these asymptotic quantities accurately reflect the behaviour of finite populations, even for only moderately sized finite populations. The model is compared and contrasted with related models previously studied in the literature. The effects of the amount of clustering present in the overall population structure and the infectious period distribution on the outcomes of the model are also explored.     

Effects of behavioral changes in a smallpox attack model

The impact of individual and community behavioral changes in response to an outbreak of a disease with high mortality is often not appreciated. Response strategies to a smallpox bioterrorist attack have focused on interventions such as isolation of infectives, contact tracing, quarantine of contacts, ring vaccination, and mass vaccination. We formulate and analyze a mathematical model in which some individuals lower their daily contact activity rates once an epidemic has been identified in a community. Transmission parameters are estimated from data and an expression is derived for the effective reproduction number. We use computer simulations to analyze the effects of behavior change alone and in combination with other control measures. We demonstrate that the spread of the disease is highly sensitive to how rapidly people reduce their contact activity rates and to the precautions that the population takes to reduce the transmission of the disease. Even gradual and mild behavioral changes can have a dramatic impact in slowing an epidemic. When behavioral changes are combined with other interventions, the epidemic is shortened and the number of smallpox cases is reduced. We conclude that for simulations of a smallpox outbreak to be useful, they must consider the impact of behavioral changes. This is especially true if the model predictions are being used to guide public health policy.     

Dynamic population epidemic models.

Most multipopulation epidemic models are of the contact distribution type, in which the locations of successive contacts are chosen independently from appropriate contact distributions. This paper is concerned with an alternative class of models, termed dynamic population epidemic models, in which infectives move among the populations and can infect only within their current population. Both the stochastic and deterministic versions of such models are considered. Their threshold behavior is analyzed in some depth, as are their final outcomes. Velocities of spread of infection are considered when the populations have a spatial structure. A criterion for finding the equivalent contact distribution epidemic for any given dynamic population epidemic is provided, enabling comparisons to be made for the velocities and final outcomes displayed by the two classes of models. The relationship between deterministic and stochastic epidemic models is also discussed briefly.     

A two-stage model for the SIR outbreak: accounting for the discrete and stochastic nature of the epidemic at the initial contamination stage

The evolution of an infectious disease outbreak in an isolated population is split into two stages: a stochastic Markov process describing the initial contamination and a linked deterministic dynamical system with random initial conditions for the continued development of the outbreak. The initial contamination stage is well approximated by the randomized SI (susceptible/infected) model. We obtain the probability density function for the early behavior of the epidemic. This provides an appropriate distribution for the initial conditions with which to describe the subsequent deterministic evolution of the system. We apply the method of matching asymptotic expansions to link the two stages. This allows us to estimate the standard deviation of the number of infectives in the developed outbreak, and the statistical characteristics of the outbreak time. The potential trajectories caused by the stochastic nature of the contamination stage show greatest divergence at the initial and fade-out stages and coincide most tightly just after the peak of the epidemic. The time to the peak of the outbreak is not strongly dependent on the initial trajectory.     

Epidemic patch models applied to pandemic influenza: Contact matrix, stochasticity, robustness of predictions

Due to the recent emergence of H5N1 virus, the modelling of pandemic influenza has become a relevant issue. Here we present an SEIR model formulated to simulate a possible outbreak in Italy, analysing its structure and, more generally, the effect of including specific details into a model. These details regard population heterogeneities, such as age and spatial distribution, as well as stochasticity, that regulates the epidemic dynamics when the number of infectives is low. We discuss and motivate the specific modelling choices made when building the model and investigate how the model details influence the predicted dynamics. Our analysis may help in deciding which elements of complexity are worth including in the design of a deterministic model for pandemic influenza, in a balance between, on the one hand, keeping the model computationally efficient and the number of parameters low and, on the other hand, maintaining the necessary realistic features.     

The basic reproduction number and the probability of extinction for a dynamic epidemic model

We consider the spread of an epidemic through a population divided into n sub-populations, in which individuals move between populations according to a Markov transition matrix Σ and infectives can only make infectious contacts with members of their current population. Expressions for the basic reproduction number, R₀, and the probability of extinction of the epidemic are derived. It is shown that in contrast to contact distribution models, the distribution of the infectious period effects both the basic reproduction number and the probability of extinction of the epidemic in the limit as the total population size N → ∞. The interactions between the infectious period distribution and the transition matrix Σ mean that it is not possible to draw general conclusions about the effects on R₀ and the probability of extinction. However, it is shown that for n = 2, the basic reproduction number, R₀, is maximised by a constant length infectious period and is decreasing in ?, the speed of movement between the two populations.     

Asymptotic behavior in a deterministic epidemic model

The effects of a periodic contact rate and of carriers are considered for a generalization of Bailey's simple epidemic model. In this model it is assumed that individuals become susceptible again as soon as they recover from the infection so that a fixed population can be divided into a class of infectives and a class of susceptibles which vary with time. If the contact rate is periodic, then the number of infectives as time approaches infinity either tends to zero or is asymptotically periodic depending on whether the total population size is less than or greater than a threshold value. The behavior for large time of the number of infectives is determined for three modifications of the model which involve carriers.     

Effect of a sharp change of the incidence function on the dynamics of a simple disease

We investigate two cases of a sharp change of incidencec functions on the dynamics of a susceptible-infective-susceptible epidemic model. In the first case, low population levels have mass action incidence, while high population levels have proportional incidence, the switch occurring when the total population reaches a certain threshold. Using a modified Dulac theorem, we prove that this system has a single equilibrium which attracts all solutions for which the disease is present and the population remains bounded. In the second case, an increase of the number of infectives leads to a mass action term being added to a standard incidence term. We show that this allows a Hopf bifurcation to occur, with periodic orbits being generated when a locally asymptotically stable equilibrium loses stability.     

Saturation recovery leads to multiple endemic equilibria and backward bifurcation

The number of patients need to be treated may exceed the carry capacity of local hospitals during the spreading of a severe infectious disease. We propose an epidemic model with saturation recovery from infective individuals to understand the effect of limited resources for treatment of infectives on the emergency disease control. It is shown that saturation recovery from infective individuals leads to vital dynamics, such as bistability and periodicity, when the basic reproduction number R₀ is less than unity. An interesting dynamical behavior of the model is a backward bifurcation which raises many new challenges to effective infection control.     

Spatial spread of an epidemic through public transportation systems with a hub

This article investigates an epidemic spreading among several locations through a transportation system, with a hub connecting these locations. Public transportation is not only a bridge through which infections travel from one location to another but also a place where infections occur since individuals are typically in close proximity to each other due to the limited space in these systems. A mathematical model is constructed to study the spread of an infectious disease through such systems. A variant of the next generation method is proposed and used to provide upper and lower bounds of the basic reproduction number for the model. Our investigation indicates that increasing transportation efficiency, and improving sanitation and ventilation of the public transportation system decrease the chance of an outbreak occurring. Moreover, discouraging unnecessary travel during an epidemic also decreases the chance of an outbreak. However, reducing travel by infectives while allowing susceptibles to travel may not be enough to avoid an outbreak.     

Transition from growth to decay of an epidemic due to lockdown

We study the transition of an epidemic from growth phase to decay of the active infections in a population when lockdown health measures are introduced to reduce the probability of disease transmission. Although in the case of uniform lockdown, a simple compartmental model would indicate instantaneous transition to decay of the epidemic, this is not the case when partially isolated active clusters remain with the potential to create a series of small outbreaks. We model this using the Gillespie stochastic simulation algorithm based on a connected set of stochastic susceptible-infected-removed/recovered networks representing the locked-down majority population (in which the reproduction number is less than 1) weakly coupled to a large set of small clusters in which the infection may propagate. We find that the presence of such active clusters can lead to slower than expected decay of the epidemic and significantly delayed onset of the decay phase. We study the relative contributions of these changes, caused by the active clusters within the population, to the additional total infected population. We also demonstrate that limiting the size of the inevitable active clusters can be efficient in reducing their impact on the overall size of the epidemic outbreak. The deceleration of the decay phase becomes apparent when the active clusters form at least 5% of the population.     

Backward bifurcation of an epidemic model with treatment

An epidemic model with a limited resource for treatment is proposed to understand the effect of the capacity for treatment. It is assumed that treatment rate is proportional to the number of infectives below the capacity and is a constant when the number of infectives is greater than the capacity. It is found that a backward bifurcation occurs if the capacity is small. It is also found that there exist bistable endemic equilibria if the capacity is low.     

A waterborne outbreak of epidemic diarrhea due to group A rotavirus in Malatya, Turkey

We characterized an outbreak of acute diarrheal disease caused by group A rotavirus that occurred during the Autumn of 2005 in Malatya City, Turkey. A total 9907 patients between 0 to 91 years old (mean age: 25.05?19.67) were included in the epidemic. The patients? data were prospectively collected and statistically analyzed. Microbiologic analyses were performed to determine the etiologic agent. Rapid onset diarrhea (98.36%), abdominal cramps (69%), fever (44.4%) and vomiting (69.6%) were the most common symptoms observed in patients. Rotavirus antigen was detected in 52.7% of the studied patients. RT-PCR analysis led to identification of Group A rotavirus as the causative agent of this epidemic. Simultaneous measurements of the drinking water samples yielded very low chlorine levels; as low as 0 to 0.05 mg/L. The outbreak investigation team indicated possible contamination of a large water depository from a water well, which supplies drinking water to two major districts of the city. Effective chlorination and blockage of the passage between the well and the water depository stopped the outbreak. This outbreak shows the high epidemic potency of rotavirus in large human populations, including all age groups, and underlines the importance of water safety in pipeline systems.     

A Stochastic Epidemic Model with Seasonal Variations in Infection Rate

In this paper we consider a modification of Bailey's stochastic model for the spread of an epidemic when there are seasonal variations in infection rate. The resulting nonlinear model is analyzed by employing the diffusion approximation technique. We have shown that for a large population the process, on suitable scaling and normalization, converges to a non-stationary Ornstein-Uhlenbeck process. Consequently the number of infectives has in the steady state a gaussian distribution.     

Resource allocation for epidemic control in metapopulations

Deployment of limited resources is an issue of major importance for decision-making in crisis events. This is especially true for large-scale outbreaks of infectious diseases. Little is known when it comes to identifying the most efficient way of deploying scarce resources for control when disease outbreaks occur in different but interconnected regions. The policy maker is frequently faced with the challenge of optimizing efficiency (e.g. minimizing the burden of infection) while accounting for social equity (e.g. equal opportunity for infected individuals to access treatment). For a large range of diseases described by a simple SIRS model, we consider strategies that should be used to minimize the discounted number of infected individuals during the course of an epidemic. We show that when faced with the dilemma of choosing between socially equitable and purely efficient strategies, the choice of the control strategy should be informed by key measurable epidemiological factors such as the basic reproductive number and the efficiency of the treatment measure. Our model provides new insights for policy makers in the optimal deployment of limited resources for control in the event of epidemic outbreaks at the landscape scale.     

Sexual dimorphism and resource allocation in male and female shrubs of Simmondsia chinensis

Summary Desert populations of the evergreen dioecious shrub Simmondsia chinensis exhibit sex-related leaf and canopy dimorphisms not present in populations from more mesic coastal environments. Leaves on female shrubs have characteristically larger sizes, greater specific weights, and greater water-holding capacity than male leaves in desert habitats. In coastal scrub environments no significant difference is present, with leaf characteristics of both sexes similar to those of desert male shrubs. Desert female shrub canopies are typically relatively open with little mutual branch shading. In male shrubs canopies are more densely branched with considerable mutual shading of branches. Female plants allocate a greater proportion of their vegetative resources to leaves than do male plants. Considering total biomass, male plants allocate 10–15% of their resources (biomass, calories, glucose-equivalents, nitrogen, phosphorus) to reproductive tissues. Female allocation is dependent on seed set. At 100% seed set females would allocate 30–40% of their resources to reproduction, while female reproductive investment would equal that of males at approximately 30% seed set. Sexual dimorphism and the associated physiological characteristics in Simmondsia act as an alternative to differential habitat selection by male and female plants. Female plants respond to limited water resources in desert areas by increasing their efficiency in allocating limited resources to reproductive structures.     

Modelling the effects of pre-exposure and post-exposure vaccines in tuberculosis control

Epidemic control strategies alter the spread of the disease in the host population. In this paper, we describe and discuss mathematical models that can be used to explore the potential of pre-exposure and post-exposure vaccines currently under development in the control of tuberculosis. A model with bacille Calmette-Guerin (BCG) vaccination for the susceptibles and treatment for the infectives is first presented. The epidemic thresholds known as the basic reproduction numbers and equilibria for the models are determined and stabilities are investigated. The reproduction numbers for the models are compared to assess the impact of the vaccines currently under development. The centre manifold theory is used to show the existence of backward bifurcation when the associated reproduction number is less than unity and that the unique endemic equilibrium is locally asymptotically stable when the associated reproduction number is greater than unity. From the study we conclude that the pre-exposure vaccine currently under development coupled with chemoprophylaxis for the latently infected and treatment of infectives is more effective when compared to the post-exposure vaccine currently under development for the latently infected coupled with treatment of the infectives.     

Micro-economics based resource allocation in Grid-Federation environment

A Grid-Federation environment is composed of a collection of autonomous and selfish distributed cluster resource managers. These selfish managers participate in Grid-Federation to share their resources. Market models could be used to motivate the self-interested participants to share their resources. In this paper, firstly, a market for resource exchange in grid federation environment is established. Then, in order that the market reaches a Walrasian equilibrium, a computationally tractable mechanism is proposed. A Walrasian equilibrium problem consists of finding a set of prices and allocations of resources in such a way that the cluster resource managers could maximize their utilities and the market clears. Market clears when the resource supply equals to the demand. We show that in a Walrasian equilibrium, the Grid Federation market reaches an efficient resource allocation.