Interpretation of the Basic and Effective Reproduction Number

Article information

J Prev Med Public Health. 2020;53(6):405-408
Publication date (electronic) : 2020 October 20
doi : https://doi.org/10.3961/jpmph.20.288
1College of Veterinary Medicine and Institute of Veterinary Science, Kangwon National University, Chuncheon, Korea
2Department of Public Health Science, Graduate School of Public Health, Seoul National University, Seoul, Korea
3Department of Preventive Medicine, Konyang University College of Medicine, Daejeon, Korea
Corresponding author: Son-Il Pak, DVM, PhD College of Veterinary Medicine and Institute of Veterinary Science, Kangwon National University, Chuncheon 24341, Korea E-mail: paksi@kangwon.ac.kr
Received 2020 June 19; Revised 2020 October 2; Accepted 2020 October 7.

Abstract

In epidemiology, the basic reproduction number (R0) is a term that describes the expected number of infections generated by 1 case in a susceptible population. At the beginning of the coronavirus disease 2019 (COVID-19) pandemic, R0 was frequently referenced by the public health community and the wider public. However, this metric is often misused or misinterpreted. Moreover, the complexity of the process of estimating R0 has caused difficulties for a substantial number of researchers. In this article, in order to increase the accessibility of this concept, we address several misconceptions related to the threshold characteristics of R0 and the effective reproduction number (Rt). Moreover, the appropriate interpretation of the metrics is discussed. R0 should be considered as a population-averaged value that pools the contact structure according to a stochastic transmission process. Furthermore, it is necessary to understand the unavoidable time lag for Rt due to the incubation period of the disease.

INTRODUCTION

As a key parameter of infectious disease epidemiology, the basic reproduction number (R0) and reproduction number (R), defined as the expected number of infected individuals from an infectious person in a totally or partially susceptible population, respectively, were frequently referenced during the early stages of the coronavirus disease 2019 (COVID-19) pandemic [1,2]. Although common misconceptions surrounding the metrics have been addressed [3], confusion persists related to the threshold characteristics of R0, R, and the effective reproduction number (Rt). To respond to infectious disease outbreaks, it is crucial to understand these metrics in more detail in order to improve communication and understanding among the fields of public health and social science. Thus, this article aims to highlight some points of caution related to R0 and Rt.

Stochastic Population-averaged Characteristics of Basic Reproduction Number and Reproduction Number

The estimation of R0 and R is largely reliant upon on the estimation method, contact structure, and epidemiological parameters such as the infectious period and incubation period [3]. In addition, it must be kept in mind that these estimates are population-averaged values during the epidemic period. Since the study population is usually not totally susceptible, it is appropriate to term these estimates as R, not R0. The estimation of these metrics involves averaging values at the study population level, which in turns requires averaging variations across regions and contact structures. Therefore, the estimates should be applied only at the study population level and not to subgroups of the study population. Specifically, these metrics can have values above 1 in a certain region or contact network, even if the estimated population-averaged value is below 1. Even more importantly, the presence of asymptomatic and undetected cases of COVID-19 may bias the estimation of R. However, since a certain proportion of infected persons would be asymptomatic at the population level, the estimated R is expected to be reliable since R is the ratio of the number of infected cases to the number of cases showing symptoms during the study period.

In epidemiological theory, R has a threshold centered on 1. However, R is not related to the intensity of transmission. Furthermore, it neither measures nor reflects the risk of infection [3]. Thus, when R is below 1, there are still newly infected cases for a certain period. In particular, in the early phase of an epidemic, the disease dynamics are strongly affected by the initial conditions [4]. Therefore, the number of infected individuals attributed to an infected person in the early phase can be above 1, even though R is less than 1 [4]. Secondly, even if R is above 1, it is possible that the disease will die out since the transmission of an infectious disease is a branching process with a negative binomial probability distribution. The distribution is defined by 2 parameters, mean and dispersion, with the former being R and the latter referring to variance. Even if R is greater than 1, if the variance is high, the probability of extinction increases. The reason for this is that under such circumstances, most of the infected individuals have a low R, whereas a few infected people have a high R; in this sense, the overall situation resembles the disease dynamics when R<1. Nonetheless, this possibility clearly decreases as R increases or dispersion decreases [5]. Whether or not a certain disease goes extinct is also related to critical community size, which is defined as the smallest population in which the disease would not go extinct without re-introduction [4]. Thus, a disease with R>1 can be eradicated if the population number is under the critical community size.

Taken together, R is a population-averaged value during an epidemic period. The literature suggests that each R and its variance differ according to the region, contact structure, and time. It is appropriate to estimate and interpret each R based on a careful consideration of its epidemiological implications.

Effective Reproduction Number

In order to assess the current dynamics of infectious disease transmission in a timely manner, Rt is commonly utilized. This metric can be used to quantify the efficacy of control measures for a disease in real time [6]. In order to estimate Rt, it is crucial to use the number of the cases, the serial interval, and time of symptom onset. As these data are commonly not publicly available during epidemics, especially for an emerging infectious disease, the absence of data can bias Rt. In particular, a surveillance system can only report the incidence of symptoms, not the incidence of infections. Thus, Rt illustrates the delayed dynamics of transmission due to time lag, including the incubation period and delayed reports. That is, the estimates reflect the dynamics of the time period when the current symptomatic cases were infected, not the current dynamics. Likewise, the estimates should be carefully interpreted if the metric is based on the reported dates [7]. For example, as shown in Figure 1, if the infected case is ideally reported at time T2 (“Report 1” in Figure 1), the Rt estimated from symptom onset will show the transmission dynamics at time T1. However, if Rt is estimated using the report data at time T3 (“Report 2” in Figure 1) due to the time lag of surveillance (“Time lag” in Figure 1), there is a high likelihood of uncertainty whether the estimated Rt describes the dynamics at time T1. Although information on the incubation period could overcome these uncertainties, there are still limitations for emerging infectious diseases [8]. Lastly, as Rt is similar to R0 in that it is a population-averaged value, it is recommended to be estimate and apply Rt for each cluster or administrative level.

Figure. 1.

Different reporting time points by the surveillance system based on the incidence of symptoms (report 1) and by the one with time lag (report 2).

Other Factors Affecting Variation

In this study, we discussed the R only in terms of the contact rate, not in terms of the probability of infection per contact and infectious period, for which variation can occur across different types of the pathogen responsible for a disease (e.g., subtypes, clades, or strains). This factor is highly likely to add uncertainty to the estimated R in that the estimates show the population-averaged characteristics of infection dynamics. However, these possible variations could be reflected through appropriate components in a mathematical model. Moreover, if cases of the same type are spatiotemporally clustered, as in previous COVID-19 clusters, the suggestions provided in this article to estimate the R for clustered or epidemiologically-related populations could increase the likelihood that the estimates will reflect the epidemiological characteristics of different types of the pathogen.

CONCLUSION

In order to understand transmission dynamics, it should be recognized that the R is a population-averaged value reflecting a stochastic transmission process. Moreover, to understand the transmission dynamics of infectious diseases, Rt, which requires data related to the serial interval and time of symptom onset, should be considered as well as R0.

Ethics Statement

This paper is a perspective, so it did not need ethical approval.

Notes

CONFLICT OF INTEREST

The authors have no conflicts of interest associated with the material presented in this paper.

FUNDING

This study was supported by a grant from the Institute of Information & Communication Technology Planning & Evaluation (IITP) funded by the Korea government (No. 2018-0-00430).

ACKNOWLEDGEMENTS

The authors would like to thank Jin A Kim (Daegu Center for Infectious Diseases Control and Prevention) and Kyung-Duk Min (Seoul National University) for their comments on the manuscript.

Notes

AUTHOR CONTRIBUTIONS

Conceptualization: JSL, SIC, SR, SIP. Funding acquisition: SIP. Writing – original draft preparation: JSL. Writing – review & editing: JSL, SIC, SR, SIP.

References

1. Ryu S, Chun BC; Korean Society of Epidemiology 2019-nCoV Task Force Team. An interim review of the epidemiological characteristics of 2019 novel coronavirus. Epidemiol Health 2020;42:e2020006.
2. Ryu S, Ali ST, Lim JS, Chun BC. Estimation of the excess COVID-19 cases in Seoul, South Korea by the students arriving from China. Int J Environ Res Public Health 2020;17(9):3113.
3. Delamater PL, Street EJ, Leslie TF, Yang YT, Jacobsen KH. Complexity of the basic reproduction number (R0). Emerg Infect Dis 2019;25(1):1–4.
4. Keeling MJ, Rohani P. Modeling infectious diseases in humans and animals Princeton: Princeton University Press; 2011. p. 62.
5. Lloyd-Smith JO, Schreiber SJ, Kopp PE, Getz WM. Superspreading and the effect of individual variation on disease emergence. Nature 2005;438(7066):355–359.
6. Cauchemez S, Boëlle PY, Thomas G, Valleron AJ. Estimating in real time the efficacy of measures to control emerging communicable diseases. Am J Epidemiol 2006;164(6):591–597.
7. Hwang J, Park H, Jung J, Kim SH, Kim N. Basic and effective reproduction numbers of COVID-19 cases in South Korea excluding Sincheonji cases. medRxiv 2020;doi: https://doi.org/10.1101/2020.03.19.20039347.
8. Thompson RN, Stockwin JE, van Gaalen RD, Polonsky JA, Kamvar ZN, Demarsh PA, et al. Improved inference of time-varying reproduction numbers during infectious disease outbreaks. Epidemics 2019;29:100356.

Article information Continued

Figure. 1.

Different reporting time points by the surveillance system based on the incidence of symptoms (report 1) and by the one with time lag (report 2).