COVID-19 vaccination hesitancy model: The impact of vaccine education on controlling the outbreak in the United States ====================================================================================================================== * Bismark Oduro * Attou Miloua * Ofosuhene O. Apenteng * Prince P. Osei ## Abstract The coronavirus outbreak continues to pose a significant challenge to human lives globally. Many efforts have been made to develop vaccines to combat this virus. However, with the arrival of the COVID-19 vaccine, there is hesitancy and a mixed reaction toward getting the vaccine. We develop a mathematical model to analyze and investigate the impacts of education on individuals hesitant to get vaccinated. The findings indicate that vaccine education can substantially minimize the daily cumulative cases and deaths of COVID-19 in the United States. The results also show that vaccine education significantly increases the number of willing susceptible individuals, and with a high vaccination rate and vaccine effectiveness, the outbreak can be controlled in the US. Key words and phrases * COVID-19 disease model * vaccination rate * vaccine hesitancy * vaccine education ## 1 Introduction The severe acute respiratory syndrome coronavirus (SARS–CoV–2) strain that caused 2019 novel coronavirus disease (COVID-19) pandemic was first observed in Wuhan, China, in December 2019 and was later declared a pandemic by the World Health Organization (WHO) on March 11, 2020 [1]. The first case of coronavirus in the United States of America was revealed on January 20, 2020, in Washington State. Since then, the number of positive cases has reached more than 30 million across the country in just over one year. COVID-19 is spread from person to person mainly through respiratory droplets released when an infected person sneezes, coughs or speaks [1, 2]. Several non-pharmaceutical measures, including face mask use, social distance, quarantining, etc., are recommended to reduce the virus’s spread. However, there is still a need for public health and clinical interventions to successfully contain the disease. Health experts agree that the best way to end the pandemic is to vaccinate most of the population [3, 4]. Currently, three COVID-19 vaccines are already known, recommended, and being used in the USA. The first is the Pfizer-BioNTech; one must complete two doses, about three weeks apart, recommended for individuals aged 16 years and older. The second vaccine available is called Moderna COVID-19 vaccine. It is also to be taken in two doses, one month or 28 days apart, introduced by a shot in the muscle of the upper arm. It is recommended for people aged 18 years and older. The third vaccine is called Johnson and Jonson’s Janssen COVID-19 vaccine, got emergency use authorized on February 27, 2021. These vaccines are intended for the prevention of coronavirus disease. They may help to mitigate the spread of the coronavirus once most individuals make an effort to be vaccinated. However, a large proportion of the American population is reluctant about the COVID-19 vaccines. Lack of information about the side effects, especially the long-term effects, time-line of the COVID vaccines production, political [5] and conspiracy theory, are among the reasons for vaccine hesitancy. The population opinion and the trust in the vaccine are of the most significant importance for appropriate coverage. A report in the American Medical Association Journal shows that skepticism toward the COVID-19 Vaccines is on the rise among Americans. A survey by the Kaiser Family Foundation shows that about 29% of health workers were hesitant to accept the vaccine [6]. In [7], a sample of 1878 adult Americans was asked different questions related to vaccination. The sample is composed of (52%) Females, (74%) Whites, (81%) non-Hispanic, (56%) married, (68%) employed full time, (77%) with a bachelor’s degree and higher. The probability of receiving the vaccine in the research was as follows: (52%) for very likely, (27%) for somewhat likely, (15%) not likely, (7%) definitely not. Vaccine hesitancy was also higher in African Americans, Hispanics, and pregnant women, and breastfeeding moms. Mathematical models are powerful tools for investigating human infectious diseases, such as COVID-19, contributing to the understanding of infections’ dynamics, and can provide valuable information for public-health policymakers [8, 9]. Numerous mathematical models have been used to provide insights into public health measures for mitigating the spread of the coronavirus pandemic [10, 11, 12, 13]. For example, Eikenberry et al. in [14] developed a mathematical model to assess the impact of mask use by the general public on the transmission dynamics of the COVID-19 pandemic. Their results showed that broad adoption of even relatively ineffective face masks might reduce community transmission of COVID-19 and decrease peak hospitalizations and deaths. Given the pervasiveness of vaccine hesitancy, we develop a compartmental model to study the impacts of education for individuals unwilling to accept the COVID-19 vaccines. We explore policy-related questions, including investigating the vaccination rate and vaccine education impact on disease dynamics in the USA. ## 2 Model formulation We consider a compartmental model for the infection’s transmission dynamics and control. With the arrival of various COVID–19 vaccines, there is a mixed reaction to get vaccinated or not. We classify the US’s total population into two subgroups: Those willing to get the vaccine and those unwilling to receive the vaccine. We further sub-divide the populations into eleven mutually exclusive compartments of willing susceptible (*S**w*), unwilling susceptible (*S**u*), willing early-exposed (*E**ew*), unwilling early-exposed (*E**eu*), willing late-exposed (*E**lw*), unwilling late-exposed (*E**lu*), willing infected (*I**w*), unwilling infected (*I**u*), willing recovered (*R**w*), unwilling recovered (*R**u*), and vaccinated population (*V*) so that the total population at time *t*, denoted by*N* (*t*) is ![Formula][1] We incorporate education for the unwilling populations at a rate *σ*. The willing populations receive the vaccine at rate *r* and progress to the vaccinated class. There is evidence that individuals exposed or infected should wait for 90 days before receiving the COVID-19 vaccine; we exclude exposed and infected individuals from the vaccination until they are recovered. Figure 1 fully illustrates the flow of populations in the various compartments; the model’s parameters and their descriptions are listed in Table 2. View this table: [Table 1:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/T1) Table 1: Description of state variables of the COVID-19 model. View this table: [Table 2:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/T2) Table 2: Description of the parameters of the COVID-19 model (1). ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/25/2021.05.21.21257612/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/F1) Figure 1: Schematic diagram of the COVID–19 model The dynamics in Figure 1 can be represented as a system of nonlinear ordinary differential equations given by ![Formula][2] where the forces of infection are given by ![Formula][3] ## 3 Results ### 3.1 Disease-free equilibrium and reproduction number The model (1) has a disease-free equilibrium (DFE) given by ![Formula][4] where *N* *∗* is the initial total population size, ![Graphic][5], and ![Graphic][6]. The next generation operator method [15, 16, 17] can be used to analyze the asymptotic stability property of the disease-free equilibrium, 𝒟. In particular, using the notation in [15, 16, 17], it follows that the associated next generation matrices, *F* and *V*, for the new infection terms and the transition terms, are given, respectively, by ![Formula][7] and, ![Formula][8] The control reproduction number is given by ![Formula][9] Where ![Formula][10] with *K*1 = *k* + *σ, K*2 = *α* + *σ, K*3 = *γ* + *δ*, and *K*4 = *γ* + *δ* + *σ*. The reproduction number is the average number of new COVID-19 cases generated by a typical infectious individual introduced into a population where a certain fraction is protected; it is a measure of contagiousness of infectious diseases. Theorem 3.1 *The disease-free equilibrium (DFE) of the model* (1) *is locally-asymptotically stable if ℛ**c* *<* 1. *When ℛ**c* *>* 1, *we will initially observe a near exponential growth of infectious cases, reach a peak, and eventually declines to an equilibrium*. ### 3.2 Parameter estimation The proposed model is fitted and validated using the USA COVID–19 daily cumulative cases and deaths from December 14, 2020, to January 20, 2021. The choice of this data is motivated by the commencement of national vaccination on December 14, 2020, in the USA. There are nine parameters underlying the model; however, three of the parameters, i.e., *γ, δ*, and *θ*, were fixed, and the rest were estimated. It is worth noting that the model parameters may not be uniquely identifiable based only on cumulative cases and deaths of USA COVID–19 data available. We addressed this identifiability problem by using an inverse modelling, sensitivity and Monte Carlo analysis method built in FME [18] package in R (see also [19]). This method analyzes mathematical models with data, performs local and global sensitivity, and Monte Carlo analysis. It addresses parameter identifiability issues and fits a model to data using existing optimization methods such as the constrained quasi-Newton method. The estimated parameter values obtained from the best model prediction are given in Table 3. The fitted models and the observed data for the coronavirus cumulative cases and deaths in the USA are displayed in Figure 2. We observe that the model fit the observed data significantly. View this table: [Table 3:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/T3) Table 3: Estimated parameter values for the model (1) using data for the US. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/25/2021.05.21.21257612/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/F2) Figure 2: The fitted model and the observed COVID–19 cumulative cases and deaths in USA from December 14, 2020 to January 20, 2021 ### 3.3 Sensitivity analysis This section uses the Latin Hypercube Sampling and Partial Rank Correlation Coefficients (PRCC) to perform sensitivity analysis [21, 22] on the model parameters. The analysis is needed to identify model parameters that significantly impact model outcomes using the reproduction number (ℛ*c*) as the response function [12, 21, 22]; that is, to determine the model robustness to parameter values. Parameters with large PRCC greater than +0.50 are strongly positively correlated with the response function. In contrast, those less than *−*0.50 are said to be largely negatively associated with the response function [21, 22]. Figure 3 displays the PRCC analysis plot of the model parameters considered. The results show that the effective contact rate (*β*) and the modification parameter (*ρ*) positively affect ℛ*c*, meaning that an increase in these parameters will increase ℛ*c*. On the other hand, the education rates for unwilling individuals (*σ*), the vaccination rate *r*, and vaccine efficacy (*θ*) have a negative effect on the *ℛ**c*, and increasing these parameters would lower the *ℛ**c*. The results further indicate that the parameters; the effective contact rate, education rate for unwilling individuals, the vaccination rate, and the vaccine efficacy mainly influence the response function (*ℛ**c*). ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/25/2021.05.21.21257612/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/F3) Figure 3: Partial rank correlation coefficients (PRCCs) showing the impact of model parameters on the reproduction number (ℛ*c*) of the model. Parameter values used are as given in Table 3. ### 3.4 Contour plots results We generated contour plots to analyze the reproduction number (ℛ*c*) of the model as a function of desired parameters as displayed Figure 4. Parameter values used for the simulation are as given in Table 3. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/25/2021.05.21.21257612/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/F4) Figure 4: Contour plot of the reproduction number (*ℛ**c*) of the model (1), as a function of: (a) Education rate (*σ*) and the ratio ![Graphic][11], (b) Vaccination rate (*r*) and the ratio ![Graphic][12], and (c) Vaccination efficacy (*θ*) and vaccination rate (*r*). Based on the contour plot results in Figure 4, the following observation inferences are made: 1. To effectively curb the outbreak, that is, reducing (ℛ *c*) to a value less than unity, the Figure 4 (a) suggests that a high education rate (*σ*) is needed. The lower the fraction of individuals hesitant to accept the vaccine, the lower the education rate required and vice versa. For example, to decrease the reproduction number below one, at least an education rate of 0.3 is needed if about 40% of the susceptible population is unwilling to vaccinate. 2. The ℛ*c* decreases as more individuals are being educated and are willing to receive the vaccine, and the vaccination rate is high, Figure 4 (b). 3. Figure 4 (c) illustrates the impact of a relatively high vaccination rate and efficacy. One can observe that a high combined effects cause the *ℛ* *c* to drop to below one. ### 3.5 Time evolution analysis, and predictions In this section, we analyze the effects of vaccine education *σ*, the vaccination rate *r*, and the efficacy of the vaccine *θ*, on the cumulative US COVID-19 cases and deaths. In addition, we provide some future predictions of the cumulative cases and deaths with their associated uncertainty. First, we begin with the effects of the education for unwilling individuals and its impact on the cumulative cases and deaths. Table 4 provides numerical description while Figure 5 graphically explains the decrease and increase in cumulative cases and deaths when the education rate is increased or decreased, respectively. For example, from Table 4, a 25% increase in education rate reduces the cumulative cases from 25441752 to 25108128 and from 425513 to 422423 for the cumulative deaths. View this table: [Table 4:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/T4) Table 4: A summary of various increase in education rate. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/25/2021.05.21.21257612/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/F5) Figure 5: Effects of increasing and decreasing the baseline education rate parameter *σ* We observe that an increase in vaccine education decreases the COVID-19 cumulative cases and deaths. On the other hand, a lower vaccine education increases the baseline cumulative cases and deaths. This suggest that, an increase of education parameter from the baseline, lowers the spread of COVID-19. It is worth noting that while the curves for the baseline education parameter and its increase or decrease overlap up until the end of the year 2020, they are separated beginning of January 2021. Second, we analyze the effects of increasing or decreasing the vaccination rate *r*, on cumulative cases and deaths numerically and graphically. The results are displayed in Table 5 and Figure 6. View this table: [Table 5:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/T5) Table 5: A summary of various increase in vaccination rate. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/25/2021.05.21.21257612/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/F6) Figure 6: Effects of increasing and decreasing the baseline of vaccination rate, *r* Note that 50% increase or decrease in the vaccination rate decreases or increases the cumulative cases and deaths significantly compared to the baseline scenario. For instance, a 50% increase in the baseline rate of vaccination reduces the cumulative cases and deaths from 25441752 and 425513 to 24192245 and 414785 respectively. On the other hand, the same reduction rate from the baseline increases the cumulative cases and deaths to 27141116 and 439270, respectively. Note also that the 50% increase of the vaccination rate is significantly below the baseline after January 15, 2021, while a 50% reduction of the rate of vaccination is significantly above the baseline curve. This suggests that as the vaccination rate increases, the cumulative cases and deaths of COVID–19 in the USA decreases. Finally, we consider the effects of the efficacy of the vaccine and report the numerical and graphical results of varying the parameter *θ*. Note that the baseline vaccine efficacy was high, and thus the maximum increment adopted in this communication stands at 5%, which represents 99.4% efficacy. The quantitative results and graphical description of various decrease and an increase in the efficacy of the vaccine is displayed in Table 6 and Figure 7. View this table: [Table 6:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/T6) Table 6: A summary of various decrease and an increase in the efficacy of COVID–19 vaccine. ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/25/2021.05.21.21257612/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/F7) Figure 7: Effects of increasing and decreasing the baseline efficacy of the COVID–19 vaccine *θ* We observe that a lower vaccine efficacy has the highest cumulative cases and deaths of COVID-19. Note that a 45% reduction of vaccine efficacy from the baseline, i.e., a 52.1% vaccine efficacy increases the cumulative cases and deaths from the baseline 25441752 to 26945280 and 425513 to 437726, respectively. Also, note that the curve of the 5% increase of the baseline efficacy of the vaccine is similar to the baseline line. We validate our model by predicting the cumulative cases and deaths every two weeks after January 20, 2021 for the USA COVID–19 outbreak; all the data are given in the format of 95% Confidence Interval (CI) of the predicted values. The uncertainty of the predicted values are characterized using the the 95% confidence intervals (CI). We compared the estimated predicted values to the actual values reported in [23] within the 95% CI, as shown in Table 7. Our model is robust and can be replicated for other countries. View this table: [Table 7:](http://medrxiv.org/content/early/2021/05/25/2021.05.21.21257612/T7) Table 7: Prediction and reported values. ## 4 Discussion and conclusions The coronavirus outbreak poses a severe threat to human lives globally. Since the beginning of the pandemic, a campaign to use non-pharmaceutical measures to prevent the virus’s spread has been underway, including a face mask, social distancing, and frequent hand washing. The emergence of vaccines sounds promising; however, there is a challenge of vaccine hesitancy. It is therefore essential to understand how educating the general public will impact the fight against the outbreak. In this paper, we use a compartmental model with vaccine education to study the dynamics of the COVID-19 infection. We classify the US’s total population into two subgroups: Those willing to accept the vaccine and those unwilling to receive the vaccine. The vaccine education is incorporated for the general public hesitant to take the vaccine. We assessed the impact of the education campaign on the control of the outbreak. First, we computed an expression for the reproduction number (ℛ*c*), a threshold that measures the contagiousness of infectious diseases. We performed a sensitivity analysis of the reproduction number. The result shows that vaccine education negatively influences the reproduction number; that is, an increase in the vaccine education implies a decrease in the ℛ*c*. Epidemiologically, when ℛ*c* *≤*1, then the transmission will fade or die out. In contrast, the infected number of people is expected to increase if ℛ*c* *>* 1. The sensitivity analysis result also shows that vaccine efficacy and vaccination rate negatively impact theℛ*c*, and raising them will reduce the ℛ*c*. Using contour plots, we further analyzed the reproduction number as a function of two independent variables, the vaccine education and the proportion of unwilling susceptible individuals ![Graphic][13]. The prospect of curtailing the outbreak is achievable with a high education rate. For instance, the result in Figure 4 (a) shows that if a high proportion of the susceptible populations is unwilling to accept the vaccine, high public education can help diminish the reproduction number and mitigate the virus. Also, vaccine education raises the willing susceptible individuals, and with a high vaccination rate and vaccine efficacy, it will be possible to control the spread of the outbreak. Next, we analyze the effects of vaccine education, vaccination rate, and the efficacy of the vaccine on the daily cumulative cases and deaths from December 14, 2020 to January 20, 2021. The results suggest that higher rates of vaccine education, vaccination, and vaccine efficacy contribute to the mitigation of the spread of the COVID-19 disease outbreak. For example, the result in Figure 5 shows that a lower vaccine education contributes to the higher projected cumulative cases and deaths. On the other hand, higher education rates help to control the outbreak. The effect of these higher rates becomes evident over time. These results suggest that a combined effect of higher vaccine education, vaccination, and vaccine efficacy rates would contribute to slow the pace of the spread of the COVID-19 outbreak and thus help the government to control the disease. The mitigation effect of higher rates of education, vaccination, and vaccine efficacy is clearly evident over time as we observe from the curves in Figures 5, 6 and 7; thus, these mitigation strategies need to be enforced consistently over time. Finally, to validate and test the robustness of our proposed mathematical model, we predicted the cumulative cases and deaths every two weeks after January 20, 2021. The reported cumulative cases and fatalities in Table 7 almost all fall within the prediction interval of 95% confidence level. The results suggest that our model is robust and is applicable for COVID–19 vaccine hesitancy models in other countries. ## Data Availability The data can be accessed via https://www.worldometers.info/coronavirus/country/us/ [https://www.worldometers.info/coronavirus/country/us/](https://www.worldometers.info/coronavirus/country/us/) ## Footnotes * miloua{at}calu.edu * ofap{at}food.dtu.dk * prinpep{at}gmail.com * Received May 21, 2021. * Revision received May 21, 2021. * Accepted May 25, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1].World Health Organization (WHO). 2. [2].Centers for Disease Control and Prevention (CDC). [https://www.cdc.gov/coronavirus/2019-ncov/index.html](https://www.cdc.gov/coronavirus/2019-ncov/index.html). 3. [3].Fontanet, A., Autran, B., Lina, B., Kieny, M. P., Karim, S. S. A., & Sridhar, D. (2021). SARS-CoV-2 variants and ending the COVID–19 pandemic. The Lancet. 4. [4].Callaghan, T., Moghtaderi, A., Lueck, J. A., Hotez, P., Strych, U., Dor, A., … & Motta, M. (2020). Correlates and disparities of intention to vaccinate against COVID–19. Social science & medicine. 5. [5].The COCONEL Group. (2020) A future vaccination campaign against COVID-19 at risk of vaccine hesitancy and politicisation. The lancet. Volume 20(7), P769–770, JULY 01, 2020. DOI:[https://doi.org/10.1016/S1473-3099(20)30426-6](https://doi.org/10.1016/S1473-3099(20)30426-6). 6. [6].Kaiser Family Foundation (KFF) COVID–19 Vaccine Monitor. December 2020, [https://www.kff.org/coronavirus-covid-19/report/kff-covid-19-vaccine-monitor-december-2020/](https://www.kff.org/coronavirus-covid-19/report/kff-covid-19-vaccine-monitor-december-2020/) 7. [7].Khubchandani J., Sharma S., Price J. H., Wiblishauser M. J., Sharma M., Webb F. J. (2021) COVID-19 Vaccination Hesitancy in the United States: A Rapid National Assessment. J Community Health. 46(2):270–277. doi: 10.1007/s10900-020-00958-x. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10900-020-00958-x&link_type=DOI) 8. [8].Bramson, H., et al. (2015). State laws, syringe exchange, and HIV among persons who inject drugs in the United States: History and effectiveness. Journal of Public Health Policy, 36(2), 212e230. 9. [9].Luboga, S., Galukande, M., Mabweijano, J., Ozgediz, D., & Jayaraman, S. (2010). Key aspects of health policy development to improve surgical services in Uganda. 10. [10].Iboi E., Richardson A., Ruffin R., Ingram D., Clark J., Hawkins J., McKinney M., Horne N., Ponder R., Denton Z., Agusto F. B., Oduro B., Akinyemi L. (2021). Impact of Public Health Education Program on the Novel Coronavirus Outbreak in the United States. Front. Public Health. doi: 10.3389/fpubh.2021.630974 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fpubh.2021.630974&link_type=DOI) 11. [11].Owusu-Mensah, I., Akinyemi, L., Oduro, B. et al. A fractional-order approach to modeling and simulations of the novel COVID–19. Adv Difference Equations, 2020, 683. 12. [12]. Ofosuhene O. Apenteng, Bismark Oduro, Isaac Owusu-Mensah. A Compart-mental Model to Investigate the Dynamics of the COVID–19 Pandemic: A case Study in Five Countries. International Journal of Biomathematics, 2020. [https://doi.org/10.1142/S1793524521500273](https://doi.org/10.1142/S1793524521500273) 13. [13].Oduro, B. and Magagula, V. M. COVID-19 intervention models: An initial aggressive treatment strategy for controlling the infection. Infectious Disease Modelling. Volume 6, 2021, Pages 351–361. 14. [14].Eikenberry et al. (2020). To mask or not to mask: Modeling the potential for face mask use by the general public to curtail the COVID–19 pandemic. Infectious disease modeling, 293–308. 15. [15].Diekmann, O., Heesterbeek, J. A. P, Britton, T. (2012). Mathematical Tools for Understanding Infectious Disease Dynamics. Kindle Edition. Princeton University Press. 16. [16].Driessche, P. V., Watmough, J. (2002). Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Mathematical Biosciences 180 29–48. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0025-5564(02)00108-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12387915&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F25%2F2021.05.21.21257612.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000179220600004&link_type=ISI) 17. [17].Heesterbeek, J. A. P. (2002). A brief history of R0 and a recipe for its calculation. Acta Biotheo 50 189—204. 18. [18].Soetaert, K., Petzoldt, T. (2010). Inverse modelling, sensitivity and Monte Carlo anal-ysis in R using package FME. Journal of Statistical Software 33 (3), 1–28. 19. [19].Apenteng, O. O., Osei, P. P., Oduro, B., Kwabla, M. P., Ismail, N. A. (2020). The impact of implementing HIV prevention policies therapy and control strategy among HIV and AIDS incidence cases in Malaysia. Infectious Disease Modelling 5 755–765. 20. [20]. C. N. Ngonghala, et al. Mathematical assessment of the impact of non-pharmaceutical interventions on curtailing the 2019 novel Coronavirus. Mathematical Biosciences. 2020. Volume 325, 108364. 21. [21].Cariboni J., Gatelli D., Liska R., and Saltelli A. The role of sensitivity analysis in ecological modelling, Ecological modelling 203, 167–182 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ecolmodel.2005.10.045&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246020300017&link_type=ISI) 22. [22].Blower S. M. and Dowlatabadi H. Sensitivity and uncertainty analysis of complex models of disease transmission: an hiv model, as an example. International Statistical Review/Revue Internationale de Statistique, 229–243 (1994). 23. [23].Worldometer (2020). [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/) Accessed: 2 June 2020. 24. [24].Wang, C., Horby, P.W., Hayden, F.G. and Gao, G.F. A novel coronavirus outbreak of global health concern. The Lancet, 395(10223), 2020. pp.470–473. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)&link_type=DOI) 25. [25]. B. Tang, X. Wang, Q. Li, N. L. Bragazzi, S. Tang, Y. Xiao, and J. Wu, “Estima-tion of the transmission risk of the 2019-nCoV and its implication for public health interventions,” Journal of Clinical Medicine 9, 462. 2020. 26. [26].John Hopkins University, “Daily state-by-state testing trends,” John Hop-kins University Coronavirus resource center (Assessed on January 5, 2021). [https://coronavirus.jhu.edu/region/united-states](https://coronavirus.jhu.edu/region/united-states). 27. [27].Food and Drug Administration (FDA) [https://www.fda.gov/vaccines-blood-biologics/vaccines/emergency-use-authorization-vaccines-explained](https://www.fda.gov/vaccines-blood-biologics/vaccines/emergency-use-authorization-vaccines-explained). [1]: /embed/graphic-1.gif [2]: /embed/graphic-5.gif [3]: /embed/graphic-6.gif [4]: /embed/graphic-7.gif [5]: /embed/inline-graphic-1.gif [6]: /embed/inline-graphic-2.gif [7]: /embed/graphic-8.gif [8]: /embed/graphic-9.gif [9]: /embed/graphic-10.gif [10]: /embed/graphic-11.gif [11]: F4/embed/inline-graphic-3.gif [12]: F4/embed/inline-graphic-4.gif [13]: /embed/inline-graphic-5.gif