Research Article  Open  Published:
Estimating the basic reproductive ratio for the Ebola outbreak in Liberia and Sierra Leone
Infectious Diseases of Povertyvolume 4, Article number: 13 (2015)
Abstract
Background
Ebola virus disease has reemerged as a major public health crisis in Africa, with isolated cases also observed globally, during the current outbreak.
Methods
To estimate the basic reproductive ratio R_{0}, which is a measure of the severity of the outbreak, we developed a SEIR (susceptibleexposedinfectedrecovered) type deterministic model, and used data from the Centers for Disease Control and Prevention (CDC), for the Ebola outbreak in Liberia and Sierra Leone. Two different data sets are available: one with raw reported data and one with corrected data (as the CDC suspects underreporting).
Results
Using a deterministic ordinary differential equation transmission model for Ebola epidemic, the basic reproductive ratio R_{0} for Liberia resulted to be 1.757 and 1.9 for corrected and uncorrected case data, respectively. For Sierra Leone, R_{0} resulted to be 1.492 and 1.362 for corrected and uncorrected case data, respectively. In each of the two cases we considered, the estimate for the basic reproductive ratio was initially greater than unity leading to an epidemic outbreak.
Conclusion
We obtained robust estimates for the value of R_{0} associated with the 2014 Ebola outbreak, and showed that there is close agreement between our estimates of R_{0}. Analysis of our model also showed that effective isolation is required, with the contact rate in isolation less than one quarter of that for the infected nonisolated population, and that the fraction of highrisk individuals must be brought to less than 10% of the overall susceptible population, in order to bring the value of R_{0} to less than 1, and hence control the outbreak.
Multilingual abstracts
Please see Additional file 1 for translations of the abstract into the six official working languages of the United Nations.
Background
Ebola virus disease (EVD), named after the Ebola River in Zaire, is known to be a highly contagious disease with a high mortality rate [1,2]. Previously known as Ebola hemorrhagic Fever, EVD has a number of different strains. Originating in Sudan and Zaire in 1976, there has been a number of over the years [1,2]. From 1976 to 2008, the total case fatality rate for EVD victims was 79% [3]. The ongoing outbreak of EVD is affecting multiple countries in Central and Western Africa [2]. Beginning in December 2013 in West Africa, precisely Guinea, the EVD outbreak spread to Sierra Leone, Liberia and Nigeria with fatality rates of 73.2%, 43.0%, 52.5% and 33.3% respectively [4]. In May 2014, the second outbreak was confirmed in nearby regions, including Sierra Leone and Liberia [3]. As of October 14th 2014, 4,555 Ebola deaths have been reported in these countries, including one death in the United States [5].
EVD is a viral infection caused by a virus of the family Filoviridae, genus Ebolavirus. There are five identified subspecies of the Ebolavirus. Four of the five subspecies: i) Ebolavirus (Zaire ebolavirus), EBOV ii) Sudan virus (Sudan ebolavirus), SUDV, iii) Bundibugyo virus (Bundibugyo ebolavirus), BDBV and iv) Ta Forest virus (Ta Forest ebolavirus), TAFV, have caused disease in humans. The fifth, called the Reston virus (Reston ebolavirus), has caused disease in nonhuman primates [2]. The primary source of Ebola virus is considered to be fruit bats of the Pteropodidae family, with monkeys, gorillas, and chimpanzees believed to further transmit the disease [4].
Ebola is transmitted through direct contact with the skin, blood or bodily fluids of an infected individual or animal and with contaminated objects [2,6]. Individuals who take care of an infected person or bury someone who has died from the disease can also acquire the virus [2]. There is evidence that health care workers and relatives may become infected after contact with a patient or patient’s bodily fluids. It must be noted that a recovered individual cannot spread the virus but the Ebola virus has been found to remain in semen for up to three months. Therefore, abstinence from sex is recommended for at least this period of time [6]. From the outbreaks since 1976, it has been observed that the Ebola virus cannot naturally transmit through air, water, or food like influenza or diarrheal diseases [2,4]. Furthermore, individuals suffering from EVD do not infect other individuals during the incubation period, which can last between two and twentyone days [4]. Common symptoms of EVD include fever, myalgia, malaise, sore throat, chest pain, red eyes, hiccups, rash, weakness, severe headaches, joint and muscle pain, diarrhea, vomiting, stomach pain, dehydration, dry and hacking cough, and loss of appetite. These symptoms typically start two days to three weeks after acquiring EVD. As the infection spreads, the body undergoes severe blood loss and coagulation abnormalities. Ultimately, the liver, kidney, and micro vascular endothelial cells (capillary walls) become infected, leading to compromise of vascular integrity. If not diagnosed and treated, death usually occurs in the second week of symptoms, and is usually caused by massive blood loss [4].
Recovery from Ebola is dependent on good supportive clinical care and the infected individual's immune response. Fortunate individuals who recover from EVD develop antibodies that last for at least 10 years [6]. These individuals may still experience weakness, fatigue, headaches, hair loss, hepatitis, sensory changes, and inflammation of organs [4].
Over the years, a few models for EVD have been studied and analyzed. Thomas E. Sutto has used an exponential fitting of data provided by the Centers of Disease Control and Prevention (CDC) and the World Health Organization (WHO) to develop formulae that best fit infection rate totals [7]. Similarly, Camacho et al. depicted the potential for large EVD outbreaks by fitting a mathematical model to time series, estimating epidemiological factors responsible for disease transmission [8]. Another significant contribution estimated parameters from daily incidence and mortality time series for the 1995 Congo Ebola outbreak [9]. Similarly, Chowell et al. used epidemic modeling to estimate the number of secondary cases generated by an index case in the absence of control interventions [10]. In another study, Chowell et al. carried out a comparative review of mathematical models of the spread and control of Ebola [11]. Legrand et al. have previously studied transmission dynamics in order to derive a relationship between the hospitalization rate and epidemic size [12].
Apart from mathematical modeling of Ebola, noteworthy nonmathematical contributions have also been made. Recently, Tambo and XiaoNong examined research, prevention, detection, and management related issues of the Ebola outbreak and reflected upon the major gaps in frontline and airport Ebola control and containment, providing possible structured opportunities to the public [13]. Lai et al. provided insight into nonmutable host cell therapeutic agents targeting different steps of the life cycle of the Ebola virus [14]. Tambo et al. proposed surveillance response systems for controlling the Ebola outbreak shedding light on the use of early warnings, critical human resources development, and methods to enhance tracking and managing challenges and urging further development in new drug discovery and vaccines [15]. Tambo also shed light on nonconventional humanitarian interventions on Ebola in another recent study [16].
None of the above mentioned models have used the aspect of two susceptible populations: highrisk and lowrisk. The model we present in this paper explains different aspects of the disease dynamics. In the following section we explore the formulation of our model. Then we explain the data fitting technique and present our results.
The purpose of our study is to estimate the basic reproductive ratio (R_{0}) for the EVD outbreak in Liberia and Sierra Leone. Using data obtained from the CDC for the period of May 1st, 2014 up until October 1st, 2014, we present a deterministic SEIR type model for the transmission dynamics of the Ebola virus to estimate R_{0}. We also present an illustration of the required disease control scenario to achieve R_{0} < 1.
Methods
Mathematical model formulation
We base our study on a deterministic ordinary differential equations (ODE) epidemic model in which the population size is divided into six mutually exclusive compartments. The total population at any time instant t, denoted by N(t), is the sum of individual populations in each compartment that includes lowrisk susceptible individuals S _{ L } (t), highrisk susceptible individuals S _{ H } (t), exposed individuals E(t), infected individuals I(t), hospitalized individuals H(t), and recovered individuals R(t) such that,
The highrisk susceptible population includes healthcare workers and providers (including all frontline workers), relatives of infected individuals, and the people involved in burial processes. The rest of the susceptible population is considered to be at a low risk of acquiring EVD.
Since there is no vertical transmission of the infection, we assume that all newborns are susceptible. The susceptible population increases at a constant recruitment rate Π (all recruited individuals are assumed to be susceptible), and p is the fraction of recruited individuals who are at a high risk of acquiring the infection. Lowrisk susceptible individuals acquire the infection at a rate λ. Furthermore, the susceptible population also decreases at the natural death rate μ. An increase in the highrisk population S _{ H } means that there are more healthcare workers and more people involved in the burial process. The rate at which infection is acquired by susceptible classes, also recognized as the force of infection, depends on the proportion of the infected and hospitalized individuals, taking into account the infectiousness of the hospitalized individuals (η). In our model, the force of infection is given by λ. The Exposed population increases after susceptible individuals acquire the infection at a rate λ from lowrisk individuals or at a rate ψ _{ H } λ from highrisk individuals respectively. The population of infected individuals generated at a rate α decreases when these individuals go to a hospital at rate τ. It also decreases because of the natural death rate μ, and the diseaseinduced death rate δ _{ I }. Infected individuals recover from disease at rate of θ _{ I }. The number of hospitalized individuals is generated when infected individuals are hospitalized at a rate τ. It diminishes when individuals recover at a rate θ _{ H }, and die naturally or due to infection at rates μ and δ _{ H }, respectively. The flow diagram of Ebola model is shown in Figure 1.
The model in this study is given by the following system of nonlinear differential equations:
where, \( \lambda =\beta \frac{\left(I+\eta H\right)}{N} \) is called the force of infection.
All model parameters are summarized in Table 1 in Section 2.4.
The basic reproductive ratio
The basic reproductive ratio (R_{0}) is the number of individuals infected by a single infected individual during the infectious period in an entirely susceptible population [17]. As our model consists of multiple classes, the next generation operator is used to calculate R_{0} [17]. Because our population consists of a hospitalized class in addition to exposed and infected classes, our final R _{ 0 } is the sum of R _{ 1 } and R _{ 2 } (resulting from the infectiousness linked to community and hospitals respectively).
We use the next generation matrix approach, as introduced by Diekmann et al., to calculate R_{0} [18]. Using differential equations associated with the exposed (E), infected (I), and hospitalized (H) compartments as stated below, we compute a function (F) for the rate of new infection terms entering, and another function (V) for the rate of transfer into and out of the exposed, infected, and hospitalized compartments by all possible means depicted in Model 1.
The matrices F (for the new infection terms) and V (of the transition terms) are given by,
where, \( {K}_3={\theta}_H+{\delta}_H+\mu, \kern1em \mathrm{and}\kern1em \varOmega \kern0.5em =\kern0.75em \frac{\Pi \left(1p\right)+{\psi}_Hp}{\mu }. \)
Reproductive ratio R_{0} is then given as
where ρ is the spectral radius (the maximum Eigen value of the matrix) and FV ^{− 1} is the next generator matrix. This leads to the following expression
Here, R _{ 1 } and R _{ 2 } reflect the continuation of infectious individuals from the community and from hospitals respectively. The epidemiological significance of the basic reproductive ratio R_{0}  which represents the average number of new cases generated by a primary infectious individual in a population where some susceptible individuals are at high risk and some infected individuals go to hospital  is that the Ebola pandemic can be effectively controlled by reducing the number of highrisk individuals and by decreasing peoples’ contact with hospitalized individuals with other individuals, be they relatives, healthcare workers, people involved in burial processes, etc. This can bring the threshold quantity (R _{ 0 }) to a value less than unity. Biologically, this implies that the Ebola pandemic can be eliminated from the population when R _{ 0 } < 1.
Data sources
The epidemic data used in this study wacollected by the WHO during the current outbreak; the data is available at http://apps.who.int/ebola/en/currentsituation/ebolasituationreport. The CDC analyzed this data and proposed that there is underreporting of the cases. They estimate that the true number of cases is 2.5 times more than the ones reported [19]. We use the raw or uncorrected data from the WHO, and then the corrected data from the CDC from May 1st, 2014 until October 1st, 2014 to estimate R_{0} for Liberia and Sierra Leone.
Parameter values
Based on previous studies, the infection period for EVD is six days. The latent period of EVD is between two and seventeen days with a mean of ten days. The generation time of the disease is about 16.6 days. So, the values are 1/θ _{ I } = 6 and 1/α = 10. The parameter η is the relative transmissibility of hospitalized individuals compared to infected individuals. The value of η < 1 would indicate that isolation in hospitals is effective. The value η > 1 would indicate ineffectiveness of transmissibility in hospitals. Such would be the situation in some developing countries. For West Africa, the value of η lies between 0.6 and 0.8. The parameter ψ _{ H } accounts for infection that arises from the individuals involved in burial processes as well as health care workers. This indicates that highrisk individuals are more likely to get infected as compared to lowrisk individuals. Its value lies between 1.5 and 2. The average time from hospitalization to recovery and from infection to recovery is five and 10 days respectively. Moreover only 45% of the infected and 50% of hospitalized individuals recover from the disease. The death rates of the infected and hospitalized individuals, δ _{ I } and δ _{ H }, are 0.10 and 0.5, respectively. Of the infected individuals, 80% are hospitalized, and on average, hospitalization lasts four to five days after getting infected. These parameter values are taken from [11,12,20,21]. The assigned values are summarized in Table 1 below. All rates are defined per day except for μ.
Estimation scheme
In order to calculate R_{0}, we use parameter values for Model 1 as stated in Table 1 in the previous section. The estimates for several of the model parameters used in model (1) have been obtained from existing studies on EVD. The effective contact rate β, which is a measure of the rate at which contact between an infected and a susceptible individual occurs, and the probability that such contact will lead to an infection, is extremely difficult to determine directly. Therefore, we adopt an indirect approach, similar to previous studies such as [22] and [23], by first finding the value of parameter β for which Model 1 has the best agreement with the epidemic data, and then using the resultant parameter values to estimate R_{0}.
Furthermore, we require knowledge of the initial conditions to be used for simulation of the trajectories of Model 1. It is possible to consider the initial conditions (S _{ L } (0), S _{ H } (0), E(0), I(0) H(0),R(0)) as model parameters, along with the effective contact rate and estimate values for all parameters. Such a technique, however, produces slightly unreliable results. This is explained by the fact that the available epidemic data is restricted to the reported cumulative case number, while the optimization scheme that we employ produces estimates for six variables.
There are thus too many degrees of freedom and the ‘bestfit’ may result in unrealistic estimates for the initial conditions used. We use the number of cases first recorded as the initial conditions and restrict ourselves to optimizing only the effective contact rate. The following initial conditions were used; S _{ H }(0) = 20000, S _{ L }(0) = one million, E(0) = 15, I(0) = 10, H(0) = 0, and R(0) = 0. This initial data is based on the fact that the total population in the region that was under threat of Ebola was one million. Among this population, 2% are at highrisk. Initially there were only a few exposed and infected individuals, and no hospitalized or recovered individuals.
In the following section we employ ordinary least squares (OLS) estimation to estimate the parameter β by minimizing the difference between predictions of Model 1 and the epidemic data. This is implemented by using the fminsearch function in the builtin optimization toolbox in MATLAB. The function fminsearch searches the minimum of a function of several variables. An initial guess of the variables is provided. Then, by using that guess, it searches for the values that are local minimizers of the function. This allows us to estimate the parameter β in order to calculate R_{0}.
Results
As described in previous section, we fitted the two different data sets, with and without correction for under reporting, for Liberia and Sierra Leone. Figure 2a shows the fit for the uncorrected Liberia data until October 1st, 2014 for which β is 0.371 and R_{0} is 1.757. Figure 2b shows the fit for the corrected Liberia data until October 1st, 2014 for which β is 0.3906 and R_{0} is 1.9.
For Sierra Leone, β is 0.361 and R_{0} is 1.492, as shown in Figure 3a for the uncorrected case. After correcting for underreporting, β becomes 0.344 and R_{0} is 1.3682, as shown in Figure 3b.
Finally, we study the variation in R_{0} with the relative risk of the highly susceptible population group and its percentage in the total susceptible population group. We note that less than 10% of the susceptible population should be in the highrisk group in order to bring R_{0} to less than 1. This is shown in Figure 4 below.
We also use our model to study the effects of isolation on R_{0}. As shown in Figure 5 we note that not only does isolation have to be very effective, reducing the infectivity to less than 0.25, but at the same time around 45% or more of the population has to be isolated in order to bring R_{0} to a value less than 1.
Discussion and conclusion
We have developed a model for the transmission dynamics of EVD, incorporating the important factor of the individuals who are at a much greater risk of contracting the disease than the general population in the affected areas. These include frontline healthcare workers, family members of EVD patients, and individuals involved in the burial process of deceased EVD patients.
Using data from the WHO and CDC, we have calculated estimates for Liberia and Sierra Leone for the ongoing EVD outbreak. The values are significantly above 1, indicating the severity of the disease. The estimated values for Liberia are consistent with published estimates for the current outbreak [11,20,24], while those for Sierra Leone are consistent with Nishiura et al., however another recent study [8] estimated that R_{0} in this country is 2.53, which is significantly greater.
We have also looked at the effect of interventions to control the outbreak. In absence of any vaccine or medication for EVD, the only control measure available is isolation. Ideally during isolation strict control should be in place so that the isolated individuals do not transmit the disease. However, in practice, there is a nonzero risk of transmission from isolated individuals. Our analysis suggests that in order for R_{0} to reduce to less than 1, the transmission rate of the isolated individuals should be less than one quarter of that for the nonisolated. This means that strict protocols should be followed at treatment facilities. Further analysis of the model also leads to the conclusion that the fraction of highrisk individuals has to be controlled and must be brought to less than 10% of the overall susceptible population in order to bring R_{0} to less than 1 and hence control the outbreak.
Our model is an attempt to capture the most important features of the transmission dynamics of EVD. As an extension of this work, optimal, timedependent strategies should be developed and advised to public health authorities in order to control the disease.
Abbreviations
 EVD:

Ebola virus disease
 CDC:

Centers for disease control and prevention
 WHO:

World Health Organization
 SEIR:

Susceptibleexposedinfectedrecovered
 ODE:

Ordinary differential equation
 OLS:

Oordinary least squares
References
 1.
Brief General History of Ebola, 29 Oct 2014 https://web.stanford.edu/group/virus/filo/history.html.
 2.
Haradhan M. The most fatal 2014 outbreak of ebola virus disease in Western Africa. Am J Epidemiol Infect Dis. 2014;2(4).
 3.
Derek G. The 2014 Ebola virus disease outbreak in West Africa. J Gene Virolog. 2014;95:1919–24.
 4.
Moran KI. What do we really fear? The epidemiological characteristics of Ebola and our preparedness. Epidemiol Health KAMJE, 2014;36.
 5.
Ebola mapping the outbreak, BBC News Africa, BBC Web 29 Oct 2014. http://www.bbc.com/news/worldafrica28755033.
 6.
Ebola (Ebola virus disease) Transmission. Cent Dis Contr Prev, 20 Nov 2014. http://www.cdc.gov/vhf/ebola/transmission/index.html.
 7.
Thomas S: Mathematical analysis of the West Africa Ebola outbreak. Res Gate. Oct 2014. doi:10.13140/2.1.2295.5527.
 8.
Camachoa A, Kucharski J, Funka S, Bremanb J, Piot P, Edmunds WJ. Potential for large outbreaks of Ebola virus disease. Epedemics. 2014;154:70–8.
 9.
Lekone PE, Finkenstdt BF. Statistical inference in a stochastic epidemic SEIR model with control intervention: Ebola as a case study. Biometrics. 2006;62(4):1170–7.
 10.
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.
 11.
Chowell G, Nishiura H. Transmission dynamics and control of Ebola virus disease (EVD): a review. BMC Med. 2014;12(196):1–12.
 12.
Legrand J, Gais RF, Boelle PY, Valleron AJ, Flahault A. Understanding the dynamics of Ebola epidemics. Epidemiol Infect. 2007;135(4):610–21.
 13.
Tambo E, XiaoNong Z. Acquired immunity and asymptomatic reservoir impact on frontline and airport ebola outbreak syndromic surveillance and response. Infect Dis Pov. 2014;3(41):1–10.
 14.
KY Lai, WYG Ng, FF Cheng: Human Ebola virus infection in West Africa: a review of available therapeutic agents that target different steps of the life cycle of Ebola virus. Infect Dis Pov. 2014;3(43).
 15.
E Tambo, EC Ugwu, JY Ngogang. Need of surveillance response systems to combat Ebola outbreaks and other emerging infectious diseases in African countries. Infect Dis Pov. 2014;3(29).
 16.
E Tambo. Nonconventional humanitarian interventions on Ebola outbreak crisis in West Africa: health, ethics and legal implications. Infect Dis Pov. 2014;3(42).
 17.
Diekmann O, Heesterbeek JAP, Metz JAJ. On the definition and the computation of the basic reproduction ratio R _{0} in models for infectious diseases. J Math Biol. 1990;35:503–22.
 18.
Heffernan JM, Smith RJ, Wahl LM. Perspectives on the basic reproductive ratio. R Soc Pub. 2005;2(4):281–93.
 19.
Meltzer NI, Atkins CY, Santibanez S, Knust B, Petersen BW, Ervin ED, et al. Estimating the future number of cases in the Ebola epidemic – Liberia and Sierra Leone. MMWR Surveill Summ. 2014;6.
 20.
Nishiura H, Chowell G. Early transmission dynmaics of eboal virus disease (EVD), West Africa, March to August 2014. Eurosurveillance. 2014;19(36).
 21.
Althaus CL. Estimating the reproduction number of Zaire ebolavirus (EBOV) during the 2014 outbreak in West Africa. PLOS Curr Outbr. 2014.
 22.
Cintron A, CastilloChavez C, Bettencourt LM, Llyod AL Banks T. The estimation of the effective reproductive Number from Disease outbreak Data. Mathematical Biosciences and Engineering. 2009;6(2).
 23.
Adnan K, Muhammad H, Imran M. Estimating the basic reproduction number for singlestrain dengue fever epidemics. J Infect Dis Pov. 2014;3(1).
 24.
David F, Edwin K, Ashleigh T: Early epidemic dynamics of the West African 2014 Ebola outbreak: estimates derived with a simple twoparameter model. PLoS. 2014.
Acknowledgements
The authors would like to acknowledge CDC for providing case data used in this research article.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
MI, AK and DA initiated and directed the research study. MI and AK implemented algorithms, conducted programming and drafted the manuscript. MN wrote introduction, collected parameters, ran simulations, analyzed results, compared results with findings of other authors, and edited the manuscript. All authors read and approved the final version of the manuscript.
Additional file
Additional file 1:
Translation of the abstract into the six official working languages of the United Nations.
Rights and permissions
About this article
Received
Accepted
Published
DOI
Keywords
 Epidemiology
 Ebola virus disease
 Transmission model
 Basic reproductive ratio