Analysis of key factors of a SARS-CoV-2 vaccination program: A mathematical modeling approach
=============================================================================================

* David Martínez-Rodríguez
* Gilberto Gonzalez-Parra
* Rafael-J. Villanueva-Micó

## Abstract

The administration of vaccines against the coronavirus disease 2019 (COVID-19) just started in early December of 2020. Currently, there are only a few approved vaccines, each with different efficacies and mechanisms of action. Moreover, vaccination programs in different regions may vary due to differences in implementation, for instance, simply the availability of the vaccine. In this article, we study the impact of the pace of vaccination and the intrinsic efficacy of the vaccine on prevalence, hospitalizations, and deaths related to the SARS-CoV-2 virus. Then we study different potential future scenarios regarding the burden of the COVID-19 pandemic in the near future. We construct a compartmental mathematical model and use computational methodologies to study these different scenarios. Thus, we are able identify some key factors to reach the aims of the vaccination programs. We use some metrics related to the outcomes of the COVID-19 pandemic in order to assess the impact of the efficacy of the vaccine and the pace of the vaccine inoculation. We found that both factors have a high impact on the outcomes. However, the rate of vaccine administration has a higher impact in reducing the burden of the COVID-19 pandemic. This result shows that health institutions need to focus in increasing the vaccine inoculation pace and create awareness in the population about the importance of COVID-19 vaccines.

## 1 Introduction

The world is facing the COVID-19 pandemic and just in the middle of December of 2020 the first vaccines are being given to a few people [42, 128, 78, 119, 103, 129, 130, 115]. This pandemic has caused more than 80 million confirmed cases and more than 1.7 million deaths [30, 136].

The SARS-CoV-2 virus causes an illness called COVID-19 that can result in severe pneumonia and death [81, 102]. The complex process of the SARS-CoV-2 spread involves several factors that are currently not very well understood [84, 112, 139, 102, 31, 32]. There are several aspects that impact the spread of the virus in the human population, such as, social behavior, age, weather variables, mutation of the virus, and immunocompetence [152, 111, 147]. Other factors may affect the spread of the SARS-CoV-2 virus but still be unknown. Regarding mutations, SARS-CoV-2 could acquire mutations with fitness advantages and immunological resistance [71]. Therefore, studying evolutionary transitions is important to ensure effectiveness of the vaccines and immunotherapeutic interventions [151, 50, 104, 157, 71]. It has been stateded that not only is the efficacy of the vaccine important, but whether a vaccine reduces infection and transmission as well as disease progression [69].

The genomic analysis suggested that the base sequence of SARS-CoV-2 is almost 80% similar to that of SARS-CoV. Further, both these viruses bind to same host cell receptor ACE-2 [6]. However, the development of vaccines for this novel SARS-CoV-2 virus took nearly a year. There are now more than 80 vaccines in development [3, 16, 69, 71, 149, 105, 135]. Vaccination programs have recently begun (December) in a few countries, and there are many uncertainties regarding the optimal implementation of these vaccination programs and the probable outcomes [1, 69, 137]. Therefore, studying the COVID-19 vaccination programs is of paramount importance. Effective vaccination helps tackle the transmission of the SARS-CoV-2 virus in the population [9, 29, 69, 72, 137, 145].

Mathematical models, statistical analyses and computational techniques are very useful tools to study different processes, including testing hypotheses and understanding how factors affect the processes. For infectious disease processes, mathematical models can be used to perform *in silico* simulations of different potential scenarios, vaccination programs, and test different strategies to slow down epidemics [92, 12, 73, 40, 2, 20, 114, 45, 5, 46]. The outcomes of the complex infectious disease processes under different scenarios are generally impossible to predict without mathematical models and computational techniques. In some cases, results of the simulations can be counter-intuitive and very interesting from a predictive point of view.

There are many articles related to the use of mathematical modeling in combination with computational and statistical techniques to study the spread of the SARS-CoV-2 virus [31, 73, 126, 40, 130, 74, 155]. Some mathematical models used the SIR (Susceptible-Infected-Recovered) mathematical model [107, 8, 110]. Previous studies have used susceptible-exposed-infected-recovered (SEIR) type models [130, 61, 153]. Other mathematical models use a curve fitting of some particular growth model to the data, and also artificial intelligence techniques have been considered for fitting models to real data related to COVID-19 [41, 118, 80].

The main advantage of mathematical models is that many different simulations can be done and this allows us to study the main driving factors of pandemics under a variety of complex scenarios [92, 12, 73, 40, 20, 114, 45, 5, 46]. However, many forecasts related to the COVID-19 pandemic disagree with each other due to many related uncertainties in key characteristics of the SARS-CoV-2 virus [117, 59, 38, 64, 125, 133, 73, 126, 40, 130, 74, 155]. Moreover, currently we are facing new strains due to mutations of the virus, which has raised questions about the efficacy of the vaccines against the mutations of the SARS-CoV-2 virus. Recently, it has been found that the SARS-CoV-2 is mutating and its transmission is more efficient [151, 98, 33, 77]. There is a growing literature about mutations of the SARS-CoV-2, but it is not clear what further mutations could occur in the near future [151, 50, 104, 157, 71].

A new vaccine campaign against the SARS-CoV-2 virus began in December of 2020 in the United States and other countries. Currently, there are only a few approved vaccines with different efficacies and mechanisms of action. Our principal aim in this article is to study the impact of the pace of vaccination and the efficacy of the vaccine on the outcome or dynamics of the incidence, prevalence, and deaths related to the SARS-CoV-2 virus [58, 70, 79, 82, 101]. This will help us to explain different potential patterns in different countries related to their vaccination programs [83, 98, 112, 51]. Despite the huge health crisis caused by the spread of the SARS-CoV-2 virus around the world, there are few studies related to the prediction of feasible scenarios in the year 2021 [49, 87, 91, 99]. In [91] the authors used a SEIR mathematical model based on differential equations to study the situation in South-Africa with respect to the number of reported cases of COVID-19. They found that a vaccine with 70% efficacy had the capacity to contain the COVID-19 outbreak but only at a very high vaccination coverage of 94.44%. In [49], the authors proposed a distribution of vaccines in time and space, which sequentially prioritizes regions with the most new cases of infection in a given time period. They used a SEIR type model that includes an extra class for mild infected individuals and with spatio-temporal effects. They found that, for a locally well-mixed population, the proposed strategy strongly reduces the number of deaths. In [87], the authors implemented a SEIR type model that includes pre-symptomatic, asymptomatic and the entire equivalent vaccinated classes. They studied the impact of different vaccination coverages, efficacy and reduction of symptoms on different metrics such as deaths, and ICU hospitalizations. The outcomes of these studies are important because they help to understand better the impact of vaccination programs, and generate optimal actions to diminish the spread of the SARS-CoV-2 virus [73, 126, 40, 55, 87, 152, 111, 147]. One factor that brings uncertainty to the outcomes and that must be taken into account is the fact that vaccines are under threat in different places and there are antivaccine movements that have gained traction with some people [1, 27].

In this article, we construct a compartmental mathematical model and use computational methodologies to study different scenarios. In particular, we include the asymptomatic carriers of the virus, who are nevertheless able to spread the virus. It has been mentioned that asymptomatic people are in some way a key contributor in the spread of the SARS-CoV-2 virus, and are real threat for the control policies [7, 62, 90, 100, 124, 43, 68, 52, 132, 31, 62]. For instance, it has been found that quantitative SARS-CoV-2 viral loads were similarly high for infected individuals with symptoms, pre-symptomatic, or asymptomatic. Moreover, it has been found 6 to 24 times more estimated infections per site with seroprevalence than with coronavirus disease 2019 (COVID-19) case report data [56]. Thus, we construct a mathematical model taking into account asymptomatic people, which have been missed in other studies. In addition, we will study scenarios with different SARS-CoV-2 virus transmission rates, which results in different effective reproduction numbers *R**t* of COVID-19 [75, 85, 14, 25].

## 2 Materials and Methods

### 2.1 Mathematical Model

We constructed a compartmental model based on differential equations that includes individuals in the susceptible, latent, infected, asymptomatic, and hospitalized stages. The mathematical model considers transitions of individuals through the aforementioned stages depending on the COVID-19 progression. In addition to the previous stages, the model incorporates vaccinated individuals that might be in analogous stages such as susceptible or asymptomatic vaccinated. Thus, in some way we can classify individuals in two disjoint groups: unvaccinated and vaccinated. We assume that unvaccinated individuals in the susceptible, latent, and asymptomatic compartments can receive the vaccine against the SARS-CoV-2 virus. On the other hand, we assume that symptomatic, recovered and hospitalized unvaccinated individuals do not receive the vaccine. The individuals can transit from the unvaccinated susceptible class to vaccinated susceptible if they get the vaccine. In an analogous way, the latent and asymptomatic unvaccinated individuals can move to the respective vaccinated compartment. It is important to mention that the model incorporates the type of vaccine that diminishes the progression to the COVID-19 disease [87, 99]. Individuals in the latent stage are not yet infectious. The individuals remain in the latent stage for a certain time which is chosen from an exponential distribution with mean time *α*. The individuals then transit into the infective symptomatic or asymptomatic stages, where they are able to spread the SARS-CoV-2 virus to other individuals. They stay in the infectious stage for a time chosen from an exponential distribution with mean time *γ*. After that, individuals in the asymptomatic stage move to the recovered stage. However, individuals in the infective symptomatic stage can move to the recovered or to the hospitalized stages, depending on the level of disease progression. Even though we assume exponential distributions, the Erlang distributions are more realistic but at the expense of more complex models and more parameters [63, 48, 116, 47, 38, 140]. Thus, many studies assume exponential distributions to avoid greater complexity in the models and in the analysis. However, in some cases exponential distributions are not far from reality. We have found that the length of stay in the hospital is not far from an exponential distribution [37]. Finally, hospitalized individuals can die due the COVID-19 disease [37, 39, 40, 156]. This last metric (or outcome) is of paramount importance [40, 138, 142, 150].

We use a mathematical model that is similar to a *SEIR*-type epidemiological model to explain the dynamics of COVID-19 spread on the human population under a vaccination program. This model has parameters that can be varied in order to study different possible scenarios. For instance, the pace of vaccination and efficacy of the vaccine can be modified. This is important since it is known that the efficacy of vaccines varies and they have different underlying mechanisms of action [69, 72, 58, 79, 9]. Moreover, different countries and regions would apply the vaccines at different rates due to a variety of factors such as availability and resources [1, 35, 69, 79, 27].

The constructed mathematical model based on differential equations is given by ![Formula][1]</img>  where *S*(*t*) denotes the number of susceptible individuals. When a susceptible and an infectious individual come into infectious contact, the susceptible individual contracts the disease and transitions to the latent compartment *E*(*t*). Individuals in compartment *E*(*t*) are infected (carry the virus) but cannot spread the virus. Compartment *I*(*t*) represents individuals who have been infected and show symptoms. These individuals are capable of infecting susceptible individuals after being in the *E*(*t*) subpopulation. The subpopulation *A*(*t*) represents the number of individuals who have been infected but are asymptomatic. These individuals are capable of infecting susceptible individuals after being in the *E*(*t*) subpopulation. The variable *H*(*t*) denotes the number of hospitalized individuals at time *t*. The compartment *D*(*t*) represents the number of deaths due to the SARS-CoV-2 virus from the beginning of the simulation period. Similarly, *SV* (*t*), *EV* (*t*), *IV* (*t*), *AV* (*t*) and *HV* (*t*) denote the analogous vaccinated population at time *t*. Individuals in the *S*(*t*), *E*(*t*) and *A*(*t*) classes are vaccinated with rates *ν p**S*, *ν p**E* and *ν p**A* respectively. These proportions are related to the their respective size populations. The model assumes that people in states *E*(*t*), *EV* (*t*), *HV* (*t*), *H*(*t*), and *R*(*t*) do not transmit the infection.

In this model we consider that COVID-19 confers immunity after recovery (currently assumed but not con-firmed), and assume that when an individual is in the latent and latent vaccinated stages the virus cannot be transmitted. In addition, we consider that once the individuals receive the vaccine inoculation they transit to the vaccination compartments. The model assumes one vaccine inoculation to cause the individual to transit to one of the vaccinated compartments. The model also considers that hospitalized individuals are not able to transmit the SARS-CoV-2 virus. This assumption is arguable, but we assume the conditions in the hospitals are safe regarding the transmission of the SARS-CoV-2 virus. We also assume that individuals in the susceptible, latent, and asymptomatic compartments are those that can be vaccinated. This assumption might sound arguable, but recently it has been mentioned that a nurse tested positive for COVID-19 more than a week after receiving Pfizer Inc’s vaccine [113]. One hypothesis to explain this observation is that the nurse was in the latent stage before being inoculated.

### 2.2 Parameter values

In this work we are interested in the impact of the vaccination rate and the efficacy of the vaccines on the infected, hospitalized, and death cases. We assume that the rates of virus transmission in asymptomatic and symptomatic individuals are constant from the beginning of the period of study i.e. when the vaccination program starts. This implicitly assumes that people would not change behavior (on average) until the vaccination program is well advanced. This is a credible assumption in the USA, and previous physical and social behavior changes can be included in the transmissibility. Many health policies and guidelines would have been implemented before the vaccination program started. In some cases it is more realistic to include time-varying transmissibility, which has been used to study other infectious diseases and in particular one closely related to the SARS-CoV-2 virus [131, 76, 73, 67]. This latter approach is more troublesome to implement since it is necessary to estimate a time-varying parameter, and identifiability issues thus arise. Even with accurate data from the past it is difficult to estimate a time-varying transmission. Moreover, in this study we can not predict how the behavior of individuals might change in the future. Thus, we take an approximation and a conservative assumption that the transmissibility would not change during the beginning of the vaccination program.

We assume that the parameters related to the COVID-19 disease progression are the same for vaccinated and unvaccinated. In addition, we consider that vaccinated individuals are not able to get the disease unless the vaccine was not effective. This aspect is not clear in the relevant scientific literature due to the different types of COVID-19 vaccines [69, 72, 58, 79, 9]. In addition, some studies have indicated that the antibody titers may decline over time in patients recovered from COVID-19, particularly in those who were asymptomatic [143]. However, we do not consider that recovered individuals can return to the susceptible stage. One reason for this is that further studies are needed to check how long the immunity lasts, and furthermore, the time horizon of this study is less than 16 months. We also consider that for this period the immunity provided by the vaccines does not diminish.

Reasons for doubting this undiminished immunity come from studies on the SARS virus. For example, in a study of 56 patients recovered from SARS it was found that the neutralizing and IgG antibodies quickly declined after 16 months and continued to decline further to a very low level after 3 years [143]. Moreover, trials of SARS vaccines also suggest that the neutralizing antibody responses may decline over time [143]. On the other hand, it is unclear whether vaccine induced antibody levels could persist and, if not, whether the long-lasting memory T cells could affect susceptibility and pathogenesis of SARS-CoV-2 infection [143, 19].

It has been mentioned that US federal officials hoped for twenty million people to get their first of two required shots by the end of 2020. However, they recently changed that goal and just over 1 million doses of vaccines had been administered (Dec. 26th 2020) [134]. Therefore, we assume as a lower bound rate *ν* for the inoculation of the vaccine a value of one million per week. This rate can be increased since it is expected that the process of the vaccine administration will be improved. However, this value of the parameter is subject to variation due to the reluctance of some people to vaccinate for because of doubtfulness about the preliminary tests of efficacy to pressure from anti-vaccination movements [27].

For the death rate of hospitalized individuals we use a variety of data from the scientific literature [40, 66, 148, 91, 99]. We used the weighted average of the probability of dying for severe and critical cases (ICU), and in addition we took into account the average length of stay in the hospital [99]. We varied in a reasonable way the death rate in order to take into account the possible uncertainty in the data.

For the asymptomatic cases and proportions we also relied on data from the scientific literature [21, 88, 40, 68, 90, 97, 96, 154]. However, the discrepancies in the relevant data are great. We chose as a conservative starting point that the percentage of infections that are asymptomatic is 50% [30]. However, for the numerical simulations we additionally considered a percentage of 40% [30, 96].

For the *β**s* parameters we assume values in the range of [0.1 − 0.5], which are values found in some studies. Currently, in USA there are several non-pharmaceutical interventions. We also assume for the numerical simulations that *β**A* *≤β**I*. This assumption is based on the uncertainty in these values as well as results from the literature that the infectiousness of asymptomatic carriers is similar or smaller to the symptomatic [21, 88, 40, 68, 121, 154]. One interesting article found that asymptomatic carriers have a higher viral load, and, taking into account that asymptomatic carriers might have more physical contacts, it is possible to assume that *β**A* *≥ β**I* [54].

### 2.3 Initial conditions for the scenarios

For the initial conditions we assume the particular situation of the USA since is one of the first countries that started a vaccination program [134, 106]. We rely on data from the scientific literature and demographics of the USA. As expected, there are some uncertainties related to data of the COVID-19 pandemic and which is usual in many epidemics. For instance, the infected reported cases have uncertainties due to many factors such sensitivity and specificity of COVID-19 tests [10, 127]. Moreover, asymptomatic cases represent a great uncertainty [31, 43, 52, 62, 90, 100, 124, 132, 7]. Taking into account these uncertainties, we set the initial conditions presented in Table 2. The total initial population *N* (0) is taken from the current USA population [15]. The birth and death rates are taken from the official website of the CDC in USA [144, 86]. All the initial vaccinated subpopulations, are set to zero since the simulations are performed at the beginning of the vaccination program. Two key initial subpopulations are those corresponding to the infected and asymptomatic, since they affect the initial dynamics of the COVID-19 pandemic under the vaccination program. We took the seven day average of the infected reported cases and then multiplied by seven days (assumed infectiousness period) and by 0.8 to obtain the initial number of symptomatic cases (assumption of the proportion of symptomatic cases in the reported cases) [17, 4, 89, 123, 23, 154, 65, 122, 44]. The percentage of asymptomatic cases in the official statistics varies for each country. In some countries it may be close to zero, since no random tests are performed. However, the detection of asymptomatic infections is possible in the case of the USA in situations in which testing is mandatory (as in some universities) or random. We approximated this value by relying on data from different studies [17, 4, 89, 123, 23, 154, 65, 122, 44]. However, in our simulations we varied the parameter through reasonable values. This variability only affects the initial conditions of some of the populations.

In the reported cases we have a subpopulation of asymptomatic cases since the testing programs take into account the entire population. We take the initial subpopulation of asymptomatic carriers as equal to the symptomatic one. This implicitly assumes that the percentage of asymptomatic cases is 50%. We use this value based on information from the CDC official website, even though it is mentioned that there is uncertainty in this percentage. We found that there is a large uncertainty in scientific literature regarding this percentage [30, 31, 43, 52, 62, 90, 100, 124, 132, 7]. Therefore, we vary it in order to deal with its uncertainty. For the initial latent subpopulation we take into account that the latent period is around 5.2 days and the latent stage includes individuals who will become either asymptomatic or symptomatic [40]. For the initial hospitalized subpopulation we take into account that hospitalized individuals spend an average of 10.4 days in the hospital and that around 4% of the symptomatic infected transit to the hospitalization stage [40, 78, 60]. For the recovered COVID-19 cases we take into account the current total of reported infected cases and the fact that a subset of the asymptomatic cases are not reported. In addition, we notice that we need to subtract the current number of infected and asymptomatic cases. This approximation gives us a plausible number of recovered cases that exceeds the reported recovered cases (*≈* 11 millions) [30, 136]. Finally, we use for the initial susceptible subpopulation the fact that initially there are not vaccinated individuals and therefore *S*(0) = *—N* (0) *−E*(0) *−I*(0) *−A*(0) *−R*(0) *− H*(0). In Table 2, we present the initial conditions for the subpopulations.

## 3 Results

In this section, we perform numerical simulations of the mathematical model (1) to analyze the impact of the vaccination rate and the efficacy of the vaccine on the dynamics of the COVID-19 pandemic. We use the parameter values of Table 1 and the initial conditions given in Table 2. We vary the values of the vaccination rate, the efficacy of the vaccine and the transmission rates in order to include a variety of scenarios that take into account the uncertainty in the aforementioned factors. We introduce some important metrics related to the outcomes of the COVID-19 pandemic in order to assess the impact of the inoculation rate and the efficacy of the vaccine.

View this table:
[Table 1:](http://medrxiv.org/content/early/2021/02/23/2021.02.19.21252095/T1)

Table 1: 
Mean values of parameters used to perform numerical simulations of the different scenarios.

View this table:
[Table 2:](http://medrxiv.org/content/early/2021/02/23/2021.02.19.21252095/T2)

Table 2: 
Initial conditions assumed for the different subpopulations using the USA current situation (Early-December).

### 3.1 Vaccination rate, efficacy and transmission rate scenarios

Here we present the results of the numerical simulations for different scenarios varying the inoculation rate, efficacy of the vaccine, percentage of infections that are asymptomatic, and the transmission rates. We consider two different plausible efficacies for the vaccines. We set the efficacy (*E*) to 80% and 94%. These values were chosen based on some results of vaccine trials and the current approved vaccines [3, 16, 69, 134, 149, 105, 135]. We could simulate scenarios with lower efficacies if we desire and based on the fact that the FDA established a minimum efficacy threshold of at least 50% [58, 99]. We also vary the inoculation rate (vaccination pace) to test different potential vaccination program scenarios [134, 106]. It is important to remark that despite the plans that health institutions make regarding vaccination, there are uncertainties present in the logistics [134, 95, 141, 106]. For instance, currently there is a significant delay in coronavirus vaccinations while hospitalizations continue to set records in the USA [106]. Therefore, here we considered two different plausible inoculation rates based on the current situation. Specifically, we chose vaccination rates of two and four millions per week. It is important to mention that even though these rates might not be 100% accurate, this approach helps to elucidate the impact of the inoculation rate on the main outcomes of the COVID-19 pandemic under a vaccination program.

Regarding the values of the SARS-CoV-2 virus transmission rate that plays an important role in the value of the effective reproduction number ℛ*t*, we chose two different values. These two values of the SARS-CoV-2 virus transmission rate between humans correspond to two different reproduction numbers ℛ. Thus, we can relate them using the following equation ![Graphic][2]</img>. It is important to remark that the effective reproduction number ℛ*t* varies over the time, and several methods have been proposed to compute it [24, 146, 120, 94, 75, 26, 131]. For instance, under certain conditions ℛ*t* = ℛ*S*(*t*)*/N*, which relates the value of the virus transmissibility *β* to the effective reproduction number [146]. It is important to remark that independently of the method that is used to compute the effective reproduction number ℛ*t* all of them show that it depends explicitly or implicitly on the value of the transmission rates (*β**s*). Therefore varying these rates implies a variation in the basic reproduction number ℛ and on the effective reproduction number ℛ*t*. Thus, we are considering different scenarios regarding the risk of becoming infected by the SARS-CoV-2.

The aim is to test the impact of the inoculation rate and efficacy of the vaccine under two different SARS-CoV-2 virus transmission rate scenarios. As it has been mentioned before, there are a lot of uncertainties in the transmission rates for different regions and they vary over the time depending on official and unofficial non-pharmaceutical interventions [40, 55, 85, 93, 120, 25, 75]. However, the approach used here allows us to understand the impact of inoculation rate and vaccine efficacy under two different transmission of the SARS-CoV-2 virus environments. Then results can be extrapolated to other environment settings.

### 3.2 Numerical simulation of scenarios

Here we present the results of the numerical simulations for different scenarios varying the inoculation rate and efficacy, under low transmission rate scenarios. Table 3 shows the peak of the number of infected (symptomatic plus asymptomatic) and hospitalized people for the different vaccine efficacies and inoculation rates. In addition, the number of deaths and recovered cases at the end of the simulation period of 365 days is shown. In this Table we considered several scenarios with two different transmission rates for *β**I*, assumed that the infectiousness of the asymptomatic individuals is the same as the symptomatic (*β**A* = *β**I*), and that the percentage of infections that are asymptomatic is 50% (*a* = 0.5). It can be seen that the impact of the vaccination rate is greater than the efficacy of the vaccine. For instance, observing the first row of Table 3 it is seen that the peak of the number of infected people is 5,765,525 under a scenario with transmission rate of *β**I* = 0.2, vaccine efficacy of *E* = 94% and a vaccination rate of two millions per week. Additionally, it is seen that this peak becomes 5,864,871 if the vaccine efficacy decreases to 80%. However, the peak of the number of infected people is 4,791,577, if under the same scenario we change the vaccination rate to four million per week. Thus, it can be seen than the impact of the vaccination rate in the number of infected people is larger than that from vaccine efficacy. We performed additional numerical simulations varying all the parameters in reasonable ranges, and the impact of the vaccination rate was always larger than the vaccine efficacy. In Figure 2, we show the peak of the number of infected and hospitalized individuals for a wide range of different vaccine efficacies and inoculation rates. In addition, the number of deaths and recovered cases are shown. The variation of all these outcomes is larger when the vaccination rate is varied. The results in Table 3 make sense since the effect of the vaccination rate and vaccine efficacy reduce the number of infected, hospitalized and deaths.

View this table:
[Table 3:](http://medrxiv.org/content/early/2021/02/23/2021.02.19.21252095/T3)

Table 3: 
Impact of the inoculation rate (Vac.) and the efficacy (ϵ) of the vaccine on the peak of the infected (I) and hospitalized (H) subpopulations. In addition, on the deaths (D) and recovered (R) cases. In these scenarios, the infectiousness of the asymptomatic individuals relative to symptomatic is 100% (βA = βI), and the percentage of infections that are asymptomatic is 50%.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/23/2021.02.19.21252095/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2021/02/23/2021.02.19.21252095/F1)

Figure 1: 
Diagram for the COVID-19 mathematical model (1). The boxes represent the subpopulation and the arrows the transition between the subpopulations. Arrows are labeled by their corresponding model parameters.

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/23/2021.02.19.21252095/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2021/02/23/2021.02.19.21252095/F2)

Figure 2: 
Impact of the inoculation rate and efficacy of the vaccine on the peak of the infected (I) and hospitalized (H) subpopulations. In addition, on the deaths (D) and recovered (R) cases. In these scenarios, the infectiousness of the asymptomatic individuals relative to symptomatic is 100% (βA = βI), and the percentage of infections that are asymptomatic is 50%. In all these scenarios the transmission rate is β = 0.2.

We present additional Tables for different scenarios varying the transmission of the asymptomatic carriers and the proportion of asymptomatic individuals. The Table 4 shows the same outcomes that we mentioned above, but we consider now that the infectiousness of the asymptomatic individuals relative to symptomatic is 75%. The impact of the vaccination rate is greater than the efficacy of the vaccine. The first row of Table 4 shows that the peak of the number of infected people is 2,932,727 when transmission rate is *β**I* = 0.2, the vaccine efficacy is *E* = 94% and the vaccination rate is two millions per week. This metric is just 2,954,850 if the vaccine efficacy decreases to 80%. This is a small change if we compare it to that when the rate of vaccination is increased to four million per week. Thus, this scenario also supports the importance of a high vaccination rate. Figure 3 shows the different outcomes for a wide range of different vaccine efficacies and inoculation rates. The variation of all these outcomes is larger when the vaccination rate is varied.

View this table:
[Table 4:](http://medrxiv.org/content/early/2021/02/23/2021.02.19.21252095/T4)

Table 4: 
Impact of the inoculation rate (Vac.) and the efficacy (ϵ) of the vaccine on the peak of the infected (I) and hospitalized (H) subpopulations. In addition, on the deaths (D) and recovered (R) cases. In these scenarios, the infectiousness of the asymptomatic individuals relative to symptomatic is 75% (βA = βI), and the percentage of infections that are asymptomatic is 50%.

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/23/2021.02.19.21252095/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2021/02/23/2021.02.19.21252095/F3)

Figure 3: 
Impact of the inoculation rate and efficacy of the vaccine on the peak of the infected (I) and hospitalized (H) subpopulations. In addition, on the deaths (D) and recovered (R) cases. In these scenarios, the infectiousness of the asymptomatic individuals relative to symptomatic is 75% (βA = βI), and the percentage of infections that are asymptomatic is 50%. In all these scenarios the transmission rate is β = 0.2.

Finally, Table 5 and Table 6, show the outcome when the percentage of infections that are asymptomatic is 40% (*a* = 0.4) and the infectiousness of the asymptomatic individuals relative to symptomatic is 100% and 75% respectively. The numerical simulation results show similar trends to the two previously studied cases. Again, that the impact of the vaccination rate is greater than the efficacy of the vaccine can be observed. This qualitative effect can be seen under a variety of scenarios regarding vaccine efficacy and vaccination rate in Figures 4 and 5, respectively. The numerical simulations include many different parameter values for the infectiousness of asymptomatic individuals, percentage of infections that are asymptomatic, efficacy of the vaccine, and the vaccination rate. Thus, uncertainty in these parameters has been considered in this study.

View this table:
[Table 5:](http://medrxiv.org/content/early/2021/02/23/2021.02.19.21252095/T5)

Table 5: 
Impact of the inoculation rate (Vac.) and the efficacy (ϵ) of the vaccine on the peak of the infected (I) and hospitalized (H) subpopulations. In addition, on the deaths (D) and recovered (R) cases. In these scenarios, the infectiousness of the asymptomatic individuals relative to symptomatic is 100% (βA = βI), and the percentage of infections that are asymptomatic is 40%.

View this table:
[Table 6:](http://medrxiv.org/content/early/2021/02/23/2021.02.19.21252095/T6)

Table 6: 
Impact of the inoculation rate (Vac.) and the efficacy (ϵ) of the vaccine on the peak of the infected (I) and hospitalized (H) subpopulations. In addition, on the deaths (D) and recovered (R) cases. In these scenarios, the infectiousness of the asymptomatic individuals relative to symptomatic is 75% (βA = βI), and the percentage of infections that are asymptomatic is 40%.

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/23/2021.02.19.21252095/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2021/02/23/2021.02.19.21252095/F4)

Figure 4: 
Impact of the inoculation rate and efficacy of the vaccine on the peak of the infected (I) and hospitalized (H) subpopulations. In addition, on the deaths (D) and recovered (R) cases. In these scenarios, the infectiousness of the asymptomatic individuals relative to symptomatic is 100% (βA = βI), and the percentage of infections that are asymptomatic is 40%. In all these scenarios the transmission rate is β = 0.2.

![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/23/2021.02.19.21252095/F5.medium.gif)

[Figure 5:](http://medrxiv.org/content/early/2021/02/23/2021.02.19.21252095/F5)

Figure 5: 
Impact of the inoculation rate and efficacy of the vaccine on the peak of the infected (I) and hospitalized (H) subpopulations. In addition, on the deaths (D) and recovered (R) cases. In these scenarios, the infectiousness of the asymptomatic individuals relative to symptomatic is 75% (βA = βI), and the percentage of infections that are asymptomatic is 40%. In all these scenarios the transmission rate is β = 0.2.

Based on the previous results we can conclude that under some plausible scenarios that the impact of the inoculation rate is more relevant to control the burden of the COVID-19 pandemic. Thus, these results suggest that health authorities should focus in increasing the inoculation rate in order to avert more infected people, hospitalizations and deaths. Our results agree with previous result under different assumptions, and with recommendations made by some scholars [99, 141].

## 4 Discussion

Currently there are authorized and recommended vaccines to prevent COVID-19 in the United States. The COVID-19 vaccination program started in early December. Depending on the specific vaccine, the people will get a second shot 3-4 weeks after the first in order to achieve the most protection against the disease caused by the SARS-CoV-2 virus [30, 106]. The vaccines against the SARS-CoV-2 virus have different efficacies and mechanisms of action [3, 16, 69, 71, 149, 105, 135]. Vaccination programs have been recently begun in other countries, using several different types of program and different outcomes can therefore be anticipated [1, 69, 137]. For instance, a vaccination program can focus first on health care workers or on elderly people [87]. However, whatever group the vaccination program targets first; there is an inoculation rate of the vaccine. In this study we propose a mathematical model to assess the impact of the vaccination programs as a function of the efficacy of the vaccine and the inoculation pace. The study of different COVID-19 vaccination programs is of paramount importance to reduce the burden of the COVID-19 pandemic. An optimal vaccination program helps to tackle the transmission of the SARS-CoV-2 virus in the population in an efficient way [9, 29, 69, 72, 137, 145]. It is important to mention that the vaccination programs in different regions or countries vary due to multiple constraints. For instance, there are vaccines that require special storage and transportation, and this affects the availability of the vaccine. Moreover, the current number of vaccine shots are not enough to vaccinate the whole world population [11, 16, 18, 53, 87].

In this article, we studied the impact of the vaccination pace and the efficacy of the vaccine on the dynamics of the COVID-19 pandemic. We studied the particular scenario of USA, but the methodology presented here can be extrapolated to other countries or regions. We were able to study different potential scenarios regarding the burden of the COVID-19 pandemic. We varied the inoculation rate, efficacy of the vaccine and the SARS-CoV-2 virus transmission rates. The constructed compartmental mathematical model allows the variation of the aforementioned factors, and using computational methodologies we obtained metrics that indicate which are the most important factors to decrease the burden of the COVID-19 pandemic.

We found that the efficacy of the vaccine and the vaccine inoculation rate have a high impact on the outcomes. However, the rate of vaccine administration has a larger impact on reducing the infected and hospitalized subpopulations. In a similar way, it has a greater impact on the number of deaths caused by the SARS-CoV-2 virus. Another important finding is that the impact of the inoculation rate and vaccine efficacy is larger for scenarios with higher SARS-CoV-2 virus transmission rates. Thus, our results suggest that health institutions need to focus in increasing the vaccine inoculation rate in the regions with higher rate of new infections. Our results are in accordance with previous recommendations made by some scholars [99, 141].

As expected from a vaccination program the results show that the benefits depend on how it is implemented and the efficacy of the vaccine. As we have mentioned, the total coverage of the population would depend on the production of the vaccine doses and the deployment of resources to execute the vaccination programs. In addition, there is a potential limitation of the total coverage due to the relctance of some parts of the population who have been influenced by doubting the science or by pressure from anti-vaccination groups [27, 53, 108, 137]. Our results also show that the impact of a COVID-19 vaccination program is highly dependent on the SARS-CoV-2 virus transmission rates and these affect the effective reproductive number ℛ*t* of the SARS-CoV-2 virus. Thus, it is important to educate the population about the importance of maintaining non-pharmaceutical control interventions such as the use of facial masks and physical distancing [30, 40, 73, 84, 110].

Furthermore, the additional benefit of a vaccine with 80% or 94% efficacy depends on the SARS-CoV-2 virus transmission rate as has been observed in this study. When we have low SARS-CoV-2 virus transmission rates (equivalently: lower effective reproductive number) the vaccine with 80% efficacy has a smaller impact on the COVID disease related metrics in comparison with a scenario with a high SARS-CoV-2 virus transmission rate (equivalently higher effective reproductive number). Thus, even with a highly effective vaccine it is important to maintain as low as possible the SARS-CoV-2 virus transmission rate to reduce the burden of the current pandemic. Moreover, if the immunity against the SARS-CoV-2 virus diminishes over the time, then lowering transmission rates is even more crucial.

The constructed compartmental model is a SEIR type but with some additional features such the compartment for asymptomatic cases. We expanded the model to include vaccinated people even if the vaccine is not effective in some subset. The SARS-CoV-2 virus spread is mainly driven by the values of the parameters, which have some uncertainty, as is usual in this type of epidemiological model. The uncertainty related to the COVID-19 pandemic is higher in comparison with other diseases such influenza due to the novelty of the SARS-CoV-2 virus. The parameter values were chosen from scientific literature. Despite the limitations of this type of mathematical model, they have been useful in many epidemics and are a classical method to deal with epidemics [57, 22, 28, 45, 34, 13, 36, 73, 12]. Some particular limitations of this study are that constant inoculation rates were used and the vaccination programs do not target any specific subpopulation. We hopefully anticipate that the inoculation rates will increase due to an increase in vaccine production and improving of the logistics. However, the vaccination programs might face several obstacles along the way. The model does not consider a subpopulation that is not willing to take the vaccine, and this has been an issue for other vaccines [27, 53, 108, 137]. Further studies are needed to extend the mathematical model for other vaccination programs. For instance, those that target first health care workers or specific age groups. This would require more parameters and therefore more uncertainty and details. In addition, our mathematical model does not consider the fact that immunity wanes. In fact, the US FDA recommends that follow-up of study participants should continue for as long as is feasible, to assess the duration of protection [58].

## 5 Conclusions

The results presented in this study show that the effectiveness of a COVID-19 vaccination program strongly depends on the vaccination rate and the efficacy of the vaccine. Moreover, the SARS-CoV-2 virus human transmission rates and consequently the effective reproductive number impact the outcome of the vaccination programs. It is important to remark that vaccination rate depends on many variables or resources such as health care facilities or logistical transportation aspects. On the other hand, the efficacy of the vaccine is out of the hands of health institutions and official entities. However, the rate of vaccine administration plays a more important role to reduce the burden of the COVID-19 pandemic. Our results show that health institutions need to focus in increasing the vaccine inoculation pace and create awareness in the population about the importance of the COVID-19 vaccines. In some countries the vaccination rate would be limited due to the availability of the vaccine. Currently, in the USA there are issues with the vaccination rate due to logistics, but not regarding availability [134, 106]. As we mentioned in the introduction, at some point there might be difficulties keeping a constant vaccination rate since a proportion of the population is not willing to be vaccinated. This topic is interesting and can be studied in the future.

The type of mathematical model, based on ordinary differential equations, used here suffers from the following limitations: exponential distributions in the transitions from one stage to another are implicitly assumed as well as homogeneous mixing in the population. Additionally, the behavior of individuals is averaged in order to avoid more complex models that in turn have their own limitations. For instance, individual agent based models have many parameters and in several cases the values of these parameters are very difficult to obtain. However, in some cases exponential transition between stages are not very far from reality. Despite the limitations of our model, we found valuable results to face the current COVID-19 pandemic. Support is given to characteristics of efficient vaccine campaigns. In particular, our study encourages governments and their health institutions to increase the pace of the vaccination in the population in order to diminish the consequences of the catastrophic COVID-19 pandemic.

## Data Availability

Data used in the article is public.

*   Received February 19, 2021.
*   Revision received February 19, 2021.
*   Accepted February 23, 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/)

## Bibliography

1.  [1]. Derrick Bary Abila et al. “We need to start thinking about promoting the demand, uptake, and equitable distribution of COVID-19 vaccines NOW!” In: Public Health in Practice (2020), p. 100063.
    
    

2.  [2]. Luis Acedo et al. “Uncertainty and sensitivity of the sexual behavior changes to the current human papillomavirus vaccination campaign in Spain”. In: Mathematical Methods in the Applied Sciences ().
    
    

3.  [3]. Dae-Gyun Ahn et al. “Current status of epidemiology, diagnosis, therapeutics, and vaccines for novel coronavirus disease 2019 (COVID-19)”. In: Journal of Microbiology and Biotechnology (2020).
    
    

4.  [4]. Manaf Al-Qahtani et al. “The prevalence of asymptomatic and symptomatic COVID-19 in a cohort of quarantined subjects”. In: International Journal of Infectious Diseases 102 (2020), pp. 285–288.
    
    

5.  [5]. V. Andreasen,  C. Viboud, and  L. Simonsen. “Epidemiologic Characterization of the 1918 Influenza Pandemic Summer Wave in Copenhagen:Implications for Pandemic Control Strategies”. In: The Journal of Infectious Diseases 197.2 (2008), pp. 270–278.
    
    

6.  [6]. Kirtikumar C Badgujar,  Vivek C Badgujar, and  Shamkant B Badgujar. “Vaccine development against coronavirus (2003 to present): An overview, recent advances, current scenario, opportunities and challenges”. In: Diabetes & Metabolic Syndrome: Clinical Research & Reviews 14.5 (2020), pp. 1361–1376.
    
    

7.  [7]. Yan Bai et al. “Presumed asymptomatic carrier transmission of COVID-19”. In: Jama 323.14 (2020), pp. 1406–1407.
    
    

8.  [8]. Georgios D Barmparis and  GP Tsironis. “Estimating the infection horizon of COVID-19 in eight countries with a data-driven approach”. In: Chaos, Solitons & Fractals (2020), p. 109842.
    
    

9.  [9]. Sarah M Bartsch et al. “Vaccine efficacy needed for a COVID-19 coronavirus vaccine to prevent or stop an epidemic as the sole intervention”. In: American journal of preventive medicine 59.4 (2020), pp. 493–503.
    
    

10. [10]. Zeno Bisoffi et al. “Sensitivity, specificity and predictive values of molecular and serological tests for COVID-19: a longitudinal study in emergency room”. In: Diagnostics 10.9 (2020), p. 669.
    
    

11. [11]. Barry R Bloom,  Glen J Nowak, and  Walter Orenstein. ““When Will We Have a Vaccine?”—Understanding Questions and Answers about Covid-19 Vaccination”. In: New England Journal of Medicine 383.23 (2020), pp. 2202–2204.
    
    

12. [12]. Fred Brauer,  Carlos Castillo-Chavez, and  Carlos Castillo-Chavez. Mathematical models in population biology and epidemiology. Vol. 40. Springer, 2001.
    
    

13. [13]. Tom Britton,  Frank Ball, and  Pieter Trapman. “A mathematical model reveals the influence of population heterogeneity on herd immunity to SARS-CoV-2”. In: Science 369.6505 (2020), pp. 846–849.
    
    

14. [14]. Shelby R Buckman et al. “Replicating and projecting the path of COVID-19 with a model-implied reproduction number”. In: Infectious Disease Modelling 5 (2020), pp. 635–651.
    
    

15. [15].United States Census Bureau. [https://www.census.gov/popclock/](https://www.census.gov/popclock/). 2020.
    
    

16. [16]. Talha Khan Burki. “The Russian vaccine for COVID-19”. In: The Lancet Respiratory Medicine 8.11 (2020), e85–e86.
    
    

17. [17]. Oyungerel Byambasuren et al. “Estimating the extent of asymptomatic COVID-19 and its potential for community transmission: systematic review and meta-analysis”. In: Official Journal of the Association of Medical Microbiology and Infectious Disease Canada 5.4 (2020), pp. 223–234.
    
    

18. [18]. Ewen Callaway. “OXFORD COVID VACCINE RESULTS PUZZLE SCIENTISTS”. In: Nature 588 (2020).
    
    

19. [19]. Wu-Chun Cao et al. “Disappearance of antibodies to SARS-associated coronavirus after recovery”. In: New England Journal of Medicine 357.11 (2007), pp. 1162–1163.
    
    

20. [20]. Ishanu Chattopadhyay et al. “Conjunction of factors triggering waves of seasonal influenza”. In: eLife 7 (2018).
    
    

21. [21]. Yi Chen et al. “The epidemiological characteristics of infection in close contacts of COVID-19 in Ningbo city”. In: Chin J Epidemiol 41.5 (2020), pp. 668–672.
    
    

22. [22]. G. Chowell et al. “Transmission dynamics of the great influenza pandemic of 1918 in Geneva, Switzerland: Assessing the effects of hypothetical interventions”. In: Journal of Theoretical Biology 241.2 (2006), pp. 193–204.
    
    

23. [23]. Candice Clarke et al. “High prevalence of asymptomatic COVID-19 infection in hemodialysis patients detected using serologic screening”. In: Journal of the American Society of Nephrology 31.9 (2020), pp. 1969–1975.
    
    

24. [24].Sebastián Contreras et al. “Real-time estimation of Rt for supporting public-health policies against COVID-19”. In: Frontiers in public health 8 (2020).
    
    

25. [25]. Aritra Das. “An approximation-based approach for periodic estimation of effective reproduction number: a tool for decision-making in the context of coronavirus disease 2019 (COVID-19) outbreak”. In: Public Health 185 (2020), pp. 199–201.
    
    

26. [26]. Paul L Delamater et al. “Complexity of the basic reproduction number (R0)”. In: Emerging infectious diseases 25.1 (2019), p. 1.
    
    

27. [27]. Terence S Dermody,  Daniel DiMaio, and  Lynn W Enquist. “Vaccine Safety, Efficacy, and Trust Take Time”. In: Annual Review of Virology 8 (2020).
    
    

28. [28]. O Diekmann,  JAP Heesterbeek, and  MG Roberts. “The construction of next-generation matrices for compartmental epidemic models”. In: Journal of the Royal Society Interface (2009), rsif20090386.
    
    

29. [29]. Ener Cagri Dinleyici et al. “Vaccines and routine immunization strategies during the COVID-19 pandemic”. In: Human vaccines & immunotherapeutics (2020), pp. 1–8.
    
    

30. [30].Centers for Disease Control and Prevention. [https://www.cdc.gov/coronavirus/2019-nCoV/index.html](https://www.cdc.gov/coronavirus/2019-nCoV/index.html). 2020.
    
    

31. [31]. Hana M Dobrovolny. “Modeling the role of asymptomatics in infection spread with application to SARS-CoV-2”. In: Plos one 15.8 (2020), e0236976.
    
    

32. [32]. Hana M Dobrovolny. “Quantifying the effect of remdesivir in rhesus macaques infected with SARS-CoV-2”. In: Virology 550 (2020), pp. 61–69.
    
    

33. [33]. Lucy van Dorp et al. “Emergence of genomic diversity and recurrent mutations in SARS-CoV-2”. In: Infection, Genetics and Evolution (2020), p. 104351.
    
    

34. [34]. P. van den Driessche and  James Watmough. “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission”. In: Mathematical Biosciences 180.12 (2002), pp. 29 –48.
    
    

35. [35]. Eve Dubé and  Noni E MacDonald. “How can a global pandemic affect vaccine hesitancy?” In: Expert Review of Vaccines 19.10 (2020), pp. 899–901.
    
    

36. [36]. Matthias Ehrhardt,  Ján Gašper, and  Soňa Kilianová. “SIR-based mathematical modeling of infectious diseases with vaccination and waning immunity”. In: Journal of Computational Science 37 (2019), p. 101027.
    
    

37. [37]. Christel Faes et al. “Time between symptom onset, hospitalisation and recovery or death: Statistical analysis of belgian covid-19 patients”. In: International Journal of Environmental Research and Public Health 17.20 (2020), p. 7560.
    
    

38. [38]. Baylor Fain and  Hana M Dobrovolny. “Initial Inoculum and the Severity of COVID-19: A Mathematical Modeling Study of the Dose-Response of SARS-CoV-2 Infections”. In: Epidemiologia 1.1 (2020), pp. 5–15.
    
    

39. [39]. Jeremy Samuel Faust and  Carlos del Rio. “Assessment of Deaths From COVID-19 and From Seasonal Influenza”. In: JAMA Internal Medicine (May 2020). ISSN: 2168-6106. DOI: 10.1001/jamainternmed.2020.2306. eprint: [https://jamanetwork.com/journals/jamainternalmedicine/articlepdf/2766121/jamainternal\_faust\_2020\_vp\_200013.p](https://jamanetwork.com/journals/jamainternalmedicine/articlepdf/2766121/jamainternal\_faust_2020_vp_200013.p). URL: [https://doi.org/10.1001/jamainternmed.2020.2306](https://doi.org/10.1001/jamainternmed.2020.2306).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamainternmed.2020.2306&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32407441&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F23%2F2021.02.19.21252095.atom) 

40. [40]. Neil M Ferguson et al. “Impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand”. In: Imperial College, London. DOI: [https://doi.org/10.25561/77482](https://doi.org/10.25561/77482) (2020).
    
    

41. [41]. Simon James Fong et al. “Finding an accurate early forecasting model from small dataset: A case of 2019-ncov novel coronavirus outbreak”. In: International Journal of Interactive Multimedia and Artificial Intelligence 6.1 (2020), pp. 132–40.
    
    

42. [42]. Monica Gandhi and  George W Rutherford. “Facial masking for Covid-19—potential for “variolation” as we await a vaccine”. In: New England Journal of Medicine 383.18 (2020), e101.
    
    

43. [43]. Monica Gandhi,  Deborah S. Yokoe, and  Diane V. Havlir. “Asymptomatic Transmission, the Achilles’ Heel of Current Strategies to Control Covid-19”. In: New England Journal of Medicine 382.22 (2020), pp. 2158–2160. DOI: 10.1056/NEJMe2009758.eprint: [https://doi.org/10.1056/NEJMe2009758](https://doi.org/10.1056/NEJMe2009758). URL: [https://doi.org/10.1056/NEJMe2009758](https://doi.org/10.1056/NEJMe2009758).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMe2009758.eprint&link_type=DOI) 

44. [44]. Navid Ghaffarzadegan and  Hazhir Rahmandad. “Simulation-based estimation of the early spread of COVID-19 in Iran: actual versus confirmed cases”. In: System Dynamics Review 36.1 (2020), pp. 101–129.
    
    

45. [45]. Gilberto González-Parra,  A. Arenas, and  Lupe Segovia Diego F. Aranda. “Modeling the epidemic waves of AH1N1/09 influenza around the world”. In: Spatial and Spatio-temporal Epidemiology 2 (2011), pp. 219– 226.
    
    

46. [46]. Gilberto González-Parra,  Abraham J Arenas, and  Benito M Chen-Charpentier. “A fractional order epidemic model for the simulation of outbreaks of influenza A (H1N1)”. In: Mathematical Methods in the Applied Sciences 37.15 (2014), pp. 2218–2226.
    
    

47. [47]. Gilberto González-Parra and  Hana M Dobrovolny. “Assessing uncertainty in A2 respiratory syncytial virus viral dynamics”. In: Computational and mathematical methods in medicine 2015 (2015).
    
    

48. [48]. Gilberto González-Parra et al. “Modeling the social obesity epidemic with stochastic networks”. In: Physica A: Statistical Mechanics and its Applications 389.17 (2010), pp. 3692–3701.
    
    

49. [49]. Jens Grauer,  Hartmut Löwen, and  Benno Liebchen. “Strategic spatiotemporal vaccine distribution increases the survival rate in an infectious disease like Covid-19”. In: Scientific reports 10.1 (2020), pp. 1–10.
    
    

50. [50]. Nathan D Grubaugh,  William P Hanage, and  Angela L Rasmussen. “Making sense of mutation: what D614G means for the COVID-19 pandemic remains unclear”. In: Cell 182.4 (2020), pp. 794–795.
    
    

51. [51]. Madhu Gupta et al. “The need for COVID-19 research in low-and middle-income countries”. In: Global Health Research and Policy 5.1 (2020), pp. 1–4.
    
    

52. [52]. Dongsheng Han et al. “COVID-19: Insight into the asymptomatic SARS-COV-2 infection and transmission”. In: International Journal of Biological Sciences 16.15 (2020), p. 2803.
    
    

53. [53]. Azizul Haque and  Anudeep B Pant. “Efforts at COVID-19 Vaccine Development: Challenges and Successes”. In: Vaccines 8.4 (2020), p. 739.
    
    

54. [54]. Imran Hasanoglu et al. “Higher viral loads in asymptomatic COVID-19 patients might be the invisible part of the iceberg”. In: Infection (2020), pp. 1–10.
    
    

55. [55]. Johannes Haushofer and C  Jessica E Metcalf. “Which interventions work best in a pandemic?” In: Science 368.6495 (2020), pp. 1063–1065.
    
    

56. [56]. Fiona P Havers et al. “Seroprevalence of antibodies to SARS-CoV-2 in 10 sites in the United States, March 23-May 12, 2020”. In: JAMA Internal Medicine (2020).
    
    

57. [57]. Herbert W Hethcote. “Mathematics of infectious diseases”. In: SIAM Review 42 (2005), pp. 599–653.
    
    

58. [58]. Susanne H Hodgson et al. “What defines an efficacious COVID-19 vaccine? A review of the challenges assessing the clinical efficacy of vaccines against SARS-CoV-2”. In: The Lancet Infectious Diseases (2020).
    
    

59. [59]. Inga Holmdahl and  Caroline Buckee. “Wrong but useful-what covid-19 epidemiologic models can and cannot tell us”. In: New England Journal of Medicine (2020).
    
    

60. [60]. Amirhoshang Hoseinpour Dehkordi et al. “Understanding epidemic data and statistics: A case study of COVID-19”. In: Journal of medical virology (2020).
    
    

61. [61]. Can Hou et al. “The effectiveness of quarantine of Wuhan city against the Corona Virus Disease 2019 (COVID-19): A well-mixed SEIR model analysis”. In: Journal of medical virology (2020).
    
    

62. [62]. Lei Huang et al. “Rapid asymptomatic transmission of COVID-19 during the incubation period demon-strating strong infectivity in a cluster of youngsters aged 16-23 years outside Wuhan and characteristics of young patients with COVID-19: a prospective contact-tracing study”. In: Journal of Infection (2020).
    
    

63. [63].Zeynep Gökçe İşlier,  Refik Güllü, and  Wolfgang Hörmann. “An exact and implementable computation of the final outbreak size distribution under Erlang distributed infectious period”. In: Mathematical Biosciences (2020), p. 108363.
    
    

64. [64]. Nicholas P Jewell,  Joseph A Lewnard, and  Britta L Jewell. Caution warranted: using the Institute for Health Metrics and Evaluation model for predicting the course of the COVID-19 pandemic. 2020.
    
    

65. [65]. Xiaoqian Jia et al. “Modeling the prevalence of asymptomatic COVID-19 infections in the Chinese mainland”. In: The Innovation 1.2 (2020), p. 100026.
    
    

66. [66]. Lindsay Kim et al. “Risk factors for intensive care unit admission and in-hospital mortality among hospi-talized adults identified through the US coronavirus disease 2019 (COVID-19)-associated hospitalization surveillance network (COVID-NET)”. In: Clinical Infectious Diseases (2020).
    
    

67. [67]. Yunhwan Kim,  Ana Vivas Barber, and  Sunmi Lee. “Modeling influenza transmission dynamics with media coverage data of the 2009 H1N1 outbreak in Korea”. In: Plos one 15.6 (2020), e0232580.
    
    

68. [68]. Ryo Kinoshita et al. “Containment, Contact Tracing and Asymptomatic Transmission of Novel Coronavirus Disease (COVID-19): A Modelling Study”. In: Journal of Clinical Medicine 9.10 (2020), p. 3125.
    
    

69. [69]. Maria Deloria Knoll and  Chizoba Wonodi. “Oxford–AstraZeneca COVID-19 vaccine efficacy”. In: The Lancet (2020).
    
    

70. [70]. Archana Koirala et al. “Vaccines for COVID-19: The current state of play”. In: Paediatric respiratory reviews 35 (2020), pp. 43–49.
    
    

71. [71]. Bette Korber et al. “Tracking changes in SARS-CoV-2 Spike: evidence that D614G increases infectivity of the COVID-19 virus”. In: Cell 182.4 (2020), pp. 812–827.
    
    

72. [72]. Philip Krause et al. “COVID-19 vaccine trials should seek worthwhile efficacy”. In: The Lancet 396.10253 (2020), pp. 741–743.
    
    

73. [73]. Adam J Kucharski et al. “Early dynamics of transmission and control of COVID-19: a mathematical modelling study”. In: The lancet infectious diseases (2020).
    
    

74. [74]. Toshikazu Kuniya. “Prediction of the Epidemic Peak of Coronavirus Disease in Japan, 2020”. In: Journal of Clinical Medicine 9.3 (2020), p. 789.
    
    

75. [75]. Alessia Lai et al. “Early phylogenetic estimate of the effective reproduction number of SARS-CoV-2”. In: Journal of medical virology 92.6 (2020), pp. 675–679.
    
    

76. [76]. Kian Boon Law et al. “Tracking the early depleting transmission dynamics of COVID-19 with a time-varying SIR model”. In: Scientific reports 10.1 (2020), pp. 1–11.
    
    

77. [77]. Qianqian Li et al. “The impact of mutations in SARS-CoV-2 spike on viral infectivity and antigenicity”. In: Cell 182.5 (2020), pp. 1284–1294.
    
    

78. [78]. Qun Li et al. “Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia”. In: New England Journal of Medicine (2020).
    
    

79. [79]. Marc Lipsitch and  Natalie E Dean. “Understanding COVID-19 vaccine efficacy”. In: Science 370.6518 (2020), pp. 763–765.
    
    

80. [80]. Dianbo Liu et al. “Real-Time Forecasting of the COVID-19 Outbreak in Chinese Provinces: Machine Learning Approach Using Novel Digital Data and Estimates From Mechanistic Models”. In: Journal of medical Internet research 22.8 (2020), e20285.
    
    

81. [81]. Andrea Lovato,  Cosimo de Filippis, and  Gino Marioni. “Upper airway symptoms in coronavirus disease 2019 (COVID-19)”. In: American Journal of Otolaryngology (2020).
    
    

82. [82]. Nicole Lurie et al. “Developing Covid-19 vaccines at pandemic speed”. In: New England Journal of Medicine 382.21 (2020), pp. 1969–1973.
    
    

83. [83]. Pooja Mahajan and  Jyotsna Kaushal. “Epidemic trend of COVID-19 transmission in India during lockdown-1 phase”. In: Journal of Community Health (2020), pp. 1–10.
    
    

84. [84]. Sandip Mandal et al. “Prudent public health intervention strategies to control the coronavirus disease 2019 transmission in India: A mathematical model-based approach”. In: The Indian journal of medical research 151.2-3 (2020), p. 190.
    
    

85. [85]. S Marimuthu et al. “Modelling of reproduction number for COVID-19 in India and high incidence states”. In: Clinical Epidemiology and Global Health (2020).
    
    

86. [86]. J. A. Martin et al. “Predictive modeling of Covid-19 data in the US: Adaptive phase-space approach”. In: National Vital Statistics System 68 (2019), pp. 1–47.
    
    

87. [87]. Laura Matrajt et al. “Vaccine optimization for COVID-19, who to vaccinate first?” In: medRxiv (2020).
    
    

88. [88]. David Mc Evoy et al. “The relative infectiousness of asymptomatic SARS-CoV-2 infected persons compared with symptomatic individuals: A rapid scoping review.” In: medRxiv (2020).
    
    

89. [89]. Nir Menachemi et al. “Population point prevalence of SARS-CoV-2 infection based on a statewide random sample—Indiana, April 25–29, 2020”. In: Morbidity and Mortality Weekly Report 69.29 (2020), p. 960.
    
    

90. [90]. Kenji Mizumoto et al. “Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020”. In: Eurosurveillance 25.10 (2020), p. 2000180.
    
    

91. [91]. Zindoga Mukandavire et al. “Quantifying early COVID-19 outbreak transmission in South Africa and exploring vaccine efficacy scenarios”. In: PloS one 15.7 (2020), e0236003.
    
    

92. [92]. James D Murray. Mathematical biology I: an introduction, Vol. 17 of interdisciplinary applied mathematics. New York, NY, USA: Springer, 2002.
    
    

93. [93]. Jiaming Na et al. “Probabilistic approximation of effective reproduction number of COVID-19 using daily death statistics”. In: Chaos, Solitons & Fractals 140 (2020), p. 110181.
    
    

94. [94]. Farid Najafi et al. “Serial interval and time-varying reproduction number estimation for COVID-19 in western Iran”. In: New microbes and new infections 36 (2020), p. 100715.
    
    

95. [95]. Roxanne Nelson. “COVID-19 disrupts vaccine delivery”. In: The Lancet Infectious Diseases 20.5 (2020), p. 546.
    
    

96. [96]. Daniel P Oran and  Eric J Topol. “Prevalence of Asymptomatic SARS-CoV-2 Infection: A Narrative Review”. In: Annals of Internal Medicine (2020).
    
    

97. [97].Esteban Orenes-Piñero et al. “Evidences of SARS-CoV-2 virus air transmission indoors using several untouched surfaces: A pilot study”. In: Science of The Total Environment 751 (), p. 142317.
    
    

98. [98]. Maria Pachetti et al. “Impact of lockdown on Covid-19 case fatality rate and viral mutations spread in 7 countries in Europe and North America”. In: Journal of Translational Medicine 18.1 (2020), pp. 1–7.
    
    

99. [99]. A David Paltiel et al. “Clinical Outcomes Of A COVID-19 Vaccine: Implementation Over Efficacy: Study examines how definitions and thresholds of vaccine efficacy, coupled with different levels of implementation effectiveness and background epidemic severity, translate into outcomes.” In: Health Affairs (2020), pp. 10– 1377.
    
    

100.[100]. Sang Woo Park et al. “The time scale of asymptomatic transmission affects estimates of epidemic potential in the COVID-19 outbreak”. In: Epidemics (2020), p. 100392.
    
    

101.[101]. Malik Peiris and  Gabriel M Leung. “What can we expect from first-generation COVID-19 vaccines?” In: The Lancet 396.10261 (2020), pp. 1467–1469.
    
    

102.[102]. Lubna Pinky and  Hana M Dobrovolny. “SARS-CoV-2 coinfections: Could influenza and the common cold be beneficial?” In: Journal of Medical Virology (2020).
    
    

103.[103]. Silvio Daniel Pitlik. “COVID-19 Compared to Other Pandemic Diseases”. In: Rambam Maimonides Medical Journal 11.3 (2020).
    
    

104.[104]. Jessica A Plante et al. “Spike mutation D614G alters SARS-CoV-2 fitness”. In: Nature (2020), pp. 1–6.
    
    

105.[105]. Fernando P Polack et al. “Safety and efficacy of the BNT162b2 mRNA covid-19 vaccine”. In: New England Journal of Medicine (2020).
    
    

106.[106].The Washington Post. [https://www.washingtonpost.com/nation/2020/12/29/coronavirus-covid-updates/](https://www.washingtonpost.com/nation/2020/12/29/coronavirus-covid-updates/). 2020.
    
    

107.[107]. Eugene B Postnikov. “Estimation of COVID-19 dynamics “on a back-of-envelope”: Does the simplest SIR model provide quantitative parameters and predictions?” In: Chaos, Solitons & Fractals 135 (2020), p. 109841.
    
    

108.[108]. Atika Qazi et al. “Analyzing situational awareness through public opinion to predict adoption of social distancing amid pandemic COVID-19”. In: Journal of medical virology (2020).
    
    

109.[109]. Pipetius Quah,  Andrew Li, and  Jason Phua. “Mortality rates of patients with COVID-19 in the intensive care unit: a systematic review of the emerging literature”. In: Critical Care 24 (2020), pp. 1–4.
    
    

110.[110]. Mohammad Mahmudur Rahman et al. “Impact of control strategies on COVID-19 pandemic and the SIR model based forecasting in Bangladesh.” In: medRxiv (2020).
    
    

111.[111]. Li Ran et al. “Risk factors of healthcare workers with corona virus disease 2019: a retrospective cohort study in a designated hospital of Wuhan in China”. In: Clinical Infectious Diseases (2020).
    
    

112.[112]. Ruy Freitas Reis et al. “Characterization of the COVID-19 pandemic and the impact of uncertainties, mitigation strategies, and underreporting of cases in South Korea, Italy, and Brazil”. In: Chaos, Solitons & Fractals (2020), p. 109888.
    
    

113.[113].Reuters. [https://reut.rs/3pBySaG](https://reut.rs/3pBySaG). 2020.
    
    

114.[114]. D. Rios-Doria and  G. Chowell. “Qualitative analysis of the level of cross-protection between epidemic waves of the 1918 − 1919 influenza pandemic”. In: Journal of Theoretical Biology 261.4 (2009), pp. 584–592.
    
    

115.[115]. Jennifer D Roberts and  Shadi O Tehrani. “Environments, behaviors, and inequalities: reflecting on the impacts of the influenza and coronavirus pandemics in the United States”. In: International journal of environmental research and public health 17.12 (2020), p. 4484.
    
    

116.[116]. Mick Roberts et al. “Nine challenges for deterministic epidemic models”. In: Epidemics 10 (2015), pp. 49– 53.
    
    

117.[117]. Weston C Roda et al. “Why is it difficult to accurately predict the COVID-19 epidemic?” In: Infectious Disease Modelling (2020).
    
    

118.[118]. Lionel Roques et al. “Using early data to estimate the actual infection fatality ratio from COVID-19 in France”. In: Biology 9.5 (2020), p. 97.
    
    

119.[119]. Steven Sanche et al. “The novel coronavirus, 2019-nCoV, is highly contagious and more infectious than initially estimated”. In: arXiv preprint arXiv:2002.03268 (2020).
    
    

120.[120].Luis Santamaría and  Joaquín Hortal. “COVID-19 effective reproduction number dropped during Spain’s nationwide dropdown, then spiked at lower-incidence regions”. In: Science of the Total Environment 751 (2020), p. 142257.
    
    

121.[121]. Joshua L Santarpia et al. “Aerosol and surface contamination of SARS-CoV-2 observed in quarantine and isolation care”. In: Scientific Reports 110.10 (2020), p. 12732.
    
    

122.[122]. Gino Sartor et al. “COVID-19 in Italy: Considerations on official data”. In: International journal of infectious diseases 98 (2020), pp. 188–190.
    
    

123.[123]. M Shakiba et al. “Prevalence of asymptomatic COVID-19 infection using a seroepidemiological survey”. In: Epidemiology & Infection (2020), pp. 1–7.
    
    

124.[124]. Siyao Shao et al. “Risk assessment of airborne transmission of COVID-19 by asymptomatic individuals under different practical settings”. In: Journal of aerosol science 151 (2020), p. 105661.
    
    

125.[125]. Matthew Sperrin,  Stuart W Grant, and  Niels Peek. “Prediction models for diagnosis and prognosis in Covid-19”. In: BMJ 369 (2020). DOI: 10.1136/bmj.m1464.eprint: [https://www.bmj.com/content/369/bmj.m1464.full.pdf](https://www.bmj.com/content/369/bmj.m1464.full.pdf). xURL: [https://www.bmj.com/content/369/bmj.m1464](https://www.bmj.com/content/369/bmj.m1464).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1136/bmj.m1464.eprint&link_type=DOI) 

126.[126]. Richard OJH Stutt et al. “A modelling framework to assess the likely effectiveness of facemasks in combination with lock-down in managing the COVID-19 pandemic”. In: Proceedings of the Royal Society A 476.2238 (2020), p. 20200376.
    
    

127.[127]. Elena Surkova,  Vladyslav Nikolayevskyy, and  Francis Drobniewski. “False-positive COVID-19 results: hidden problems and costs”. In: The Lancet Respiratory Medicine 8.12 (2020), pp. 1167–1168.
    
    

128.[128]. Kirtimaan Syal. “COVID-19: herd immunity and convalescent plasma transfer therapy”. In: Journal of Medical Virology (2020).
    
    

129.[129]. Keke Tang,  Yining Huang, and  Meilian Chen. “Novel Coronavirus 2019 (Covid-19) epidemic scale estimation: topological network-based infection dynamic model”. In: medRxiv (2020).
    
    

130.[130]. Zhou Tang,  Xianbin Li, and  Houqiang Li. “Prediction of New Coronavirus Infection Based on a Modified SEIR Model”. In: medRxiv (2020).
    
    

131.[131].IHME COVID-19 Forecasting Team. “Modeling COVID-19 scenarios for the United States”. In: Nature Medicine (2020).
    
    

132.[132]. Samuel Cota Teixeira. “Mild and asymptomatic cases of COVID-19 are potential threat for faecal–oral transmission”. In: Brazilian Journal of Infectious Diseases 24.4 (2020), pp. 368–368.
    
    

133.[133]. Diana M Thomas et al. “A primer on COVID-19 Mathematical Models”. In: Obesity (2020).
    
    

134.[134].New York Times. U.S. inoculations are moving more slowly than expected, with 1 million vaccinated so far. [https://nyti.ms/370L7ur](https://nyti.ms/370L7ur). 2020.
    
    

135.[135]. John S Tregoning et al. “Vaccines for COVID-19”. In: Clinical & Experimental Immunology 202.2 (2020), pp. 162–192.
    
    

136.[136].Johns Hopkins University and Medicine. [https://coronavirus.jhu.edu](https://coronavirus.jhu.edu). 2020.
    
    

137.[137]. Pierre Verger and Eve Dubé. “Restoring confidence in vaccines in the COVID-19 era”. In: Expert Review of Vaccines 0.0 (2020). PMID: 32940574, pp. 1–3. DOI:10.1080/14760584.2020.825945.eprint: [https://doi.org/10.1080/14760584.2020.1825945](https://doi.org/10.1080/14760584.2020.1825945).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/14760584.2020.825945.eprint&link_type=DOI) 

138.[138]. Froma Walsh. “Loss and resilience in the time of COVID-19: Meaning making, hope, and transcendence”. In: Family process 59.3 (2020), pp. 898–911.
    
    

139.[139]. Xutong Wang et al. “Impact of social distancing measures on coronavirus disease healthcare demand, central Texas, USA”. In: Emerging Infectious Diseases 26.10 (2020), p. 2361.
    
    

140.[140]. Helen J Wearing,  Pejman Rohani, and  Matt J Keeling. “Appropriate models for the management of infectious diseases”. In: PLoS Med 2.7 (2005), e174.
    
    

141.[141]. Rebecca L Weintraub et al. “COVID-19 Vaccine To Vaccination: Why Leaders Must Invest In Delivery Strategies Now: Analysis describe lessons learned from past pandemics and vaccine campaigns about the path to successful vaccine delivery for COVID-19.” In: Health Affairs (2020), pp. 10–1377.
    
    

142.[142]. Bárbara Bowoniuk Wiegand and  Solange do Carmo Bowoniuk Wiegand. “Losses in times of a pandemic: Bioethical perceptions about mourning caused by COVID-19 in Brazil”. In: Voices in Bioethics 6 (2020).
    
    

143.[143]. Shengli Xia et al. “Effect of an inactivated vaccine against SARS-CoV-2 on safety and immunogenicity outcomes: interim analysis of 2 randomized clinical trials”. In: Jama 324.10 (2020), pp. 951–960.
    
    

144.[144]. JQ Xu et al. “Mortality in the United States, 2018.” In: NCHS Data Brief 355 (2020).
    
    

145.[145]. Gavin Yamey et al. “Ensuring global access to COVID-19 vaccines”. In: The Lancet 395.10234 (2020), pp. 1405–1406.
    
    

146.[146]. Ping Yan and  Gerardo Chowell. “Beyond the Initial Phase: Compartment Models for Disease Transmission”. In: Quantitative Methods for Investigating Infectious Disease Outbreaks. Springer, 2019, pp. 135–182.
    
    

147.[147]. HY Yang and  GC Duan. “Analysis on the epidemic factors for the corona virus disease”. In: Zhonghua yu Fang yi xue za zhi [Chinese Journal of Preventive Medicine] 54 (2020), E021–E021.
    
    

148.[148]. Baligh R Yehia et al. “Association of race with mortality among patients hospitalized with coronavirus disease 2019 (COVID-19) at 92 US hospitals”. In: JAMA network open 3.8 (2020), e2018039–e2018039.
    
    

149.[149]. Jingyou Yu et al. “DNA vaccine protection against SARS-CoV-2 in rhesus macaques”. In: Science (2020).
    
    

150.[150]. Yusen Zhai and  Xue Du. “Loss and grief amidst COVID-19: A path to adaptation and resilience”. In: Brain, Behavior, and Immunity (2020).
    
    

151.[151]. Lizhou Zhang et al. “The D614G mutation in the SARS-CoV-2 spike protein reduces S1 shedding and increases infectivity”. In: bioRxiv (2020). DOI: 10.1101/2020.06.12.148726.eprint: [https://www.biorxiv.org/content/early/2020/06/12/2020.06.12.148726.full.pdf](https://www.biorxiv.org/content/early/2020/06/12/2020.06.12.148726.full.pdf). URL: [https://www.biorxiv.org/content/early/2020/06/12/2020.06.12.148726](https://www.biorxiv.org/content/early/2020/06/12/2020.06.12.148726).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2020.06.12.148726.eprint&link_type=DOI) 

152.[152]. Xiaonan Zhang et al. “Viral and host factors related to the clinical outcome of COVID-19”. In: Nature 583.7816 (2020), pp. 437–440.
    
    

153.[153]. Yong Zhang et al. “Applicability of time fractional derivative models for simulating the dynamics and mitigation scenarios of COVID-19”. In: Chaos, Solitons & Fractals (2020), p. 109959.
    
    

154.[154]. Hong-jun Zhao et al. “COVID-19: asymptomatic carrier transmission is an underestimated problem”. In: Epidemiology & Infection (2020), pp. 1–7.
    
    

155.[155]. Linhao Zhong et al. “Early Prediction of the 2019 Novel Coronavirus Outbreak in the Mainland China based on Simple Mathematical Model”. In: IEEE Access (2020).
    
    

156.[156]. Fei Zhou et al. “Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study”. In: The Lancet (2020).
    
    

157.[157]. Wentao Zhu et al. “Beta-and Novel Delta-Coronaviruses Are Identified from Wild Animals in the Qinghai-Tibetan Plateau, China”. In: Virologica Sinica (2020), pp. 1–10.

 [1]: /embed/graphic-1.gif
 [2]: /embed/inline-graphic-1.gif