S-I-R Model of EpidemicsPart 1Basic Model and ExamplesRevised September 22, 20051. IntroductionüDescription of the ModelIn this notebook, we develop in detail the standard S-I-R model for epidemics. For the integration of the nonlinear differential equations, we use the package DynPac. Although some familiarity with DynPac is assumed, brief descriptions of some DynPac commands are given as they are used. The integrations in this notebook also can be done easily with the Mathematica function NDSolve.Much of our presentation is tied to the specific example given in section 3, an influenza epidemic in a British Boarding School, and for that example we follow the treatment given by J.D. Murray (Mathematical Biology I, An Introduction, p. 325-326, Springer-Verlag, 2002). The S-I-R model was introduced by W.O. Kermack and A.G. McKendrick ("A Contribution to the Mathematical Theory of Epidemics," Proc. Roy. Soc. London A 115, 700-721, 1927), and has played a major role in mathematical epidemiology. A summary of the model and its uses is given by Murray. In the model, a population is divided into three groups: the susceptibles S, the infectives I, and the recovered R, with numbers s, i, and r respectively. The total population is n = s + i + r . (1)The susceptibles are those who are not infected and not immune, the infectives are those who are infected and can transmit the disease, and the recovered are those who have been infected and are immune (recovered or dead). In the simple form of the model considered here, we assume that the recovered are permanently immune. We also assume that the whole event is of sufficiently short duration that we may ignore natural births and deaths during the epidemic. Finally, we ignore any subdivisions of the population by age, sex, mobility, or other factors, although such distinctions are obviously of importance. üDifferential Equations and Parameters of the ModelThe basic differential equations are as follows: dsÅÅÅÅÅÅÅÅdt= -asi , H2Lsir1.nb 1diÅÅÅÅÅÅÅÅdt= asi - bi , H3LdrÅÅÅÅÅÅÅÅdt= bi , H4LThese equations describe the transitions of individuals from S to I to R. By adding the three equations, we show easily that the total population n is constant. We will call the parameter a the transmissivity and the parameter b the recovery rate. The term asi is a standard kinetic term, based on the idea that the number of encounters per unit time between susceptibles and infectives will be proportional to the numbers of each. The transmissivity a is determined by both the encounter frequency and the efficiency with which the disease is transmitted per encounter. In later notebooks, we will look more closely at the problem of estimating a. We will see later that there are advantages to formulating the problem in terms of population fractions rather than total population. For example, it is obvious that if a given population is placed into a much larger area, the encounter rate will decrease and a will decrease. A proper formulation in terms of fractions will give a transmissivity independent of total population in the first approximation. For the present, we continue with the formulation in terms of total populations. We will not be concerned further in this notebook with the estimation of a from first principles. The interpretation and estimation of b is more straightforward. We now show that b-1is the average duration of the infection. Consider the special case of no new infections, so that equation (3) reduces to di/dt = -bi. Then i(t) = i(0)e- bt. The number of infectives recovering in time (t, t + dt) is |di| = |(di/dt)dt | = i(t)bdt = i(0)be- bt dt, and this is the number of infectives with infection durations in the range (t, t + dt). The average duration tdur is then the average of t with this as a weighting function, hencetdur=Ÿ0¶ i H0L be-bt t dtÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅŸ0¶ i H0L be- bt dt=1ÅÅÅÅÅÅb. H5LWe may treat the three equations given by (2) - (4) as a set of two equations in s and i (equations (2) and (3)). After these equations are solved for s and i, we may calculate r from equation (1). üConditions for an EpidemicIn epidemiology, the word "epidemic" has a technical meaning: it is a situation in which the number of infectives increases from the initial value. The condition for an epidemic follows directly from equation (3): di/dt will be positive whenever s > b/a . (6) Thus for given a and b there is a critical number of susceptibles for an epidemic, and if s0 is the initial number of susceptibles, then the condition for an epidemic is s0> b/a. There is a simple interpretation of equation (6). A given infective will, on the average, be infectious for a time 1/b. The number of susceptibles infected by one infective per unit time is as. Hence the total number of infections produced by one infective is as/b, and for an event to qualify as an epidemic, this number must exceed 1, which gives condition (6). The parameter appearing here, R0= as0/b , (7)is called the reproduction ratio of the epidemic.It is somewhat surprising that the number of initial infectives does not play a role in determining whether there is an epidemic. Of course other aspects of the epidemic will depend on the number of initial infectives -- especially the time of the peak of the epidemic. The peak of the epidemic is the time of maximum number of infectives. We see from equation (3) that this will occur when s = b/a -- that is, when the reproduction ratio is 1. sir1.nb 2It is somewhat surprising that the number of initial infectives does not play a role in determining whether there is an epidemic. Of course other aspects of the epidemic will depend on the number of initial infectives -- especially the time of the peak of the epidemic. The peak of the epidemic is the time of maximum number of infectives. We see from equation (3) that this will occur when s = b/a -- that is, when the reproduction ratio is 1. üEquilibrium and StabilityIt is easy to show that the equilibrium states of equations (2) - (4) are those states in which i = 0, and s = s*, where s* is any positive constant. Thus the equilibria are non-isolated. A formal stability analysis by linearization yields two eigenvalues: 0, and as* - b. If the condition for an epidemic is satisfied, this second