Estimation of basic reproduction number of the Middle East respiratory syndrome coronavirus (MERSCoV) during the outbreak in South Korea, 2015
 HyukJun Chang^{1}Email author
Received: 31 January 2017
Accepted: 5 June 2017
Published: 13 June 2017
Abstract
Background
In South Korea, an outbreak of Middle East respiratory syndrome (MERS) occurred in 2015. It was the second largest MERS outbreak. As a result of the outbreak in South Korea, 186 infections were reported, and 36 patients died. At least 16,693 people were isolated with suspicious symptoms. This paper estimates the basic reproduction number of the MERS coronavirus (CoV), using data on the spread of MERS in South Korea.
Methods
The basic reproduction number of an epidemic is defined as the average number of secondary cases that an infected subject produces over its infectious period in a susceptible and uninfected population. To estimate the basic reproduction number of the MERSCoV, we employ data from the 2015 South Korea MERS outbreak and the susceptibleinfectedremoved (SIR) model, a mathematical model that uses a set of ordinary differential equations (ODEs).
Results
We fit the model to the epidemic data of the South Korea outbreak minimizing the sum of the squared errors to identify model parameters. Also we derive the basic reproductive number as the terms of the parameters of the SIR model. Then we determine the basic reproduction number of the MERSCoV in South Korea in 2015 as 8.0977. It is worth comparing with the basic reproductive number of the 2014 Ebola outbreak in West Africa including Guinea, Sierra Leone, and Liberia, which had values of 1.5–2.5.
Conclusions
There was no intervention to control the infection in the early phase of the outbreak, thus the data used here provide the best conditions to evaluate the epidemic characteristics of MERS, such as the basic reproduction number. An evaluation of basic reproduction number using epidemic data could be problematic if there are stochastic fluctuations in the early phase of the outbreak, or if the report is not accurate and there is bias in the data. Such problems are not relevant to this study because the data used here were precisely reported and verified by Korea Hospital Association.
Keywords
Basic reproduction number Susceptibleinfectedremoved (SIR) model Middle East respiratory syndrome (MERS) MERS coronavirus (MERSCoV)Background
The Middle East respiratory syndrome (MERS) is caused by a coronavirus (CoV), the MERSCoV. In Saudi Arabia, the first case of the disease was reported in 2012 [1]. The first case of MERS in the Republic of Korea was identified on 20 May 2015 [2]. A significant outbreak of MERS occurred in South Korea and lasted for almost three months, from May to July 2015.
The 2015 MERS spread in South Korea is the second largest outbreak recorded to date [3]. As a result of the outbreak in South Korea, 186 infections were reported, and 36 patients died. At least 16,693 people were isolated with suspicious symptoms [4]. This paper evaluates the basic reproduction number of MERSCoV, using data from the 2015 South Korea outbreak.
The basic reproduction number (generally denoted as \(R_0\)) of an epidemic is defined as the average number of secondary cases that an infected subject produces over its infectious period in a susceptible and uninfected population [5–11]. It can estimate the growth rate of an infectious disease at the early stage of the outbreak, when most individuals are susceptible [12]. The basic reproduction number of an epidemic is useful for determining whether an outbreak of the disease will occur or not [13], and for analyzing epidemic properties of the disease further [14].
Note that \(R_0\) is referred to as the basic reproductive (or reproduction) number (or ratio). The basic reproductive (or reproduction) rate is incorrect nomenclature because \(R_0\) is a dimensionless number that is not related to any physical quantity corresponding to rate.
Based on [11], the reasons to estimate the basic reproduction number of an epidemic are summarized as follows: First, we can relatively evaluate the risk of the corresponding epidemic using \(R_0\). In other words, we can compare the infectivity of the epidemic with others, already familiar to us. Second, the reproduction number can be evaluated multiple times, e.g., before (and after) an infection control measure intervention. To this end, it is needed to distinguish the reproduction number after control intervention from the basic reproductive number \(R_0\), which is estimated before the intervention. Now we refer to the reproduction number after control intervention as the effective reproduction number, denoted by \(R_{\text {eff}}\). Then we can compare \(R_{\text {eff}}\) with \(R_0\), and we can evaluate the efficacy of a control measure quantitatively based on \(R_{\text {eff}}\). By doing so, eventually we can determine how to apply control intervention to reduce \(R_{\text {eff}}\) to less than one. If \(R_{\text {eff}} < 1\), it is concluded that the control intervention works effectively, and that the outbreak will eventually be controlled by reducing the reproduction number to less than the threshold level, 1.
In this paper, to estimate the basic reproduction number of the MERSCoV, we employ data from the 2015 South Korea MERS outbreak and the susceptibleinfectedremoved (SIR) model [13, 15, 16], a mathematical model that uses a set of ordinary differential equations (ODEs). Because the availability of epidemic data is limited, we usually employ nonstructured deterministic models to evaluate \(R_0\) [11].
Based on the data reported from the 2015 South Korea outbreak of MERS, we evaluate the basic reproductive number of the virus, MERSCoV. We fit the model to epidemic data from the South Korea outbreak, and identify model parameters and the basic reproduction number, \(R_0\). Note that other epidemiological parameters, such as incubation period and serial interval, have been discussed in [17] for the outbreak.
Preliminary work relating to this paper is presented in [18]. This paper includes an analysis of the derivation of the basic reproduction number, careful screening of the reported data, and a sophisticated approach using the sum of the squared errors to evaluate the basic reproduction number precisely.
A number of papers and books have been dedicated to the study of \(R_0\) for other infectious diseases. A few of them are as follows: see [19] for the \(R_0\) of severe acute respiratory syndrome(SARS), [20] for the \(R_0\) of influenza, [21, 22] for the \(R_0\) of Ebola, and [23] for the \(R_0\) of malaria.
Although some literature on the study of \(R_0\) of the MERSCoV has been reported such as [12, 24], the MERS outbreak in South Korea is unique [25]. MERS spread almost naturally without any intervention in the early stage, and the Korean government did not respond appropriately [3, 25]. The list of medical facilities involved was not even announced to public. Ironically, it is the ideal condition to fit a mathematical model to the clinical epidemic data and to evaluate epidemic properties of the MERSCoV, including the basic reproduction number.
The paper is organized as follows: In “Methods” section, a mathematical model, which comprises a set of ordinary differential equations, is introduced, and an estimation method is discussed for the basic reproductive number. “Results” section describes the evaluation result, providing further discussion. “Conclusions” section concludes the paper by suggesting additional research.
Methods
In this section, we briefly discuss the method employed in this paper, including a definition of the basic reproduction number, and we introduce the SIR model.
Basic reproduction number
The basic reproduction number is defined as the number of secondary cases that one infected primary subject causes on average in an uninfected and totally susceptible population, over the infectious period [8–10]. Based on this definition, we obtain a mathematical description of \(R_0\) via the socalled ‘survival function’ [7, 11].
Formula (1) is derived from the \(R_0\) definition and can thus be used for any mathematical model, not just models given by ODEs. However, it requires explicit expressions for \(F(\cdot )\) and \(b(\cdot )\), which are functions of time. This paper employs the SIR model described by ODEs that is introduced in the next section. Because the model does not provide explicit descriptions for \(F(\cdot )\) and \(b(\cdot )\), we use an alternative expression for \(R_0\), derived from the SIR model.
SIR model
State variables, parameters  Descriptions 

S  Number of susceptible subjects 
I  Number of infected subjects 
R  Number of removed subjects 
\(\beta\)  Disease transmission rate 
\(\nu\)  Removed rate 
\(\tau\)  Transmissibility of the infection 
\(\kappa\)  Number of transmittable contacts by infected patient per unit time 
Any term related to birth and death in the population that is not caused by MERS is not included in the model (2). The dynamics of the disease (e.g., infection or recovery) is assumed to be significantly faster than that of birth and death in the population. Generally, epidemic models such as SIR do not include birth and death because zero net change of the population is assumed. If we model an infectious disease with comparatively slow dynamics (e.g., an endemic disease), we must consider dynamic terms describing birth and death.
Results
In Table 2, the history of the MERSCoV spread status is presented. The Ministry of Health and Welfare, Korea officially announced the data. “Infected” represents the accumulated number of infected patients. “Deceased” is the number of dead subjects. “Recovered” is the number of individuals returning to healthy status. All entries in the table are as of the “Date”. The number of infected patients includes both removed and recovered patients.
Table 2 shows the MERSCoV spread data for only the initial phase of the outbreak, i.e., from 20 May 2015 to 12 June 2015.
On 7 June 2015, the South Korean government disclosed to the public the list of all hospitals exposed to MERSCoV, with the dates and duration of exposure [4]. This is the first intervention of the government to control the spread. Before this date, there was no control action that could affect estimation of the basic reproduction number of MERSCoV.
Accumulated MERSCoV patients in Korea, 2015 [4]
Date  Infected  Deceased  Recovered 

20 May  2  0  0 
21 May  3  0  0 
22 May  3  0  0 
23 May  3  0  0 
24 May  3  0  0 
25 May  3  0  0 
26 May  5  0  0 
27 May  5  0  0 
28 May  7  0  0 
29 May  13  0  0 
30 May  15  0  0 
31 May  18  0  0 
1 June  25  1  0 
2 June  30  1  0 
3 June  30  3  0 
4 June  36  4  0 
5 June  42  5  1 
6 June  64  5  1 
7 June  87  5  1 
8 June  95  7  2 
9 June  108  7  3 
10 June  122  9  4 
11 June  126  10  7 
12 June  138  13  9 
We search for the parameter pair (\(\kappa \tau\), \(\nu\)) such that can respond appropriately with the data in Table 2. To evaluate how closely system response is fitted to the data in Table 2, we employ a quantitative measure, the sum of the squared errors. Once we obtain the optimal values for the parameters with respect to this measure, we can estimate the basic reproduction number as described in the “Methods” section.
To compare the quantitative measures for each pair of parameters, we consider the plane, i.e., 2dimensional space, of the parameters, (\(\kappa \tau\), \(\nu\)). We can obtain a surface in 3dimensional space by plotting the corresponding measure as the value along the third axis.
The function \(f_E(\cdot , \cdot )\) describes the sum of squared error between the outbreak data of Table 2 and the simulation result of model (5) with the initial condition (6) and the function arguments. Thus, if we find the parameters minimizing the function (7), then these parameters can be considered to correspond to the case of Table 2.
We determine the \(R_0\) of the MERSCoV in South Korea in 2015 as 8.0977 (i.e., 0.21153/0.026122). It is worth comparing with \(R_0\) of the 2014 Ebola outbreak, which had values of 1.5–2.5 [21].
Conclusions
In this paper, we evaluated the basic reproduction number of the MERSCoV outbreak that occurred in 2015 in South Korea, using officially reported data. We employed a mathematical dynamic model, the SIR model. We first fit the response of the SIR model to the epidemic curve data reported from the MERS outbreak. Then, we identified the system parameters of the model to estimate the basic reproduction number.
Because there was no intervention to control the infection in the early phase of the outbreak, the data used here provide the best conditions to evaluate the epidemic characteristics of MERS, such as the basic reproduction number. An evaluation of \(R_0\) using epidemic data could be problematic if there are stochastic fluctuations in the early phase of the outbreak, or if the report is not accurate and there is bias in the data [11]. Such problems are not relevant to this study because the data used here were precisely reported and verified by [4].
We conclude this paper with the following discussion on future work to overcome the limitations of research, derived from assumptions in the paper.
Further research direction

Behind the SIR model (2), there are several strong assumptions, one of which is a zero latent period, i.e., the incubation period is zero. This implies that a patient becomes infectious immediately after infection. However, the incubation phase occurs during the course of the MERS outbreak. To address this weak point, in future work we could consider the 4dimensional SEIR (i.e., susceptibleexposedinfectiousremoved) model, which has been employed in [28] to study Ebola epidemic model. The additional state in the SEIR model can help us deal with the latent period.

In this paper, we considered the epidemic curve data in [4] only from the early stage of the 2015 MERS outbreak in South Korea, where there was no intervention to control the spread. Accordingly, we evaluated \(R_0\) based on the data. In future work, we will also consider the epidemic data in [4] from the later (or closing) stage of the MERS spread in South Korea in 2015, so we can estimate the effective production number (i.e., \(R_\mathrm{eff}\)), which is the production number resulting from interventions, such as education, quarantine, and the tracing of contacts by infected patients. By doing so, we can evaluate the effectiveness of each control measure on the spread of the infectious disease[29]. Eventually, such evaluation could help us improve public health policy.
Declarations
Acknowledgements
A preliminary form of this research was presented at the 16th International Conference on Control, Automation and Systems, ICCAS 2016.
Competing interests
The author declares that no competing interests.
Availability of data and materials
All data generated or analysed during this study are included in this published article.
Funding
This work was supported by National Research Foundation of Korea  Grant funded by the Korean Government (NRF2014R1A1A1003056).
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
 World Health Organization: Fact Sheet 401: Middle East respiratory syndrome coronavirus (MERSCoV). 2015. http://www.who.int/mediacentre/factsheets/merscov/en/. Accessed 8 May 2017.
 Su S, Wong G, Liu Y, Gao GF, Li S, Bi Y. MERS in South Korea and China: a potential outbreak threat? Lancet. 2015;385(9985):2349–50.View ArticleGoogle Scholar
 Hyun MH, Park Y. Re: Why the panic? South Korea’s MERS response questioned. BMJ. 2015. doi:https://doi.org/10.1136/bmj.h3403.View ArticleGoogle Scholar
 Korea Hospital Association MERS. White Paper by Korea Hospital Association. Korea Hospital Association, Seoul; 2015.Google Scholar
 Anderson RM, May RM. Infectious diseases of humans: dynamics and control. New York: Oxfor University Press; 1992.Google Scholar
 Diekmann O, Heesterbeek JAP, Metz JAJ. On the definition and the computation of the basic reproduction ratio \({R}_0\) in models for infectious diseases in heterogeneous populations. J Math Biol. 1990;28(4):365–82.MathSciNetView ArticleMATHGoogle Scholar
 Heesterbeek JAP, Dietz K. The concept of \(R_0\) in epidemic theory. Stat Neerl. 1996;50(1):89–110.MathSciNetView ArticleMATHGoogle Scholar
 Diekmann O, Heesterbeek JAP. Mathematical epidemiology of infectious diseases: model building. Analysis and interpretation. New York: Wiley; 2000.MATHGoogle Scholar
 van den Driessche P, Watmough J. Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002;180(1–2):29–48.MathSciNetView ArticleMATHGoogle Scholar
 Heesterbeek JAP. A brief history of \(R_0\) and a recipe for its calculation. Acta Biotheor. 2002;50(3):189–204.View ArticleGoogle Scholar
 Heffernan JM, Smith RJ, Wahl LM. Perspectives on the basic reproductive ratio. J R Soc Interface. 2005;2(4):281–93.View ArticleGoogle Scholar
 Chowell G, Blumberg S, Simonsen L, Miller MA, Viboud C. Synthesizing data and models for the spread of MERSCoV, 2013: key role of index cases and hospital transmission. Epidemics. 2014;9:40–51.View ArticleGoogle Scholar
 DiStefano J III. Dynamic systems biology modeling and simulation. London: Academic Press; 2015.Google Scholar
 Yashima K, Sasaki A. Spotting epidemic keystones by \(R_0\) sensitivity analysis: highrisk stations in the Tokyo metropolitan area. PLOS ONE. 2016;11(9):1–19.View ArticleGoogle Scholar
 Anderson RM. Discussion: the KermackMcKendrick epidemic threshold theorem. Bull Math Biol. 1991;53(1–2):1–32.View ArticleGoogle Scholar
 Weiss H. The SIR model and the foundations of public health. MATerials MATemàtics. 2013;3(17):1–17.Google Scholar
 Cowling BJ, Park M, Fang VJ, Wu P, Leung GM, Wu JT. Preliminary epidemiologic assessment of MERSCoV outbreak in South Korea, May to June 2015. Euro Surveill. 2015;20(25):1–7.View ArticleGoogle Scholar
 Chang H. Evaluation of the basic reproduction number of MERSCoV during the 2015 outbreak in South Korea. In: Proceedings of 16th international conference on control, automation and systems. 2016. p. 981–984Google Scholar
 Wallinga J, Teunis P. Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures. Am J Epidemiol. 2004;160(6):509–16.View ArticleGoogle Scholar
 Mills CE, Robins JM, Lipsitch M. Transmissibility of 1918 pandemic influenza. Nature. 2004;432:904–6.View ArticleGoogle Scholar
 Althaus CL. Estimating the reproduction number of ebola virus (EBOV) during the 2014 outbreak in west africa. PLOS Curr Outbreaks. 2014;6:1–5.Google Scholar
 Chowell G, Hengartner NW, CastilloChavez C, Fenimore PW, Hyman JM. The basic reproductive number of Ebola and the effects of public health measures: the cases of Congo and Uganda. J Theor Biol. 2004;229(1):119–26.MathSciNetView ArticleGoogle Scholar
 Smith DL, McKenzie FE, Snow RW, Hay SI. Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biol. 2007;5(3):1–12.View ArticleGoogle Scholar
 Ejima K, Aihara K, Nishiura H. Probabilistic differential diagnosis of middle east respiratory syndrome (MERS) using the time from immigration to illness onset among imported cases. J Theor Biol. 2014;346:47–53.View ArticleGoogle Scholar
 Ha KM. A lesson learned from the MERS outbreak in South Korea in 2015. J Hosp Infect. 2016;92(3):232–4.View ArticleGoogle Scholar
 Assiri A, McGeer A, Perl TM, Price CS, Al Rabeeah AA, Cummings DAT, Alabdullatif ZN, Assad M, Almulhim A, Makhdoom H, Madani H, Alhakeem R, AlTawfiq JA, Cotten M, Watson SJ, Kellam P, Zumla AI, Memish ZA. Hospital outbreak of middle east respiratory syndrome coronavirus. N Engl J Med. 2013;369(5):407–16.View ArticleGoogle Scholar
 Area I, Losada J, Ndaïrou F, Nieto JJ, Tcheutia DD. Mathematical modeling of 2014 Ebola outbreak. Math Methods Appl Sci. 2015. doi:https://doi.org/10.1002/mma.3794.
 Area I, Batarfi H, Losada J, Nieto JJ, Shammakh W, Torres Á. On a fractional order Ebola epidemic model. Adv Differ Equ. 2015;2015(1):278.MathSciNetView ArticleMATHGoogle Scholar
 Nishiura H, Ejima K, Mizumoto K, Nakaoka S, Inaba H, Imoto S, Yamaguchi R, Saito M. Costeffective length and timing of school closure during an influenza pandemic depend on the severity. Theor Biol Med Model. 2014;11(1):5.View ArticleGoogle Scholar