Why estimating population-based case fatality rates during epidemics may be misleading ====================================================================================== * Lucas Böttcher * Mingtao Xia * Tom Chou ## Abstract Different ways of calculating mortality ratios during epidemics can yield widely different results, particularly during the COVID-19 pandemic. We formulate both a survival probability model and an associated infection duration-dependent SIR model to define individual- and population-based estimates of dynamic mortality ratios. The key parameters that affect the dynamics of the different mortality estimates are the incubation period and the length of time individuals were infected before confirmation of infection. We stress that none of these ratios are accurately represented by the often misinterpreted case fatality ratio (CFR), the number of deaths to date divided by the total number of infected cases to date. Using available data on the recent SARS-CoV-2 outbreaks and simple assumptions, we estimate and compare the different dynamic mortality ratios and highlight their differences. Informed by our modeling, we propose a more systematic method to determine mortality ratios during epidemic outbreaks and discuss sensitivity to confounding effects and errors in the data. ## 1. INTRODUCTION The mortality ratio is a key metric describing the severity of a viral disease. It changes in time and can be measured in a number of ways during an epidemic. One common metric is the case fatality ratio (CFR), found by dividing the total number of deaths to date, *D*(*t*), by the total number of all cases to date *N*(*t*) [1, 2]. In the recent outbreaks of SARS-CoV-2, the CFR has been estimated from aggregated population data. We show examples of CFR curves in Fig. 1 and in the *Supplemental Information* (SI). As of March 26, 2020, the global CFR(*t*) = *D*(*t*)*/N*(*t*) = 21, 306*/*472, 762 ≈ 4.5% [3], while CFRs in individual regions vary significantly. Clearly, this estimate would correspond to an actual mortality ratio only if *all* the remaining unresolved individuals recover. However, some of these unresolved cases will lead to death and thus to a gradual increase of the estimated mortality ratio over time. Despite the underestimation of this type of population-based measurement, it is still commonly being used by various health officials and is often inconsistently defined as deaths/(deaths + recovereds) even though this difference has been clearly distinguished [4]. ![FIG. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.26.20044693/F1.medium.gif) [FIG. 1.](http://medrxiv.org/content/early/2020/03/30/2020.03.26.20044693/F1) FIG. 1. Mortality-ratio estimates. (a) Evolution of the cumulative number of infected (red), death (black), and recovered (green) cases. The size of the circles indicates the number of cases in the respective compartments on a certain day. (b–c) Estimates of mortality ratios (see Eqs. (8) and (14)) of SARS-CoV-2 infections in China and Italy. The “delayed” mortality-ratio estimate CFRd corresponds to the number of deaths to date divided by total number of cases at time *t* − *τ*res. Many studies use CFRd, although this metric underestimates the individual-based mortality (defined below). Another population-based mortality ratio is *M*p(*t*), the number of deaths divided by the sum of death and recovered cases, up to time *t*. The data are based on Ref. [5]. During the severe 2003 acute respiratory syndrome (SARS) outbreak in Hong Kong, the World Health Organization (WHO) also used the aforementioned estimate to obtain a CFR of 4.5% while the final values approached 17.0% [6, 7]. For the ongoing SARS-CoV-2 outbreaks, analyses by WHO and other institutions still use the CFR = *D*(*t*)*/N*(*t*) metric (see Table I). Since actual mortality probabilities are important measures for assessing the risks associated with epidemic outbreaks, typical underestimations by CFRs may lead to insufficient countermeasures and a more severe epidemic [8, 9]. View this table: [TABLE I.](http://medrxiv.org/content/early/2020/03/30/2020.03.26.20044693/T1) TABLE I. Different CFR estimates of COVID-19. An unambiguous, physiologically-based definition of mortality ratio is the probability that a single, newly infected individual will eventually die of the disease. If there are sufficient individual-level or cohort data, these probabilities can be further stratified according to patient age, gender, health condition, etc. [10]. The mortality ratio or probability of death should be an *intrinsic* property of the virus and the infected individual, depending on age, health, access to health care, etc. This intrinsic probability ought not to be *directly* dependent on the population-level dynamics of infected and recovered individuals. Thus, it can in principle be framed by a model for the survival probability of a single infected individual. Whether this individual infects others does not directly affect his probability of eventually dying. In section II A, we derive a model describing the probability *M*1(*t*) that an infected individual dies or recovers before time *t*. Importantly, these models incorporate the duration of infection (including an incubation period) before a patient tests positive at time *t* = 0. However, as mentioned above, the CFR and other mortality measures are typically reported based on population data. Do these population-based measures, including CFR, provide reasonable measures of the probability of death of an individual? In section II B, we describe how mortality ratios are defined within population-level models, specifically, a disease duration-structured SIR model. We will show that population-based estimates are typically not a meaningful measure of mortality, but that under simplifying assumptions, the mortality ratio *M*p(*t*) is more closely related to the number of deaths to date divided by the number of dead plus the number of *recovered* individuals to date [4]. In the simplest approximation, the mortality ratio is currently (as of March 26, 2020) 21, 306*/*(21, 306+114, 749) ≈ 15.7% [3], significantly higher than the March 26, 2020 CFR ≈ 4.5% estimate. We use the same estimates for the rate parameters in our individual and population models to compute the different mortality ratios. Note that in general, both the individual mortality probability *M*1(*t*) and the population-based estimates *M*p(*t*) depend on the time of measurement *t*. By critically analyzing these estimates and another *ad hoc* “delayed” ratio CFRd, we illustrate and interpret the differences among these measures and discuss how changes or uncertainty in the data affect them. In section III, we summarize our results and identify a correction factor to transform population-level mortality estimates into individual mortality probabilities. ## II. RESULTS ### A. Intrinsic individual mortality rate Consider an individual that, at the time of positive testing (*t* = 0), had been infected for a duration *τ*1. A “survival” probability density can be defined such that *P*(*τ, t*|*τ*1)d*τ* is the probability that the patient is still alive and infected (not recovered) at time *t >* 0 and has been infected for a duration between *τ* = *t* + *τ*1 and *τ* + d*τ*. Since *τ*1 is unknown, it must be estimated or averaged over. The individual survival probability evolves according to ![Formula][1] where the death and recovery rates, *µ*(*τ, t*| *τ*1) and *c*(*τ, t*| *τ*1), depend explicitly on the duration of infection at time *t* and implicitly on patient health and age *a* [15]. They may also depend explicitly on time *t* to reflect changes in clinical policy or available health care. For example, enhanced medical care may decrease the death rate *µ*, giving the individual’s intrinsic physiological processes a chance to cure the patient. These intrinsic individual-based death and recovery rates do not directly depend on population-level viral transmission. Equation (1), assuming an initial condition of one particular individual who has been infected for time *τ*1 at the time of positive test, can be solved using the method of characteristics. From the solution *P*(*τ, t*| *τ*1) one can derive the probabilities of death and recovery by time *t* as ![Formula][2] The probability that an individual died before time *t*, conditioned on resolution (either death or recovery), is then defined as ![Formula][3] where we have explicitly indicated the dependence on the duration of the infection prior to confirmation of infection. These formulae also depend on all other relevant patient attributes such as age, accessibility to health care, etc. In the long-time limit, when resolution has occurred (*P*d(∞) + *P*r(∞) = 1), the individual mortality ratio is simply *M*1(∞) = *P*d(∞). This result relies only on intrinsic individual rate parameters and is completely independent of disease transmission at the population level. In order to capture the dependence of death and recovery rates on the time an individual has been infected, we propose a constant recovery rate *c* and a simple piece-wise constant death rate *µ*(*τ* |*τ*1) that is not explicitly a function of time *t*: ![Formula][4] where *τ*inc is the incubation-time parameter, the time after infection during which an individual remains asymptomatic. During this incubation period, the patient has zero death rate but can recover by clearing the virus. In other words, some patients fully recover without ever developing serious symptoms. For coronavirus infections, the incubation period appears to be highly variable with a mean of *τ*inc ≈ 6.4 days [17]. We can estimate *µ*1 and *c* using individual patient data where 19 patients (outside Hubei) had been tracked from the date on which their first symptoms occurred until the disease resolved [16]. Two out of 19 patients died, on average, 20.5 days after first symptoms occurred and the mean recovery time of the remaining 17 patients is 16.8 days. We show the recovery-time distribution in Fig. 2(a). Since we know that the mortality ratio in this dataset is 2*/*19, we can determine the dependence between *µ*1 and *c* according to *µ*1*/*(*µ*1 + *c*) ≈ 2*/*19 (or *c/µ*1 ≈ 8.5). The constant recovery and after-incubation period death rates [18] are estimated to be ![Formula][5] ![FIG. 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.26.20044693/F2.medium.gif) [FIG. 2.](http://medrxiv.org/content/early/2020/03/30/2020.03.26.20044693/F2) FIG. 2. Individual mortality ratio. (a) Recovery time after first symptoms occurred based on individual data of 17 patients [16]. (b) Death- and recovery rates as defined in Eq. (4). The death rate *µ*(*τ*1) approaches *µ*1 for *τ*1 *> τ*inc, where *τ*inc is the incubation period and *τ*1 is the time the patient has been infected before first being tested positive. (c) The individual mortality ratio *M*1(*t* |*τ*1) for *τ*inc = 6.4 days at different values of *τ*1. Note that the individual death probability *P*d(*t* |*τ*1) and *M*1(*t*| *τ*1) are nonzero only after *t > τ*inc− *τ*1. (d) The asymptotic individual mortality ratio *M*1(∞) (see Eq. (3)) as a function of *τ*1. Using these numbers, the recovery and death rate functions *c*(*τ, t*| *τ*1) and *µ*(*τ* |*τ*1) are plotted as functions of *τ* in Fig. 2(b). We show the evolution of *M*1(*t* |*τ*1) at different values of *τ*1 in Fig. 2(c). The corresponding long-time limit *M*1() is readily apparent in Fig. 2(d): for *τ*1 ≥ *τ*inc, *M*1(∞) = *µ*1*/*(*µ*1 +*c*) ≈ 0.105, while *M*1(∞) *< µ*1*/*(*µ*1 +*c*) when *τ*1 *< τ*inc. The smaller expected mortality associated with early identification of infection arises from the remaining incubation time during which the patient has a chance to recover without possibility of death. When conditioned on testing positive at or after the incubation period, the patient immediately suffers a positive death rate, increasing his *M*1(∞). Finally, in order to infer *M*1 (and also indirectly *µ* and *c*) during an outbreak, a number of statistical issues must be considered. First, if the outbreak is ongoing, there may not be sufficient long-time cohort data. Second, *τ*1 is unknown. Since testing typically occurs at the onset of symptoms, most positive patients will have been infected a few days earlier. The uncertainty in *τ*1 can be represented by a probability density *ρ*(*τ*1) for the individual. The expected mortality can then be constructed as an average over *ρ*(*τ*1): ![Formula][6] where ![Graphic][7] and ![Graphic][8] are the the *τ*1-averaged probabilities death and cure probabilities. Note that this averaging is different from the population-level averaging ![Graphic][9], which would describe the average of mortality ratios over a population with heterogeneous initial durations *τ*1. Some properties of the distribution *ρ*(*τ*1) can be inferred from the behavior of patients. Before symptoms arise, only very few patients will know they have been infected, seek medical care, and get their case confirmed (*i*.*e*., *ρ*(*τ*1) ≈ 0 for *τ*1 ≈ 0). The majority of patients will contact hospitals/doctors when they have been infected for a duration of *τ*inc. The distribution *ρ*(*τ*1) thus reaches its maximum near or shortly after *τ*inc. Since patients are most likely to test positive after experiencing symptoms, we choose a gamma distribution ![Formula][10] with shape parameter *n* = 8 and rate parameter *γ* = 1.25*/*day so that the mean *n/γ* is equal to *τ*inc = 6.4. Upon using the rates in Eqs. (4) and averaging over *ρ*(*τ*1), we derived expressions for ![Graphic][11], and ![Graphic][12] which are explicitly given in the SI. Using the values in Eq. (5) we find an expected individual mortality ratio ![Graphic][13] (which are subsequently plotted in Fig. 3) and its asymptotic value ![Graphic][14]. Of course, it is also possible to account for more complex time-dependent forms of *c* and *µ*1 [19], but we will primarily use Eqs. (4) in our subsequent analyses. ![FIG. 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.26.20044693/F3.medium.gif) [FIG. 3.](http://medrxiv.org/content/early/2020/03/30/2020.03.26.20044693/F3) FIG. 3. Population-level mortality-ratio estimates. Outbreak evolution and mortality ratios without containment measures (a,c) and with quarantine (b,d). The curves are based on numerical solutions of Eqs. (9) using the initial condition *I*(*τ*, 0) = *ρ*(*τ* ; 8, 1.25) (see Eq. (7)). The death and recovery rates are defined in Eqs. (4) and (5). We use a constant infection rate *β*1*S*(0) = 0.158/day, which we estimated from the basic reproduction number of SARS-CoV-2 [17]. To model quarantine effects, we set *β*1 = 0 for *t >* 50. We show the mortality-ratio estimates ![Graphic][15] and ![Graphic][16] (see Eq. (14)) and CFRd(*t, τ*res) (see Eqs. (8), (11), (12), and (14)). In the next section, we define population-based estimates for mortality ratios, *M*p(*t*), and explore how they can be computed using SIR-type models. By comparing *M*1(*t*) to *M*p(*t*), we gain insight into whether population-based metrics are good proxies for individual mortality ratios. We will outline the mathematical differences and additional errors that confound population-level estimates. ### B. Infection duration-dependent SIR model While individual mortalities can be estimated by tracking many individuals from infection to recovery or death, oftentimes, the available data are not resolved at the individual level and only total populations are given. Typically, one has the total number of cases accumulated up to time *t, N*(*t*), the number of deaths to date *D*(*t*), and the number of cured/recovered patients to date *R*(*t*) (see Fig. 1). The CFR is simply *D*(*t*)*/N*(*t*). Note that *N*(*t*) includes unresolved cases and that *N*(*t*) ≥ *R*(*t*)+ *D*(*t*). Resolution (death or recovery) of all patients, *N* (∞) = *R*(∞)+ *D*(∞), occurs only well after the epidemic passes. A variant of the CFR commonly used in the literature [1, 2] is the delayed CFR ![Formula][17] where *τ*res is a corresponding time lag that accounts for the duration from the day when first symptoms occurred to the day of cure/death. Many estimates of the COVID-19 mortality ratio assume that *τ*res = 0 [1, 2] and thus underestimate the number of death cases *D*(*t*) that result from a certain number of infected individuals. Similar underestimations using CFRd have been reported in previous epidemic outbreaks of SARS [4, 6] and Ebola [20]. Alternatively, a simple and interpretable population-level mortality ratio is *M*p(*t*) = *D*(*t*)*/*(*R*(*t*) + *D*(*t*)), the death ratio of all *resolved* cases. To provide a concrete model for *D*(*t*) and *R*(*t*), and hence *M*p(*t*), we will use a variant of the standard infection duration-dependent susceptible-infected-recovered (SIR)-type model described by [21] ![Formula][18] and ![Graphic][19], where *S*(*t*) is the number of susceptibles, *I*(*τ, t*) is density of individuals at time *t* who have been infected for time *τ*, and *R*(*t*) is the number of recovered individuals. The rate at which an individual infected for time *τ* at time *t* transmits the infection to a susceptible is denoted by *β*(*τ, t*)*S*(*t*). Note that the equation for *I*(*τ, t*) is identical to the equation for the survival probability described by Eq. (1). It is also equivalent to McKendrick age-structured models [22, 23] [24]. Infection of susceptibles is described by the boundary condition ![Formula][20] which is similar to that used in age-structured models to represent birth [22]. Finally, we use an initial condition consistent with the infection duration density given by Eq. (7): *I*(*τ*, 0) = *ρ*(*τ* ; *n* = 8, *γ* = 1.25). Note that Eq. (10) assumes that all newly infected individuals are immediately identified; *i*.*e*., these newly infected individuals start with *τ*1 = 0. After solving for the infected population density, the total number of deaths and recoveries to date can be found via ![Formula][21] The corresponding total number of cases *N*(*t*) in Eq. (8) is ![Formula][22] In the definitions of *D*(*t*), *R*(*t*), and *N*(*t*), we account for all possible death and recovery cases to date (see SI) and that newly infected individuals are immediately identified. We use these case numbers as approximations of the reported case numbers to study the evolution of mortality-ratio estimates. Mortality ratios based on these numbers underestimate the actual individual mortality *M*1 (see section II A) since they involve individuals that have been infected for different durations *τ*, particularly recently infected individuals who have not yet died. An alternative way to compute populations is to exclude the newly infecteds and consider only the initial cohort. The corresponding populations in this case are defined as ![Formula][23] Since *D*1(*t*) and *R*1(*t*) do not include infecteds with *τ < t*, they exclude the effect of newly infected individuals, but may yield more accurate mortality-ratios as they are based on an initial cohort of individuals in the distant past. The infections that occur after *t* = 0 contribute only to *I*(*τ < t, t*); thus, *D*1(*t*) and *R*1(*t*) do not depend on the transmission rate *β* or the number of susceptibles *S*(*t*). Note that all the populations derived above implicitly average over *ρ*(*τ*1; *n, γ*) for the first cohort of identified infecteds (but not subsequent infecteds). Moreover, the population density *I*(*τ t, t*) follows the same equation as ![Graphic][24] provided the same *ρ*(*τ*1; *n, γ*) is used in their respective calculations. The two different ways of partitioning populations (Eqs. (11) and (13)) lead to two different population-level mortality ratios ![Formula][25] Since the populations *D*(*t*) and *R*(*t*), and hence ![Graphic][26], depend on disease transmission through *β* (*τ*; *t*) and *S*(*t*), we expect ![Graphic][27] to carry a different interpretation from M1(*t*) and ![Graphic][28]. In the special case in which *µ* and *c* are constants, the time-integrated populations ![Graphic][29] and ![Graphic][30] factor out of ![Graphic][31] and ![Graphic][32], rendering them time-independent and ![Formula][33] Thus, only in the special time-homogeneous case do both population-based mortality ratios become *independent* of the population (and transmission *β*) and coincide with the individual death probability. To illustrate the differences between ![Graphic][34], and CFRd(*t, τ*res) in more general cases, we use the simple death and cure rate functions given by Eqs. (4) in solving Eqs. (1) and (9). For *β*(*τ, t*) in Eq. (10), we account for incubation effects by neglecting transmission during the asymptomatic incubation period (*τ* ≤ *τ*inc) and assume ![Formula][35] We use the estimated basic reproductive number ℛ = *β*1*S*(0)*/*(*µ*1 + *c*) ≈ 2.91 [17] to fix *β*1*S*(0) = (*µ*1 + *c*) ℛ ≈ 0.158*/*day. We also first assume that the susceptible population does not change appreciably before quarantine and set *S*(*t*) = *S*(0). Thus, we only need to solve for *I*(*τ, t*) in Eqs. (9) and (10). We solve Eqs. (9) and (10) numerically (see the *Methods* section for further details) and use these numerical solutions to compute *D*0,1(*t*), *R*0,1(*t*), and *N*0,1(*t*) (see Fig. 3(a) and (b)), which are then used in Eqs. (14) and CFRd(*t* − *τ*res). To determine a realistic value of the time lag *τ*res, we use data on death/recovery periods of 36 tracked patients [16] and find that patients recover/die, on average, *τ*res = 16.5 days after first symptoms occurred. We show in Figs. 3(c) and (d) that ![Graphic][36] approaches the individual mortality ratio ![Graphic][37] of section II A. This occurs because the model for *P*(*τ, t*) and *I*(*τ, t*) are equivalent and we assumed that the initial distribution of *τ* for both quantities are given by *ρ*(*τ* ; 8, 1.25). However, the population-level mortality ratios CFRd(*t, τ*res) and ![Graphic][38] also take into account recently infected individuals who may recover before symptoms. This difference yields different mortality ratios because newly infecteds are implicitly assumed to be detected immediately and all have *τ*1 = 0. Thus, the underlying infection-time distribution is not the same as that used to compute ![Graphic][39] (see SI for further details). The mortality ratios CFRd(*t, τ*res) and ![Graphic][40] should not be used to quantify the individual mortality probability of individuals who tested positive after their incubation period. During the course of an outbreak, the measures CFRd(*t, τ*res) and ![Graphic][41] are subject to another confounding influence. Since *D*(*t*), *R*(*t*), and *N*(*t*) do not change with the same rates at the same time, these population-level mortality estimates only reach their steady state after sufficiently long times (see Fig. 3(c) and (d)). To summarize, we described two confounding factors that complicate the direct use of population-level mortality ratio to estimate individual mortality probabilities. First, infection-time distributions *ρ*(*τ* ; *n, γ*) that are meaningful on an individual level may not correspond to those in population-level data. Second, population-level mortality ratios are often time-dependent and most informative only in the steady state after the outbreak stopped. The evolution of the mortality ratios in Fig. 3 qualitatively resembles the behavior of the mortality-ratio estimates in Fig. 1. As shown in Fig. 1, the population-based estimates for coronavirus varies, decreasing in time for China but fluctuating for Italy. These changes could result from changing practices in data collecting, or from explicitly time-inhomogeneous parameters *µ*(*τ, t*), *c*(*τ, t*), and/or *β*(*τ, t*). Although population-level quarantining does not directly affect the individual mortality *M*1(*t* | *τ*1) or ![Graphic][42], it can be easily incorporated into the SIR-type population dynamics equations through changes in *β*(*τ, t*)*S*(*t*). For example, we have set *S*(*t > t*q) = 0 to represent implementation of a quarantine after *t*q = 50 days of the outbreak. After *t*q = 50 days, no new infections occur and the estimates CFRd(*t, τ*res) and ![Graphic][43] start converging immediately towards their steady-state values (see Fig. 3(d)). Since the number of deaths decreases after the implementation of quarantine measures, the delayed CFRd(*t, τ*res = 17) is first decreasing until *t* = *t*q + *τ*res = 67. For *t >* 67, the CFRd(*t, τ*res = 17) measures no new cases and is thus equal to the CFR. ## III. DISCUSSION AND SUMMARY After an outbreak, it is important to assess the severity of the disease by estimating its mortality and other disease characteristics. Assuming accurate data, the often-used CFR and delayed CFR typically underestimate the true, final death ratio. For example, during the SARS outbreaks in Hong Kong, the WHO first estimated the fatality rate to 2.5% (March 30, 2003) whereas the final estimates reached values of about 17.0% (June 30, 2003) [7]. Standard metrics like the CFR are seen to be easily confounded by and sensitive to uncertainty in intrinsic disease parameters such as the incubation period and the time *τ*1 a patient had been infected before clinical confirmation of infection. For the recent COVID-19 outbreaks, CFR-based measures may still provide reasonable estimates of the actual mortality across different age classes due to a counter-acting error in the numbers of unreported mild-symptom cases. Here, we stress that more mechanistically meaningful and interpretable metrics can be defined and be as easily estimated from data as CFRs. Our proposed mortality ratios for viral epidemics are defined in terms of (i) individual survival probabilities and (ii) population ratios using numbers of deaths and recovered individuals. Both of these measures are based on the within-host evolution of the disease, and in the case of ![Graphic][44], the population-level transmission dynamics. Thus, these metrics directly incorporate key parameters operating on the weeks or months timescale, the incubation time *τ*inc and time of prior infection *τ*, through the solution of age-structured PDEs. Among the metrics we describe, ![Graphic][45] is structurally closest to ![Graphic][46] in that both are independent of transmission *β* since new infections are not considered. Both of these converge after an incubation time *τ*inc to a value smaller than or equal to *µ*1*/*(*µ*1 + *c*). The most accurate estimates of *M*1 can be obtained if we keep track of the fate of cohorts that were infected within a small time window in the past. By following only these individuals, one can track how many of them died as a function of time. As more cases arise, one should stratify them according to estimated *τ* to gather improving statistics for *M*1(∞). These data should also be collated according to the other central factor in COVID-19 mortality: patient age. With the further spread of SARS-CoV-2 in different countries, data on more individual cases of death and recovery can be more easily stratified by age, health condition, and other individual characteristics. Using identical initial infection time distributions *ρ*(*τ*1; *n, γ*) (see Eq. (7)), the long-time limit of ![Graphic][47] approaches the individual mortality ![Graphic][48] (see Eq. (6)). Besides accurate cohort data, for which at present there are few for coronavirus, cumulative population data has been used to estimate the mortality ratio. The metrics ![Graphic][49] and CFR(*t*) are based on these aggregate populations but implicitly depend on new infections and the transmission rate *β*. Despite this confounding factor, ![Graphic][50] and CFRd(*t, τ*res) approach ![Graphic][51] as *t* → ∞, where ![Graphic][52] is the probability that no recovery occurred during the incubation time *τ*inc. Based on these results, we can establish the following connection between the different mortality ratios for initial infection times with distribution *ρ*(*τ*1; *n, γ*) and mean *τ* = *n/γ*: ![Formula][53] According to Eq. (17), population-level mortality estimates (*e*.*g*., CFR and ![Graphic][54]) can be transformed, at least approximately, into individual mortality probabilities using the correction factor ![Graphic][55] with ![Graphic][56]. Besides the mathematical differences between *M*1(*t*) and ![Graphic][57], CFR, estimating ![Graphic][58] and CFR(*t*) from aggregate populations implicitly incorporate a number of confounding factors that lead to variability in these estimates. In Fig. 4, we plot the population-level mortality-ratio estimates ![Graphic][59] against the CFR for different regions and observe large variations and very little correlation between countries [25]. As of March 26, 2020, the value of ![Graphic][60] in Italy is almost 45% and can increase further if the current conditions (e.g., treatment methods, age group proportion of infecteds, etc.) do not change. Differences between the mortality ratios in China and Italy (see Figs. 1(b) and (c)) might be a result of varying medical treatment strategies, different practices in data collecting (e.g., post-mortem testing), and differences in the age demographics between the countries. ![FIG. 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.26.20044693/F4.medium.gif) [FIG. 4.](http://medrxiv.org/content/early/2020/03/30/2020.03.26.20044693/F4) FIG. 4. Region-dependence of COVID-19 mortality-ratio estimates. Mortality-ratio estimates of COVID-19 in different regions (see Eqs. (8) and 14 (*τ*res = 0)). We used data on the cumulative number of cases, recoveries, and deaths in Ref. [5] as of March 24, 2020. The marker sizes indicate the population of the corresponding countries. The metrics ![Graphic][61] and and CFR are largely uncorrelated with correlation coefficient 0.33. In general, even if the cohort initially tested was only a fraction of the total infected population, tracking ![Graphic][62] or ![Graphic][63] of this cohort still provides an accurate estimation of the mortality rate. However, the newly infecteds that contribute to CFR and ![Graphic][64] at later times may not all be tested or may be tested at different times after they were infected. A reported/tested fraction *f <* 1 *would not* directly affect the CFRs or mortality ratios if the unreported/untested individuals die and recover in the same proportion as the tested infecteds. Undertesting will overestimate true CFR or mortality rates if the untested infecteds are less likely to die than the tested infecteds. In other words, if the untested (presumably because they were mildly or asymptomatic) population predominantly recovers instead of dying, the actual CFR and mortality ratios would be significantly lower than those based on tested individuals. If untested infecteds do not die, the asymptotic mortality of all infected individuals ![Graphic][65] (see the SI). Current estimates show that only a minority of SARS-CoV-2 infections are reported (e.g., *f* 14% in China before January 23, 2020) [26]. Besides under-reporting, the delay in transmission after becoming infected will also affect ![Graphic][66]. Although we have assumed that transmission occurs only after the incubation period when symptoms arise, there is evidence of asymptomatic transmission of coronavirus [26, 27]. Asymptomatic transmission can be modeled by setting *β*(*τ*) *>* 0 even for *τ < τ*inc. An undelayed transmission in a nonquarantine scenario causes relatively more new infecteds who have not had the chance to die yet, leading to a *smaller* mortality ratio ![Graphic][67]. Within our SIR model, delaying transmission reduces the number of infected individuals and deaths at any given time but *increases* the measured mortality ratio ![Graphic][68]. Without quarantine, the asymptotic values ![Graphic][69] and CFR(∞) will also change as a result of changing the transmission latency period, as shown in the SI. With perfect quarantining, the asymptote ![Graphic][70] is eventually determined by a cohort that does not include new infections and is thus independent of the transmission delay. In this work, we have explicitly defined a number of interpretable mathematical metrics that represent the risk of death. By rigorously defining these metrics, we are able to reveal the inherent assumptions and factors that affect their estimation. Within survival probability and SIR-type models, we explicitly illustrate how physiologically important parameters such as incubation time, death rate, cure rate, and transmissibility influence the metrics. We also discussed how statistical factors such as time of testing after infection (*τ*1) and testing ratio (*f*) affect our estimates. Given the uncertainty in the testing fraction, we conclude that *M*1(*t*) and ![Graphic][71] is best interpreted as approximately the mortality probability conditioned on being tested positive. In practice, these are probably also good estimates of mortality of patients conditioned on showing symptoms. In addition to our metrics and mathematical models, we emphasize the importance of curating individual cohort data. These data are more directly related to the probability of death *M*1(*t*) and are subject to the fewest confounding factors and statistical uncertainty. ## METHODS To numerically solve Eqs. (9) and (10), we used a uniform discretization *τ**k* = *k*Δ*τ, k* = 0, 1, …, *K*. A backward difference operator [*I*(*τ**k*, *t*) − *I*(*τ**k*−1, *t*)] */*(Δ*τ*) is used to approximate ∂*τ* *I*(*τ, t*) and a predictor-corrector Euler scheme is used to advance time [28]. Setting the cut-off *I*(Δ*τ, t*) ≡ 0 and *I*(*K*Δ*τ, t*) ≡ 0, the resulting discretized equations for the full SIR model are ![Formula][72] where ![Graphic][73] is the initial predicted guess, and the last term proportional to *δ**k*,0 encodes the boundary condition Eq. (10). Note that we use ![Graphic][74] to indicate the numerical evaluation of ![Graphic][75]. Quadrature methods such as Simpson’s rule and the trapezoidal rule can be used to approximate the integral more efficiently. The total deaths, recovereds, and infecteds at time *t* are found by ![Formula][76] with analogous expressions for *D*1(*m*Δ*t*) and *R*1(*m*Δ*t*). To obtain a stable integration scheme, the time steps Δ*t* and Δ*τ* have to satisfy Δ*t*/(2Δ*τ*) < 1. In all of our numerical computations, we thus set Δ*t* = 0.002, Δ*τ* = 0.02, and *K* = 104. In the SI, we show additional plots of the magnitude of *I*(*τ, t*) in the *t* − *τ* plane. ## Data Availability The datasets that we used in this study are stored in the publicly accessible repositories of Refs. [3, 5, 16]. ## DATA AVAILABILITY The datasets that we used in this study are stored in the publicly accessible repositories of Refs. [3, 5, 16]. ## COMPETING INTERESTS The authors declare no competing interests. ## AUTHOR CONTRIBUTIONS LB and TC developed the analyses and wrote the manuscript. LB analyzed data and MX performed the numerical computations. ## SUPPLEMENTARY INFORMATION ### A. Additional examples of mortality-ratio evolutions ![FIG. S1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.26.20044693/F5.medium.gif) [FIG. S1.](http://medrxiv.org/content/early/2020/03/30/2020.03.26.20044693/F5) FIG. S1. Mortality ratio estimates. Estimates of mortality ratios (see Eqs. (8) and (14) in the main text) of SARS-CoV-2 infections in different countries. The case fatality ratio, CFR, corresponds to the number of deaths to date divided by the total number of cases to date. Another population-based mortality ratio is ![Graphic][77], the number of deaths divided by the sum of deaths and recovereds, up to time *t*. The data are derived from Ref. [5]. In Fig. S1, we show additional examples of mortality-ratio estimates for Iran, South Korea, Spain, Germany, Switzerland, and the United Kingdom. As in Fig. 1 in the main text, we observe that, by definition, the population-based mortality ratio ![Graphic][78] is significantly larger than the corresponding CFR in all cases. ### B. Solutions for *τ*1-averaged probabilities Using the method of characteristics, we find the formal solution to Eq. (1): ![Formula][79] which can be used to construct the death and cure probabilities ![Formula][80] If we now invoke the functional forms of *µ* and *c* given in Eq. (4), we find explicitly ![Formula][81] and ![Formula][82] ![FIG. S2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.26.20044693/F6.medium.gif) [FIG. S2.](http://medrxiv.org/content/early/2020/03/30/2020.03.26.20044693/F6) FIG. S2. Phase plot for *P*(*τ > t, t*) and *I*(*τ > t, t*). The regions delineating different forms for the solution (Eq. (S5)). Here, we have included an incubation time *τ*inc before which no death occurs. The solution for ![Graphic][83] or *I*(*τ, t*) in the *τ < t* region must be self-consistently solved using the boundary condition Eq. (10). At any fixed time, the integral of *I*(*τ, t*) over *t < τ* ≤ ∞ captures only the initial population, excludes newly infecteds, and is used to compute *D*1(*t*), *R*1(*t*), and ![Graphic][84]. To compute *D*(*t*), *R*(*t*), and ![Graphic][85], we integrate across all infecteds (including the integral over *t > τ >* ≥ 0 shown in magenta). Finally, we can also find the *τ*1-averaged probabilities for *τ* ≥ *t* by weighting over *ρ*(*τ*1; *n, γ*). For example, ![Formula][86] ![FIG. S3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.26.20044693/F7.medium.gif) [FIG. S3.](http://medrxiv.org/content/early/2020/03/30/2020.03.26.20044693/F7) FIG. S3. Density plots of *I*(*τ, t*) in the *t*− *τ* plane. Numerical solution of the equation for *I*(*τ, t*) in Eqs. (9) under the assumption of a fixed susceptible size *β*1*S* = 0.158/day. (a) The density without quarantine monotonically grows with time *t* in the region *τ < t* as an unlimited number of susceptibles continually produces infecteds. (b) With quarantining after *t*q = 50 days, we set *β*1*S* = 0 for *t > t*q, which shuts off new infections. Both plots were generated using the same initial density *ρ*(*τ*1) defined in Eq. (7). In both cases, the density *I*(*τ > t*) is identical to *P*(*τ > t*) if the same *ρ*(*τ*1) is used and is independent of disease transmission, susceptible dynamics, etc. (c-d) Probability-density functions (PDFs) of the number of infected *I*(*τ, t*) for *t* = 0, 60 (b) without and (c) with quarantine. The blue solid line corresponds to the initial distribution *ρ*(*τ* ; *n* = 8, *γ* = 1.25) (see Eq. (7)). These solutions hold for the different regions shown in the phase plot of Fig. S2 and are equivalent to those for *I*(*τ > t, t*). Corresponding expressions for ![Graphic][87] and ![Graphic][88] can be found and used to construct ![Graphic][89]. Fig. S3(a) shows the magnitude of *I*(*τ, t*) in the *t* − *τ* plane when we set *S*(*t*) = *S* constant (so that the first equation in Eq. (18) does not apply) such that *β*1*SS* ≈ 0.158/day. In this case, the epidemic continues to grow in time, but the mortality rates ![Graphic][90] nonetheless converge as *t*. In Fig. S3(b), we set *β*1*S* = 0 for *t > t*q to model strict quarantining after *t*q = 50 days. We observe no new infections after the onset of strict quarantine measures. In both cases (quarantine and no quarantine), we use *ρ*(*τ* ; *n* = 8, *γ* = 1.25) (see Eq. (7) in the main text) to describe the initial distribution of infection times *τ*. As time progresses, more of the distribution of *τ* moves towards smaller values until quarantine measures take effect (see Fig. S3(c) and (d)). ### C. Effects of undertesting Note that *I*(*τ, t*) in the SIR equations determines the dynamics of the actual infected population. However, (i) typically only a fraction *f* of the total number of infecteds might be tested and confirmed positive and (ii) the testing of newly infecteds may also be delayed by a distribution *ρ*(*τ* ; *n, γ*). If positive tests represent only a fraction *f* of the total infected population, and the confirmation of newly infecteds occurs immediately, the known infected density is given by *I**∗*(*τ, t*) = *f I*(*τ, t*) where *I*(*τ, t*) is the true total infected population. If testing of newly infecteds occurs after a distribution *ρ*(*τ* ; *n, γ*) of infection times, ![Graphic][91] In our development of ![Graphic][92] and CFRd(*t, τ*res) in the manuscript, we assumed the entire infected population was tested and confirmed. Thus, ![Graphic][93] and CFRd(*t, τ*res) were computed using *f* = 1 and more accurately represent the mortality ratios of the population *conditioned* on being tested positive. ![FIG. S4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.26.20044693/F8.medium.gif) [FIG. S4.](http://medrxiv.org/content/early/2020/03/30/2020.03.26.20044693/F8) FIG. S4. Fractional testing. An example of fractional testing in which a fixed fraction *f* of the real total infected population is assumed to be tested. The remaining 1 − *f* proportion of infecteds are untested. Equivalently, if the total tested fraction has unit population, then the total population of the untested pool is 1*/f*− 1. (a) At short times after an outbreak, the known tested infected population has not yet resolved and is composed of deaths (gray), recovereds (green), and infecteds (red). We assume that the untested fraction of infecteds (red) have mild or no symptoms, do not die, and can only recover (green). (b) At longer times, the infecteds further resolve. The *M*p(*t*) and CFR metrics that are based on only the tested fraction will overestimate the true mortality fraction of all infected cases. To estimate the mortality ratio of the population conditioned simply on being infected, we have to estimate the larger number of recovereds that went untested. As shown in Fig. S4, the untested recovered fraction can be estimated by assuming that the death rate for the untested infecteds is zero and by writing an SIR model without death for the untested pool of infected ![Formula][94] where *I*(*τ, t*) = *I**∗*(*τ, t*) + *I*u(*τ, t*). The true mortality ratio is then straightforwardly defined by, for example, ![Formula][95] Where ![Formula][96] with analogous expressions for ![Graphic][97], and ![Graphic][98]. At long times, after resolution of all infecteds, the untested recovered population is ![Formula][99] which yields the asymptotic true ratio ![Graphic][100] as described in the Discussion and Summary. In this simple rescaling to account for untested populations, we have assumed that all deaths come from the tested pool and that the recovery rate *c* is the same in the tested and untested pools. ### D. Influence of different transmission rates In Fig. 3 of the main text, we observe that the population-level mortality ratio ![Graphic][101] approaches a plateau during the initial exponential growth phase of an epidemic (*i*.*e*., for *S*(*t*) ≈ *S*). If the number of new infections decreases (*e*.*g*., due to quarantine measures), ![Graphic][102] starts growing until it reaches its asymptotic value ![Graphic][103] Interestingly, the pre-asymptotic values of ![Graphic][104] are smaller for larger infection rates *β*1 (see Fig. S5(a)). This counter-intuitive effect arises because larger values of *β*1 generate relatively larger numbers of new infected which have a lower chance of dying before *τ*inc (see Eq. (4) in the main text). A similar effect occurs for non-delayed transmission (*i*.*e*., *τ**β* ≈ 0). ![FIG. S5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.26.20044693/F9.medium.gif) [FIG. S5.](http://medrxiv.org/content/early/2020/03/30/2020.03.26.20044693/F9) FIG. S5. Population-level mortality for different infection rates. (a) The population-level mortality ratio ![Graphic][105] for different values of *β*1 and an incubation time of *τ*inc = 6.4 days. In the initial exponential growth phase of the epidemic (*i*.*e*., *S*(*t*) ≈ *S*), larger infection rates *β*1 lead to smaller values of ![Graphic][106]. (b) We observe a similar effect for non-delayed transmissions (*i*.*e*., *τ**β* ≈ 0). As long as *S*(*t*) ≈ *S*, smaller transmission delays *τ**β* lead to larger relative numbers of new infections and smaller ![Graphic][107]. As the transmission delay decreases, more secondary cases will result from one infection, leading to smaller values of ![Graphic][108] in the initial exponential growth phase of an epidemic (see Fig. S5(b)). ## ACKNOWLEDGEMENTS LB acknowledges financial support from the SNF Early Postdoc.Mobility fellowship on “Multispecies interacting stochastic systems in biology”. The authors also acknowledge financial support from the Army Research Office (W911NF-18-1-0345), the NIH (R01HL146552), and the National Science Foundation (DMS-1814364). ## Footnotes * * lucasb{at}ethz.ch * † xiamingtao97{at}g.ucla.edu * ‡ tomchou{at}ucla.edu * Received March 26, 2020. * Revision received March 26, 2020. * Accepted March 30, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1]. Z. Xu, L. Shi, Y. Wang, J. Zhang, L. Huang, C. Zhang, S. Liu, P. Zhao, H. Liu, L. Zhu, et al., Lancet Resp. Med. (2020). 2. [2]. Z. Wu and J. M. McGoogan, JAMA. 3. [3].“COVID-19 statistics,” [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/), accessed: 2020-02-26. 4. [4]. A. Ghani, C. Donnelly, D. Cox, J. Griffin, C. Fraser, T. Lam, L. Ho, W. Chan, R. Anderson, A. Hedley, et al., Am. J. Epidemiology 162, 479 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwi230&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16076827&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.26.20044693.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231363800012&link_type=ISI) 5. [5]. E. Dong, H. Du, and L. Gardner, Lancet Inf. Dis. (2020). 6. [6]. P. S. Yip, K. Lam, E. H. Lau, P.-H. Chau, K. W. Tsang, and A. Chao, J. Roy. Stat. Soc. A 168, 233 (2005). 7. [7]. P. S. Yip, E. H. Lau, K. Lam, and R. M. Huggins, Am. J. Epidemiol. 161, 700 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwi088&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15781959&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.26.20044693.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000227954600011&link_type=ISI) 8. [8]. L. Böttcher, O. Woolley-Meza, N. A. Araújo, H. J. Herrmann, and D. Helbing, Sci. Rep. 5, 1 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/srep08026&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26082330&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.26.20044693.atom) 9. [9]. L. Böttcher, O. Woolley-Meza, E. Goles, D. Helbing, and H. J. Herrmann, Phys. Rev. E 93, 042315 (2016). 10. [10]. J. Riou, A. Hauser, M. J. Counotte, and C. L. Althaus, medRxiv (2020), doi:10.1101/2020.03.04.20031104. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMy4wNC4yMDAzMTEwNHYzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMzAvMjAyMC4wMy4yNi4yMDA0NDY5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 11. [11]. Z. Xu, S. Li, S. Tian, H. Li, and L.-q. Kong, The Lancet (2020). 12. [12]. E. Mahase, “Coronavirus: COVID-19 has killed more people than SARS and MERS combined, despite lower case fatality rate,” (2020). 13. [13].W. H. Organization, “WHO Director-General’s opening remarks at the media briefing on COVID-19 - 24 February 2020,” [https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19\---|24-february-2020](https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19\---|24-february-2020) (2020), accessed: 2020-02-28. 14. [14].W. H. Organization, “WHO Director-General’s opening remarks at the media briefing on COVID-19 - 3 March 2020,” [https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19\---|3-march-2020](https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19\---|3-march-2020) (2020), accessed: 2020-03-05. 15. [15].The timescale of infections are much smaller than the aging process. In other words, the aging process across the disease time scale is negligible and the age-dependent transport terms ∂P/∂a can be neglected. 16. [16].nCoV 2019DataWorking Group, “Epidemiological Data fromthe nCoV-2019 Out-break:Early Descriptionsfrom Publicly Available Data,” [http://virological.org/t/epidemiological-data-from-the-ncov-2019-outbreak-early-descriptions-from-publicly-available-data/337](http://virological.org/t/epidemiological-data-from-the-ncov-2019-outbreak-early-descriptions-from-publicly-available-data/337) (2020), waccessed: 2020-02-26. 17. [17]. C.-C. Lai, T.-P. Shih, W.-C. Ko, H.-J. Tang, and P.-R. Hsueh, Int. J. Antimicrob. Agents, 105924 (2020). 18. [18]. M. J. Keeling and P. Rohani, Modeling infectious diseases in humans and animals (Princeton University Press, 2011). 19. [19]. L. Böttcher and N. Antulov-Fantulin, arXiv preprint arxiv:2002.11765 (2020). 20. [20]. K. E. Atkins, N. S. Wenzel, M. Ndeffo-Mbah, F. L. Altice, J. P. Townsend, and A. P. Galvani, BMJ 350, h1115 (2015). [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNTAvbWFyMTZfMy9oMTExNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAzLzMwLzIwMjAuMDMuMjYuMjAwNDQ2OTMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 21. [21]. J. T. Wu, K. Leung, M. Bushman, N. Kishore, R. Niehus, P. M. de Salazar, B. J. Cowling, M. Lipsitch, and G. M. Leung, Nature Medicine (2020), doi:10.1038/s41591-020-0822-7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0822-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32284616&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.26.20044693.atom) 22. [22]. A. G. McKendrick, Proc. Edinburgh Math. Soc. 44, 98 (1926). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0013091500034428&link_type=DOI) 23. [23]. T. Chou and C. D. Greenman, Journal of Statistical Physics 164, 49 (2016. PMCID: PMC3894939). 24. [24].In both the individual model (Eq. (1)) and population models (Eq. (9)), the death and cure rates are insensitive to changes in age a over the time scale of epidemic if it occurs over, say, only 1-2 years. In this limit, we consider only infection duration dependence in the population dynamics. 25. [25]. M. Battegay, R. Kuehl, S. Tschudin-Sutter, H. H. Hirsch, A. F. Widmer, and R. A. Neher, Swiss Medical Weekly 150 (2020). 26. [26]. R. Li, S. Pei, B. Chen, Y. Song, T. Zhang, W. Yang, and J. Shaman, Science (2020). 27. [27]. C. Rothe, M. Schunk, P. Sothmann, G. Bretzel, G. Froeschl, C. Wallrauch, T. Zimmer, V. Thiel, C. Janke, W. Guggemos, M. Seilmaier, C. Drosten, P. Vollmar, K. Zwirglmaier, S. Zange, R. Wlfel, and M. Hoelscher, New England Journal of Medicine 382, 970 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2001468&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.26.20044693.atom) 28. [28]. W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes: The Art of Scientific Computing (Cambridge University Press, 2007). [1]: /embed/graphic-3.gif [2]: /embed/graphic-4.gif [3]: /embed/graphic-5.gif [4]: /embed/graphic-6.gif [5]: /embed/graphic-7.gif [6]: /embed/graphic-9.gif [7]: /embed/inline-graphic-1.gif [8]: /embed/inline-graphic-2.gif [9]: /embed/inline-graphic-3.gif [10]: /embed/graphic-10.gif [11]: /embed/inline-graphic-4.gif [12]: /embed/inline-graphic-5.gif [13]: /embed/inline-graphic-6.gif [14]: /embed/inline-graphic-7.gif [15]: F3/embed/inline-graphic-8.gif [16]: F3/embed/inline-graphic-9.gif [17]: /embed/graphic-12.gif [18]: /embed/graphic-13.gif [19]: /embed/inline-graphic-10.gif [20]: /embed/graphic-14.gif [21]: /embed/graphic-15.gif [22]: /embed/graphic-16.gif [23]: /embed/graphic-17.gif [24]: /embed/inline-graphic-11.gif [25]: /embed/graphic-18.gif [26]: /embed/inline-graphic-12.gif [27]: /embed/inline-graphic-13.gif [28]: /embed/inline-graphic-14.gif [29]: /embed/inline-graphic-15.gif [30]: /embed/inline-graphic-16.gif [31]: /embed/inline-graphic-17.gif [32]: /embed/inline-graphic-18.gif [33]: /embed/graphic-19.gif [34]: /embed/inline-graphic-19.gif [35]: /embed/graphic-20.gif [36]: /embed/inline-graphic-20.gif [37]: /embed/inline-graphic-21.gif [38]: /embed/inline-graphic-22.gif [39]: /embed/inline-graphic-23.gif [40]: /embed/inline-graphic-24.gif [41]: /embed/inline-graphic-25.gif [42]: /embed/inline-graphic-26.gif [43]: /embed/inline-graphic-27.gif [44]: /embed/inline-graphic-28.gif [45]: /embed/inline-graphic-29.gif [46]: /embed/inline-graphic-30.gif [47]: /embed/inline-graphic-31.gif [48]: /embed/inline-graphic-32.gif [49]: /embed/inline-graphic-33.gif [50]: /embed/inline-graphic-34.gif [51]: /embed/inline-graphic-35.gif [52]: /embed/inline-graphic-36.gif [53]: /embed/graphic-21.gif [54]: /embed/inline-graphic-37.gif [55]: /embed/inline-graphic-38.gif [56]: /embed/inline-graphic-39.gif [57]: /embed/inline-graphic-40.gif [58]: /embed/inline-graphic-41.gif [59]: /embed/inline-graphic-42.gif [60]: /embed/inline-graphic-43.gif [61]: F4/embed/inline-graphic-44.gif [62]: /embed/inline-graphic-45.gif [63]: /embed/inline-graphic-46.gif [64]: /embed/inline-graphic-47.gif [65]: /embed/inline-graphic-48.gif [66]: /embed/inline-graphic-49.gif [67]: /embed/inline-graphic-50.gif [68]: /embed/inline-graphic-51.gif [69]: /embed/inline-graphic-52.gif [70]: /embed/inline-graphic-53.gif [71]: /embed/inline-graphic-54.gif [72]: /embed/graphic-23.gif [73]: /embed/inline-graphic-55.gif [74]: /embed/inline-graphic-56.gif [75]: /embed/inline-graphic-57.gif [76]: /embed/graphic-24.gif [77]: F5/embed/inline-graphic-58.gif [78]: /embed/inline-graphic-59.gif [79]: /embed/graphic-26.gif [80]: /embed/graphic-27.gif [81]: /embed/graphic-28.gif [82]: /embed/graphic-29.gif [83]: F6/embed/inline-graphic-60.gif [84]: F6/embed/inline-graphic-61.gif [85]: F6/embed/inline-graphic-62.gif [86]: /embed/graphic-31.gif [87]: /embed/inline-graphic-63.gif [88]: /embed/inline-graphic-64.gif [89]: /embed/inline-graphic-65.gif [90]: /embed/inline-graphic-66.gif [91]: /embed/inline-graphic-67.gif [92]: /embed/inline-graphic-68.gif [93]: /embed/inline-graphic-69.gif [94]: /embed/graphic-34.gif [95]: /embed/graphic-35.gif [96]: /embed/graphic-36.gif [97]: /embed/inline-graphic-70.gif [98]: /embed/inline-graphic-71.gif [99]: /embed/graphic-37.gif [100]: /embed/inline-graphic-72.gif [101]: /embed/inline-graphic-73.gif [102]: /embed/inline-graphic-74.gif [103]: /embed/inline-graphic-75.gif [104]: /embed/inline-graphic-76.gif [105]: F9/embed/inline-graphic-77.gif [106]: F9/embed/inline-graphic-78.gif [107]: F9/embed/inline-graphic-79.gif [108]: /embed/inline-graphic-80.gif