Journal of The Royal Society Interface
Restricted access Research articles

The construction of next-generation matrices for compartmental epidemic models

O. Diekmann

O. Diekmann

Department of Mathematics, Utrecht University, Budapestlaan 6, 3584 CD, Utrecht, The Netherlands

Google Scholar

Find this author on PubMed

,
J. A. P. Heesterbeek

J. A. P. Heesterbeek

Faculty of Veterinary Medicine, Utrecht University, Yalelaan 7, 3584 CL, Utrecht, The Netherlands

[email protected]

Google Scholar

Find this author on PubMed

and
M. G. Roberts

M. G. Roberts

Centre for Mathematical Biology, Institute of Information & Mathematical Sciences, Massey University, Private Bag 102 904, North Shore Mail Centre, Auckland, New Zealand

Google Scholar

Find this author on PubMed

Published:https://doi.org/10.1098/rsif.2009.0386

    The basic reproduction number ℛ0 is arguably the most important quantity in infectious disease epidemiology. The next-generation matrix (NGM) is the natural basis for the definition and calculation of ℛ0 where finitely many different categories of individuals are recognized. We clear up confusion that has been around in the literature concerning the construction of this matrix, specifically for the most frequently used so-called compartmental models. We present a detailed easy recipe for the construction of the NGM from basic ingredients derived directly from the specifications of the model. We show that two related matrices exist which we define to be the NGM with large domain and the NGM with small domain. The three matrices together reflect the range of possibilities encountered in the literature for the characterization of ℛ0. We show how they are connected and how their construction follows from the basic model ingredients, and establish that they have the same non-zero eigenvalues, the largest of which is the basic reproduction number ℛ0. Although we present formal recipes based on linear algebra, we encourage the construction of the NGM by way of direct epidemiological reasoning, using the clear interpretation of the elements of the NGM and of the model ingredients. We present a selection of examples as a practical guide to our methods. In the appendix we present an elementary but complete proof that ℛ0 defined as the dominant eigenvalue of the NGM for compartmental systems and the Malthusian parameter r, the real-time exponential growth rate in the early phase of an outbreak, are connected by the properties that ℛ0 > 1 if and only if r > 0, and ℛ0 = 1 if and only if r = 0.

    Footnotes

    All authors contributed equally to the study, the order is alphabetical.

    References