Assessing required SARS-CoV-2 blanket testing rates for possible control of the outbreak in the epicentre Lusaka province of Zambia with consideration for asymptomatic individuals: a simple mathematical modelling study ========================================================================================================================================================================================================================== * Joseph Sichone * Musalula Sinkala * Mervis Kikonko * Sody M. Munsaka * Martin Simuunza ## Abstract **Introduction** The novel Coronavirus disease (COVID-19), caused by the severe acute respiratory syndrome coronavirus - 2 (SARS-CoV-2), in Africa is characterised by a more substantial proportion of asymptomatic (or mildly symptomatic) individuals thought to be playing a role in the spread of the infection. The exact proportion and degree of infectiousness of asymptomatic individuals remains unclear. Studies however indicate that their management is crucial for control of SARS-CoV-2 transmission. **Methodology** We developed a simplified deterministic susceptible-exposed-infectious-removed (SEIR) mathematical model to assess the effect of active isolation of SARS-CoV-2 infected but asymptomatic individuals through blanket testing for control of the outbreak in Lusaka Province of Zambia. Here we modelled two scenarios; (1) assuming asymptomatic individuals comprised 70% of all COVID-19 cases and (2) asymptomatic individuals comprised only 50% of the cases. For contrast, the model was assessed first under the assumption that asymptomatic individuals are equally as infectious as symptomatic individuals and then secondly, and more likely, assuming asymptomatic individuals are only half as infectious as symptomatic individuals. **Results** For the model assuming 70% asymptomatic cases, a minimum sustained blanket testing rate of ≥ 7911 tests/100000 population was sufficient to control the outbreak if asymptomatic individuals are only half as infectious while if equal infectiousness was assumed then a testing rate of ≥ 10028 tests/ 100000 population would be required. For 50% asymptomatic, minimum blanket testing rates of ≥ 4540 tests/ 100000 population was sufficient to control the outbreak at both assumed levels of infectiousness for asymptomatic individuals relative to symptomatic individuals. **Discussion and conclusion** Our model predicts that the current testing rates of ≈ 150/100,000 population are inadequate to control transmission of SARS-Cov-2 in Lusaka. Active isolation of COVID-19 cases including asymptomatic individuals through blanket testing can be used as a possible measure for control of the SARS-Cov-2 transmission in Lusaka, Zambia. Key words * SARS-Cov-2 * asymptomatic transmission * deterministic model * blanket testing ## Introduction Since the first reported case on 31 December 2019 in China, the current pandemic of the novel Coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus - 2 (SARS-CoV-2), has killed over 327,700 and infected over 4.9 million people globally by 21st May, 2020 [52]. The disease is characterised by a more substantial proportion of asymptomatic (or mildly symptomatic) individuals thought to be playing an “in-dismissible” role in the spread of the infection [2,3,4,24,25,33,64]. The degree of infectiousness of the asymptomatic individuals remains unclear although recent data from the World Health Organisation (WHO) and others suggest that relative to individuals with symptomatic infections, those with asymptomatic infection are considerably less infectious. [57,61,62]. Nonetheless, various studies, including modelling studies, have suggested that asymptomatic individuals comprise about 50 - 80% of all COVID-19 cases and may be responsible for as much as 40 – 73% of new infections [3,4,53,55,56]. Therefore, their management is critical for control of the disease [3,4,33]. To give a striking example, identification and isolation of asymptomatic people through blanket testing helped eliminate the virus in a completely isolated village of about 3000 people in northern Italy which saw the number of people with COVID-19 symptoms fall by over 90% within ten days [3]. The current study applied a simple mathematical modelling approach to explore the effect of increased blanket testing rates as a possible measure to capture and isolate asymptomatic individuals and control the COVID-19 outbreak in Lusaka Province of Zambia which is the epicenter of the outbreak in Zambia since the first recorded case in the country on 18th March 2020. [44]. Recent studies have modelled the spread and expected burden of the COVID-19 outbreak in Africa and Zambia and explored the effects of various control measures such as applying different levels of physical distancing and shielding in the population [42,43,51,72]. Although this provides vital information to guide policy for Zambia, some interventions may not be easy to monitor in practice. Additionally, while such interventions have already been instituted, cases continue to rise in Zambia and other African countries [20,43,44,50]. Assessment of more COVID-19 control options through mathematical modelling based on the known epidemiology of the disease would therefore serve to supplement current information on the possible management of the outbreak in Zambia. ## Materials and Methods ### Study Area Lusaka Province is the smallest and highly urbanised province in Zambia (83.5% urbanisation) with seven districts over an area of about 23,490 km2 [45]. It is one of the most densely populated provinces in the country [46,47]. It has a total population of about 3,308,438 - Density: 140.8/km2 (2019) [45]. The provincial capital, which is also the capital of Zambia, is the highly-populated Lusaka district with the latest population estimated at 2,627,700 – Density: 6,288/km2 (2019) [45]. Lusaka is a busy corporate and commercial hub of Zambia and an outlet to the rest of the world with the busy Kenneth Kaunda International Airport. It is therefore no surprise that the first recorded COVID-19 cases in Zambia occurred in Lusaka district as imported cases [44]. The province also shares borders with neighboring Zimbabwe and Mozambique [47]. ### Lusaka Province COVID-19 Outbreak On 18th March 2020, Zambia recorded the first confirmed cases of COVID-19 in Lusaka from two residents who had previously travelled to France [44]. By 21st May 2020, Lusaka had recorded about 291 confirmed cases and 6 deaths [44]. From the onset of the outbreak in Lusaka, the Zambian Ministry of Health and other stakeholders had implemented required preparedness and disease control measures including mandatory physical distancing, surveillance and case notification, heightened sanitation and handwashing in public places, closure of some public institutions as well as general sensitisation [44]. The main laboratory testing approach employed has been targeted testing based on prescribed case definitions and contact tracing to optimise positivity rates (Zambia COVID-19 situation reports No. 1-14) [44]. A few mass testing campaigns have also been conducted in some areas of Lusaka province as by 21st May 2020 e.g. Rhodes park (about 1,190 persons tested), Chirundu, (about 1,000 persons tested), and Kafue (undetermined) (Zambia COVID-19 situation reports No. 29,44-46,62-64) [44]. ### The Model The spread of SARS-Cov-2 in Lusaka province was modelled through a simplified deterministic susceptible-exposed-infectious-removed (SEIR) compartmental mathematical model as shown in Fig. 1. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/14/2020.07.12.20152124/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/07/14/2020.07.12.20152124/F1) Figure 1. The SEIR model for spread of COVID-19 in Lusaka province. The solid arrows represent flow between compartments and dotted arrows represent additive contribution relationship to a compartment. Infectious individuals are split into symptomatic and asymptomatic individuals. The model both directly and indirectly incorporated the current mitigation measures in place to attempt to predict the trajectory of the outbreak in Lusaka accurately. We started by denoting the infection states as total number of susceptible *S(t)*, exposed *E(t)*, infectious *I(t)*, removed *R(t)* and dead persons *D(t)* at any given time *(t)* in the population of size *N*. For our analysis, the total population size was assumed to be constant and demographics of natural birth and deaths rates were considered negligible [17,70]. Table 1 shows the average values of parameters used and equations in set 1 and set 2 give the systems of ordinary differential equations (ODEs) describing the flow of individuals in the model. The ODEs in set-1 describe spread of infection under the assumption of equal infectiousness for both asymptomatic and symptomatic individuals while the ODEs in set-2 describe the spread of infection assuming asymptomatic individuals are only half as infectious as symptomatic individuals. For the ODEs in set-1, we assume that susceptible individuals are infected at the rate ![Graphic][1] where *I**T*(*t*) is the total number of active infectious individuals in the population both symptomatic *I**S*(*t*) and asymptomatic *I**A*(*t*). On the other hand, for the ODEs in set-2, the susceptible individuals are infected at the rate ![Graphic][2] with the probability of infection from asymptomatic individuals only half that of symptomatic individuals. In both equations, after infection occurs the susceptible individuals are exposed *E*(*t*) and enter an incubation period of about 4.8 days before they become infectious *I*(*t*). Once they have become infectious, they may belong to either one of two classes of infectious individuals; symptomatic infectious (*I**S*), or asymptomatic infectious (*I**A*) determined by the fraction for symptomatic persons ɛ. Based on current interventions in Lusaka, the symptomatic infectious individuals are effectively managed in that once a person becomes symptomatic in the community, they are identified through targeted testing under the current testing criteria and quarantined. However, in our model we assumed an average delay of about 2.5 days in the identification of these symptomatic individuals because sometimes people delay in seeking medical attention when symptoms are unclear, may be in self-denial, plus have an extra 24hrs turnaround time delay in laboratory results [44]. The symptomatic infectious individuals *I**S**Q*(*t*) are therefore quarantined at a rate *π* (1/diagnosis delay). It is taken that during this delay period these symptomatic individuals can infect others but are no longer able to do so once isolated in quarantine. The quarantined infectious individuals recover at rate *γ* and are removed as *I**S**R*(*t*). On the other hand, in our model the asymptomatic (or mildly symptomatic) infectious individuals (*I**A*) are taken to be generally unnoticed in the community but also recover at a rate *γ* in which time they can infect susceptible individuals before they become removed as *I**A**R*(*t*) and no longer infectious. Note that some asymptomatic individuals develop symptoms much later in their infection but this does not substantially affect our model because at that time they would still be removed (*I**A**R*) if they become diagnosed and quarantined through targeted testing. In this study, we assumed that the current mostly targeted testing for COVID-19 in Lusaka province is restrictive and probably missing some asymptomatic individuals [53]. Therefore, a parameter (*θ*) was introduced in the model which describes blanket testing applied as tests per 100000 populations used to identify and isolate all infectious individuals in the community (symptomatic and asymptomatic) through sustained random mass testing. The total removed individuals for the model *R*(*t*) are given as *I**S**R*(*t*) + *I**A**R*(*t*) while the total confirmed cases are given as *I**S**Q*(*t*) + *I**S**R*(*t*) of which a fraction *α* (Case fatality rate -CFR) are recorded dead *D*(*t*). View this table: [Table 1.](http://medrxiv.org/content/early/2020/07/14/2020.07.12.20152124/T1) Table 1. Parameters used for the SEIR model for spread of COVID-19 in Lusaka Province ### Model optimisation and simulation Model optimisation and simulation was done using Vensim PLE systems dynamics modelling software for Windows (version-7) [48]. This was done for two scenarios of 70% and a modest 50% assumed proportion of asymptomatic infectious individuals in the population. Data from the first three months of the outbreak in Lusaka as given in the Zambia COVID-19 situation reports No. 1-64 [44] was used to configure the model and optimise parameters. However, due to presence of imported cases in the early days and the considerable variations in recorded cases between some days (probably influenced by variations in availability of testing kits), only data from 10th April 2020 to 16th May 2020 was used. This is because this period had more consistent data and by then community infections had been established [44]. The model initial conditions were estimated from the available data as follows: *S (0)* = *N*-92, *E (0)* = 24, *I**S* *(0)* = 9, *I**S*Q *(0)* = 38, *I**S**R (0)* = 0, *I**A* *(0)* = 21, *I**A**R (0)* = 0, *R (0)* = 0, and *D (0)* = 0. With other parameter values fixed, the model was calibrated to the cumulative number of confirmed cases over time. This was done by adjusting values of R until the best model fit was achieved (since R was expectedly affected by the current mitigation measures). R estimates the average number of secondary infections arising from a single infectious individual in a naive population. Model fit was statistically evaluated using Spearman correlation and chi-square goodness of fit test at a significance level of 0.05. For the calibration, *θ* was approximated at 150 tests/100000 populations – reflective of the mass testing rates achieved in this period [44]. After calibration, the model simulation was extended to 631 days (10th April 2020 – 31st December 2021) to predict the spread of the outbreak in Lusaka under the current transmission rate. To assess which value of *θ* would sufficiently flatten the curves of both the total number of active infectious individuals and the cumulative number of confirmed cases over time (as a key indicator of control of the outbreak), several iterations of this simulation were then performed using increasingly higher values for *θ*. ### Ethical Considerations No ethical issues were encountered as no human or animal subjects were used in this study and cases were anonymous. ## Results The model had significant fit to outbreak data under all the assessed conditions and therefore could be used for the general purpose of analysing the outbreak under all these general scenarios. Figure 2 shows results of model optimisation and fit to outbreak data for both the 70% and 50% asymptomatic scenarios assuming equal infectiousness for asymptomatic and symptomatic individuals (denoted as A1 and A2). Figure 3 shows the model fit results when asymptomatic individuals are only half as infectious as symptomatic individuals (B1 and B2). All R values were within the range estimated for COVID-19 [1,2,12,16,17,21,27,28,32,33,37]. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/14/2020.07.12.20152124/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/07/14/2020.07.12.20152124/F2) Figure 2. Model optimisation and fit to outbreak data for both the 70% (A1) and 50% (A2) proportion asymptomatic individual’s scenarios under the assumption of equal infectiousness for asymptomatic and symptomatic individuals. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/14/2020.07.12.20152124/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/07/14/2020.07.12.20152124/F3) Figure 3. Model optimisation and fit to outbreak data for both the 70% (B1) and 50% (B2) proportion asymptomatic individual’s scenarios under the assumption that asymptomatic individuals are only half as infectious as symptomatic individuals At both assessed levels of infectiousness for asymptomatic individuals, Table 2 shows that the model rightly predicts slightly lower expected peak values for the cumulative number of confirmed cases if there is more (70%) of asymptomatic individuals in the community who are poorly detected. However, this more substantial proportion of asymptomatic individuals (70%) is associated with higher expected peak values of both total active infectious individuals and total quarantined individuals. It is also observed that the outbreak is expected to peak much earlier for all variables in Table 2 for the 70% asymptomatic individual’s scenario compared to only 50% asymptomatic. Note that the model predictions for the Lusaka COVID-19 outbreak given in Table 2 are however subject to the effectiveness of the containment policies in Lusaka province over time. Figures 4 – 7 and Table 3 give the results of the effect of increasing blanket testing rates (θ) on control of the outbreak (flattening the curve for the expected cumulative number of confirmed cases and expected total number of active infectious individuals over time) for the assessed model scenarios. If asymptomatic individuals make up 70% of all COVID-19 cases and we assume equal infectiousness for asymptomatic and symptomatic individuals, then a minimum blanket testing rate of about ≥10028/100000 would be sufficient to flatten the curve for both the expected cumulative number of confirmed cases (Fig4 A1.1) and the total number of active infectious individuals (Fig4 A1.2) by more than 90% (Table 3). However, if the asymptomatic individuals only make up 50% of all cases and are also as infectious as symptomatic individuals, then a minimum blanket testing rate of ≥4540/100000 would be sufficient to flatten the curve for both the expected cumulative number of confirmed cases (Fig5 A2.1) and the total number of active infectious individuals (Fig5 A2.2) by more than 90% (Table 3). Alternatively, if asymptomatic individuals make up 70% of all COVID-19 cases but are only half as infectious as symptomatic individuals, then a minimum blanket testing rate of about ≥ 7911/100000 would be sufficient to flatten the curve for both the expected cumulative number of confirmed cases (Fig6 B1.1) and the total number of active infectious individuals (Fig6 B1.2) by more than 90% (Table 3). However, if the asymptomatic individuals only make up 50% of all cases and are only half as infectious relative to the symptomatic individuals, then a minimum blanket testing rate of ≥ 4383/100000 would be sufficient to flatten the curve for both the expected cumulative number of confirmed cases (Fig7 B2.1) and the total number of active infectious individuals (Fig7 B2.2) by more than 90% (Table 3) which is comparable to the result for A2.1 and A2.2. View this table: [Table 2.](http://medrxiv.org/content/early/2020/07/14/2020.07.12.20152124/T2) Table 2. The model predicted key outbreak variables for the Lusaka COVID-19 outbreak for period 10th April 2020 – 31st December 2021 under the different model assumptions. View this table: [Table 3.](http://medrxiv.org/content/early/2020/07/14/2020.07.12.20152124/T3) Table 3. Approximated percentage reduction in the predicted cumulative number of confirmed cases and total number of active infectious individuals for higher values of θ (blanket testing rate). ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/14/2020.07.12.20152124/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/07/14/2020.07.12.20152124/F4) Figure 4. The estimated reduction in the cumulative number of confirmed cases (A1.1) and the total number of active infectious individuals (A1.2) at different values of θ (blanket testing rate) for the 70% asymptomatic scenario assuming equal infectiousness for asymptomatic and symptomatic individuals. ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/14/2020.07.12.20152124/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2020/07/14/2020.07.12.20152124/F5) Figure 5. The estimated reduction in the cumulative number of confirmed cases (A2.1) and the total number of active infectious individuals (A2.2) at different values of θ (blanket testing rate) for the 50% asymptomatic scenario assuming equal infectiousness for asymptomatic and symptomatic individuals. ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/14/2020.07.12.20152124/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2020/07/14/2020.07.12.20152124/F6) Figure 6. The estimated reduction in the cumulative number of confirmed cases (B1.1) and the total number of active infectious individuals (B1.2) at different values of θ (blanket testing rate) for the 70% asymptomatic scenario assuming asymptomatic individuals are only half as infectious as symptomatic individuals. ![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/14/2020.07.12.20152124/F7.medium.gif) [Figure 7.](http://medrxiv.org/content/early/2020/07/14/2020.07.12.20152124/F7) Figure 7. The estimated reduction in the cumulative number of confirmed cases (B2.1) and the total number of active infectious individuals (B2.2) at different values of θ (blanket testing rate) for the 50% asymptomatic scenario assuming asymptomatic individuals are only half as infectious as symptomatic individuals ## Discussion In this study, we developed a simple deterministic model to forecast the spread of infection and assess required blanket testing rates for the control of the novel SARS-Cov-2 outbreak in Lusaka province, Zambia, with specific consideration for asymptomatic infectious individuals. In the early stages of novel infectious diseases, stochastic models are used due to uncertainty in parameter values. However, deterministic models which require average values can be used after the pandemic has progressed as is the case in this study and other studies that have used deterministic models to study COVID-19 [17,64]. In the context of Africa, earlier modelling studies have given gloomy predictions of the pandemic for the continent in case of failed early containment [42,43,50,51,53] with Zambia projected to have more than 2.8 million total infections [53] or between 4.8 million and 5.7 million total symptomatic cases [43] in the first year with estimated outbreak peaks as early as July 2020 [42,43, 53]. Further, the model in [43] also predicted that Zambia would record between 2.5 million – 4.4 million total number of symptomatic cases over a 12 months even under various levels of physical distancing and shielding interventions applied at country level. However, early implementation of the outbreak containment measures in Zambia, and most African countries, appears to be substantially averting the given predictions [53,44]. This may be because Africa was on high alert and prepared for the virus given the weaker health systems on the continent hence early containment of the disease significantly slowed the pandemic. Further, other studies have pointed to favourable bio-socioecological factors on the African continent leading to the lower rate of transmission in the region with fewer expected cases and deaths compared to other countries such as USA and Italy [53,64,65]. With limited incorporation of current containment measures implemented in Zambia, the highest number of cumulative confirmed cases predicted for Lusaka province in our model was 772,161 cases by September 2021 (Table 2). This lower projection is plausible given that it is an estimate for Lusaka province alone and it conforms to the relatively slow progression of the outbreak in Zambia compared to earlier predictions [44,53,64,65]. Further, our model predicts that the outbreak in Lusaka would spread much faster and peak earlier if asymptomatic individuals make up 70% of all COVID-19 cases (earliest peaks October - November 2020 for total active infectious individuals). This is because a higher proportion of poorly detected asymptomatic individuals inflated the population of active infectious individuals in the population (*I**S* + *I**A* = *I**T*) thereby increasing the infection rates in the community; ![Graphic][3] (ODEs in set-1) or ![Graphic][4] (ODEs in set-2), resulting in a faster progression of the outbreak. The higher proportion of asymptomatic individuals was, however, associated with less overall cumulative number of confirmed cases since asymptomatic individuals are generally less diagnosed under targeted testing approaches such as those employed in Lusaka province [53]. In our model, we refer as asymptomatic all “silent spreaders” of COVID-19 which includes all those that have become infectious but never develop symptoms (asymptomatic or mildly symptomatic) and those that have become infectious but only develop symptoms much later in their infection (pre-symptomatic spreaders) [53,55-59]. Studies however agree that collectively these asymptomatic individuals (or silent spreaders) are less infectious compared to symptomatic individuals [57,61,62]. This is most likely because their viral shedding rate (how much virus an infected person releases) may be limited without the symptoms of coughing and sneezing which produces infectious respiratory droplets even if these asymptomatic individuals may have similar viral load, duration of viral shedding, and contact rates in the community as symptomatic individuals [2,60]. Therefore, we suspect that transmission of the infection from asymptomatic individuals may mostly rely on the other possible routes of transmission for COVID-19 including contact with fomites, unwitting close intimate facial contact, and possibly even fecal-oral route [53,62,63]. Even with a lowered ability to transmit the infection, lack of detection and awareness of the carrier state may cause the asymptomatic individuals to be playing an important role in the spread of COVID-19 and sustaining the outbreak [2,3,4,24,25,33]. Active isolation of all infected individuals (symptomatic and asymptomatic) through blanket testing therefore offers a possible solution to this challenge [3,4,53]. In our study, we assessed the effect of various levels of blanket testing for control of the COVID-19 outbreak in Lusaka province under different assumed scenarios of the proportion and infectiousness of asymptomatic individuals with results tabulated in Figures 4 – 7 and Table 3. Out of all the assessed scenarios, we nominate scenario B1 to represent the most likely epidemiological dynamics for COVID-19 in the population where asymptomatic individuals make up 70% of cases but are only half as infectious as symptomatic individuals. Therefore, for this scenario, a blanket testing rate of about ≥ 7911 tests /100000 (≈ 80/1000) would be sufficient to control the outbreak in Lusaka. To compare with other possible outbreak control strategies, note that the model fitted value of R for the selected scenario B1 was 3.53 (Figure 3 B1). Therefore, in case of available vaccine for COVID-19, a minimum protective vaccination rate of about 1 − 1/*R* (≈70%) would be required to control the outbreak in the population; assuming vaccination provides full protection [66-69]. This proportionately translates into roughly 70000 vaccinations /100000 populations which is much higher than the estimated required blanket testing intervention. The effect of blanket testing is that it acts to directly reduce the size of the total number of active infectious individuals which leads to a disproportionately higher reduction in the infection rate in the population. This is achieved through active isolation of even the otherwise largely un-detected asymptomatic individuals in the population who possibly participate in the transmission of the infection [2,3,4,24]. In the case of Africa, control of asymptomatic infections may be even more pertinent than other regions of the world. This is because Africa has been found to have a younger population and with lower personal vulnerabilities to severe infection compared to Europe, USA, and other regions leading to most of the infections in Africa being asymptomatic [53,64,65]. Therefore, more outbreak control strategies with specific consideration of the role of asymptomatic individuals will be required in Africa. ## Conclusion Active isolation of SARS-Cov-2 asymptomatic infectious individuals through blanket testing has been shown to have indeed the potential to control the outbreak in Lusaka province of Zambia. However, this can only be achieved at a rate of approximately ≥ 7911 tests /100000 (≈ 80/1000) following the most likely epidemiological characteristics of the COVID-19 outbreak in Lusaka. ## Limitations As is the case in creating infectious disease models, appropriateness of the model design to incorporate containment measures implemented in Lusaka and quality of data used were some of the challenges faced in this study [70,71]. Also, limited outbreak data was used due to significant irregularities in the reported outbreak data. Cases were not reported in real-time and there is a possibility of under-reporting due to low testing capacities in the earlier stages of the outbreak in Zambia [44,54]. Additionally, the model was fitted to cumulative confirmed cases which could still introduce uncertainties [54]. Further modelling studies using more refined models and more outbreak data should, therefore, be conducted to study the COVID-19 outbreak in Lusaka. ## Data Availability The data underlying the results presented in the study are available from the Zambia MOH/ZNPH/WHO, | New Coronavirus Disease of 2019 (COVID-19) Situation Reports [1-64]. March-May 2020. Accessed from http://znphi.co.zm/news/situation-reports-new-coronavirus-COVID-19-sitreps/. [http://znphi.co.zm/news/situation-reports-new-coronavirus-COVID-19-sitreps/](http://znphi.co.zm/news/situation-reports-new-coronavirus-COVID-19-sitreps/). ## Abbreviations COVID-19 : Novel Coronavirus disease SARS-CoV-2 : Severe acute respiratory syndrome Coronavirus 2 CFR : Case fatality rate SEIR : susceptible-exposed-infectious-removed compartmental mathematical model ODEs : Ordinary differential equations * Received July 12, 2020. * Revision received July 12, 2020. * Accepted July 14, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1).Yi Y, Lagniton PN, Ye S, Li E, Xu RH. COVID-19: what has been learned and to be learned about the novel coronavirus disease. Int J Biol Sci. 2020;16(10):1753. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7150/ijbs.45134&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F14%2F2020.07.12.20152124.atom) 2. 2).Singhai T. A review of the coronavirus disease-2019. Indian J Pediatr. 2020;87:281–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12098-020-03263-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F14%2F2020.07.12.20152124.atom) 3. 3).Day M. Covid-19: identifying and isolating asymptomatic people helped eliminate virus in Italian village. BMJ: (Online). 2020 Mar 23;368. 4. 4).Day M. Covid-19: four fifths of cases are asymptomatic, China figures indicate. BMJ 2020; 369: m1375 [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE4OiIzNjkvYXByMDJfMTEvbTEzNzUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNy8xNC8yMDIwLjA3LjEyLjIwMTUyMTI0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 5. 5).Linton NM, Kobayashi T, Yang Y, Hayashi K, Akhmetzhanov AR, Jung SM, Yuan B, Kinoshita R, Nishiura H. Incubation period and other epidemiological characteristics of 2019 novel coronavirus infections with right truncation: a statistical analysis of publicly available case data. J Clin Med. 2020 Feb;9(2):538 6. 6).Hellewell J, Abbott S, Gimma A, Bosse NI, Jarvis CI, Russell TW, Munday JD, Kucharski AJ, Edmunds WJ. Centre for the Mathematical Modelling of Infectious Diseases COVID- 19 Working Group. Funk S, Eggo RM. Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts. Lancet Glob Health. 2020 Apr; 8 (4): e488–e496. doi: 10.1016/S2214-109X(20)30074-7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2214-109X(20)30074-7&link_type=DOI) 7. 7).Huang Y, Yang L, Dai H, Tian F, Chen K. Epidemic situation and forecasting of COVID- 19 in and outside China. Bull World Health Organ. 2020 Mar 16;10. 8. 8).Lin Q, Zhao S, Gao D, Lou Y, Yang S, Musa SS, Wang MH, Cai Y, Wang W, Yang L, He D. A conceptual model for the outbreak of Coronavirus disease 2019 (COVID-19) in Wuhan, China with individual reaction and governmental action. Int J Infect Dis. 2020 Mar 4. 9. 9).Hamzah FB, Lau C, Nazri H, Ligot DV, Lee G, Tan CL. CoronaTracker: worldwide COVID-19 outbreak data analysis and prediction. Bull World Health Organ. 2020 Mar;1:32. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2471/BLT.20.255695&link_type=DOI) 10. 10).Kucharski AJ, Russell TW, Diamond C, Liu Y, Edmunds J, Funk S, Eggo RM, Sun F, Jit M, Munday JD, Davies N. Early dynamics of transmission and control of COVID-19: a mathematical modelling study. Lancet Infect Dis. 2020 Mar 11. 11. 11).Du Z, Xu X, Wu Y, Wang L, Cowling BJ, Meyers LA. Serial interval of COVID-19 among publicly reported confirmed cases. Emerg Infect Dis. 2020 Jun;26(6):1341. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F14%2F2020.07.12.20152124.atom) 12. 12).Anastassopoulou C, Russo L, Tsakris A, Siettos C. Data-based analysis, modelling and forecasting of the novel coronavirus (2019-Ncov) outbreak. medRxiv. 2020 Jan 1. 13. 13).Singh R, Adhikari R. Age-structured impact of social distancing on the COVID-19 epidemic in India. arXiv preprint arXiv:2003.12055. 2020 Mar 26. 14. 14).Peng L, Yang W, Zhang D, Zhuge C, Hong L. Epidemic analysis of COVID-19 in China by dynamical modeling. arXiv preprint arXiv:2002.06563. 2020 Feb 16. 15. 15).Giordano G, Blanchini F, Bruno R, Colaneri P, Di Filippo A, Di Matteo A, Colaneri M. Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy. Nat Med. 2020 Apr 22:1–6. 16. 16).Choi S, Ki M. Estimating the reproductive number and the outbreak size of COVID-19 in Korea. Epidemiol Health. 2020 Mar 12;42. 17. 17).Prem K, Liu Y, Russell TW, Kucharski AJ, Eggo RM, Davies N, Flasche S, Clifford S, Pearson CA, Munday JD, Abbott S. The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study. Lancet Public Health. 2020 Mar 25. 18. 18).Bastola A, Sah R, Rodriguez-Morales AJ, Lal BK, Jha R, Ojha HC, Shrestha B, Chu DK, Poon LL, Costello A, Morita K. The first 2019 novel coronavirus case in Nepal. Lancet Infect Dis. 2020 Mar 1;20(3):279–80. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30067-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F14%2F2020.07.12.20152124.atom) 19. 19).Wang D, Hu B, Hu C, Zhu F, Liu X, Zhang J, Wang B, Xiang H, Cheng Z, Xiong Y, Zhao Y. Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus– infected pneumonia in Wuhan, China. JAMA. 2020 Mar 17;323(11):1061–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.1585&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F14%2F2020.07.12.20152124.atom) 20. 20).Ritchie H, Ortiz-Ospina E, Beltekian D, Mathieu E, Hasell J, Macdonald B, Giattino C, Roser M. Our World in data: Coronavirus Pandemic (COVID-19) – the data. 2020. Available from: [https://ourworldindata.org/coronavirus-data](https://ourworldindata.org/coronavirus-data) 21. 21).Rothan HA, Byrareddy SN. The epidemiology and pathogenesis of coronavirus disease (COVID-19) outbreak. J Autoimmun. 2020 Feb 26:102433. 22. 22).Wu P, Hao X, Lau EH, Wong JY, Leung KS, Wu JT, Cowling BJ, Leung GM. Real-time tentative assessment of the epidemiological characteristics of novel coronavirus infections in Wuhan, China, as at 22 January 2020. Euro Surveill. 2020 Jan 23;25(3):2000044. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2807/1560-7917.ES.2020.25.3.2000044&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31992388&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F14%2F2020.07.12.20152124.atom) 23. 23).Rodriguez-Morales AJ, Cardona-Ospina JA, Gutiérrez-Ocampo E, Villamizar-Peña R, Holguin-Rivera Y, Escalera-Antezana JP, Alvarado-Arnez LE, Bonilla-Aldana DK, Franco-Paredes C, Henao-Martinez AF, Paniz-Mondolfi A. Clinical, laboratory and imaging features of COVID-19: A systematic review and meta-analysis. Travel Med Infect Dis. 2020 Mar 13:101623. 24. 24).Weiss P, Murdoch DR. Clinical course and mortality risk of severe COVID-19. The lancet. 2020 Mar 28;395(10229):1014–5. 25. 25).Lai CC, Liu YH, Wang CY, Wang YH, Hsueh SC, Yen MY, Ko WC, Hsueh PR. Asymptomatic carrier state, acute respiratory disease, and pneumonia due to severe acute respiratory syndrome coronavirus 2 (SARSCoV-2): facts and myths. J Microbiol Immunol Infect. 2020 Mar 4. 26. 26).Yliu ZF. ZXiang JWang YSong BGu X et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. The lancet. 2020;395:1054–62. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/SO140-6736(20)30566-3&link_type=DOI) 27. 27).Velavan TP, Meyer CG. Mild versus severe COVID-19: laboratory markers. Int J Infect Dis. 2020 Apr 25. 28. 28).Lai CC, Shih TP, Ko WC, Tang HJ, Hsueh PR. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and corona virus disease-2019 (COVID-19): the epidemic and the challenges. Int J Antimicrob Agents. 2020 Feb 17:105924. 29. 29).Surveillances V. The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19)—China, 2020. China CDC Weekly. 2020;2(8):113–22. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.46234/ccdcw2020.032&link_type=DOI) 30. 30).Tuite AR, Fisman DN, Greer AL. Mathematical modelling of COVID-19 transmission and mitigation strategies in the population of Ontario, Canada. CMAJ. 2020 May 11;192(19):E497–505. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY21haiI7czo1OiJyZXNpZCI7czoxMToiMTkyLzE5L0U0OTciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNy8xNC8yMDIwLjA3LjEyLjIwMTUyMTI0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 31. 31).Nishiura H, Linton NM, Akhmetzhanov AR. Serial interval of novel coronavirus (COVID- 19) infections. Int J Infect Dis. 2020 Mar 4. 32. 32).Zhang Y, Xu J, Li H, Cao B. A novel coronavirus (COVID-19) outbreak: a call for action. Chest. 2020 Apr 1;157(4):e99–101. 33. 33).Wang Y, Wang Y, Chen Y, Qin Q. Unique epidemiological and clinical features of the emerging 2019 novel coronavirus pneumonia (COVID-19) implicate special control measures. J Med Virol. 2020 Jun;92(6):568–76. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F14%2F2020.07.12.20152124.atom) 34. 34).Perlman S. Another decade, another coronavirus. NEJM. 2020. [CrossRef] [PubMed] 35. 35).Chan JF, Yuan S, Kok KH, To KK, Chu H, Yang J, Xing F, Liu J, Yip CC, Poon RW, Tsoi HW. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. The Lancet. 2020 Feb 15;395(10223):514–23. 36. 36).Suxin W, Yi X, Wei F, Yu Z, Boqun L, Yanjun H, Chunhui L, Daoqiu H, Qiuyan S, Yan X, Xia H. Clinical features and treatment of COVID-19 patients in northeast Chongqing. [https://onlinelibrary.wiley.com/doi/pdfdirect/10.1002/](https://onlinelibrary.wiley.com/doi/pdfdirect/10.1002/). J Med Virol. p25783. 2020. 37. 37).Liu Y, Gayle AA, Wilder-Smith A, Rocklöv J. The reproductive number of COVID-19 is higher compared to SARS coronavirus. J Travel Med. 2020 Mar 13. 38. 38).Zhu N, Zhang D, Wang W. China Novel Coronavirus Investigating and Research Team. A novel coronavirus from patients with pneumonia in China, 2019 [published January 24, 2020]. NEJM. 39. 39).Baud D, Qi X, Nielsen-Saines K, Musso D, Pomar L, Favre G. Real estimates of mortality following COVID-19 infection. Lancet Infect Dis. 2020 Mar 12. 40. 40).Lauer SA, Grantz KH, Bi Q, Jones FK, Zheng Q, Meredith HR, Azman AS, Reich NG, Lessler J. The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: estimation and application. Ann Intern Med. 2020 May 5;172(9):577–82. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-0504&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32150748&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F14%2F2020.07.12.20152124.atom) 41. 41).Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, Zhang L, Fan G, Xu J, Gu X, Cheng Z. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. The Lancet. 2020 Feb 15;395(10223):497–506. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0146736(20)30183-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F14%2F2020.07.12.20152124.atom) 42. 42).MOBS Laboratory, Northeastern University, Boston. Center for Inference and Dynamics of Infectious Diseases (CIDID). Fred Hutchinson Cancer Research Center, Seattle.: Modeling local transmission and the burden of the COVID-19 epidemic in African countries. Report 2.0, v1.0. April 12, 2020. 43. 43).Nightingale ES, Munday JD, Medley G, Gibbs HP, Abbott S, Houben RMGJ, et al. Modelling projections for COVID-19 epidemic in Zambia. LSHTM CMMID COVID19 Working Group (2020). 44. 44).MOH/ZNPH/WHO, Situation Reports | New Coronavirus Disease of 2019 (COVID-19) Sitreps [1-64]. March-May 2020. Accessed from [http://znphi.co.zm/news/situation-reports-new-coronavirus-COVID-19-sitreps/](http://znphi.co.zm/news/situation-reports-new-coronavirus-COVID-19-sitreps/). 45. 45).City population. Population statistics for countries, administrative divisions, cities, urban areas and agglomerations – interactive maps and charts. Lusaka Province in Zambia. 2020. Available from: [https://www.citypopulation.de/php/zambia-admin.php?adm1id=05](https://www.citypopulation.de/php/zambia-admin.php?adm1id=05). 46. 46).ZAMBIA DATA PORTAL. Central Statistical Office. Lusaka. 2020. Available from: [https://zambia.opendataforafrica.org/apps/atlas/Lusaka](https://zambia.opendataforafrica.org/apps/atlas/Lusaka) 47. 47).Wikipedia. The Free Encyclopedia. Lusaka Province. Available from: 2020: [https://en.wikipedia.org/wiki/Lusaka\_Province](https://en.wikipedia.org/wiki/Lusaka_Province) 48. 48).Vensim. VENTANA systems.inc. 2020. Available from: [https://vensim.com/vensim-personal-learning-edition/](https://vensim.com/vensim-personal-learning-edition/) 49. 49).Ritchie H, Roser M. What do we know about the risk of dying from COVID-19? Our World in Data. March 25, 2020. Available from: [https://ourworldindata.org/covid-mortality-risk](https://ourworldindata.org/covid-mortality-risk) 50. 50).Pearson CA, Van Schalkwyk C, Foss A, O’Reilly K, Pulliam J. Projection of early spread of COVID-19 in Africa as of 25 March 2020. CMMID. 2020. 51. 51).Pearson CA, Van Schalkwyk C, Foss AM, O’Reilly KM, Pulliam JR, CMMID COVID-19 working group. Projected early spread of COVID-19 in Africa through 1 June 2020. Euro Surveill. 2020 May 7;25(18):2000543. 52. 52).WHO Coronavirus disease (COVID-19) Situation Reports – 123. 2020. Available from: [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports) 53. 53).Cabore JW, Karamagi HC, Kipruto H, Asamani JA, Droti B, Seydi AB, Titi-Ofei R, Impouma B, Yao M, Yoti Z, Zawaira F. The potential effects of widespread community transmission of SARS-CoV-2 infection in the World Health Organization African Region: a predictive model. BMJ Glob Health. 2020 May 1;5(5):e002647. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiYm1qZ2giO3M6NToicmVzaWQiO3M6MTE6IjUvNS9lMDAyNjQ3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMTQvMjAyMC4wNy4xMi4yMDE1MjEyNC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 54. 54).Lu FS, Nguyen A, Link N, Santillana M. Estimating the prevalence of COVID-19 in the United States: three complementary approaches. April 23, 2020 ([https://www.medrxiv.org/content/10.1101/2020.04.18](https://www.medrxiv.org/content/10.1101/2020.04.18) 55. 55).Xia W, Liao J, Li C, Li Y, Qian X, Sun X, Xu H, Mahai G, Zhao X, Shi L, Liu J. Transmission of corona virus disease 2019 during the incubation period may lead to a quarantine loophole. MedRxiv. 2020 Jan 1. 56. 56).Ganyani T, Kremer C, Chen D, Torneri A, Faes C, Wallinga J, Hens N. Estimating the generation interval for COVID-19 based on symptom onset data. MedRxiv. 2020 Jan 1. 57. 57).Ferretti L, Wymant C, Kendall M, Zhao L, Nurtay A, Abeler-Dörner L, Parker M, Bonsall D, Fraser C. Quantifying SARS-CoV-2 transmission suggests epidemic control with digital contact tracing. Science. 2020 May 8;368(6491). 58. 58).Mizumoto K, Kagaya K, Zarebski A, Chowell G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. ES. 2020 Mar;25(10). DOI: 10.2807/1560-7917. Euro Surveill.2020.25.10.2000180. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2807/1560-7917&link_type=DOI) 59. 59).Arons MM, Hatfield KM, Reddy SC, Kimball A, James A, Jacobs JR, Taylor J, Spicer K, Bardossy AC, Oakley LP, Tanwar S. Presymptomatic SARS-CoV-2 infections and transmission in a skilled nursing facility. NEJM. 2020 Apr 24. 60. 60).Lavezzo E, Franchin E, Ciavarella C, Cuomo-Dannenburg G, Barzon L, Del Vecchio C, Rossi L, Manganelli R, Loregian A, Navarin N, Abate D. Suppression of COVID-19 outbreak in the municipality of Vo, Italy. medRxiv. 2020 Jan 1. 61. 61).Schraer R. Coronavirus: Asymptomatic transmission still an ‘open question’. BBC NEWS. 9 June 2020. Available from: [https://www.bbc.com/news/health-52977940](https://www.bbc.com/news/health-52977940). 62. 62).WHO Coronavirus disease (COVID-19) Situation Reports – 73. Available from: [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports) 63. 63).Cabore JW, Karamagi HC, Kipruto H, Asamani JA, Droti B, Seydi AB, Titi-Ofei R, Impouma B, Yao M, Yoti Z, Zawaira F. The potential effects of widespread community transmission of SARS-CoV-2 infection in the World Health Organization African Region: a predictive model. BMJ Glob Health. 2020 May 1;5(5):e002647. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiYm1qZ2giO3M6NToicmVzaWQiO3M6MTE6IjUvNS9lMDAyNjQ3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMTQvMjAyMC4wNy4xMi4yMDE1MjEyNC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 64. 64).Munjita SM, Samutela M, Ndashe K, Munsaka SM. Immunity, parasites, genetics and sex hormones: contributors to mild inflammatory responses in COVID-19?. Pan Afr Med J. 2020 May 15;35(36). 65. 65).Diop BZ, Ngom M, Biyong CP, Biyong JN. The relatively young and rural population may limit the spread and severity of Covid-19 in Africa: a modelling study. BMJ Glob Health. 2020 May 1;5(5):e002699. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiYm1qZ2giO3M6NToicmVzaWQiO3M6MTE6IjUvNS9lMDAyNjk5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMTQvMjAyMC4wNy4xMi4yMDE1MjEyNC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 66. 66).Dietz K. The estimation of the basic reproduction number for infectious diseases. Stat Methods Med Res. 1993 Mar;2(1):23–41. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/096228029300200103&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8261248&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F14%2F2020.07.12.20152124.atom) 67. 67).van den Driessche P. Reproduction numbers of infectious disease models. Infect Dis Model. 2017 Aug 1;2(3):288–303. 68. 68).Mukandavire Z, Liao S, Wang J, Gaff H, Smith DL, Morris JG. Estimating the reproductive numbers for the 2008–2009 cholera outbreaks in Zimbabwe. Proc Natl Acad Sci. 2011 May 24;108(21):8767–72. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTA4LzIxLzg3NjciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNy8xNC8yMDIwLjA3LjEyLjIwMTUyMTI0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 69. 69).Mukandavire Z, Smith DL, Morris Jr JG. Cholera in Haiti: reproductive numbers and vaccination coverage estimates. Scientific reports. 2013 Jan 10;3:997. 70. 70).Y Li M. An introduction to mathematical modeling of infectious diseases. Spriger; 2018. 71. 71).Fung IC. Cholera transmission dynamic models for public health practitioners. Emerg Themes Epidemiol. 2014 Dec 1;11(1):1. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1742-7622-11-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24520853&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F14%2F2020.07.12.20152124.atom) 72. 72).Sinkala M, Nkhoma P, Zulu M, Kafita D, Tembo R, Daka V. The COVID-19 Pandemic in Africa: Predictions using the SIR Model Indicate the Cases are Falling. medRxiv. 2020 Jan 1. [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif