Mathematical Medicine and Biology Advance Access originally published online on April 28, 2006
Mathematical Medicine and Biology 2006 23(3):231-254; doi:10.1093/imammb/dql011
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Bifurcations of an epidemic model with non-linear incidence and infection-dependent removal rate
Institute for Biodiagnostics, National Research Council Canada, Winnipeg, Manitoba, Canada R3B 1Y6
** Email: Murray.Alexander{at}nrc-cnrc.gc.ca
An epidemic model with a generalized non-linear incidence is extended to incorporate the effect of an infection-dependent removal strategy, which is defined as a function of the number of infected individuals. It is assumed that the removal rate decreases from a maximum capacity for removing infected individuals as their number increases. The existence and stability of the associated equilibria are analysed, and the basic reproductive number (
[graphic: see PDF]
) is formulated. It is shown that
[graphic: see PDF]
is independent of the functional form of the incidence, but depends on the removal rate. Normal forms are derived to show the different types of bifurcation the model undergoes, including transcritical, generalized Hopf (Bautin), saddle-node and BogdanovTakens. A degenerate Hopf bifurcation at the Bautin point, where the first Lyapunov coefficient vanishes, is discussed. Sotomayor's theorem is applied to establish a saddle-node bifurcation at the turning point of backward bifurcation. The BogdanovTakens normal form is derived, from which the local bifurcation curve for a family of homoclinic orbits is formulated. Bifurcation diagrams and numerical simulations, using parameter values estimated for some infectious diseases, are also presented to provide more intuition to the theoretical findings. The results show that sufficiently increasing the removal rate can reduce
[graphic: see PDF]
below a subthreshold domain, which leads to disease eradication.
Keywords: epidemic models; non-linear incidence; removal rate; Hopf bifurcation; saddle-node bifurcation; BogdanovTakens bifurcation
Received on 31 May 2005. revised on 24 March 2006. accepted on 31 March 2006.