Visualizing the invisible: The effect of asymptomatic transmission on the outbreak dynamics of COVID-19 ======================================================================================================= * Mathias Peirlinck * Kevin Linka * Francisco Sahli Costabal * Jay Bhattacharya * Eran Bendavid * John P.A. Ioannidis * Ellen Kuhl ## Abstract Understanding the outbreak dynamics of the COVID-19 pandemic has important implications for successful containment and mitigation strategies. Recent studies suggest that the population prevalence of SARS-CoV-2 antibodies, a proxy for the number of asymptomatic cases, could be an order of magnitude larger than expected from the number of reported symptomatic cases. Knowing the precise prevalence and contagiousness of asymptomatic transmission is critical to estimate the overall dimension and pandemic potential of COVID-19. However, at this stage, the effect of the asymptomatic population, its size, and its outbreak dynamics remain largely unknown. Here we use reported symptomatic case data in conjunction with antibody seroprevalence studies, a mathematical epidemiology model, and a Bayesian framework to infer the epidemiological characteristics of COVID-19. Our model computes, in real time, the time-varying contact rate of the outbreak, and projects the temporal evolution and credible intervals of the effective reproduction number and the symptomatic, asymptomatic, and recovered populations. Our study quantifies the sensitivity of the outbreak dynamics of COVID-19 to three parameters: the effective reproduction number, the ratio between the symptomatic and asymptomatic populations, and the infectious periods of both groups. For nine distinct locations, our model estimates the fraction of the population that has been infected and recovered by Jun 15, 2020 to 24.15% (95% CI: 20.48%-28.14%) for Heinsberg (NRW, Germany), 2.40% (95% CI: 2.09%-2.76%) for Ada County (ID, USA), 46.19% (95% CI: 45.81%-46.60%) for New York City (NY, USA), 11.26% (95% CI: 7.21%-16.03%) for Santa Clara County (CA, USA), 3.09% (95% CI: 2.27%-4.03%) for Denmark, 12.35% (95% CI: 10.03%-15.18%) for Geneva Canton (Switzerland), 5.24% (95% CI: 4.84%-5.70%) for the Netherlands, 1.53% (95% CI: 0.76%-2.62%) for Rio Grande do Sul (Brazil), and 5.32% (95% CI: 4.77%-5.93%) for Belgium. Our method traces the initial outbreak date in Santa Clara County back to January 20, 2020 (95% CI: December 29, 2019 - February 13, 2020). Our results could significantly change our understanding and management of the COVID-19 pandemic: A large asymptomatic population will make isolation, containment, and tracing of individual cases challenging. Instead, managing community transmission through increasing population awareness, promoting physical distancing, and encouraging behavioral changes could become more relevant. Keywords * COVID-19 * epidemiology * asymptomatic transmission * machine learning * uncertainty quantification * Bayesian inference ## 1. Motivation Since its outbreak in December 2019, the COVID-19 pandemic has rapidly swept across the globe and is now affecting 188 countries with more than 24 million cases reported worldwide [10]. In the early stages of a pandemic, doctors, researchers, and political decision makers mainly focus on symptomatic individuals that come for testing and address those who require the most urgent medical attention [15]. In the more advanced stages, the interest shifts towards mildly symptomatic and asymptomatic individuals who–by definition–are difficult to trace and likely to retain normal social and travel patterns [40]. In this manuscript, we collectively use the term “asymptomatic” for undetected individuals who either have mild symptoms that are not directly associated with COVID-19 or display no symptoms at all. Recent antibody seroprevalence studies suggests that the number of asymptomatic COVID-19 cases outnumbers the symptomatic cases by an order of magnitude or more [3, 4, 9, 12, 14, 18, 24, 32, 34, 46, 58, 59, 64, 65, 68, 66, 67, 71, 72, 74, 76, 78, 81]. Estimating the prevalence and contagiousness of these asymptomatic cases is critical since it will change our understanding of the overall dimension and the pandemic potential of COVID-19 [16]. Yet, at this stage, the effect of the asymptomatic population, its size, and its outbreak dynamics remain largely unknown [**?**]. The first evidence of asymptomatic individuals in a family cluster of three was reported in late January, where one individual was mildly symptomatic and two remained asymptomatic, with normal lymphocyte counts and chest computer tomography images, but positive quantitative reverse transcription polymerase chain reaction tests [53]. As of mid June, more than 50 studies have reported an asymptomatic population, twenty-three of them with a sample size of at least 500 [30], with a median undercount of 20 across all studies, suggesting that only one in twenty COVID-19 cases is noticed and reported. These studies are based on polymerase chain reaction or antibody seroprevalence tests in different subgroups of the population, at different locations, at different points in time [3, 7, 71]. To no surprise, the reported undercount varies hugely, ranging from 3.5 and 5.0 in Luxembourg [68] and Germany [71] to 543 and 627 in Iran [65] and Japan [9] respectively. Most of these studies are currently only available on preprint servers, but an increasing number is now passing peer review, including a study of 1402 individuals in Wuhan City with an undercount of 22.1 [81], a study of 400 health care workers in London with an undercount of 35.0 [77], a community spreading study of 131 patients with influenza-like symptoms in Los Angeles with an undercount of 100.0 [70], and a seroprevalence study in Los Angeles county with an undercount of 43.5 [69]. The reported trend across all studies is strikingly consistent: A much larger number of individuals displays antibody prevalence than we would expect from the reported symptomatic case numbers. Knowing the exact dimension of the asymptomatic population is critical for two reasons: first, to truly estimate the severity of the outbreak, e.g., hospitalization or mortality rates [16], and second, to reliably predict the success of surveillance and control strategies, e.g., contact tracing or vaccination [19]. While there is a pressing need to better understand the prevalence of asymptomatic transmission, it is also becoming increasingly clear that it will likely take a long time until we can, with full confidence, deliver reliable measurements of this asymptomatic group [6]. In the meantime, mathematical modeling can provide valuable insight into the tentative outbreak dynamics and outbreak control of COVID-19 for varying asymptomatic scenarios [40]. Many classical epidemiology models base their predictions on compartment models in which individuals pass through different stages as they experience the disease [33]. A popular model to simulate the outbreak dynamics of COVID-19 is the SEIR model [15], which is made up of four compartments for the susceptible, exposed, infectious, and recovered populations [2]. Here, to explicitly account for the asymptomatic population, we introduce an SEIIR model, which further divides the infectious population into symptomatic and asymptomatic groups [35]. Similar models have recently been used to study the general role of asymptomatic carriers in disease transmission [54] and to illustrate how asymptomatic individuals have facilitated the rapid spread of COVID-19 throughout China [40], South Korea [75], and Italy [20]. While it is tempting-and easily possible-to introduce many more sub-populations into the model, for example a pre-symptomatic, hospitalized, or mortality group [57], here, we focus on the simplest possible model that allows us to explore the role of the asymptomatic population throughout the COVID-19 pandemic. To systematically probe different scenarios, we combine this deterministic SEIIR model with a dynamic effective reproduction number, machine learning, and uncertainty quantification to learn the reproduction number-in real time-and quantify and propagate uncertainties in the symptomatic-to-asymptomatic ratio and in the initial exposed and infectious populations. We show that this not only allows us to visualize the dynamics and uncertainties of the dynamic contact rate, the effective reproduction number, and the symptomatic, asymptomatic, and recovered populations, but also to estimate the initial date of the outbreak. ## 2. Methods ### 2.1. Epidemiology modeling #### SEIIR epidemiology model We model the epidemiology of COVID-19 using an SEIIR model with five compartments, the susceptible, exposed, symptomatic infectious, asymptomatic infectious, and recovered populations. Figure 1 illustrates our SEIIR model, which is governed by a set of five ordinary differential equations, ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/29/2020.05.23.20111419/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/F1) Figure 1: SEIIR epidemiology model. The SEIIR model contains five compartments for the susceptible, exposed, symptomatic infectious, asymptomatic infectious, and recovered populations. The transition rates between the compartments, *β*, *α*, and *γ* are inverses of the contact period *B =* 1/*β*, the latent period *A* = 1*/α*, and the infectious period *C =* 1/*γ*. The symptomatic and asymptomatic groups have the same latent period *A*, but they can have individual contact periods *B*s = 1/β*s* and *B*a = 1/*β*a and individual infectious periods *C*s = 1/*γ*s and *C*a = 1/*γ*a. The fractions of the symptomatic and asymptomatic subgroups of the infectious population are *v*s and *v*a. We assume that the infection either goes through the symptomatic or the asymptomatic path, but not both for one individual. ![Formula][1] where the fractions of all five populations add up to one, *S + E + Is + Ia + R =* 1. We assume that both the symptomatic group *I*s and the asymptomatic group *I*a can generate new infections. We introduce these two groups by fractions *v*s and *v*a of the subjects transferring from the exposed group *E*. We postulate that the two infectious groups *I*s and *I*a have the same latent period *A* = 1/*α*, but can have individual contact periods *B*s = 1/*β*s and *B*a = 1/*β*a to mimic their different community spreading, and individual infectious periods *C*s = 1/*γ*s and *C*a = 1/*γ*a to mimic their different likelihood of isolation. From the infectious fractions, we can derive the overall contact and infectious rates *β* and *γ* from their individual symptomatic and asymptomatic counterparts, *β*s, *β*a, *γ*s, and *γ*a, ![Formula][2] Similarly, we can express the overall contact and infectious periods *B* and *C* in terms of their symptomatic and asymptomatic counterparts, *B*s, *B*a, *C*s, and *C*a, ![Formula][3] Naturally, the different dynamics for the symptomatic and asymptomatic groups also affect the basic reproduction number *R*, the number of new infections caused by a single one individual in an otherwise uninfected, susceptible population, ![Formula][4] For a large asymptomatic group *v*a → 1, the basic reproduction number approaches the ratio between the infectious and contact periods of the asymptomatic population, *R* → *C*a/*B*a, which could be significantly larger than the basic reproduction number for the symptomatic group, *R* = *C*s/*B*s, that we generally see reported in the literature. To characterize the effect of changes in social behavior and other interventions that may affect contact, we assume that the contact rate *β*(*t*) can vary as a function of time [43], but is the same for the symptomatic and asymptomatic groups, ![Formula][5] This introduces a time-varying effective reproduction number *R*(*t*), which is an important real time characteristic of the current outbreak dynamics. #### SEIR epidemiology model For the special case when the dynamics of the symptomatic and asymptomatic groups are similar, i.e., *β*s = *β*a = *β* and *γ*s = *γ*a = *γ*, we can translate the SEIIR model (1) into the classical SEIR model (6) with four compartments, the susceptible, exposed, infectious, and recovered populations [42, 44, 82, 83]. For this special case, we can back-calculate the symptomatic and asymptomatic groups from equation (6.3) as *I*s = *v*s *I* and *I*a = *v*a *I*. Figure 2 illustrates the SEIR model, which is governed by a set of four ordinary differential equations [28], ![Formula][6] ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/29/2020.05.23.20111419/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/F2) Figure 2: SEIR epidemiology model. The SEIR model contains four compartments for the susceptible, exposed, infectious, and recovered populations. The transition rates between the compartments, *β*, *α*, and *γ* are inverses of the contact period *B* = 1/*β*, the latent period *A* = 1/*α*, and the infectious period *C* = 1 /*γ*. If the transition rates are similar for the symptomatic and asymptomatic groups, the SEIIR model simplifies to the SEIR model with *I*s = *v*s *I* and *I*a = *v*a*I*. #### Discrete SEIIR epidemiology model We focus on the SEIIR model with five ordinary differential equations (1), which we discretize in time using a finite difference approximation ![Graphic][7] of the evolution of all five populations (o) = *S*, *E*, *I*s, *I*a, *R*. Here Δ*t* = *t*n+1 − *t*n denotes the discrete time increment, usually Δ*t* = 1 day, and (o)n+1 and (o)n denote the populations of the new and previous time steps. We apply an explicit time integration scheme to obtain the following discrete system of equations, ![Formula][8] At each location, we begin our simulation on day *t*, the day of the beginning of the local lockdown, with initial conditions, *S*, *E*, *I*s0, *I*a0, *R*. We further assume that *S* = 1 − *E* − *I*s0 − *I*a0, *E* = [*A/C*s] [*I*s0 + *I*a0] [13], and *V*s *I*a0 = [1 − *v*s] *I*s0 to express all five initial conditions in terms of the initial symptomatic infectious population *I*s0, such that *S* = 1 − [1 + *A/C*s]/*v*s *I*s0, *E* = [*A/C*s]/*v*s *I*s0, *I*a0 = [1 − *v*s]/*v*s *I*s0, and *R* = 0. From the solution of the time discrete SEIIR system (7), for each point in time, we calculate the of detected population, ![Formula][9] as the discrete sum of the symptomatic infectious and recovered populations *I*s,n+1 and *R*s,n+1 with ![Graphic][10], or, equivalently, as the explicit update of the previous detected population *D*n with the new symptomatic infectious population, *v*s *α E*n, that transitions from the exposed to the symptomatic infectious group. We compare the daily simulated detected population *D*(*t*) against the daily reported detected population ![Graphic][11], which we obtain by scaling the daily number of detected cases by the total population *N* for each location [5, 11, 27, 48, 49, 73]. ### 2.2. Seroprevalence studies We consider seroprevalence studies up to mid June 2020 where the representation ratio, the ratio between the amount of seroprevalence samples and the representative population, is larger than 0.02%. Subsequently, we draw the daily number of confirmed reported cases for Heinsberg, Ada County, New York City, Santa Clara County, Denmark, Geneva Canton, Netherlands, Rio Grande do Sul, and Belgium from online data repositories [5, 11, 27, 48, 49, 73]. The symptomatic fraction *v*s follows from the ratio between the reported confirmed cases for the studied region and the amount of seroprevalence-estimated cases on the last day of the seroprevalence study. Table 1 summarizes the study date, number of samples, representation ratio, population, and symptomatic fraction of each location. View this table: [Table 1:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/T1) Table 1: Seroprevalence studies with a representative population larger than 0.02%, location, study date, number of samples, representation ratio, population, and symptomatic fraction vs. ### 2.3. Bayesian learning To calibrate our SEIIR model (1), we compare the model output, the simulated detected population *D*(*t*) for a given parameter set ![Graphic][12], against the data, the reported detected population ![Graphic][13] for each location [5, 11, 27, 48, 49, 73]. Our time-discrete SEIIR model (7) uses the following set of parameters, ![Formula][14] the latency period *A*, the symptomatic and asymptomatic infectious periods *C*s and *C*a, the symptomatic fraction *v*s, the dynamic contact rate *β*(*t*), the initial exposed and infectious populations *E*, *I*s0, and *I*a0, and the likelihood width o. To reduce the set of unknowns, we fix the latency period to *A* = 2.5 days and the symptomatic infectious period to *C*s = 6.5 days [26, 37, 39, 62]. Since the asymptomatic infectious period *C*a is unreported, we study three cases with *C*a = 0.5,1.0,2.0 *C*s resulting in infectious periods of *C*a = 3.25, 6.5, and 13.0 days. We further express the initial exposed and asymptomatic infectious populations in terms of the symptomatic infectious population as *E* = [*A*/*C*s]/*v*s *I*s0 [13] and *I*a0 = [1 − *v*s]/*v*s *I*s0. This results in the following reduced set of parameters, ![Formula][15] including the symptomatic fraction *v*s, the dynamic contact rate *β*(*t*), the initial symptomatic infectious population *I*s0, and the likelihood witdth *σ*. #### Bayesian inference We use Bayesian inference to estimate the set of parameters ![Graphic][16] so that the statistics of the output of the model *D*(*t*) agree with that of the data ![Graphic][17] [41, 51]. We relate the prior, likelihood, and posterior using Bayes’ theorem [52], ![Formula][18] where ![Graphic][19] is the prior, i.e., the probability distribution of the model parameters, ![Graphic][20] is the likelihood, i.e., the conditional probability of the data when the parameter set is fixed, ![Graphic][21] is the posterior, i.e., the conditional probability of the parameters for given data ![Graphic][22], and ![Graphic][23] is the evidence. #### Priors For the prior probability distributions ![Graphic][24], we assume weakly informative priors for our model parameters, the symptomatic fraction *v*s, the dynamic contact rate *β*(*t*), the initial symptomatic infectious population *I*s0, and the likelihood witdth *σ*. For the symptomatic fraction *v*s, we adopt the normal distributions from the individual local antibody seroprevalence studies of all nine locations [3,4, 14, 58, 64, 66, 67, 71, 72]. Specifically, we set ![Graphic][25], where *μ*s and *σ*s are the mean and standard deviation of the symptomatic fraction based on the ratio of the reported and seroprevalence-estimated numbers of confirmed cases for each location in Table 1. For the dynamic contact rate *β*(*t*), we postulate that its logarithm follows a Gaussian random walk, ![Graphic][26]. Specifically, this Gaussian random walk assumes a sequence of inter-dependent Gaussian distributions where the contact rate *β*(*t*n+1) on day *t*n+1 depends the contact rate *β*(*t*n) of the previous day *t*n with the initial condition, ![Graphic][27], and the daily update, ![Graphic][28], where *μβ* is the mean initial contact rate, *σβ* is the daily step width, and *μβ* is the overall drift between the initial value log(*β*(*t*)) and the final value log(*β*(*t*n*+*1)). For the initial symptomatic infectious population *I*s0, we postulate a log-normal distribution, ![Graphic][29], where we select a weakly informative prior for the mean ![Graphic][30], as the reported detected population on day *t* and a standard deviation of *σ*s *=* 1. Table 2 summarizes our prior distributions and the SEIIR model parameters. View this table: [Table 2:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/T2) Table 2: Prior distributions and SEIIR model parameters. #### Likelihood For the likelihood ![Graphic][31], we introduce a likelihood function ![Graphic][32] that evaluates the proximity between the model output *D*(*t*), the simulated detected population for a given parameter set ![Graphic][33], and the data ![Graphic][34], the relative reported detected population in probabilistic terms, ![Formula][35] We choose a Student’s t-distribution because it resembles a Gaussian distribution around the mean with heavy tails, which make the method more robust with respect to outliers and reporting noise [8, 36]. Here *σ* represents the likelihood width for which we postulate a halfCauchy distribution, see Table 2. The likelihood, ![Formula][36] is simply the product of all likelihood functions ![Graphic][37] evaluated at each time *ti* for *i* = 0,1,…, n. #### Posteriors With the prior ![Graphic][38] and the likelihood ![Graphic][39], we estimate the posterior ![Graphic][40] using Bayes’ theorem (11) [31]. Since we cannot describe the posterior distribution over the model parameters ![Graphic][41] analytically, we adopt approximate-inference techniques to calibrate our model on the available data. We use the NO-U-Turn sampler (NUTS) [29], which is a type of Hamiltonian Monte Carlo algorithm as implemented in PyMC3 [61]. We use four chains. The first 4 times 500 samples serve to tune the sampler and are later discarded. The subsequent 4 times 1000 samples define the posterior distribution of the parameters ![Graphic][42]. From the converged posterior distribution ![Graphic][43], we sample multiple combinations of parameters ![Graphic][44] that describe the time evolution of detected population ![Graphic][45]. Using these posterior samples, we quantify the uncertainty of each parameter. As such, each parameter set ![Graphic][46] provides a set of values for the symptomatic fraction *v*s, the initial symptomatic populations *I*s0, and the dynamic contact rate *β*(*t*). From these values, we quantify the effective reproductive number *R*(*t*) using equation (5) and the time evolution of the susceptible, exposed, symptomatic infectious, asymptomatic infectious and recovered populations, *S*(*t*), *E*(*t*), *I*s(*t*), *I*a(*t*), and *R*(*t*) using equations (7) and report their values with the associated 95% credible interval [50]. #### Hierarchical modeling To estimate the asymptomatic infectious period *C*a, we create a hierarchical model and analyze the confirmed case data from all nine locations in Table 1. To prevent overfitting, we postulate that the initial reproduction number *R* is the same for all nine locations and remains constant during a two-week time window prior to lockdown. This implies that the contact rate *β* is static and constant. We assume a normal distribution for the basic reproduction number ![Graphic][47] with mean *μR* = 3.87 and standard deviation *σR* *=* 0.44 [17]. This defines the initial contact rate *β*0i at each location, i = 1,…, 9, as *β*0i = *R* (*v*si/*C*s + (1 − *v*si)/*C*ai). We create a hyper-distribution for the overall asymptomatic infectious period, ![Graphic][48], with mean *μC*a = 6.5 and standard deviation ![Graphic][49], and assume that each local asymptomatic infectious period ![Graphic][50] is drawn from this hyperdistribution. All other parameters are similar to the previous section. Table 3 summarizes the prior distributions and the SEIIR model parameters of our hierarchical model to infer the asymptomatic infectious period *C*a. View this table: [Table 3:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/T3) Table 3: Prior distributions and SEIIR model parameters to hierarchically infer the asymptomatic infectious period *C*a. To estimate the hierarchical asymptomatic infectious period *C*a, we use Bayesian inference. Following equation (11), we relate the prior, likelihood, and posterior, with the only difference that the simulated and reported detected populations *D*i(*t*) and ![Graphic][51] are now vectors that contain the detected populations of all *i =* 1,…, 9 locations. The likelihood ![Graphic][52] becomes the product of the all likelihood functions ![Graphic][53] evaluated for a period of n = 14 days prior to lockdown, ![Formula][54] For each location, we estimate the posterior ![Graphic][55] using the NO-U-Turn sampler (NUTS) [29] in PyMC3 [61]. From the converged posterior distribution, we sample multiple ![Graphic][56] combinations to quantify the evolution of the reproduction number *R* and contact rate *β* at each location, and the associated susceptible, exposed, symptomatic and asymptomatic infectious, and recovered populations with their associated 95% credible intervals. #### Outbreak date For each sample from the posterior distribution, we use the infered initial exposed and asymptomatic infectious populations *E* and *I*a0 to estimate the date of the very first COVID-19 case in Santa Clara County [63]. Specifically, for each parameter set, we create an SEIIR model and assume that the outbreak begins with one single asymptomatic infectious individual. We fix the latency and symptomatic infectious periods to *A =* 2.5 days and *C*s = 6.5 days and impose the hierarchically estimated infectious period *C*a from the previous section. For each posterior sample for the exposed, symptomatic infectious, and asymptomatic infectious population size on the lockdown date, March 16, 2020, we use the Nelder-Mead optimization method [21] to find the most probable outbreak origin date. Specifically, we solve the SEIIR model forward in time using an explicit time integration, starting from various start dates before March 16, 2020, and iteratively minimize the difference between the computed exposed, symptomatic, and asymptomatic infectious populations and the sample’s actual exposed, symptomatic, and asymptomatic infectious populations. We concomitantly fit a static contact rate *β*, which is bounded between zero and the posterior sample’s estimated contact rate *β* at the beginning of the simulation *t*. Repeating this process for each sample of the Bayesian inference generates a distribution of possible origin dates. From this distribution, we compute the most probable origin date and its uncertainty. ## 3. Results ### 3.1. Outbreak dynamics of COVID-19 in Santa Clara County Figure 3 illustrates the outbreak dynamics of COVID-19 in Santa Clara County. To explore the effects of the asymptomatic infectious period *C*a on the initial effective reproduction number, we start the simulations on March 6, 2020, 10 days before the local lockdown date. The black dots highlight the reported detected population ![Graphic][57] from this day forward. Based on these data points, we learn the posterior distributions of our SEIIR model parameters for fixed latent and symptomatic infectious periods *A* = 2.5 days and *C*s = 6.5 days, and for three asymptomatic infectious periods, *C*a = 3.25, 6.5, and 13.0 days, from left to right. The gray and green-purple regions highlights the 95% credible intervals on the confirmed cases *D*(*t*), top row, and the contact rate *β*(*t*), bottom row, based on the reported cases D(t), while taking into account uncertainties on the symptomatic infectious fraction *v*s, and the initial exposed and infectious populations *E*, *I*s0, and *I*a0. The red, orange, and green histograms on the middle row display the inferred initial exposed and infectious populations, *E*, *I*s0, and *I*a0, for the three different asymptomatic infectious periods. The graphs confirm that our dynamic SEIIR epidemiology model is capable of correctly capturing the gradual flattening and deflattening of the curve of confirmed cases in agreement with the decrease in new cases reported in Santa Clara County, top row. The consistent downward trend of the contact rate *β*(*t*) after the lockdown date (March 16, 2020) quantifies the efficiency of public health interventions. The different magnitudes in the contact rate highlight the effect of the three different asymptomatic infectious periods *C*a: For larger asymptomatic infectious periods *C*a, from left to right, to explain the same number of confirmed cases *D*(*t*) *= I*s(*t*) + *R*s(*t*), the contact rate *β*(*t*) has to decrease. On March 6, 2020, the mean contact rate *β*(*t*) was 0.688 (95% CI: 0.612 - 0.764) for an infectious period of *C*a=3.25 days, 0.529 (95% CI: 0.464 - 0.593) for *C*a= 6.5 days, and 0.463 (95% CI: 0.403 - 0.523) for *C*a =13.0 days. By March 16, 2020, the day Santa Clara County announced the first county-wide shelter-in-place order in the entire United States, these contact rates *β*(*t*) were 0.491 (95% CI: 0.462 - 0.523) for an infectious period of *C*a=3.25 days, 0.328 (95% CI: 0.305 - 0.352) for *C*a= 6.5 days, and 0.252 (95% CI: 0.234 - 0.271) for *C*a=13.0 days. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/29/2020.05.23.20111419/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/F3) Figure 3: Outbreak dynamics of COVID-19 in Santa Clara County. The simulation learns the time-varying contact rate *β*(*t*) for fixed latent and symptomatic infectious periods *A* = 2.5 days and *C*s = 6.5 days, and for three asymptomatic infectious periods *C*a = 3.25 days, 6.5 days, and 13.0 days (from left to right). Computed and reported confirmed cases in Santa Clara County, *D*(*t*) = *I*s(*t*) + *R*s(*t*) and ![Graphic][58] (top), initial exposed and infectious populations, *E*, *I*s0, and *I*a0 (middle), and dynamic contact rate, *β*(*t*) (bottom). The gray and green-blue regions highlight the 95% credible intervals on the confirmed cases *D*(*t*) (top) and the contact rate *β*(*t*) (bottom) based on the reported cases ![Graphic][59], while taking into account uncertainties on the fraction of the symptomatic infectious population *v*s = *I*s /*I*, and the initial exposed and infectious populations *E*, *I*s0, and *I*a0. ### 3.2. Effect of asymptomatic transmission ofCOVID-19 in Santa Clara County Figure 4 visualizes the effect of asymptomatic transmission in Santa Clara County. The simulation learns the dynamic contact rate *R*(*t*), and with it the time-varying effective reproduction number *R*(*t*), top row, for three asymptomatic infectious periods *C*a = 3.25 days, 6.5 days, and 13.0 days, from left to right. The effective reproduction number *R*(*t*) follows a similar downward trend as the contact rate *β*(*t*). For larger asymptomatic infectious periods *C*a, from left to right, since *R*(t) *= C*s*β*(*t*)*/*[*v*s + *v*a *C*s/*C*a], as *C*a increases, *C*s/*C*a decreases, and *R*(*t*) increases. Since *R*(*t*) represents the number of new infections from a single case, a decrease below *R*(*t*) < 1 implies that a single infectious individual infects less than one new individual, which indicates that the outbreak decays. The dashed vertical lines indicate the date *R*(*t*) = 1 during which one infectious individual, either symptomatic or asymptomatic, infects on average one other individual. For an asymptomatic infectious period of *C*a=3.25 days, it took until March 28 before Santa Clara County managed to get *R*(*t*) below 1 for the first time after the outbreak. For *C*a=6.5 days, this only occurred by April 1 and for *C*a=13.0 days, this occurred on April 8, 2020. This confirms our intuition that, the larger the asymptomatic infectious period *C*a, for example because asymptomatic individuals will not isolate as strictly as symptomatic individuals, the higher the effective reproduction number *R*(*t*), and the more difficult it will be to control *R*(*t*) by public health interventions. For each of the three cases, the symptomatic infectious, asymptomatic infectious, and recovered population, are shown in the bottom row. For larger asymptomatic infectious periods *C*a, from left to right, the total infectious population *I* increases and its maximum occurs later in time. Specifically, the maximum infectious population since March 6, 2020 amounts to 0.70% (95% CI: 0.43%-0.97%) on March 28, 2020 for *C*a = 3.25 days, 1.23% (95% CI: 0.72%-1.75%) on April 2, 2020 for *C*a = 6.5 days, and 2.10% (95% CI: 1.25%-2.94%) on April 7, 2020 for *C*a = 13.0 days. For larger asymptomatic infectious periods *C*a, from left to right, the recovered population *R* decreases. Specifically, on June 15, 2020, the recovered population *R* amounts to 10.85% (95% CI: 7.06%-16.09%) for an infectious period of *C*a =3.25 days, 10.20% (95% CI: 6.07%-14.90%) for *C*a =6.5 days, and 9.90% (95% CI: 6.07%-13.93%) for *C*a =13.0 days. Similarly, and important when considering different exit strategies, the total infectious population, *I* = *I*s + *I*a, on June 15, 2020 is estimated to 0.39% (95% CI: 0.24%-0.58%) for *C*a =3.25 days, 0.68% (95% CI: 0.40%-1.01%) for *C*a =6.5 days, and 1.25% (95% CI: 0.73%-1.77%) for *C*a =13.0 days. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/29/2020.05.23.20111419/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/F4) Figure 4: Effect of asymptomatic transmission of COVID-19 in Santa Clara County. The simulation learns the time-varying contact rate *β*(*t*), and with it the time-varying effective reproduction number *R*(*t*), for fixed latent and symptomatic infectious periods *A =* 2.5 days and *C*s = 6.5 days, and for three asymptomatic infectious periods *C*a = 3.25 days, 6.5 days, and 13.0 days (from left to right). The downward trend of the effective reproduction number *R*(*t*) reflects the efficiency of public health interventions (top row). The dashed vertical lines mark the critical time period during which the effective reproductive reproduction number fluctuates around *R*(*t*) *=* 1. The simulation predicts the symptomatic infectious, asymptomatic infectious, and recovered populations *I*s, *I*a, and *R* (bottom row). The colored regions highlight the 95% credible interval for uncertainties in the number of confirmed cases D, the fraction of the symptomatic infectious population *v*s = *I*s*/I*, the initial exposed population *E* and the initial infectious populations *I*s0 and *I*a0. ### 3.3. Outbreak dynamics of COVID-19 worldwide Figure 5 illustrates the hierarchical asymptomatic infectious period *C*a estimation for nine locations that reported COVID-19 antibody prevalence in a representative sample of the population: Heinsberg (NRW, Germany), Ada County (ID, USA), New York City (NY, USA), Santa Clara County (CA, USA), Denmark, Geneva Canton (Switzerland), Netherlands, Rio Grande do Sul (Brazil) and Belgium. This estimate assumes an initial fixed reproduction number *R* = 3.87 (95%CI: 3.01-4.66) during a two-week window before lockdown which results in an asymptomatic infectious period of *C*a = 5.76 (95%CI: 3.59-8.09) days. With this value, we simulate the outbreak dynamics of COVID-19 in all nine locations. Figure 6 illustrates the learnt effective reproduction number *R*(*t*), and the symptomatic and asymptomatic infectious populations *I*s and *I*a, and the recovered population *R* for all nine locations. Here, we assume fixed latent and infectious periods of *A* = 2.5 days, *C*s = 6.5 days, and the hierarchical asymptomatic infectious period *C*a = 5.76 (95%CI: 3.59-8.09) days from Figure 5. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/29/2020.05.23.20111419/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/F5) Figure 5: Hierarchical asymptomatic infectious period *C*a estimation. Assuming an initial fixed reproduction number *R* = 3.87 (95%CI: 3.01-4.66) [17], the simulation generates histograms of the asymptomatic infectious period *C*a for each location based on the location-specific symptomatic fraction *v*s. The black dots and grey regions represent the reported and simulated detected cases ![Graphic][60] and *D*(*t*) respectively. The hierarchical hyperdistribution for the asymptomatic infectious period results in *C*a = 5.76 (95%CI: 3.59-8.09) days, right histogram. For all nine locations, the calculated metrics display the fine balance between the dynamics of the effective reproduction number and the control of the epidemic outbreak. The downward trend of the effective reproductive number *R*(*t*) quantifies how fast each location managed to control the spreading of COVID-19. The dashed vertical line indicates the first time each location managed to lower the effective reproduction below *R*(*t*) = 1 after lockdown. This critical transition occurred on March 20 for Heinsberg, April 4 for Ada County, April 28 for New York City, March 30 for Santa Clara County, April 8 for Denmark, April 2 for Geneva, April 14 for the Netherlands, June 8 for Rio Grande do Sul and April 11 for Belgium. Based on our simulations, the maximum infectious population size amounted to 3.54% (95% CI: 2.96%-4.11%) on March 19 in Heinsberg, 0.47% (95% CI: 0.41%-0.55%) on May 5 in Ada County, 6.11% (95% CI: 5.99%-6.21%) on April 13 in New York City, 1.11% (95% CI: 0.69%-1.65%) on March 29 in Santa Clara County, 0.38% (95% CI: 0.27%-0.50%) on April 9 in Denmark, 2.10% (95% CI: 1.65%-2.69%) on April 2 in Geneva, 0.62% (95% CI: 0.56%-0.68%) on April 15 in the Netherlands, 0.28% (95% CI: 0.14%-0.48%) on June 8 in Rio Grande do Sul, and 0.70% (95% CI: 0.62%-0.78%) on April 11 in Belgium. On Jun 15, 2020, the estimated recovered population reached 24.15% (95% CI: 20.48%-28.14%) in Heinsberg, NRW, Germany 2.40% (95% CI: 2.09%-2.76%) in Ada County, ID, USA 46.19% (95% CI: 45.81%-46.60%) in New York City, NY, USA 11.26% (95% CI: 7.21%-16.03%) in Santa Clara County, CA, USA 3.09% (95% CI: 2.27%-4.03%) in Denmark 12.35% (95% CI: 10.03%-15.18%) in Geneva Canton, Switzerland 5.24% (95% CI: 4.84%-5.70%) in Netherlands 1.53% (95% CI: 0.76%-2.62%) in Rio Grande do Sul, Brazil 5.32% (95% CI: 4.77%-5.93%) in Belgium. ### 3.4. Estimating the outbreak date Figure 7 shows the estimated outbreak date of COVID-19 in Santa Clara County. For fixed latent and symptomatic infectious periods *A* = 2.5 days and *C*s = 6.5 days, and for the hierarchically estimated asymptomatic infectious periods *C*a = 5.76 (95%CI: 3.59-8.09) days, the graph highlights the estimated date of the first COVID-19 case in the county. Based on the reported case data from March 16, 2020 onward, and taking into account uncertainty on the fraction of the symptomatic infectious population *v*s, on the initial exposed population *E*, and on the initial symptomatic and asymptomatic infectious populations *I*s0 and *I*a0, we systematically backtracked the date of the first undetected infectious individual. Our results suggest that the first case of COVID-19 in Santa Clara County dates back to January 20, 2020 (95% CI: December 29, 2019 - February 13, 2020). ## 4. Discussion A key question in understanding the outbreak dynamics of COVID-19 is the dimension of the asymptomatic population and its role in disease transmission. Throughout the past months, dozens of studies have been initiated to quantify the fraction of the general population that displays antibody prevalence but did not report symptoms of COVID-19. Here we assume that this subgroup of the population has been infected with the novel coronavirus, but has remained asymptomatic, or only displayed mild symptoms that were not directly reported in the context of COVID-19. We collectively map this subgroup into an asymptomatic population and additively decompose the total infectious population, *I* = *I*s + *I*a, into a symptomatic group *I*s and an asymptomatic group *I*a. We parameterize this decomposition in terms of a single scalar valued parameter, the symptomatic fraction *v*s. Within this paradigm, we can conceptually distinguish two scenarios: the special case for which both subgroups display identical contact rates *β* latent periods *A*, and infectious periods *C*, and the general case for which these transition dynamics are different. ### For comparable dynamics, the size of the asymptomatic population does not affect overall outbreak dynamics For the special case in which both subgroups display identical contact rates *β*, latent rates *α*, and infectious rates *γ* [79], our study shows that the overall outbreak dynamics can be represented by the classical SEIR model [28] using equations (6). Importantly, however, since the reported case data only reflect the symptomatic infectious and recovered groups *I*s and *R*s, the true infectious and recovered populations *I* = *I*s/*v*s and *R = R*s/*v*s could be about an order of magnitude larger than the SEIR model predictions. From an individual’s perspective, a smaller symptomatic group *v*s, or equivalently, a larger asymptomatic group *v*a = [1 − *v*s], could have a personal effect on the likelihood of being unknowingly exposed to the virus, especially for high-risk populations: A larger asymptomatic fraction *v*a would translate into an increased risk of community transmisson and would complicate outbreak control [16]. From a health care perspective, however, the special case with comparable transition dynamics would not pose a threat to the health care system since the overall outbreak dynamics would remain unchanged, independent of the fraction *v*a of the asymptomatic population: A larger asymptomatic fraction would simply imply that a larger fraction of the population has already been exposed to the virus-without experiencing significant symptoms-and that the true hospitalization and mortality rates would be much lower than the reported rates [30]. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/29/2020.05.23.20111419/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/F6) Figure 6: Outbreak dynamics of COVID-19 worldwide. Dynamic effective reproduction number R(t) and symptomatic, asymptomatic, and recovered populations at all nine locations. The simulation learns the time-varying contact rate *β*(*t*), and with it the time-varying effective reproduction number *R*(*t*), to predict the symptomatic infectious, asymptomatic infectious, and recovered populations *I*s, *I*a, and *R*, for fixed latent and infectious periods *A* = 2.5 days, *C*s = 6.5 days, the hierarchical asymptomatic infectious period *C*a = 5.76 (95%CI: 3.59-8.09) days from Figure 5. The dashed vertical lines mark the the first time each location managed to lower the effective reproduction below *R*(*t*) = 1 after lockdown. The colored regions highlight the 95% credible interval for the effective reproductive number *R*(*t*) (top), the symptomatic and asymptomatic populations *I*s and *I*a, and the recovered population *R* (bottom plots), for uncertainties in the number of confirmed cases *D*, the fraction of the symptomatic infectious population *v*s, the initial exposed population *E*, and the initial infectious populations *I*s0 and *I*a0. ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/29/2020.05.23.20111419/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/F7) Figure 7: Estimating the outbreak date of COVID-19 in Santa Clara County varying asymptomatic infectious periods *C*a. Estimated date of the first COVID-19 case in Santa Clara County for fixed latent and symptomatic infectious periods *A =* 2.5 days and *C*s = 6.5 days, and for the hierarchical asymptomatic infectious period *C*a = 5.76 (95%CI: 3.59-8.09) days from Figure 5. The colored regions in the main plot highlight the 95% credible interval for the time evolution of the exposed and asymptomatic infectious populations *E* and *I*a estimated based on the reported cases ![Graphic][61] from March 16, 2020 onward and taking into account uncertainties on the fraction of the symptomatic infectious population *v*s = *I*s/*I*, and the exposed and asymptomatic infectious populations *E* and *I*a0 on March 16, 2020 (right plot). The bottom plot histogram shows the distribution of the most probable origin dates to January 20, 2020 (95% CI: December 29, 2019 - February 13, 2020). ### For different dynamics, the overall outbreak dynamics depend on both size and infectiousness of the asymptomatic group For the general case in which the transition rates for the symptomatic and asymptomatic groups are different, the overall outbreak dynamics of COVID-19 become more unpredictable, since little is known about the dynamics of the asymptomatic population [54]. To study the effects of different dynamics between the symptomatic and asymptomatic groups, we decided to collectively represent a lower infectivity of the asymptomatic population through a smaller infectious period *C*a < *C*s and a lack of early isolation of the asymptomatic population through a larger infectious period *C*a > *C*s, while, for simplicity, keeping the latent period A and contact rate *β* similar across both groups [39]. Our study shows that the overall reproduction number, *R*(*t*) *=* [*C*a *C*s]/[*v*s *C*a + *v*a *C*s] *β*(*t*), and with it the outbreak dynamics, depend critically on the fractions of the symptomatic and asymptomatic populations *v*s and *v*a and on the ratio of the two infectious periods *C*s and *C*a. To illustrate these effects, we report the results for three different scenarios where the asymptomatic group is half as infectious, *C*a = 0.5 *C*s, equally infectious, *C*a = *C*s, and twice as infectious *C*a = 2.0 *C*s as the symptomatic group for Santa Clara County. The second case, the middle column in Figures 3 and 4, corresponds to the special case with comparable dynamics and similar parameters. Our learnt asymptomatic infectious periods of *C*a = 5.76 days in Figure 5 suggest that *C*a is smaller than the symptomatic infectious period of *C*s = 6.5 days and that the asymptomatic population is slightly less infectious as the symptomatic population. This can be a combined effect of being less contagious but, at the same time, having more contacts, since asymptomatic individuals, by definition, do not realize that they are spreading the disease [38]. ### Dynamic contact rates are a metric for the efficiency of public health interventions Classical SEIR epidemiology models with static parameters are well suited to model outbreak dynamics under unconstrained conditions and predict how the susceptible, exposed, infectious, and recovered populations converge freely toward the endemic equilibrium [28]. However, they cannot capture changes in disease dynamics and fail to converge towards a temporary equilibrium before the entire population has become sufficiently immune to prevent further spreading [57]. To address this limitation, we introduce a time-dependent contact rate *β*(*t*), which we learn dynamically from the reported case data. Figures 3 and 6 demonstrate that our approach can successfully identify a dynamic rate that not only decreases monotonically, but is also capable of reproducing local contact fluctuations [43]. With this dynamic contact rate, our model can capture the characteristic S-shaped COVID-19 case curve that plateaus before a large fraction of the population has been affected by the disease, resembling a Gompertz function. It also naturally captures a potential regrowth of the contact rate *β*(*t*), and with it the effects of a second wave [45]. Previous studies have inferred discrete date points at which the contact rates vary [8] or used sliding windows over the amount of novel reported infections [55] to motivate dynamic contact rates. As such, our framework provides a model-based method for statistical inference of virus transmissibility: It naturally learns the most probable contact rate from the changing time evolution of new confirmed cases and concomitantly quantifies the uncertainty on that estimation. ### The dynamics of the asymptomatic population affect the effective reproduction number Our analysis in equations (4) and (5) and our simulations in Figure 4 illustrate how asymptomatic transmission affects the effective reproduction number, and with it the outbreak dynamics of COVID-19. Our results show that, the larger the infectious period *C*a of the asymptomatic group, the larger the initial effective reproduction number, *R*(*t*) = *C*s*β*(*t*)/[*v*s + *v*a *C*s/*C*a], and the later the drop of *R*(*t*) below the critical value of one. A recent study analyzed the dynamics of the asymptomatic population in three consecutive windows of two weeks during the early outbreak in China [40]. The study found relatively constant latency and infectious periods *A* and *C*, similar to our assumption, and a decrease in the contact rate *β* = 1.12,0.52,0.35 days−1 and in the effective reproduction number *R*(*t*) = 2.38,1.34,0.98, which is consistent with our results. However, rather than assuming constant outbreak parameters within pre-defined time windows, our study learns the effective reproduction number dynamically, in real time, from the available data. Figures 4 and 6 demonstrate that we can successfully learn the critical time window until *R*(*t*) drops below 1, which, in Santa Clara County, took till March 28 for *C*a =3.25 days, till April 1 for *C*a = 6.5 days and till April 8 for *C*a = 13.0 days. Our findings are consistent with the observation that the basic reproduction number will be over-estimated if the asymptomatic group has a shorter generation interval, and underestimated if it has a longer generation interval than the symptomatic group [54]. Naturally, these differences are less pronounced under current conditions where the effective reproduction number is low and the entire population has been sheltering in place. It will be interesting to see if the effects of asymptomatic transmission become more visible as we gradually relax the current constraints and allow all individuals to move around and interact with others more freely. Seasonality, effects of different temperature and humidity, and other unknown factors may also influence the extent of transmission. ### Estimates of the infectious asymptomatic population may vary, but general trends are similar Throughout the past months, an increasing number of researchers around the globe have started to characterize the size of the asymptomatic population to better understand the outbreak dynamics of COVID-19 [30]. The reported undercount varies largely with values on the order of one estimated by the World Health Organization [80] and values on the order of ten reported by the Center for Disease Control and Prevention [6]. Two major challenges drive the interest in these studies: estimating the severity of the outbreak, e.g., hospitalization and mortality rates [16], and predicting the success of surveillance and control efforts, e.g., contact tracing or vaccination [19]. This is especially challenging now-in almost complete lockdown-when the differences in transmission dynamics between the symptomatic and asymptomatic populations are small and difficult to quantify. However, as Figure 4 suggests, these transmission dynamics can have a significant effect on the size of the asymptomatic population: For infectious periods of *C*a = 0.5, 1.0, and 2.0 *C*s, the maximum infectious population in Santa Clara County varies from 0.70% and 1.23% to 2.10%. Interestingly, not only the sum of the infectious and recovered populations, but also the uncertainty of their prediction, remain relatively insensitive to variations in the infectious period. To explore whether this is a universal trend, we perform the same analysis for nine different locations at which COVID-19 antibody prevalence was measured in a representative sample of the population, Heinsberg (NRW, Germany) [71], Ada County (ID, USA) [4], New York City (NY, USA) [58], Santa Clara County (CA, USA) [3], Denmark [14], Geneva Canton (Switzerland) [72], Netherlands [67], Rio Grande do Sul (Brasil) [66], and Belgium [64]. The fraction of the symptomatic population in these nine locations is *v*s = 20.00%, 7.90%, 5.76%, 1.77%, 6.95%, 10.34%, 17.31%, 8.10%, and 10.21% respectively, broadly representing the range of reported symptomatic versus estimated total cases worldwide [3, 4, 9, 12, 14, 18, 24, 32, 34, 46, 58, 59, 64, 65, 68, 66, 67, 71, 72, 74, 76, 78, 81]. Of the nine locations we analyzed here, Heinsberg tested IgG and IgA, Ada County tested IgG, New York City tested IgG, Santa Clara County tested IgG and IgM, Denmark tested IgG and IgM, Geneva tested IgG, the Netherlands tested IgG, IgM and IgA, Rio Grande do Sul tested IgG and IgM, and Belgium tested IgG. While we did include reported uncertainty on the seroprevalence data, seroprevalence would likely have been higher if all locations had tested for all three antibodies. Despite these differences, the effective reproduction numbers *R*(*t*) and the infectious and recovered populations *I*s, *I*a, and *R* in Figure 6 display remarkably similar trends: In most locations, the effective reproduction number *R*(*t*) drops rapidly to values below one within a window of about three weeks after the lockdown date. However, the maximum infectious population, a value that is closely monitored by hospitals and health care systems, varies significantly ranging from 0.28% and 0.38% in Rio Grande do Sul and Denmark respectively to 3.54% and 6.11% in Heinsberg and New York City respectively. This is consistent with the reported ‘superspreader’ events in these last two locations. An effect that we do not explicitly address is that immune response not only results COVID-19 antibodies (humoral response), but also from innate and cellular immunity [23]. While it is difficult to measure the effects of the unreported asymptomatic group directly, and discriminate it precisely from innate and cellular immunity, mathematical models can provide valuable insight into how this population modulates the outbreak dynamics and the potential of successful outbreak control [40]. ### Simulations provide a window into the outbreak date Santa Clara County was home to the first individual who died with COVID-19 in the United States. Although this happened as early as February 6, the case remained unnoticed until April 22 [1]. The unexpected new finding suggests that the new coronavirus was circulating in the Bay Area as early as January. The estimated uncertainty on the exposed, symptomatic infectious, and asymptomatic infectious populations of our model allows us to estimate the initial outbreak date dates back to January 20 (Figure 7). This back-calculated early outbreak date is in line with our intuition that COVID-19 is often present in a population long before the first official case is reported. Interestingly, our analysis comes to this conclusion purely based on a local serology antibody study [3] and the number of reported cases after lockdown [63]. ### Limitations Our approach naturally builds in and learns several levels of uncertainties. By design, this allows us to estimate sensitivities and credible intervals for a number of important model parameters and discover important features and trends. Nevertheless, it has a few limitations, some of them by design, some simply limited by the current availability of data: First, our current SEIIR model assumes a similar contact rate *β*(*t*) for symptomatic and asymptomatic individuals. While we can easily adjust this in the model by defining individual symptomatic and asymptomatic rates *βs*(*t*) and >*β*a(*t*), we currently do not have data on the temporal evolution of the hidden asymptomatic infectious population *I*a(*t*) and longitudinal large population antibody follow-up studies would be needed to appropriately calibrate *β*a(*t*). Second, the ratio between the symptomatic and asymptomatic populations *v*s: *v*a can vary over time, especially, as we have shown, if both groups display notably different dynamics, in our model represented through *C*s and *C*a. Since this can have serious effects on the overall reproduction number *R*(*t*), and with it on required outbreak control strategies, it seems critical to perform more tests and learn the dynamics of the fractions *v*s(*t*) and *v*a(*t*) of both groups. Third, and this is not only true for our specific model, but for COVID-19 forecasts in general, all predictions can be sensitive to the amount of testing in time. As such, they crucially rely on testing policies and testing capacities. We expect to see a significant increase in the symptomatic-to-asymptomatic, or rather detected- to-undetected, ratio as we move towards systematically testing larger fractions of the population and more and more people who have no symptoms at all. The intensity of testing increases in most locations during our simulation period. For example, in Santa Clara County, testing was extremely limited until early April, increased substantially in the first three weeks of April, and even more after. Including limited testing and more undocumented cases during the early outbreak would shift the case distribution towards earlier days, and predict an even earlier outbreak date. When longitudinal antibody studies become available, additional methodologies can be developed to correct for this limitation which include a dynamic evolution for the symptomatic/asymptomatic fractions. Fourth, while we have included uncertainty in the seroprevalence data, the nine locations we analyzed here tested different types of antibodies and had different sampling procedures. Seroprevalence could have been higher if all locations had tested for the same three antibodies and data may differ depending on biases introduced by the sampling procedure. Finally, our current model does not explicitly account for innate and cellular immunity. If the fraction of the population with innate and cellular immunity is substantially high, we would anticipate a smaller susceptible population and a larger and earlier protective immunity overall. These, and other limitations related to the availability of information, can be easily addressed and embedded in our model and will naturally receive more clarification as studies and data become available in the coming months. ## 5. Conclusions The rapid and devastating development of the COVID-19 pandemic has raised many open questions about its outbreak dynamics and unsuccessful outbreak control. From an outbreak management standpoint-in the absence of effective vaccination and treatment-the two most successful strategies in controlling an infectious disease are isolating infectious individuals and tracing and quarantining their contacts. Both critically rely on a rapid identification of infections, typically through clinical symptoms. Recent antibody prevalence studies could explain why these strategies have largely failed in containing the COVID-19 pandemic: Increasing evidence suggests that the number of unreported asymptomatic cases could outnumber the reported symptomatic cases by an order of magnitude or more. Mathematical modeling, in conjunction with reported symptomatic case data, antibody seroprevalence studies, and machine learning allows us to infer, in real time, the epidemiology characteristics of COVID-19. We can now visualize the invisible asymptomatic population, estimate its role in disease transmission, and quantify the confidence in these predictions. A better understanding of asymptomatic transmission will help us evaluate strategies to manage the impact of COVID-19 on both our economy and our health care system. A large asymptomatic population is associated with a high risk of community spread and could require a conscious shift from containment to mitigation induced by behavior changes. Our study suggests that, until vaccination and treatment become available, increasing population awareness, encouraging increased hygiene, mandating the use of face masks, restricting travel, and promoting physical distancing could be the most successful strategies to manage the impact of COVID-19 on both our economy and our health care system. ## Data Availability All data are available on public servers and are cited in the references of the manuscript ## Acknowledgements This work was supported by a Stanford Bio-X IIP Seed Grant (M.P. and E.K.), by the NIH Grant U01 HL119578 (E.K.), by a DAAD Fellowship (K.L.), and by the Stanford COVID-19 Seroprevalence Study Fund (J.B. and E.B.). ## Appendix To visualize the effect of asymptomatic transmission on the outbreak dynamics of COVID-19, we perform a sensitivity analysis with respect to the two most relevant parameters associated with asymptomatic transmission, the symptomatic fraction *v*s and the asymptomatic infectious period *C*a. To highlight the effect of both parameters, we visualize the exposed, symptomatic and asymptomatic infectious, and recovered populations of our SEIIR model, ![Formula][62] We assume that the symptomatic and asymptomatic contact rates are equal, *β = β*s *= β*a, and that the differences between symptomatic and asymptomatic transmission manifest themselves exclusively in the infectious periods *C*s = 1/*γ*s and *C*a = 1/*γ*a. For illustrative purposes, rather than using a dynamic contact rate of Gaussian random walk type as we have proposed in the main body of this manuscript, we adopt a smooth, monotonically decreasing dynamic contact rate of hyperbolic tangent type [43], ![Formula][63] where *β* is the basic contact rate, *β*t is the reduced contact rate, *t** is the adaptation time, and *T* is the adaptation speed. To mimic realistic outbreak dynamics, we adopt the basic contact rate *β* *=* 0.65 /days, the reduced contact rate (*β*t = 0.10 /days, the adaptation time *t** = 18.61 days, and the adaptation speed *T* = 10.82 /days from the mean values across the 27 countries of the European union [43]. Similar to the main body of this manuscript, we select a latent period of *A* = 1/*α* = 2.5 days and a symptomatic infectious period of *C*a = 1/*γ*a = 6.5 days. For the sensitivity analysis, we vary the symptomatic fraction *v*s = [0.0,…, 1.0] at fixed *C*a = 3.25 days and *C*a = 13 days and the asymptomatic infectious period *C*a = [3.25,…, 13.00] days at fixed *v*s = 0.2 and *v*s = 0.8. For a symptomatic fraction of *v*s = 0.2, the undercount is four, i.e., for every reported case, there are four unreported cases. For a symptomatic fraction of *v*s = 0.8, the undercount is one fourth, i.e., for every four reported cases, there is one unreported case. Figure 8 highlights the effect of the symptomatic fraction *v*s. The dark red, orange, green, and blue curves correspond to the special case of the regular SEIR model (6) without asymptomatic transmission and thus *v*s = 1.0 and *v*a = 0.0. The lighter curves show the effect of decreasing the symptomatic fraction towards *v*s = 0.0 and *v*a = 1.0. Interestingly, the effect of asymptomatic transmission on the overall outbreak dynamics varies not only with the symptomatic fraction, but is also highly sensitive to the ratio of the symptomatic and asymptomatic infectious periods. For the smaller asymptomatic infectious period of *C*a = 0.5 *C*s, the outbreak dynamics are governed by the orange symptomatic population *I*s and the blue recovered population *R* decreases with increasing asymptomatic transmission *v*a. For the larger asymptomatic infectious period of *C*a = 2 *C*s, the outbreak dynamics are governed by the green asymptomatic population *I*a and the blue recovered population *R* increases with increasing asymptomatic transmission. ![Figure 8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/29/2020.05.23.20111419/F8.medium.gif) [Figure 8:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/F8) Figure 8: SEIIR epidemiology model. Sensitivity with respect to the symptomatic fraction *v*s. Darkest curves correspond to the special case of no asymptomatic transmission, *v*s = 1.0, lightest curves to only asymptomatic transmission, *v*s = 0.0. For a smaller asymptomatic infectious period of *C*a = 0.5 *C*s, the outbreak dynamics are governed by the symptomatic population *I*s and the recovered population *R* decreases with increasing asymptomatic transmission *v*a (left). For a larger asymptomatic infectious period of *C*a = 2 *C*s, the outbreak dynamics are governed by the asymptomatic population *I*a and the recovered population increases with increasing asymptomatic transmission (right). Latent period *A =* 2.5 days, symptomatic infectious period *C*s = 6.5 days, asymptomatic infectious period *C*a = 3.25 days (left) and *C*a = 13 days (right), and contact rate ![Graphic][64] with *β* = 0.65 /days, *β*t = 0.10 /days, *t** = 18.61 days, and *T* = 10.82 days. Figure 9 shows the effect of the asymptomatic infectious period *C*a. The darkest red, orange, green, and blue curves correspond to a larger asymptomatic infectious period of *C*a = 2 *C*s, the lightest curves to a smaller period *C*a = 0.5 *C*s. The middle curves correspond to the special case for which the asymptomatic and symptomatic populations display the same dynamics, *C*a = *C*s, such that the overall recovered population becomes independent of the symptomatic fraction. This implies that middle blue recovered populations of both graphs are identical and also identical to the dark blue curves of Figure 8. For the larger symptomatic fraction of *v*s = 0.8, the outbreak dynamics are governed by the orange symptomatic population *I*s and the blue recovered population *R* is relatively insensitive to the asymptomatic infectious period *C*a. For the smaller symptomatic fraction of *v*s = 0.2, the outbreak dynamics are governed by the green asymptomatic population *I*a and the blue recovered population *R* increases markedly with an increasing asymptomatic infectious period *C*a. ![Figure 9:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/29/2020.05.23.20111419/F9.medium.gif) [Figure 9:](http://medrxiv.org/content/early/2020/08/29/2020.05.23.20111419/F9) Figure 9: SEIIR epidemiology model. Sensitivity with respect to the asymptomatic infectious period *C*a. Darkest curves correspond to a larger asymptomatic infectious period of *C*a = 2 *C*s, lightest curves to a smaller period *C*a = 0.5 *C*s. For a larger symptomatic fraction of *v*s = 0.8, the outbreak dynamics are governed by the symptomatic population *I*s and the recovered population *R* is relatively insensitive to the asymptomatic infectious period *C*a (left). For a smaller symptomatic fraction of *v*s = 0.2, the outbreak dynamics are governed by the asymptomatic population *I*a and the recovered population *R* increases markedly with an increasing asymptomatic infectious period *C*a (right). Latent period *A* = 2.5 days, symptomatic infectious period *C*s = 6.5 days, asymptomatic infectious period *C*a = 3.25 days (left) and *C*a = 13 days (right), and contact rate ![Graphic][65] with *β* = 0.65 /days, *βt* = 0.10 /days, *t** = 18.61 days, and *T* = 10.82 days. Taken together, this sensitivity analysis shows that the outbreak dynamics of COVID-19 are highly sensitive to asymptomatic transmission. This sensitivity becomes more pronounced the more the asymptomatic infectious period differs from its symptomatic counterpart. This difference could be a natural result of asymptomatic individuals isolating themselves less since they are unaware that they can infect others. Knowing the precise dynamics of asymptomatic transmission is therefore critical for a successful management of COVID-19. ## Footnotes * Appendix added * Received May 23, 2020. * Revision received August 28, 2020. * Accepted August 29, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. [1]. E. Allday, M. Kawahara. First known U.S. coronavirus death occurred on Feb. 6 in Santa Clara County. San Francisco Chronicle, April 22, 2020. [https://www.sfchronicle.com/health/article/First-known-U-S-coronavirus-death-occurred-on-15217316.php](https://www.sfchronicle.com/health/article/First-known-U-S-coronavirus-death-occurred-on-15217316.php) assessed: Jun 15, 2020. 2. [2]. J. L. Aron, I. B. Schwartz. Seasonality and period-doubling bifurcation in an epidemic model. Journal of Theoretical Biology 110 (1984) 665-679. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0022-5193(84)80150-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6521486&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1984TT57500010&link_type=ISI) 3. [3]. E. Bendavic, B. Mulaney, N. Sood, S. Shah, E. Ling, R. Bromley-Dulfano, C. Lai, Z. Weissberg, R. Saavedra-Walker, J. Tedrow, D. Tversky, A. Bogan, T. Kupiec, D. Eichner, R. Gupta, J.P.A. Ioannidis, J. Bhattacharya. COVID-19 antibody seroprevalence in Santa Clara County, California. medRxiv doi:10.1101/2020.04.14.20062463. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNC4xNC4yMDA2MjQ2M3YyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 4. [4]. A. Bryan, G. Pepper, M.H. Wener, S.L. Fink, C. Morishima, A. Chaudhary, K.R. Jerome, P.C. Mathias, A.L. Greninger Performance Characteristics of the Abbott Architect SARS-CoV-2 IgG Assay and Seroprevalence in Boise, Idaho J Clin Microbiol. 2020 doi:10.1128/JCM.00941-20 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjE0OiI1OC84L2UwMDk0MS0yMCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA4LzI5LzIwMjAuMDUuMjMuMjAxMTE0MTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 5. [5].Painel de casos de doenca pelo coronavirus 2019 (COVID-19) no Brasil pelo Ministerio da Saude. [www.covid.saude.gov.br/](http://www.covid.saude.gov.br/) assessed: Jun 15, 2020. 6. [6].Center for Disease Control and Prevention. Large-scale Geographic Seroprevalence Surveys. [https://www.cdc.gov/coronavirus/2019-ncov/cases-updates/geographic-seroprevalence-surveys.html](https://www.cdc.gov/coronavirus/2019-ncov/cases-updates/geographic-seroprevalence-surveys.html) assessed: Aug 24, 2020. 7. [7]. A.M. Cuomo. Amid ongoing COVID-19 pandemic, governor cuomo announces results of completed antibody testing study. [https://www.governor.ny.gov/news/amid-ongoing-covid-19-pandemic-governor-cuomo-announces-results-completed-antibody-testing](https://www.governor.ny.gov/news/amid-ongoing-covid-19-pandemic-governor-cuomo-announces-results-completed-antibody-testing). assessed: Jun 15, 2020. 8. [8]. J. Dehning, J. Zierenberg, F.P Spitzner, M. Wibral, J.P. Neto, M. Wilczek, V. Priesemann. Inferring change points in the spread of COVID-19 reveals the effectiveness of interventions. Science (2020) doi:10.1126/science.abb9789. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNjkvNjUwMC9lYWJiOTc4OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA4LzI5LzIwMjAuMDUuMjMuMjAxMTE0MTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 9. [9]. A. Doi, K. Iwata, H. Kuroda, T. Hasuike, A. Kanda, T. Nagao, H. Nishioka, K. Tomii, T. Morimoto, Y. Kihara. Seroprevalence of novel coronavirus disease (COVID-19) in Kobe, Japan. medRxiv doi:10/1101/2020.04.26.20079822. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10/1101/2020.04.26.20079822&link_type=DOI) 10. [10]. E. Dong, L. Gardner. An interactive web-based dashboard to track COVID-19 in real time. Lancet Infectious Disease (2020) doi:10.1016/S1473-3099(20)30120-1. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30120-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32087114&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 11. [11].European Centre for Disease Prevention and Control. 2020. Situation update worldwide. [www.ecdc.europa.eu/en/geographical-distribution-2019-ncov-cases](http://www.ecdc.europa.eu/en/geographical-distribution-2019-ncov-cases) assessed: Jun 15, 2020. 12. [12]. M. Emmenegger, E. De Cecco, D. Lamparter, R.P.B. Jacquat, D. Ebner, M.M. Schneider, I. Condado Morales, D. Schneider, B. Dogancay, J. Guo, A. Wiedmer, J. Domange, M. Imeri, R. Moos, C. Zografou, C. Trevisan, A. Gonzalez-Guerra, A. Carrella, I.L. Dubach, C.K. Xu, G. Meisl, V. Kosmoliaptsis, T. Malinauskas, N. Burgess-Brown, R. Owens, J. Mongkolsapaya, S. Hatch, G.R. Screaton, K. Schubert, J.D. Huck, F. Liu, F. Pojer, K. Lau, D. Hacker, E. Probst-Mueller, C. Cervia, J. Nilsson, O. Boyman, L. Saleh, K. Spanaus, A. von Eckardstein, D.J. Schaer, N. Ban, C. Tsai, J. Marino, G.F.X. Schertler, J. Gottschalk, B.M. Frey, R. Reimann, S. Hornemann, A.M. Ring, T.P.J. Knowles, I. Xenarios, D.I. Stuart, A. Aguzzi Population-wide evolution of SARS-CoV-2 immunity tracked by a ternary immunoassay. medRxiv 2020, doi: 10.1101/2020.05.31.20118554 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4zMS4yMDExODU1NHY0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 13. [13]. R. Engbert, F.R. Drepper. Chance and chaos in population biology-Models of recurrent epidemics and food chain dynamics. Chaos, Solutions & Fractals 4 (1994) 1147-1169. 14. [14]. C. Erikstrup, C.E. Hother, O.B.V Pedersen, K. Molbak, R.L. Skov, D.K. Holm, S. Saekmose, A.C. Nilsson, P.T. Brooks, J.K. Boldsen, C. Mikkelsen, M. Gybel-Brask, E. Sorensen, K.M. Dinh, S. Mikkelsen, B.K. Moller, T. Haunstrup, L. Harritshoj, B.A. Jensen, H. Hjal- grim, S.T. Lillevang, H. Ullum Estimation of SARS-CoV-2 infection fatality rate by real-time antibody screening of blood donors medRxiv doi:10.1101/2020.04.24.20075291. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNC4yNC4yMDA3NTI5MXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 15. [15]. Y. Fang, Y. Nie, M. Penny. Transmission dynamics of the COVID-19 outbreak and effectiveness of government interventions. J. Med. Virol. (2020) 1-15. 16. [16]. A.S. Fauci, H.C. Lane, R.R. Redfield. Covid-19-Navigating the uncharted. New England Journal of Medicine 382 (2020) 1268-1269. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMe2002387&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32109011&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 17. [17]. S. Flaxman, S. Mishra, A. Gandy, H. Unwin, H. Coupland, T. Mellan, H. Zhu, T. Berah, J. Eaton, P. Perez Guzman, N. Schmit, L. Cilloni, K. Ainslie, M. Baguelin, I. Blake, A. Boonyasiri, O. Boyd, L. Cattarino, C. Ciavarella, L. Cooper, Z. Cucunuba Perez, G. Cuomo-Dannenburg, A. Dighe, A. Djaafara, I. Dorigatti, S. Van Elsland, R. Fitzjohn, H. Fu, K. Gaythorpe, L. Geidelberg, N. Grassly, W. Green, T. Hallett, A. Hamlet, W. Hinsley, B. Jeffrey, D. Jorgensen, E. Knock, D. Laydon, G. Nedjati Gilani, P. Nouvellet, K. Parag, I. Siveroni, H. Thompson, R. Verity, E. Volz, C. Walters, H. Wang, Y. Wang, O. Watson, P. Winskill, X. Xi, C. Whittaker, P. Walker, A. Ghani, C. Donnelly, S. Riley, L. Okell, M. Vollmer, N. Ferguson, S. Bhatt Report 13: Estimating the number of infections and the impact of non-pharmaceutical interventions on COVID-19 in 11 European countries Spiral Imperial College London 2020, doi: 10.25561/77731 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.25561/77731&link_type=DOI) 18. [18]. A., Fontanet, L. Tondeur, Y. Madec, R. Grant, C. Besombes, N. Jolly, S. Fernandes Pellerin M.N. Ungeheuer, I. Cailleau, L. Kuhmel, S. Temmam, C. Huon, K.Y. Chen, B. Crescenzo, S. Munier, C. Demeret, L. Grzelak, I. Staropoli, T. Bruel, P. Gallian, S. Cauchemez, S. van der Werf, O. Schwartz, M. Eloit, B. Hoen. Cluster of COVID-19 in northern France: A retrospective closed cohort study. medRxiv doi:10.1101/2020.04.20071134. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2020.04.20071134&link_type=DOI) 19. [19]. C. Fraser, S. Riley, R.M. Anderson, N.M. Ferguson. Factors that make an infectious disease outbreak controllable. Proceedings of the National Academy of Sciences 101 (2004) 6146-6151. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTAxLzE2LzYxNDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOC8yOS8yMDIwLjA1LjIzLjIwMTExNDE5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 20. [20]. M. Gatto, E. Bertuzzo, L. Mari, S. Miccoli, L. Carraro, R. Casagrandi, A. Rinaldo. Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures. Proceedings of the National Academy of Sciences; in press, doi:10.1073/pnas.2004978117. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzE5LzEwNDg0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 21. [21]. F. Gao, L. Han. Implementing the Nelder-Mead simplex algorithm with adaptive parameters. Computational Optimization and Applications 51(2012)259-277. 22. [22]. A. Gelman,, J. Hill. Data Analysis using Regression and Multilevel/Hierarchical Models. Cambridge University Press, 2006. 23. [23]. A. Grifoni, D. Weiskopf, S.I. Ramirez, J. Mateus, J.M. Dan, C. Rydyznski Moderbacher, S.A. Rawlings, A. Sutherland, L. Premku- mar R.S. Jadi, D. Marrama, A.M. de Silva, A. Frazier, A. Carlin, J.A. Greenbaum, B. Peters. F. Krammer. D.M. Smith, S. Crotty, A. Sette. Targets of T cell responses to SARS-CoV-2 coronavirus in humans with COVID-19 disease and unexposed individuals. Cell, doi: 10.1016/j.cell.2020.05.015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.05.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32473127&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 24. [24]. H. Hallal, F. Hartwig, B. Horta, G.D. Victora, M. Silveira, C. Struchiner, L.P. Vidaletti, N. Neumann, L.S. Pellanda, O.A. Dellagostin, M.N. Burattini, A.M. Menezes, F.C. Barros, A.J. Barros, C.G. Victora Remarkable variability in SARS-CoV-2 antibodies across Brazilian 1 regions: nationwide serological household survey in 27 states. medRxiv 2020, doi: 10.1101/2020.05.30.20117531 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4zMC4yMDExNzUzMXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 25. [25]. F.P. Havers, C. Reed, T. Lim, J.M. Montgomery, J.D. Klena, A.J. Hall, A.M. Fry, D.L. Cannon, C.-F. Chiang, A. Gibbons, I. Krapiunaya, M. Morales-Betoulle, K. Roguski, M.A.U. Rasheed, B. Freeman, S. Lester, L. Mills, D.S. Carroll, S.M. Owen, J.A. Johnson, V. Semenova, C. Blackmore, D. Blog, S.J. Chai, A. Dunn, J. Hand, S. Jain, S. Lindquist, R. Lynfield, S. Pritchard, T. Sokol, L. Sosa, G. Turabelidze, S.M. Watkins, J. Wiesman, R.W. Williams, S. Yendell, J. Schiffer, N.J. Thornburg. Seroprevalence of Antibodies to SARS-CoV-2 in 10 Sites in the United States, March 23-May 12, 2020. JAMA Intern Med (2020) doi:10.1001/jamainternmed.2020.4130. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamainternmed.2020.4130&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32692365&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 26. [26]. X. He, E.H. Lau, P. Wu, X. Deng, J. Wang, X. Hao, Y.C. Lau, J.Y. Wong, Y. Guan, X. Tan, X. Mo. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nature Medicine (2020), 26(5), pp.672-675. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-3012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 27. [27].Heinsberg 2020. Aktuelles aus dem Kreishaus. Coronavirus im Kreis Heinsberg. [www.kreis-heinsberg.de/aktuelles](http://www.kreis-heinsberg.de/aktuelles). assessed: Jun 15, 2020. 28. [28]. H. W. Hethcote. The mathematics of infectious diseases. SIAM Review 42 (2000) 599-653. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1137/s0036144500371907&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00016567&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 29. [29]. M.D. Hoffman, A. Gelman. The No-U-Turn sampler: adaptively setting path lengths in Hamiltonian Monte Carlo. Journal of Machine Learning Research (2014), 15(1), 1593-1623. 30. [30]. J.P.A. Ioannidis. The infection fatality rate of COVID-19 inferred from seroprevalence data. medRxiv. doi:10.1101/2020.05.13.20101253. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4xMy4yMDEwMTI1M3YzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 31. [31]. P.K. Jha, L. Cao, J.T. Oden. Bayesian-based predictions of COVID-19 evolution in Texas using multispecies mixture-theoretic continuum models. Computational Mechanics (2020) in press. 32. [32]. I. Jerkovic, T. Ljubic, Z. Basic, I. Kruzic, N. Kunac, J. Bezic, A. Vuko, A. Markotic, S. Andjelinovic SARS-CoV-2 antibody seroprevalence in industry workers in Split-Dalmatia and Sibenik-Knin County, Croatia. medRxiv 2020, doi: 10.1101/2020.05.11.20095158 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4xMS4yMDA5NTE1OHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 33. [33]. W. O. Kermack, G. McKendrick. Contributions to the mathematical theory of epidemics, Part I. Proceedings of the Royal Society London Series A 115 (1927) 700-721. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspa.1927.0118&link_type=DOI) 34. [34]. V. Kraehling, M. Kern, S. Halwe, H. Mueller, C.P.C. Rohde, M. Savini, M. Schmidt, J. Wilhelm, S. Becker, S. Ciesek, R. Gottschalk Epidemiological study to detect active SARS-CoV-2 infections and seropositive persons in a selected cohort of employees in the Frankfurt am Main metropolitan area. medRxiv 2020, doi: 10.1101/2020.05.20.20107730 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4yMC4yMDEwNzczMHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 35. [35]. E. Kuhl. Data-driven modeling of COVID-19 - Lessons learned. Extreme Mechanics Letters 2020) doi:10.1016/j.eml.2020.100921. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.eml.2020.100921&link_type=DOI) 36. [36]. K.L. Lange, R.J.A. Little, M.G. Taylor. Robust statistical modeling using the t distribution. Journal of the American Statistical Association 84 (1989) 881-896. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/2290063&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1989CD46300006&link_type=ISI) 37. [37]. S. A. Lauer, K. H. Grantz, Q. Bi, F. K. Jones, Q. Zheng, H. R. Meredith, A. S. Azman, N. G. Reich, J. Lessler. The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: estimation and application. Annals of Internal Medicine (2020) doi:10.7326/M20-0504. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-0504&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32150748&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 38. [38]. S. Lee, T. Kim, E. Lee, C. Lee, H. Kim, H. Rhee, S.Y. Park, H. Son, S. Yu, J.W. Park, E.J. Choo, S. Park, M. Loeb, T.H. Kim. Clinical course and molecular viral shedding among asymptomatic and symptomatic patients With SARS-CoV-2 infection in a community treatment center in the republic of Korea. JAMA. Internal Medicine (2020) doi:10.1001/jamainternmed.2020.3862. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamainternmed.2020.3862&link_type=DOI) 39. [39]. Q. Li, X. Guan, P. Wu, X. Wang, L. Zhou, Y. Tong, R. Ren, K.S.M. Leung, E.H.Y. Lau, J.Y. Wong, X. Xing, N. Xiang, Y. Wu, C. Li, Q. Chen, D. Li, T. Liu, J. Zhao, M. Liu W. Tu, C. Chen, L. Jin, R. Yang, Q. Wang, S. Zhou, R. Wang, H. Liu, Y. Luo, Y. Liu, G. Shao, H. Li, Z. Tao, Y. Yang, Z. Deng, B. Liu, Z. Ma, Y. Zhang, G. Shi, T.T.Y. Lam, J.T. Wu, G.F. Gao, B.J. Cowling, B. Yang, G.M. Leung, Z. Feng. Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia. New England Journal of Medicine 382 (2020) 1199-1207. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001316&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 40. [40]. R. Li, S. Pei, B. Chen, Y. Song, T. Zhang, W. Yang, J. Shaman. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV2). Science 368 (2020) 489-493. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5MC80ODkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOC8yOS8yMDIwLjA1LjIzLjIwMTExNDE5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 41. [41]. E.A.B.F. Lima, J.T. Oden, B. Wohlmuth, A. Shahmoradi, D.A. Hormuth II, T.E. Yankeelov, L. Scaraboso, T. Horger. Selection and validation of predictive models of radiation effects on tumor growth based on noninvasive imaging data. Computer Methods in Applied Mechanics and Engineering 327 (2017) 277-305. 42. [42]. K. Linka, M. Peirlinck, F. Sahli Costabal, E. Kuhl. Outbreak dynamics of COVID-19 in Europe and the effect of travel restrictions. Computer Methods in Biomechanics and Biomedical Engineering (2020) 23:710-717. 43. [43]. K. Linka, M. Peirlinck, E. Kuhl. The reproduction number of COVID-19 and its correlation with public health interventions. Computational Mechanics (2020) online first. doi:10.1007/s00466-020-01880-8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00466-020-01880-8&link_type=DOI) 44. [44]. K. Linka, P. Rahman, A. Goriely, E. Kuhl. Is it safe to lift COVID-19 travel restrictions? The Newfoundland story. Computational Mechanics (2020) online first. doi:10.1007/s00466-020-01899-x. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00466-020-01899-x&link_type=DOI) 45. [45]. K. Linka, A. Goriely, E. Kuhl. Global and local mobility as a barometer for COVID-19 dynamics. medRxiv (2020) doi:10.1101/2020.06.13.20130658. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNi4xMy4yMDEzMDY1OHYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 46. [46]. D. Ng, G. Goldgof, B. Shy, A. Levine, J. Balcerek, S.P. Bapat, J. Prostko, M. Rodgers, K. Coller, S. Pearce, S. Franz, L. Du, M. Stone, S. Pillai, A. Sotomayor-Gonzalez, V. Servellita, C. Sanchez-San Martin, A. Granados, D.R. Glasner, L.M. Han, K. Truong, N. Akagi, D.N. Nguyen, N. Neumann, D. Qazi, E. Hsu, W. Gu, Y.A. Santos, B. Custer, V. Green, P. Williamson, N.K. Hills, C.M. Lu, J.D. Whitman, S. Stramer, C. Wang, K. Reyes, J. Hakim, K. Sujishi, F. Alazzeh, L. Pharm, C. Oon, S. Miller, T. Kurtz, J. Hackett, G. Simmons, M.P. Busch, C.Y. Chiu SARS-CoV-2 seroprevalence and neutralizing activity in donor and patient blood from the San Francisco Bay Area. medRxiv 2020, doi: 10.1101/2020.05.19.20107482 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4xOS4yMDEwNzQ4MnYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 47. [47]. N.B. Noll, I. Aksamentov, V. Druelle, A. Badenhorst, B. Ronzani, G. Jefferies, J. Albert, R. Neher. COVID-19 Scenarios: an interactive tool to explore the spread and associated morbidity and mortality of SARS-CoV-2. medRxiv doi: 10.1101/2020.05.05.20091363 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4wNS4yMDA5MTM2M3YyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 48. [48].New York City Health Department. NYC Coronavirus Disease 2019 (COVID-19) Data [www.github.com/nychealth/coronavirus-data](http://www.github.com/nychealth/coronavirus-data), assessed: Jun 15, 2020. 49. [49].New York Times. An ongoing repository of data on coronavirus cases and deaths in the U.S. [www.github.com/nytimes/covid-19-data](http://www.github.com/nytimes/covid-19-data), assessed: Jun 15, 2020. 50. [50]. J.T. Oden, R. Moser, O. Ghattas. Computer predictions with quantified uncertainty I. SIAM News (2010) 43:9. 51. [51]. J.T. Oden, I. Babuska, D. Faghihi. Predictive computational science: Computer predictions in the presence of uncertainty. Encyclopedia of Computational Mechanics (2017) doi:10.1002/9781119176817.ecm2101. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/9781119176817.ecm2101&link_type=DOI) 52. [52]. J.T. Oden. Adaptive multiscale predictive modelling. Acta Numerica (2018) 353-450. 53. [53]. X. Pan, D. Chen, Y. Xia, X, Wu, T. Li, X. Ou, L. Zhou, J. Liu. Asymptomatic cases in a family cluster with SARS-CoV-2 infection. Lancet Infectious Diseases 20 (2020) 410-411. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30114-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32087116&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 54. [54]. S.W. Park, D.M. Cornforth, J. Dushoff, J.W. Weitz. The time scale of asymptomatic transmission affects estimates of epidemic potential in the COVID-19 outbreak. medRxiv 2020.03.09.20033514, doi 10.1101/2020.03.09.20033514. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMy4wOS4yMDAzMzUxNHYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 55. [55]. S.W. Park, K. Sun, C. Viboud, B.T. Grenfell, J. Dushoff. Potential roles of social distancing in mitigating the spread of coronavirus disease 2019 (COVID-19) in South Korea. medRxiv 2020.03.27.20045815, doi 10.1101/2020.03.27.20045815 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMy4yNy4yMDA0NTgxNXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 56. [56]. M. Peirlinck, F. Sahli Costabal, K.L. Sack, J.S. Choy, G.S. Kassab, J.M. Guccione, M. De Beule, P. Segers, E. Kuhl. Using machine learning to characterize heart failure across the scales. Biomechanics and Modeling in Mechanobiology 18 (2019) 1987-2001. 57. [57]. M. Peirlinck, K. Linka, F. Sahli Costabal, E. Kuhl. Outbreak dynamics of COVID-19 in China and the United States. Biomechanics and Modeling in Mechanobiology (2020) in press, doi:10.1101/2020.04.06.20055863. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNC4wNi4yMDA1NTg2M3YyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 58. [58]. J. Reifer, N. Hayum, B. Heszkel, I. Klagsbald, V.A. Streva SARS-CoV-2 IgG Antibody Responses in New York City medRxiv doi:10.1101/2020.2005.2023.20111427. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2020.2005.2023.20111427&link_type=DOI) 59. [59]. E.S. Rosenberg, J.M. Tesoriero, E.M. Rosenthal, R. Chung, M.A. Barranco, L.M. Styer, M.M. Parker, S.J. Leung, J.M. Morne, D. Greene, D. R. Holtgrave, D. Hoefer, J. Kumar, T. Udo, B. Hutton, H.A. Zucker Cumulative incidence and diagnosis of SARS-CoV-2 infection in New York. medRxiv 2020, doi: 10.1101/2020.05.25.20113050 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4yNS4yMDExMzA1MHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 60. [60]. F. Sahli Costabal, K. Matsuno, J. Yao, P. Perdikaris, E. Kuhl. Machine learning in drug development: Characterizing the effect of 30 drugs on the QT interval using Gaussian process regression, sensitivity analysis, and uncertainty quantification. Computer Methods in Applied Mechanics and Engineering 348 (2019) 313-333. 61. [61]. J. Salvatier, T.V. Wiecki, C. Fonnesbeck. Probabilistic programming in Python using PyMC3. Peer Journal Computational Science 2 (2016) e55. 62. [62]. S. Sanche, Y.T. Lin, C. Xu, R. Romero-Severson, N. Hengartner, R. Ke. High contagiousness and rapid spread of severe acute respiratory syndrome coronavirus 2. Emerging Infectious Disease (2020) doi:10.3201/eid2607.200282. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2607.200282&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32255761&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 63. [63].Santa Clara County COVID-19 Cases and Hospitalizations Dashboard. [www.sccgov.org](http://www.sccgov.org) assessed: Jun 15, 2020. 64. [64].Sciensano: National Public Health Research Institute of Belgium Weekly Epidemiological Bulleting from May 29,2020 [https://covid-19.sciensano.be](https://covid-19.sciensano.be) assessed: Jun 15, 2020. 65. [65]. M. Shakiba, S. Nazari, F. Mehrabian, S.M. Rezvani, Z. Ghasempour, A. Heidarzadeh. Seroprevalence of COVID-19 virus infection in Guilan province, Iran. medRxiv doi:10.1101/2020.04.26.20079244. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNC4yNi4yMDA3OTI0NHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 66. [66]. M. Silveira, A. Barros, B. Horta, L. Pellanda, G. Victora, O. Dellagostin, C. Struchiner, M. Burattini, A. Valim, E. Berlezi, J. Mesa, M.L. Ikeda, M. Mesenburg, M. Mantesso, M. Dall’Agnol, R. Bittencourt, F. Hartwig, A.M. Menezes, F.C. Barros, P. Hallal, C.G. Victora Repeated population-based surveys of antibodies against SARS-CoV-2 in Southern Brazil medRxiv doi: 10.1101/2020.05.01.20087205 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4wMS4yMDA4NzIwNXYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 67. [67]. E. Slot, B.M. Hogema, C. B.E.M. Reusken, J.H. Reimerink, M. Molier, H.M. Karregat, J. Ijist, V.M.J. Novotny, R.A.W. van Lier, H.L. Zaaijer. Herd immunity is not a realistic exit strategy during a COVID-19 outbreak. Research Square 2020. [https://dx.doi.org/10.21203/rs3.rs-25862/v1](https://dx.doi.org/10.21203/rs3.rs-25862/v1). 68. [68]. C.J. Snoeck, M. Vaillant, T. Abdelrahman, V.P. Satagopam, J.D. Turner, K. Beaumont, C.P.C. Gomes, J.V. Fritz, V.E. Schroder, A. Kaysen, L. Pavelka, L. Stute, G. Ramos Meyers, L. Pauly, M. Hansen, C. Pauly, G.A. Aguayo, M. Perquin, A. Hanff, S. Ghosh, M. Gantenbein, L. Huiart, M. Ollert, R. Kruger Prevalence of SARS-CoV-2 infection in the Luxembourgish population: the CON-VINCE study. medRxiv doi:10.1101/2020.05.11.20092916 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4xMS4yMDA5MjkxNnYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 69. [69]. N. Sood, P. Simon, P. Ebner, D. Eichner, J. Reynolds, E. Bendavid, J. Bhattacharya. Seroprevalence of SARS-CoV-2-specific antibodies among adults in Los Angeles County, California, on April 10-11, 2020. Journal of the American Medical Association. doi:10.1001/jama.2020.8279. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.8279&link_type=DOI) 70. [70]. B. Spellberg, M. Haddix, R. Lee, S. Butler-Wu, P. Holtom, H. Yee, P. Gounder. Community prevalence of SARS-CoV-2 among patients with influenzalike illnesses presenting to a Los Angeles medical center in March 2020. Journal of the American Medical Association. doi:10.1001/jama.2020.4958. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.4958&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32232421&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 71. [71]. H. Streeck, B. Schulte, B.M. Kummerer, E. Richter, T. Holler, C. Fuhrmann, E. Bartok, R. Dolscheid, M. Berger, L. Wessendorf, M. Eschbach-Bludau, A. Kellings, A. Schwaiger, M. Coenen, P. Hoffmann, B. Stoffel-Wagner, M.M. Nothen, A.M. Eis-Hubinger, M. Exner2, R.M. Schmithausen, M. Schmid, G. Hartmann. Infection fatality rate of SARS-CoV-2 infection in a German community with a superspreading event. medRxiv doi:10.1101/2020.05.04.20090076. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4wNC4yMDA5MDA3NnYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 72. [72]. S. Stringhini, A. Wisniak, G. Piumatti, A.S. Azman, S.A. Lauer, H. Baysson, D. De Ridder, D. Petrovic, S. Schrempft, K. Marcus, I. Arm- Vernez, S. Yerly, O. Keiser, S. Hurst, K. Posfay-Barbe, D. Trono, D. Pittet, L. Getaz, F. Chappuis, I. Eckerle, N. Vuilleumier, B. Meyer, A. Flahault, L. Kaiser, I. Guessous Repeated seroprevalence of anti-SARS-CoV-2 IgG antibodies in a population-based sample from Geneva, Switzerland medRxiv doi:10.1101/2020.05.02.20088898 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4wMi4yMDA4ODg5OHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 73. [73].Swiss Federal Office of Public Health FOPH COVID19 case numbers Cantons of Switzerland and Principality of Liechtenstein. [www.github.com/openZH/covid19/](http://www.github.com/openZH/covid19/), assessed: Jun 15, 2020. 74. [74]. M. Takita, T. Matsumura, K. Yamamoto, E. Yamashita, K. Hosoda, T. Hamaki, E. Kusumi, Regional Difference in Seroprevalence of SARS- CoV-2 in Tokyo: Results from the community point-of-care antibody testing. medRxiv 2020, doi: 10.1101/2020.06.03.20121020 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNi4wMy4yMDEyMTAyMHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 75. [75]. B. Tang, F. Xia, N.L. Bragazzi, Z. McCarthy, X. Wang, S. He, X. Sun, S. Tang, Y. Xiao, J. Wu. Lessons drawn from China and South Korea for managing COVID-19 epidemic: insights from a comparative modeling study. Bulletin of the World Health Organization. doi: [http://dx.doi.org/10.2471/BLT.20.257238](http://dx.doi.org/10.2471/BLT.20.257238). 76. [76]. C. Thompson, N. Grayson, R.S. Paton, J. Lourenco, B.S. Penman, L. Lee. Neutralising antibodies to SARS coronavirus 2 in Scottish blood donors - a pilot study of the value of serology to determine population exposure. medRxiv doi:10.1101/2020.04.13.20060467. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNC4xMy4yMDA2MDQ2N3YyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 77. [77]. T.A. Treibel, C. Manisty, M. Burton, A. McKnight, J. Lambourne, J.B. Augusto, X. Couto-Parada, T. Cutino-MOguel, M. Noursadeghi, J.C. Moon. COVID-19: PCR screening of asympotmatic healthcare workers at London hospital. Lancet. doi:10.1016/S0140-6736(20)31100-4. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)31100-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 78. [78]. L. Valenti, A. Bergna, S. Pelusi, F. Facciotti, A. Lai, M. Tarkowski, A. Berzuini, F. Caprioli, L. Santoro, G. Baselli, C. Della Ventura, E. Erba, S. Bosari, M. Galli, G. Zehender, D. Prati SARS-CoV-2 seroprevalence trends in healthy blood donors during the COVID-19 Milan outbreak. medRxiv 2020, doi: 10.1101/2020.05.11.20098442 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4xMS4yMDA5ODQ0MnYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjkvMjAyMC4wNS4yMy4yMDExMTQxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 79. [79]. A. Viguerie, G. Lorenzo, F. Auricchio, D. Baroli, T.J.R. Hughes, A. Patton, A. Reali, T.E. Yankeelov, A. Veneziani. Simulating the spread of COVID-19 via a spatially-resolved SEIRD model with heterogeneous diffusion. Oden Institute Report 20-09, Austin. 80. [80].World Health Organization. Transmission of SARS-CoV-2: implications for infection prevention precautions. [https://www.who.int/news-room/commentaries/detail/transmission-of-sars-cov-2-implications-for-infection-prevention-precautions](https://www.who.int/news-room/commentaries/detail/transmission-of-sars-cov-2-implications-for-infection-prevention-precautions); assessed: July 9, 2020. 81. [81]. X. Wu, B. Fu, L. Chen, Y. Feng. Serological tests facilitate identification of asymptomatic SARS-CoV-2 infection in Wuhan, China. Journal of Medical Virology. doi:10.1002/jmv.25904. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.25904&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32311142&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F29%2F2020.05.23.20111419.atom) 82. [82]. T.I. Zohdi. An agent-based computational framework for simulation of competing hostile planet-wide populations. Computer Methods in Applied Mechanics and Engineering (2017) 314: 513-526. 83. [83]. T.I. Zohdi. An agent-based computational framework for simulation of global pandemic and social response on Planet X. Computational Mechanics (2020) in press. [1]: /embed/graphic-2.gif [2]: /embed/graphic-3.gif [3]: /embed/graphic-4.gif [4]: /embed/graphic-5.gif [5]: /embed/graphic-6.gif [6]: /embed/graphic-7.gif [7]: /embed/inline-graphic-1.gif [8]: /embed/graphic-9.gif [9]: /embed/graphic-10.gif [10]: /embed/inline-graphic-2.gif [11]: /embed/inline-graphic-3.gif [12]: /embed/inline-graphic-4.gif [13]: /embed/inline-graphic-5.gif [14]: /embed/graphic-12.gif [15]: /embed/graphic-13.gif [16]: /embed/inline-graphic-6.gif [17]: /embed/inline-graphic-7.gif [18]: /embed/graphic-14.gif [19]: /embed/inline-graphic-8.gif [20]: /embed/inline-graphic-9.gif [21]: /embed/inline-graphic-10.gif [22]: /embed/inline-graphic-11.gif [23]: /embed/inline-graphic-12.gif [24]: /embed/inline-graphic-13.gif [25]: /embed/inline-graphic-14.gif [26]: /embed/inline-graphic-15.gif [27]: /embed/inline-graphic-16.gif [28]: /embed/inline-graphic-17.gif [29]: /embed/inline-graphic-18.gif [30]: /embed/inline-graphic-19.gif [31]: /embed/inline-graphic-20.gif [32]: /embed/inline-graphic-21.gif [33]: /embed/inline-graphic-22.gif [34]: /embed/inline-graphic-23.gif [35]: /embed/graphic-16.gif [36]: /embed/graphic-17.gif [37]: /embed/inline-graphic-24.gif [38]: /embed/inline-graphic-25.gif [39]: /embed/inline-graphic-26.gif [40]: /embed/inline-graphic-27.gif [41]: /embed/inline-graphic-28.gif [42]: /embed/inline-graphic-29.gif [43]: /embed/inline-graphic-30.gif [44]: /embed/inline-graphic-31.gif [45]: /embed/inline-graphic-32.gif [46]: /embed/inline-graphic-33.gif [47]: /embed/inline-graphic-34.gif [48]: /embed/inline-graphic-35.gif [49]: /embed/inline-graphic-36.gif [50]: /embed/inline-graphic-37.gif [51]: /embed/inline-graphic-38.gif [52]: /embed/inline-graphic-39.gif [53]: /embed/inline-graphic-40.gif [54]: /embed/graphic-19.gif [55]: /embed/inline-graphic-41.gif [56]: /embed/inline-graphic-42.gif [57]: /embed/inline-graphic-43.gif [58]: F3/embed/inline-graphic-44.gif [59]: F3/embed/inline-graphic-45.gif [60]: F5/embed/inline-graphic-46.gif [61]: F7/embed/inline-graphic-47.gif [62]: /embed/graphic-25.gif [63]: /embed/graphic-26.gif [64]: F8/embed/inline-graphic-48.gif [65]: F9/embed/inline-graphic-49.gif