Skip to main content
Advertisement
  • Loading metrics

Estimation of country-level basic reproductive ratios for novel Coronavirus (SARS-CoV-2/COVID-19) using synthetic contact matrices

  • Joe Hilton ,

    Roles Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing

    j.hilton@warwick.ac.uk

    Affiliations School of Life Sciences, University of Warwick, Coventry, United Kingdom, Zeeman Institue (SBIDER), University of Warwick, Coventry, United Kingdom, Mathematics Institute, University of Warwick, Coventry, United Kingdom

  • Matt J. Keeling

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing

    Affiliations School of Life Sciences, University of Warwick, Coventry, United Kingdom, Zeeman Institue (SBIDER), University of Warwick, Coventry, United Kingdom, Mathematics Institute, University of Warwick, Coventry, United Kingdom

Abstract

The 2019-2020 pandemic of atypical pneumonia (COVID-19) caused by the virus SARS-CoV-2 has spread globally and has the potential to infect large numbers of people in every country. Estimating the country-specific basic reproductive ratio is a vital first step in public-health planning. The basic reproductive ratio (R0) is determined by both the nature of pathogen and the network of human contacts through which the disease can spread, which is itself dependent on population age structure and household composition. Here we introduce a transmission model combining age-stratified contact frequencies with age-dependent susceptibility, probability of clinical symptoms, and transmission from asymptomatic (or mild) cases, which we use to estimate the country-specific basic reproductive ratio of COVID-19 for 152 countries. Using early outbreak data from China and a synthetic contact matrix, we estimate an age-stratified transmission structure which can then be extrapolated to 151 other countries for which synthetic contact matrices also exist. This defines a set of country-specific transmission structures from which we can calculate the basic reproductive ratio for each country. Our predicted R0 is critically sensitive to the intensity of transmission from asymptomatic cases; with low asymptomatic transmission the highest values are predicted across Eastern Europe and Japan and the lowest across Africa, Central America and South-Western Asia. This pattern is largely driven by the ratio of children to older adults in each country and the observed propensity of clinical cases in the elderly. If asymptomatic cases have comparable transmission to detected cases, the pattern is reversed. Our results demonstrate the importance of age-specific heterogeneities going beyond contact structure to the spread of COVID-19. These heterogeneities give COVID-19 the capacity to spread particularly quickly in countries with older populations, and that intensive control measures are likely to be necessary to impede its progress in these countries.

Author summary

Over 100 countries have reported laboratory-confirmed cases of atypical pneumonia caused by 2019 novel coronavirus (COVID-19). Cases are largely reported in older age groups, suggesting a strong age-dependent component to either transmission or the probability of developing symptoms and thus being detected. We introduce a mathematical model for COVID-19 transmission in which contact behaviour, susceptibility, detection probability, and transmission from undetected cases all vary with age. We fit our model to epidemiological data from the outbreak in China for the special case where asymptomatic transmission is negligible, and compare it to a null model where only contact behaviour varies with age. Our fitted model suggests that contacts involving older individuals are particularly likely to generate new detected cases, intensifying the spread of infection in countries with older populations. We estimate the basic reproductive ratio (a measure of a pathogen’s capacity for spread) of COVID-19 in 152 countries under both models, and find that estimates of the basic reproductive ratio are highly dependent on the assumed underlying transmission structure; our more complex model predicts higher values in Japan and much of Europe and lower values in much of Africa, in comparison to the contact frequency-based model where this pattern is reversed.

Introduction

The ongoing epidemic of atpyical pneumonia (COVID-19) emerged in Wuhan, China in late 2019, with the novel coronavirus SARS-CoV-2 identified as its causative agent [1]. It exhibited an early capacity for global spread [2], and, as of May 29th 2020, it is present on all inhabited continents and has spread to 216 countries and territories [3]. A surprising feature of the outbreak in all countries is that the distribution of detected cases is characterised by large numbers of cases in older individuals and fewer in younger individuals, with particularly low numbers in under-15s [1, 4]. Unfortunately, it remains unclear what role these younger age-groups play in onward transmission.

The basic reproductive ratio, commonly denoted R0, is a key epidemiological variable in any outbreak. It is defined as the expected number of infectious cases generated by a single average case in an entirely susceptible population [5, 6]. It is therefore an important measure of an infection’s capacity to spread, with sustained transmission possible only when R0 > 1, and can be used to calculate both the expected size of an outbreak and the threshold level of vaccination necessary for eradication [5]. Early estimates of R0 for COVID-19 in China have mostly clustered between 2 and 4 [1, 2, 711]. However, these estimated values of R0 may not necessarily reflect the intensity of spread in countries other than China; this is because the basic reproductive ratio is dependent on both pathogen and host population characteristics [12].

Age-structured infectious disease models have been hugely successful in describing the dynamics of a number of infectious diseases, including endemic infections and novel outbreaks [1317]. These models typically use a so-called “who acquires infection from whom” paradigm, with transmission rates between individuals in different age classes encoded either in a single matrix or in a set of matrices specific to transmission events in different locations [5, 13]. It is often assumed that age-structured transmission rates are directly proportional to the frequency and duration of age-structured contacts, measured from social contact studies [15, 16, 18]. Such studies typically find that contact intensities are highest among younger (school-)age groups [1822], so that if age-structured transmission is determined entirely by contact frequency, we expect infection to be concentrated in these younger age groups. The relatively low levels of clinically presenting COVID-19 cases in children and teenagers therefore suggests that transmission is driven by more complex age-stratified heterogeneities beyond social contact patterns.

Although studies of social contacts are increasing in number [18, 2027], they tend to be time consuming to conduct and therefore have been limited to a number of exemplar countries. Prem et al. (2017) extrapolated these contact matrices to 152 countries based on a range of social and demographic data [19]. This generates a publicly available set of country-specific contact matrices whose entries correspond to the expected total number of age-stratified contacts per day for individuals aggregated into 5-year age classes [19]. One study has already used these estimated contact patterns to model COVID-19 in China, although the potential for other age-specific heterogeneities was not addressed [28].

In this paper, we define a model of the early age-structured transmission dynamics in China which incorporates the country-specific social contact matrix (estimated by Prem et al.[19]) as well as components defining age-specific susceptibility profiles, probabilities of developing symptoms, and transmission rates from subclinical cases. Our model formulation is sufficient flexible to account for a range of different assumptions relating to the drivers of the observed age-structured heterogeneity; we estimate parameters to match age-structured data from the outbreak in China based on two different sets of assumptions which significantly simplify the model:

  1. A null model where transmission is purely related to contact pattern. This would be the default assumption for any novel pathogen, but for COVID-19 relates to when all infected individuals are equally infectious independent of symptoms.
  2. A model where there is negligible transmission from asymptomatic infections, such that younger age-groups contribute less to overall transmission.

We then combine our fitted transmission parameters with Prem et al.’s estimates of the contact matrices in the other 151 countries to generate an estimate of an age-structured transmission structure for each of those countries. The basic reproductive ratio can then be estimated from these transmission structures to give us an indication of the spreading potential of a COVID-19 outbreak in each of these countries.

Materials and methods

We develop an age-structured transmission model based upon an age-structured contact matrix and allowing for age-stratified variation in susceptibility, probability of clinical symptoms, and reduction in transmission by asymptomatic (and therefore undetected) cases relative to clinically symptomatic (and hence detectable) ones. Full details of the model structure are given in S1 Appendix.

We parameterise the contact matrices for each country using the age-dependent contact rates estimated by Prem et al. [19]. The age-stratified susceptibility, symptomatic probability and asymptomatic transmission profiles need to be inferred from data. Given the large number of parameters involved (three times the number of age classes), we only carry out fits for a few specific assumptions, corresponding to different factors which can drive age-stratified heterogeneity. An important assumption throughout this work is that these age-stratified profiles are the same from population to population; we feel this is justified since they are likely to be driven by physiological mechanisms which are unlikely to be population-specific. In particular, we assume that detection is closely correlated with showing symptoms, so that the age-structured distribution of reported cases reflects the age-structured distribution of symptomatic cases.

For reference, we first consider a simple null model in which all age groups are equally susceptible and all infected individuals (both symptomatic and asymptomatic) transmit at the same rate. That is, we assume that the next-generation matrix (which describes the number of new infections in each age-group generated by a case in a given age-group [5]) is directly proportional to the estimated contact matrix for China; such that contact patterns within the community are the only source of age-structured heterogeneity. The scaling of the contact matrix is performed such that the next-generation matrix generates a basic reproductive ratio of 2.4, in line with current literature. This scaling is applied to all country-specific contact matrices to generate their corresponding transmission matrices and associated basic reproductive ratio. (To match this model to the reported age-distribution of cases in China [4] the age-dependent probability of symptoms can be modified, which in turn determines the pattern of detection although leaves the epidemiology untouched).

We contrast the null model with one in which transmission from asymptomatic cases is negligible, and observed age-structured hetereogeneity is driven by a combination of age-dependent susceptibility and probability of symptoms. The latter has already been identified as of particularly importance, with secondary attack rates for children potentially similar to those among adults, despite the much smaller number of clinically presenting cases [29]. In this model age-structured heterogeneity in susceptibility profile or in the symptomatic probability have identical effects, and so we can fit a single set of parameters capturing both possible effects. Again, we fit this simplified model using epidemiological data from China [4] and the estimated China-level contact matrix [19], and scale the result such that the next-generation matrix generates a basic reproductive ratio of 2.4. The age-dependent parameters and scaling can then be combined with the estimated country-level contact matrices from the other countries in Prem et al.’s study [19] to produce age-stratified next-generation matrices for these countries, and hence an estimate of the country-specific basic reproductive ratio. We explain this calculation in detail in S1 Appendix.

All of the code used to implement our model and generate the results in the next section is available at https://github.com/JBHilton/hilton-keeling-estimating-R0.

Results

In what follows we use an estimate of R0 = 2.4 for the basic reproductive ratio of COVID-19 in China, consistent with values estimated in the literature so far (R0 = 2.2 [1], R0 = 2.3 − 2.6 [11], R0 = 2 − 2.7, [10], R0 = 2.35 [7], R0 = 3.11 [9]), although we note that other estimates of R0 in China can be accommodated by a linear rescaling of all the predicted values.

In Fig 1(a) we present the results of the null-model in which susceptibility is age-independent and all infectious individuals (both symptomatic and asymptomatic) have the same transmission rate(which can be matched to the age-structured case data from China by allowing the probability of displaying symptoms and therefore being detected to vary with age). We plot the estimated R0 for each of the 152 countries included in Prem et al.’s study on a map of the world, with map colours grading from blue (low values of R0) through yellow (intermediate values) to red (high values of R0). Fig 1(b) contrasts this with the estimated basic reproductive ratios based on the age distribution of the first 44,672 confirmed cases in China, as reported by China’s Novel Coronavirus Penumonia Emergency Response Epidemiology Team in China CDC Weekly [4]. Our numerical estimates of the basic reproductive ratio for each of the 152 countries are given in S1 Appendix.

thumbnail
Fig 1. Basic reproductive ratio by country.

(a) Estimated basic reproductive ratio for each country assuming contact structure only; (b) estimated basic reproductive ratio for each country based on the China CDC case data [4]. Gray countries are those not included in Prem et al.’s study [19].

https://doi.org/10.1371/journal.pcbi.1008031.g001

To interpret these findings, in Fig 2(a) we plot the two sets of basic reproductive ratios shown in Fig 1, contrasting the null-model (x-axis) with the model in which there is negligible transmission from asymptomatic infections (y-axis). The assumptions underpinning the latter model lead to an increase in the amount of variation in basic reproductive ratio by country, and in particular can generate much larger basic reproductive ratios. For the null model (or when both symptomatic and asymptomatic infections are equally infectious), the variation in basic reproductive ratio is driven by the variation in average intensity of contacts by country. In contrast, when infection is driven by the age-dependent pattern of symptomatic cases, contact patterns involving members of highly symptomatic age classes becomes particularly important—generating a core-group within the population. For COVID-19, this means that we see higher basic reproductive ratios in countries with older populations, since this generally leads to contacts involving older individuals being more common. In countries with comparatively younger populations, contacts involving older individuals are far less common and so the capacity of the infection to spread is reduced relative to the null model.

thumbnail
Fig 2. Effect of using age-specific susceptibility/symptomatic probability and underlying population pyramids.

(a) Basic reproductive ratio estimates based on age-specific susceptibility or symptomatic probability estimated from China CDC Weekly data versus estimates without age-dependent susceptibility or symptomatic probability. (b) Population pyramids for Niger, China and Italy—China being our reference case (both R0 values equal to 2.4), Niger having the highest R0 in the null model, and Italy having the second highest R0 in the age-specific susceptibility model. The highest R0 attained in the age-specific susceptibility/symptomatic probability model is in Monaco, but since Monaco’s small population is likely to make it an outlier we focus on Italy as an extreme case in the main cluster of ratios. Germany is also labelled in Figure (a); although it has a comparatively small R0 under both sets of assumptions, the proportional change from 1.22 in the null model to 1.99 based on the China CDC data is almost as dramatic as that seen for Italy (2.44 to 4.18). Data from [30].

https://doi.org/10.1371/journal.pcbi.1008031.g002

This principle, that the underlying age-structure of the population drives the estimated basic reproductive ratio, is illustrated in Fig 2(b). Here we compare population pyramids of Niger (which is predicted to have the highest reproductive ratio under the null model) and Italy (which is predicted to have one of the highest reproductive ratios under the second model) in comparison to China (which is the foundation of the parameter values). The population pyramid of Niger is dominated by young children; China has a relatively stable age-structure although there are more individuals in 30-54 age classes than in younger age-groups; the pyramid for Italy shows even fewer children and substantial proportions into older age-classes. We therefore postulate that it is the interaction between the population pyramid and the age-structured probability of symptoms (and hence significant transmission) that largely drives the scaling of the basic reproductive ratio.

This relationship between the population pyramid and the basic reproductive ratio is seen to hold for all countries investigated (Fig 1). In Fig 1(a), we observe that in the null model transmission is generally low in many European countries as well as in South Korea and Japan, and high in many African countries, consistent with the differences in daily number of contacts predicted by Prem et al driven by the proportion of children. However, Fig 1(b) shows that when age-specific susceptibility or symptomatic probability is taken into account, the pattern of infectious potential by country is generally reversed. We then expect to see higher transmission in Eastern Europe (including Italy which had the largest number of cases in Europe in mid March 2020) and Japan, and reduced transmission across Africa, central America, the Middle East and India.

Discussion

Here we have developed a flexible model for age-dependent transmission of SARS-CoV-2 with four forms of heterogeneity: an age-structured contact matrix dependent on the behaviour of the host population; age-dependent susceptibility; age-dependent symptomatic probability; and transmission profiles dependent on physiological response to infection (symptomatic vs asymptomatic). By exploiting a previously-estimated synthetic contact matrix and age-stratified data, and by focusing on the specific case where undetected transmission is negligible, we were able to estimate these age-dependent profiles based on the first 44,672 cases in China [4]. We then combined these estimated profiles with estimates of age-stratified contacts in 151 other countries to generate transmission matrices for these countries from which we can also estimate the scale of basic reproductive ratios in each country relative to China.

We explored two competing models. The null model is purely based on the frequency of age-stratified contacts, and would be the default assumption for any novel pandemic. Given the flexible way in which our model is constructed, allowing an age-dependent probability of displaying symptoms yet symptom independent transmission (asymptomatic and symptomatic infections transmit equally) leads to the same basic reproductive ratio. The second model assumes that detected cases represent a random sample of symptomatic infections and therefore provide an unbiased measure. It is also assumed that the only symptomatic infections are responsible for the overwhelming majority of onward transmission. This allows us to reliably match to the early outbreak data in China [4]. Moreover, the second model allows us to account for age-structure either through age-dependent susceptibility or age-dependent probability of symptoms, or any mixture of the two. We demonstrated that taking such age-specific factors into account results in substantially different predictions of transmission intensity by country relative to a purely contact-based null model; countries with older populations are at substantially higher risk than countries with younger populations.

There are two main limitations to our predictions. The first is the accuracy of the estimated contact matrices; although there are known issues (as discussed in [19]) they remain our best estimate of age-structured contacts to date across many regions of the world. Unfortunately, not all countries have an associated mixing matrix, as many countries (predominantly in Africa) do not have the underlying demographic data recorded that supports the generation of this matrix. This is clearly problematic since mortality from infectious disease is disproportionately concentrated in these countries [31]. Secondly, to estimate the basic reproductive ratios mapped in Fig 1(b) we chose to focus specifically on the case where transmission from undetected cases was negligible. This assumption makes it possible to analytically determine age-structured parameters, and reduces the need to separate the effects of susceptibility and probability of displaying symptoms. If asymptomatic infections transmit as strongly as symptomatic cases, then the country-specific reproductive ratio is expected to be closer to Fig 1(a).

Finally, it is worth stressing that these projections only capture the early phase of the outbreak in the absence of controls. Non-pharmaceutical interventions (contact tracing, self isolation and movement controls) can substantially reduce the infection’s reproductive ratio, with their effectiveness heavily dependent on the pre-intervention basic reproductive ratio, the proportion of contacts traced, and the timing of isolation and movement controls [8, 29, 3234].

Supporting information

S1 Appendix. Model description, calculations, tables of R0 by country.

https://doi.org/10.1371/journal.pcbi.1008031.s001

(PDF)

References

  1. 1. Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, et al. Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia. New England Journal of Medicine. 2020.
  2. 2. Wu JT, Leung K, Leung GM. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study. The Lancet. 2020;395(10225):689–697.
  3. 3. World Health Organisation. Coronavirus disease (COVID-19) Pandemic; 2020. https://web.archive.org/web/20200529095724/https://www.who.int/emergencies/diseases/novel-coronavirus-2019.
  4. 4. Novel Coronavirus Pneumonia Emergency Response Epidemiology Team. The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) in China. China CDC Weekly. 2020;41(2):145.
  5. 5. Keeling MJ, Rohani P. Modeling infectious diseases in humans and animals. Princeton University Press; 2011.
  6. 6. Diekmann O, Heesterbeek JAP, Metz JA. On the definition and the computation of the basic reproduction ratio R 0 in models for infectious diseases in heterogeneous populations. Journal of mathematical biology. 1990;28(4):365–382. pmid:2117040
  7. 7. Kucharski AJ, Russell TW, Diamond C, Liu Y, Edmunds J, Funk S, et al. Early dynamics of transmission and control of COVID-19: a mathematical modelling study. The Lancet Infectious Diseases. 2020.
  8. 8. Tang B, Xia F, Tang S, Bragazzi NL, Li Q, Sun X, et al. The Evolution of Quarantined and Suspected Cases Determines the Final Trend of the 2019-nCoV Epidemics Based on Multi-Source Data Analyses. Available at SSRN 3537099. 2020.
  9. 9. Read J, Bridgen J, Cummings D, Ho A, Jewell C. Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions. medRxiv. 2020; p. 2020.01.23.20018549.
  10. 10. Abbott S, Hellewell J, Munday J, Funk S. The transmissibility of novel Coronavirus in the early stages of the 2019-20 outbreak in Wuhan: Exploring initial point- source exposure sizes and durations using scenario analysis. Wellcome Open Research. 2020;5:17. pmid:32322691
  11. 11. Liu T, Hu J, Kang M, Lin L, Zhong H, Xiao J, et al. Transmission dynamics of 2019 novel coronavirus (2019-nCoV). bioRxiv. 2020.
  12. 12. Anderson RM, Anderson B, May RM. Infectious diseases of humans: dynamics and control. Oxford university press; 1992.
  13. 13. Schenzle D. An age-structured model of pre-and post-vaccination measles transmission. Mathematical Medicine and Biology: A Journal of the IMA. 1984;1(2):169–191.
  14. 14. Cauchemez S, Valleron AJ, Boelle PY, Flahault A, Ferguson NM. Estimating the impact of school closure on influenza transmission from Sentinel data. Nature. 2008;452(7188):750–754. pmid:18401408
  15. 15. Wallinga J, van Boven M, Lipsitch M. Optimizing infectious disease interventions during an emerging epidemic. Proceedings of the National Academy of Sciences. 2010;107(2):923–928.
  16. 16. Geard N, Glass K, McCaw JM, McBryde ES, Korb KB, Keeling MJ, et al. The effects of demographic change on disease transmission and vaccine impact in a household structured population. Epidemics. 2015;13:56–64. pmid:26616042
  17. 17. Kucharski AJ, Edmunds WJ. Characterizing the transmission potential of zoonotic infections from minor outbreaks. PLoS computational biology. 2015;11(4). pmid:25860289
  18. 18. Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS medicine. 2008;5(3). pmid:18366252
  19. 19. Prem K, Cook AR, Jit M. Projecting social contact matrices in 152 countries using contact surveys and demographic data. PLoS computational biology. 2017;13(9):e1005697. pmid:28898249
  20. 20. Read JM, Lessler J, Riley S, Wang S, Tan LJ, Kwok KO, et al. Social mixing patterns in rural and urban areas of southern China. Proceedings of the Royal Society B: Biological Sciences. 2014;281(1785):20140268. pmid:24789897
  21. 21. Danon L, House TA, Read JM, Keeling MJ. Social encounter networks: collective properties and disease transmission. Journal of The Royal Society Interface. 2012;9(76):2826–2833.
  22. 22. Kiti MC, Kinyanjui TM, Koech DC, Munywoki PK, Medley GF, Nokes DJ. Quantifying age-related rates of social contact using diaries in a rural coastal population of Kenya. PloS one. 2014;9(8). pmid:25127257
  23. 23. Edmunds WJ, O’callaghan C, Nokes D. Who mixes with whom? A method to determine the contact patterns of adults that may lead to the spread of airborne infections. Proceedings of the Royal Society of London Series B: Biological Sciences. 1997;264(1384):949–957. pmid:9263464
  24. 24. Horby P, Thai PQ, Hens N, Yen NTT, Mai LQ, Thoang DD, et al. Social contact patterns in Vietnam and implications for the control of infectious diseases. PloS one. 2011;6(2). pmid:21347264
  25. 25. Grijalva CG, Goeyvaerts N, Verastegui H, Edwards KM, Gil AI, Lanata CF, et al. A household-based study of contact networks relevant for the spread of infectious diseases in the highlands of Peru. PloS one. 2015;10(3).
  26. 26. Ajelli M, Litvinova M. Estimating contact patterns relevant to the spread of infectious diseases in Russia. Journal of Theoretical Biology. 2017;419: 1–7.
  27. 27. Johnstone-Robertson SP, Mark D, Morrow C, Middelkoop K, Chiswell M, Aquino LD, et al. Social mixing patterns within a South African township community: implications for respiratory disease transmission and control. American journal of epidemiology. 2011;174(11):1246–1255. pmid:22071585
  28. 28. Prem K, Liu Y, Russell T, Kucharski AJ, Eggo RM, Davies N, et al. The effect of control strategies that reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China. medRxiv. 2020.
  29. 29. Bi Q, Wu Y, Mei S, Ye C, Zou X, Zhang Z, et al. Epidemiology and Transmission of COVID-19 in Shenzhen China: Analysis of 391 cases and 1,286 of their close contacts. medRxiv. 2020.
  30. 30. Population Pyramids of the World from 1950 to 2010; 2019. http://www.PopulationPyramid.net.
  31. 31. World Health Organisation. Global Health Estimates 2016 Summary Tables: Deaths by cause, age and sex, by World Bank Income group, 2000-2015; 2018. https://web.archive.org/web/20200125140038/https://www.who.int/healthinfo/global_burden_disease/estimates/en/.
  32. 32. Hellewell J, Abbott S, Gimma A, Bosse NI, Jarvis CI, Russell TW, et al. Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts. The Lancet Global Health. 2020.
  33. 33. Stoecklin SB, Rolland P, Silue Y, Mailles A, Campese C, Simondon A, et al. First cases of coronavirus disease 2019 (COVID-19) in France: surveillance, investigations and control measures, January 2020. Eurosurveillance. 2020;25(6).
  34. 34. Anderson RM, Heesterbeek H, Klinkenberg D, Hollingsworth TD. How will country-based mitigation measures influence the course of the COVID-19 epidemic? The Lancet. 2020;395(10228):931–934.