Kermack–McKendrick theory

Kermack–McKendrick theory is a mathematical hypothesis about how infectious diseases spread through a population. Building on the research of Ronald Ross and others, A. G. McKendrick and W. O. Kermack published their theory in a set of three articles from 1927, 1932, and 1933. While Kermack—McKendrick theory was indeed the source of SIR models and their relatives, Kermack and McKendrick were thinking of a more subtle and empirically useful problem than the simple compartmental models discussed here. The text is somewhat difficult to read, compared to modern papers, but the important feature is it was a model where the age-of-infection affected the transmission and removal rates.

Because of their seminal importance to the field of theoretical epidemiology, these articles were republished in the Bulletin of Mathematical Biology in 1991.^[1][2][3]

Epidemic model (1927)

In its initial form, Kermack—McKendrick theory is a compartmental differential-equation model that structures the infectioned population in terms of age-of-infection, while using simple compartments for people who are susceptible (S) and recovered/removed (R). Specified initial conditions would change over time according to

$\frac{dS}{dt} = - \lambda S,$

$\frac{di}{dt} + \frac{di}{da} = \delta(a) \lambda S- \gamma(a) i,$

$\frac{dR}{dt} = \int_{0}^\infty \gamma(a) i(a,t) \, da ,$

where $\delta(a)$ is a Dirac delta-function and the infection pressure

$\lambda = \int_0^\infty \beta(a) i(a,t) \, da .$

Only in the special case when the removal rate $\gamma(a)$ and the transmission rate $\beta(a)$ are constant for all ages does the substitution $I(t) = \int_0^\infty i(a,t) \, da$ transform their theory into the simple SIR model. This basic model only accounts for infection and removal events, which are sufficient to describe a simple epidemic, including the threshold condition necessary for an epidemic to start, but can not explain endemic disease transmission or recurring epidemics.

Endemic disease (1932, 1933)

In their subsequent articles, Kermack and McKendrick extended their theory to allow for birth, migration, and death, as well as imperfect immunity. In modern notation, their model can be represented as

$\frac{dS}{dt} = b_0 + b_S S + b_I I + b_R R - \lambda S - m_S S,$

$\frac{di}{dt} + \frac{di}{da} = \delta(a) \lambda ( S + \sigma R ) - \gamma(a) i - \mu(a) i - m_i(a) i,$

$I(t) = \int_0^\infty i(a,t) \, da$

$\frac{dR}{dt} = \int_0^\infty \gamma(a) i(a,t) \, da - \sigma \lambda R - m_R R,$

where $b_0$ is the immigration rate of susceptibles, b_j is the per-capita birth rate for state j, m_j is the per-capita mortality rate of individuals in state j, $\sigma$ is the relative-risk of infection to recovered individuals who are partially immune, and the infection pressure

$\lambda = \int_0^\infty \beta(a) i(a,t) \, da .$

Kermack and McKendrick were able to show that it admits an stationary solution where disease is endemic, as long as the supply of susceptible individuals is sufficiently large. This model is difficult to analyze in its full generality, and a number of open questions remain regarding its dynamics.

References

1. Kermack, W; McKendrick, A (1991). "Contributions to the mathematical theory of epidemics—I". Bulletin of Mathematical Biology 53 (1–2): 33–55. doi:10.1016/S0092-8240(05)80040-0. PMID 2059741.
2. Kermack, W; McKendrick, A (1991). "Contributions to the mathematical theory of epidemics—II. The problem of endemicity". Bulletin of Mathematical Biology 53 (1–2): 57–87. doi:10.1016/S0092-8240(05)80041-2. PMID 2059742.
3. Kermack, W; McKendrick, A (1991). "Contributions to the mathematical theory of epidemics—III. Further studies of the problem of endemicity". Bulletin of Mathematical Biology 53 (1–2): 89–118. doi:10.1016/S0092-8240(05)80042-4. PMID 2059743.