Stochastic dynamical behavior of SIRS epidemic models with random perturbation

In this paper, we consider a stochastic SIRS model with parameter perturbation, which is a standard technique in modeling population dynamics. In our model, the disease transmission coefficient and the removal rates are all affected by noise. We show that the stochastic model has a unique positive solution as is essential in any population model. Then we establish conditions for extinction or persistence of the infectious disease. When the infective part is forced to expire, the susceptible part converges weakly to an inverse-gamma distribution with explicit shape and scale parameters. In case of persistence, by new stochastic Lyapunov functions, we show the ergodic property and positive recurrence of the stochastic model. We also derive the an estimate for the mean of the stationary distribution. The analytical results are all verified by computer simulations, including examples based on experiments in laboratory populations of mice.

[1]  Ingemar Nåsell,et al.  Stochastic models of some endemic infections. , 2002, Mathematical biosciences.

[2]  Yasuhiro Takeuchi,et al.  Global stability of an SIR epidemic model with time delays , 1995, Journal of mathematical biology.

[3]  R. Khasminskii.  Stochastic Stability of Differential Equations , 1980 .

[4]  R. May,et al.  Population biology of infectious diseases: Part II , 1979, Nature.

[5]  Wanbiao Ma,et al.  Permanence of an SIR epidemic model with distributed time delays , 2001 .

[6]  S. Levin,et al.  Dynamical behavior of epidemiological models with nonlinear incidence rates , 1987, Journal of mathematical biology.

[7]  Liangjian Hu,et al.  A Stochastic Differential Equation SIS Epidemic Model , 2011, SIAM J. Appl. Math..

[8]  Xuerong Mao,et al.  Stochastic differential equations and their applications , 1997 .

[9]  M. Pascual,et al.  Stochastic amplification in epidemics , 2007, Journal of The Royal Society Interface.

[10]  Rui Xu,et al.  Global stability of a delayed SEIRS epidemic model with saturation incidence rate , 2010 .

[11]  Yasuhiro Takeuchi,et al.  Global asymptotic stability of an SIR epidemic model with distributed time delay , 2001 .

[12]  Pasquale Vetro,et al.  Stability of a stochastic SIR system , 2005 .

[13]  Y. Kutoyants.  Statistical Inference for Ergodic Diffusion Processes , 2004 .

[14]  Sebastian Walcher,et al.  Exclusion and persistence in deterministic and stochastic chemostat models , 2005 .

[15]  Yasuhiro Takeuchi,et al.  Convergence results in SIR epidemic models with varying population sizes , 1997 .

[16]  R. May,et al.  Population Biology of Infectious Diseases , 1982, Dahlem Workshop Reports.

[17]  Qiuying Lu,et al.  Stability of SIRS system with random perturbations , 2009 .

[18]  Yasuhiro Takeuchi,et al.  Permanence of a delayed SIR epidemic model with density dependent birth rate , 2007 .

[19]  N. Ling.  The Mathematical Theory of Infectious Diseases and its applications , 1978 .

[20]  H. P. Hudson,et al.  An application of the theory of probabilities to the study of a priori pathometry.—Part I , 1917 .

[21]  Graeme C. Wake,et al.  Lyapunov functions and global stability for SIR, SIRS, and SIS epidemiological models , 2002, Appl. Math. Lett..

[22]  W. O. Kermack,et al.  A contribution to the mathematical theory of epidemics , 1927 .

[23]  A J McKane,et al.  Predator-prey cycles from resonant amplification of demographic stochasticity. , 2005, Physical review letters.

[24]  X. Mao,et al.  Stochastic Differential Equations and Applications , 1998 .

[25]  Gang George Yin,et al.  Asymptotic Properties of Hybrid Diffusion Systems , 2007, SIAM J. Control. Optim..

[26]  R. May,et al.  Infectious Diseases of Humans: Dynamics and Control , 1991, Annals of Internal Medicine.

[27]  P. Driessche,et al.  A disease transmission model in a nonconstant population , 1993, Journal of mathematical biology.

[28]  Linda J. S. Allen,et al.  Construction of Equivalent Stochastic Differential Equation Models , 2008 .

[29]  Priscilla E. Greenwood,et al.  Sustained oscillations for density dependent Markov processes , 2011, Journal of mathematical biology.

[30]  Mei Song,et al.  Global stability of an SIR epidemicmodel with time delay , 2004, Appl. Math. Lett..

[31]  Herbert W. Hethcote,et al.  An epidemiological model with a delay and a nonlinear incidence rate , 1989, Journal of mathematical biology.

[32]  Y. Chow.  Local Convergence of Martingales and the Law of Large Numbers , 1965 .

[33]  H. Hethcote,et al.  Integral equation models for endemic infectious diseases , 1980, Journal of mathematical biology.

[34]  G. Serio,et al.  A generalization of the Kermack-McKendrick deterministic epidemic model☆ , 1978 .

[35]  Alberto d’Onofrio,et al.  Bifurcation Thresholds in an SIR Model with Information-Dependent Vaccination , 2007 .

[36]  Guoting Chen,et al.  STABILITY OF STOCHASTIC DELAYED SIR MODEL , 2009 .

[37]  K. Cooke,et al.  Vertically Transmitted Diseases: Models and Dynamics , 1993 .

[38]  Y. Iwasa,et al.  Influence of nonlinear incidence rates upon the behavior of SIRS epidemiological models , 1986, Journal of mathematical biology.

[39]  Gopal K. Basak,et al.  Stability in Distribution for a Class of Singular Diffusions , 1992 .