Stability analysis and optimal vaccination of an SIR epidemic model

Almost all mathematical models of diseases start from the same basic premise: the population can be subdivided into a set of distinct classes dependent upon experience with respect to the relevant disease. Most of these models classify individuals as either a susceptible individual S, infected individual I or recovered individual R. This is called the susceptible-infected-recovered (SIR) model. In this paper, we describe an SIR epidemic model with three components; S, I and R. We describe our study of stability analysis theory to find the equilibria for the model. Next in order to achieve control of the disease, we consider a control problem relative to the SIR model. A percentage of the susceptible populations is vaccinated in this model. We show that an optimal control exists for the control problem and describe numerical simulations using the Runge-Kutta fourth order procedure. Finally, we describe a real example showing the efficiency of this optimal control.     

The role of sexually abstained groups in two-sex demographic and epidemic logistic models with non-linear mortality

We describe several gender structured population models governed by logistic growth with non-linear death rate. We extend these models to include groups of people isolated from sexual activity and individuals exposed to a mild and long-lasting sexually transmitted disease, i.e. without disease-induced mortality and recovery. The transmission of the disease is modeled through formation/separation of heterosexual couples assuming that one infected individual automatically infects his/her partner. We are interested in how the non-reproductive class may change the demographic tendencies in the general population and whether they can curb the growth of the infected group while keeping the healthy one at acceptable levels. A comparison of the equilibrium total population size in the presence and the absence of the isolated class is also provided.     

Optimal harvesting and optimal vaccination

Two optimization problems are considered: Harvesting from a structured population with maximal gain subject to the condition of non-extinction, and vaccinating a population with prescribed reduction of the reproduction number of the disease at minimal costs. It is shown that these problems have a similar structure and can be treated by the same mathematical approach. The optimal solutions have a 'two-window' structure: Optimal harvesting and vaccination strategies or policies are concentrated on one or two preferred age classes. The results are first shown for a linear age structure problem and for an epidemic situation at the uninfected state (minimize costs for a given reduction of the reproduction number) and then extended to populations structured by size, to harvesting at Gurtin-MacCamy equilibria and to vaccination at infected equilibria.     

Estimation of effective vaccination rate: pertussis in New Zealand as a case study

In some cases vaccination is unreliable. For example vaccination against pertussis has comparatively high level of primary and secondary failures. To evaluate efficiency of vaccination we introduce the idea of effective vaccination rate and suggest an approach to estimate it. We consider pertussis in New Zealand as a case study. The results indicate that the level of immunity failure for pertussis is considerably higher than was anticipated.     

Progression age enhanced backward bifurcation in an epidemic model with super-infection

 We consider a model for a disease with a progressing and a quiescent exposed class and variable susceptibility to super-infection. The model exhibits backward bifurcations under certain conditions, which allow for both stable and unstable endemic states when the basic reproduction number is smaller than one. Received: 11 October 2001 / Revised version: 17 September 2002 / Published online: 17 January 2003 Present address: Department of Biological Statistics and Computational Biology, 434 Warren Hall, Cornell University, Ithaca, NY 14853-7801 This author was visiting Arizona State University when most of the research was done. Research partially supported by NSF grant DMS-0137687. This author's research was partially supported by NSF grant DMS-9706787. Key words or phrases: Backward bifurcation – Multiple endemic equilibria – Alternating stability – Break-point density – Super-infection – Dose-dependent latent period – Progressive and quiescent latent stages – Progression age structure – Threshold type disease activation – Operator semigroups – Hille-Yosida operators – Dynamical systems – Persistence – Global compact attractor     

Mathematical analysis of a cholera model with public health interventions

Cholera, an acute gastro-intestinal infection and a waterborne disease continues to emerge in developing countries and remains an important global health challenge. We formulate a mathematical model that captures some essential dynamics of cholera transmission to study the impact of public health educational campaigns, vaccination and treatment as control strategies in curtailing the disease. The education-induced, vaccination-induced and treatment-induced reproductive numbers R(E), R(V), R(T) respectively and the combined reproductive number R(C) are compared with the basic reproduction number R(0) to assess the possible community benefits of these control measures. A Lyapunov functional approach is also used to analyse the stability of the equilibrium points. We perform sensitivity analysis on the key parameters that drive the disease dynamics in order to determine their relative importance to disease transmission and prevalence. Graphical representations are provided to qualitatively support the analytical results.     

Mathematical analysis of the global dynamics of a model for HTLV-I infection and ATL progression

Mathematical analysis is carried out that completely determines the global dynamics of a mathematical model for the transmission of human T-cell lymphotropic virus I (HTLV-I) infection and the development of adult T-cell leukemia (ATL). HTLV-I infection of healthy CD4(+) T cells takes place through cell-to-cell contact with infected T cells. The infected T cells can remain latent and harbor virus for several years before virus production occurs. Actively infected T cells can infect other T cells and can convert to ATL cells, whose growth is assumed to follow a classical logistic growth function. Our analysis establishes that the global dynamics of T cells are completely determined by a basic reproduction number R(0). If R(0)< or =1, infected T cells always die out. If R(0)>1, HTLV-I infection becomes chronic, and a unique endemic equilibrium is globally stable in the interior of the feasible region. We also show that the equilibrium level of ATL-cell proliferation is higher when the HTLV-I infection of T cells is chronic than when it is acute.     

A theoretical study on mathematical modelling of an infectious disease with application of optimal control

In this paper, we propose and analyze an epidemic problem which can be controlled by vaccination as well as treatment. In the first part of our analysis we study the dynamical behavior of the system with fixed control for both vaccination and treatment. Basic reproduction number is obtained in all possible cases and it is observed that the simultaneous use of vaccination and treatment control is the most favorable case to prevent the disease from being epidemic. In the second part, we take the controls as time dependent and obtain the optimal control strategy to minimize both the infected populations and the associated costs. All the analytical results are verified by simulation works. Some important conclusions are given at the end of the paper.     

Cluster formation for multi-strain infections with cross-immunity

Many infectious diseases exist in several pathogenic variants, or strains, which interact via cross-immunity. It is observed that strains tend to self-organise into groups, or clusters. The aim of this paper is to investigate cluster formation. Computations demonstrate that clustering is independent of the model used, and is an intrinsic feature of the strain system itself. We observe that an ordered strain system, if it is sufficiently complex, admits several cluster structures of different types. Appearance of a particular cluster structure depends on levels of cross-immunity and, in some cases, on initial conditions. Clusters, once formed, are stable, and behave remarkably regularly (in contrast to the generally chaotic behaviour of the strains themselves). In general, clustering is a type of self-organisation having many features in common with pattern formation.     

Oscillations in a patchy environment disease model

For a single patch SIRS model with a period of immunity of fixed length, recruitment-death demographics, disease related deaths and mass action incidence, the basic reproduction number R(0) is identified. It is shown that the disease-free equilibrium is globally asymptotically stable if R(0)<1. For R(0)>1, local stability of the endemic equilibrium and Hopf bifurcation analysis about this equilibrium are carried out. Moreover, a practical numerical approach to locate the bifurcation values for a characteristic equation with delay-dependent coefficients is provided. For a two patch SIRS model with travel, it is shown that there are several threshold quantities determining its dynamic behavior and that travel can reduce oscillations in both patches; travel may enhance oscillations in both patches; or travel can switch oscillations from one patch to another.     

Stability analysis and optimal control of an SIR epidemic model with vaccination

This paper focuses on the study of a nonlinear mathematical SIR epidemic model with a vaccination program. We have discussed the existence and the stability of both the disease free and endemic equilibrium. Vaccine induced reproduction number is determined and the impact of vaccination in reducing the vaccine induced reproduction number is discussed. Then to achieve control of the disease, a control problem is formulated and it is shown that an optimal control exists for our model. The optimality system is derived and solved numerically using the Runge-Kutta fourth order procedure.     

Global stability of an SIR epidemic model with information dependent vaccination

We study the global behavior of a non-linear susceptible-infectious-removed (SIR)-like epidemic model with a non-bilinear feedback mechanism, which describes the influence of information, and of information-related delays, on a vaccination campaign. We upgrade the stability analysis performed in d’Onofrio et al. [A. d’Onofrio, P. Manfredi, E. Salinelli, Vaccinating behavior, information, and the dynamics of SIR vaccine preventable diseases, Theor. Popul. Biol. 71 (2007) 301] and, at same time, give a special example of application of the geometric method for global stability, due to Li and Muldowney. Numerical investigations are provided to show how the stability properties depend on the interplay between some relevant parameters of the model.     

Sustainability and optimal control of an exploited prey predator system through provision of alternative food to predator

In the present paper, we develop a simple two species prey-predator model in which the predator is partially coupled with alternative prey. The aim is to study the consequences of providing additional food to the predator as well as the effects of harvesting efforts applied to both the species. It is observed that the provision of alternative food to predator is not always beneficial to the system. A complete picture of the long run dynamics of the system is discussed based on the effort pair as control parameters. Optimal augmentations of prey and predator biomass at final time have been investigated by optimal control theory. Also the short and large time effects of the application of optimal control have been discussed. Finally, some numerical illustrations are given to verify our analytical results with the help of different sets of parameters.     

An SIRVS epidemic model with pulse vaccination strategy

The aim of this paper is to analyze an SIRVS epidemic model in which pulse vaccination strategy (PVS) is included. We are interested in finding the basic reproductive number of the model which determine whether or not the disease dies out. The global attractivity of the disease-free periodic solution (DFPS for short) is obtained when the basic reproductive number is less than unity. The disease is permanent when the basic reproductive number is greater than unity, i.e., the epidemic will turn out to endemic. Our results indicate that the disease will go to extinction when the vaccination rate reaches some critical value.     

An age-structured epidemic model of rotavirus with vaccination

The recent approval of a rotavirus vaccine in Mexico motivates this study on the potential impact of the use of such a vaccine on rotavirus prevention and control. An age-structured model that describes the rotavirus transmission dynamics of infections is introduced. Conditions that guarantee the local and global stability analysis of the disease-free steady state distribution as well as the existence of an endemic steady state distribution are established. The impact of maternal antibodies on the implementation of vaccine is evaluated. Model results are used to identify optimal age-dependent vaccination strategies. A convergent numerical scheme for the model is introduced but not implemented. This paper is dedicated to Prof. K. P. Hadeler, who continues to push the frontier of knowledge in mathematical biology.     

A mathematical model for indirectly transmitted diseases

We consider a mathematical model for the indirect transmission via a contaminated environment of a microparasite between two spatially distributed host populations having non-coincident spatial domains. The parasite is benign in a first population and lethal in the second one. Global existence results are given for the resulting reaction-diffusion system coupled with an ordinary differential equation. Then, invasion and persistence of the parasite are studied. A simplified model for the transmission of a hantavirus from bank vole to human populations is then analysed.     

Pulse vaccination strategy in the SIR epidemic model

Theoretical results show that the measles ‘pulse’ vaccination strategy can be distinguished from the conventional strategies in leading to disease eradication at relatively low values of vaccination. Using the SIR epidemic model we showed that under a planned pulse vaccination regime the system converges to a stable solution with the number of infectious individuals equal to zero. We showed that pulse vaccination leads to epidemics eradication if certain conditions regarding the magnitude of vaccination proportion and on the period of the pulses are adhered to. Our theoretical results are confirmed by numerical simulations. The introduction of seasonal variation into the basic SIR model leads to periodic and chaotic dynamics of epidemics. We showed that under seasonal variation, in spite of the complex dynamics of the system, pulse vaccination still leads to epidemic eradication. We derived the conditions for epidemic eradication under various constraints and showed their dependence on the parameters of the epidemic. We compared effectiveness and cost of constant, pulse and mixed vaccination policies.     

Stability analysis for HIV infection delay model with protease inhibitor

In this article, we considered a model of HIV-1 infection with a protease inhibitor therapy and three delays. The frequency of the bifurcating periodic solution as well as the threshold value is approximated numerically using realistic parameter. The estimated threshold value is realistic and the frequency of the oscillations is consistent with that of the observed viral blips.     

Moment closure and the stochastic logistic model

The quasi-stationary distribution of the stochastic logistic model is studied in the parameter region where its body is approximately normal. Improved asymptotic approximations of its first three cumulants are derived. It is shown that the same results can be derived with the aid of the moment closure method. This indicates that the moment closure method leads to expressions for the cumulants that are asymptotic approximations of the cumulants of the quasi-stationary distribution.     

A discrete-time model with vaccination for a measles epidemic.

A discrete-time, age-independent SIR-type epidemic model is formulated and analyzed. The effects of vaccination are also included in the model. Three mathematically important properties are verified for the model: solutions are nonnegative, the population size is time-invariant, and the epidemic concludes with all individuals either remaining susceptible or becoming immune (a property typical of SIR models). The model is applied to a measles epidemic on a university campus. The simulated results are in good agreement with the actual data if it is assumed that the population mixes nonhomogeneously. The results of the simulations indicate that a rate of immunity greater than 98% may be required to prevent an epidemic in a university population. The model has applications to other contagious diseases of SIR type. Furthermore, the simulated results of the model can easily be compared to data, and the effects of a vaccination program can be examined.