The incidence-based dynamic reproduction index: accurate determination, diagnostic sensitivity, and predictive power ==================================================================================================================== * Robert N. J. Conradt * Stephan Herminghaus ## Abstract Two methods of calculating the reproduction index from daily new infection data are considered, one by using the generation time *t**G* as a shift (R*G*), and an incidence-based method directly derived from the differential equation system of an SIR epidemic dynamics model (R*I*). While the former is shown to have few in common with the true reproduction index, we find that the latter provides a sensitive detection device for intervention effects and other events affecting the epidemic, making it well-suited for diagnostic purposes in policy making. Furthermore, we introduce a similar quantity, ![Graphic][1], which can be calculated directly from R*G*. It shows largely the same behaviour as R*I*, with less fine structure. However, it is accurate in particular in the vicinity of R = 1, where accuracy is important for the corrrect prediction of epidemic dynamics. We introduce an entirely new, self-consistent method to derive, from both quantities, an improved ![Graphic][2] which is both accurate and contains the details of the epidemic spreading dynamics. Hence we obtain R accurately from data on daily new infections (incidence) alone. Moreover, by using R*I* instead of R*G* in plots of R versus incidence, orbital trajectories of epidemic waves become visible in a particularly insightful way, demonstrating that the widespread use of only incidence as a diagniostic tool is clearly inappropriate. PACS numbers: ## I. INTRODUCTION The outbreak of the illness COVID-19, caused by the SARS-CoV-2 virus, has resulted in a pandemic with unprecedented impact on societies all over Earth. Mitigation measures included complete lockdowns of societal life, with severe social, economic, and individual consequences [1, 2]. The dramatically varying success [3–6] of the interventions owed in part to cultural differences [3], but also to only limited understanding of infection spreading dynamics and a severe lack of established methods in epidemic state diagnosis and prediction. Improvement of this general situation is imperative, in particular as similar events are expected to strike more often in the future [7–9]. In search for optimized strategies, a two-fold view must be adopted. One the one hand, one needs to understand, in retrospect, which interventions have had what effect on the epidemic spreading dynamics, in order to properly design future interventions. This requires sensitive diagnosis tools for assessing the state of the epidemic on a (if possible) daily basis. On the other hand, tools are needed for predicting the future of epidemic dynamics as reliably as possible if conditions are known. Aside from extensive simulation, this requires careful analysis of data, such as the number of infected citizens [2, 10]. Here we discuss the system in terms of an SIR model [11, 12], referring to the number of susceptible (S), infected (I), and recovered (R) individuals, respectively, in a population of *N* citizens. Here we identify with *R* all those who are neither susceptible nor infected (*R* = *N − S − I*), which includes those who are deceased. While the ratio of deceased vs. recovered individuals is of greatest societal concern, it can be disregarded here, as we will solely discuss prevention measures addressing the spreading of the disease. We define *I* as the number of individuals who carry sufficient viral load to be contagious. They are assumed to remain in this state for an average duration *τ*. Although the viral load changes with time during the illness, contagion can be sufficiently well described by this simple picture for our purposes [13]. The spreading dynamics of an epidemic can then be described by a set of two equations [12], ![Formula][3] where *∂**t* is the derivative with respect to time, while *s*(*t*) = *S/N* and *i*(*t*) = *I/N* are the fractions of susceptible and infected individuals in the population, respectively. The constant *c* is the average number of new infections a single infected individual would cause per unit time in an otherwise infection-free (but susceptible) population. It is accessible to interventions such as closing schools, wearing facial masks etc., but this shall not concern us here, as we focus solely on methods to detect the current stage of an epidemic and to predict its near future development, at given *c*. ## II. THE REPRODUCTION INDEX The base reproduction index, R, is related to *c* via ![Formula][4] It denotes the total average number of individuals newly infected by a single infected one under the above conditions. Since the probability of infection is directly proportional to the fraction of susceptibles, we have ![Formula][5] for the dynamic reproduction index, R(*t*) [12]. The latter is of tremendous importance for assessing the current status of an epidemic. If R = 1, *I* stays constant, but when R *>* 1, each infected individual causes more than one new infection on average, such that *I*(*t*) increases exponentially. It is thus of major interest to determine R from epidemiological data as accurately as possible, in particular in the vicinity of unity. By combining eq. (3) with the first eq. (1), we obtain ![Formula][6] This relies only on quantities which can be derived from data usually available from the health care system. *∂**t**S* represents the number of new infections per unit time (daily incidence) and can be considered known accurately. *τ* is known from clinical experience with the disease, and *I*(*t*) (and hence *i*) can be estimated once *τ* is known. We will see below that it is in particular the dynamic variations of R which yield considerable insight into the infection process. Since infection data are discrete data collected on a daily basis [14], we will now write down a discrete version of eq. (4). The daily incidence will be called ![Graphic][7]. *τ* as well as *t* will henceforth be expressed in units of days, and treated as discrete variables. We then may be tempted to simply write ![Graphic][8]. However, we must be aware that for rather general infrastructural reasons, the reporting efficiency of infection numbers varies characteristically, e.g., on weekends. We therefore should provide for suitable averaging. Hence we write ![Formula][9] where *l* and *k* are parameters determining the intervals over which ![Graphic][10] and *I* are being averaged, respectively. Since the typical variability of data reflects the sequence of seven weekdays, it appears reasonable to average ![Graphic][11] over seven consecutive days. If we furthermore average the (less variable) number of infected individuals over a period *k* = *τ*, we obtain ![Formula][12] for the incidence-based dynamic reproduction index. While ![Graphic][13] is known precisely, *I* can only be estimated based, among others, on *τ*. However, this has only minor effects on the accuracy of R*I* and its dynamic variations, as *I* enters only as an average over the duration *τ*. The merits of R*I* show up clearly when compared to other definitions of reproduction indices which are cur-rently used in epidemic data based diagnostics of the infection dynamics. Since we will later use data from Germany in our analysis, we refer to what is issued by the Robert Koch Institute (RKI) in Germany as the “reproduction index”. It is based on the idea of calculating the ratio of incidence data, taken on two successive instants, separated by a delay time *t**G* [15]. The latter is called the generation time and represents the average time interval between an infection and a subsequent “successful” transmission of the infection to a third person. The RKI uses *t**G* = 4 days. The definition of this generation-time based reproduction index is then [16, 17] ![Formula][14] Its mathematical meaning becomes clearer in a continuous formulation, ![Formula][15] This can be written as ![Formula][16] where we have truncated the Taylor expansion after the first term. This reveals that R*G* is directly related to the logarithmic derivative of the daily incidence, |*∂**t**s*(*t*) |, with some delay equal to *t**G*. Hence whenever the daily incidence happens to vary exponentially, *∂**t**s*(*t*) ∝exp (R−1)*t*, R*G* can indeed be interpreted as a reproduction index. At any other time, however, when this is not the case, the use of R*G* as a reproduction index lacks mathematical foundation. Note furthermore that while *c* enters directly in R (as given by eqs. (3) and (4)), it cancels out in all expressions for R*G*. Hence public measures affecting *c* will readily show up in R*I*, which has been defined according to eq. (4), but not in R*G*. ## III. APPLICATION TO EPIDEMIC DATA From data obtained in Germany during the SARS-CoV2 pandemic in 2020 and 2021 [18], let us now calculate R*I* (*t*) and R*G*(*t*) by means of eqs. (6) and (7), respectively. The result is displayed in the top panel of Fig. 1, exposing the remarkable differences between the two quantities. There is a strong tendency of R*G*(*t*) to stay closer to unity than R*I* (*t*), which reflects the dynamics more pronouncedly. Events like the strong increase of R*I* (*t*) up to a value of 2.3 at the end of October 2020 (around day 245) hardly show up in R*G*(*t*). ![FIG. 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/04/16/2022.04.11.22273599/F1.medium.gif) [FIG. 1:](http://medrxiv.org/content/early/2022/04/16/2022.04.11.22273599/F1) FIG. 1: *Top panel:* the two reproduction indices R*G*(*t*) (dotted) and R*I* (*t*) (solid) as obtained form infection data in Germany during the SARS-CoV2 pandemic. The arrows indicate certain events and interventions (see table I) as discussed in the text. The grey shaded curve at the bottom indicates the prevalence (number of infected). the corresponding grey scale bar to the right corresponds to one half million people. *Bottom panel:* R*I* (*t*) as obtained from R*G*(*t*) through eq. (14) for different values of *τ* (black curves) and numerically from infection data (grey, same as black curve in top panel). The vertical arrows are located where the black curves reach unity, and coincide well with the extrema of the prevalence. In the interval *t* ∈ [75, 225], the prevalence curve has been scaled by a factor of 20 for visibility. For a more detailed discussion, a number of important events are listed in table I. After day 63, which corresponds to May 5, 2020, we see a sharp increase of R*I*. At the end of the first wave during March and April, the German Chancellor and the Conference of prime ministers of the Länder had decided to relax public life to almost normal conditions. Hence stores, restaurants, cultural institutions, and museums were opened. Because there was only little change in incidence (and prevalence) during the rest of spring and summer, it went unnoticed that the reproduction index was undergoing strong changes. That these were much less pronounced in R*G* (which was used) than in R*I* made their detection particularly difficult. View this table: [TABLE I:](http://medrxiv.org/content/early/2022/04/16/2022.04.11.22273599/T1) TABLE I: A number of events marked in the top panel of Fig. 1. A few days after a strong rise in R*G* was noticed, a major disease outbreak was reported on day 109 at the Tönnies sloughterhouse site near Guütersloh, among the large number of loan workers living at the site. But since prevalence remained low after that, the further rise of R was not noticed, or taken seriously, certainly in part due to the noisiness of the data. After the incidence then rose very sharply in the fall of 2020, the German government decided on what became known as a “lockdown light”, starting from November 2nd (day 245). This was associated with a hope of possibly easing restrictions for Christmas. Because the pronounced decrease of the reproduction index is not reflected by R*G*, and incidence did not seem to decrease, it was not noticed that the situation was actually relaxing during November. Consequently, a hard lockdown (including widespread closures of shops, businesses, schools, etc.) was announced on December 1st (day 275), which would start on December 16th (day 290). In particular, the relaxations for Christmas that had been previously promised were withdrawn. This announcement led to a short-lived but sharp increase in infections, as many people squeezed through the two-weeks bottleneck for their Christmas shopping. The subsequent decline (presumably due to the start of the Christmas school holidays) abruptly terminated on new year’s day (day 347), when many people had visited relatives and friends. A similar feature appears on day 400 at the Easter holidays, for similar reasons. Clearly, most (if not all) of these features are only poorly (if at all) discernible in R*G*. Nevertheless, it may seem that some features in R*G*(*t*) appear a little earlier than corresponding features in R*I* (*t*). This is apparent most clearly from the points where unity is crossed, which for R*G*(*t*) lie significantly before those for R*I* (*t*). This could be interpreted as R*G* being better suited for forecast purposes than R*I*, at it discloses the same information at an earlier time. However, this turns out to be a delusive mathematical artifact. As we see from eq. (9), with help form eq. (1), ![Graphic][17] can be written as ![Formula][18] By means of eqs. (4) and (1), the numerator of the second term can be expressed as ![Formula][19] Consequently, we have ![Formula][20] In other words, the deviation of ![Graphic][21] from unity is composed of a term proportional to the deviation of R(*t*) from unity and the logarithmic time derivative of R(*t*). Hence R*G* has a much more complex structure than R, exhibiting additional features (and additional noisiness) from the time derivative of R. It hence cannot be interpreted in terms of a true reproduction index. We have seen above that R*I* is a powerful diagnosis tool, as it reacts sensitively to events and interventions in society. Nevertheless, it comes with its downsides. From the second wave (around day 300) we see that the absolute magnitude of R*I* cannot be accurate, as the number of infections rises considerably shortly before day 300, while R*I* is clearly smaller than unity. This may be attributed to the fact that R*I* depends on *τ* (cf. eqs. (6)), which cancels out for R*G*, as it is obvious form eq. (7). The inherent problem is that *τ* may not only be known with poor accuracy, but it can also undergo gradual changes during the epidemic. An obvious cause may be the gradual appearance of mutations, which often result in shifts in the clinical picture, possibly including changes in the duration of the illness, hence in *τ*. It will only rarely be possible to receive regular reliable data on *τ*. ## IV. SELF-CONSISTENT CORRECTION OF R It turns out that the independence of R*G* from *τ* may be exploited here by inverting eq. (12), in order to calculate R*I* from R*G*. If we denote by ![Graphic][22] the deviation of ![Graphic][23] from unity, we can rewrite eq. (12) (by multiplying with *τ* and dividing by R) into ![Formula][24] This is a linear differential equation in *τ/*R and can be solved by means of the method of variation of parameters. The result is ![Formula][25] where ![Formula][26] One finds that ![Graphic][27] depends on *τ* only weakly, because its role as a prefactor im eq. (14) and its appearance in *p*(*t*) cancel each other to a large extent. In the bottom panel of Fig. 1, ![Graphic][28] thus obtained is plotted for three different values of *τ* as the black curves. The grey curve is R*I* as in the top panel. Clearly, all of the more prominent features of R*I* are reproduced. There are two main differences between ![Graphic][29] and R*I* (*t*). First, there is a vertical shift which varies with time only very slowly. Second, much finer details are visible in R*I*. The striking feature of ![Graphic][30] is that close to its transitions through unity, there is almost no sensitivity to *τ*. Hence should *τ* vary over the course of the epidemic by, e.g., as much as 40 percent (as between the dashed and the dotted curve), the shape of ![Graphic][31] would not change much. In particular, the points where it hits unity do not change their position appreciably. As the vertical arrows show, these points are very close to the extrema of the number of infected people, as one would correctly expect for the reproduction index. Hence what we display in the bottom panel of Fig. 1 may well be called the best of both worlds. In R*I* (*t*) we see very fine details which allow to identify the effects of social events, and to assess the effectiveness of public interventions in retrospect. In ![Graphic][32], which we derived from R*G*(*t*) by means of eq. (14), we see less detail, but obtain a more accurate estimate of the reproduction index. This allows for more reliable predictions of near-future epidemic dynamics, as ![Graphic][33] is particularly accurate close to unity. The difference between R*I* and R*G* is presumably due to a (slowly) varying *τ*. If we explicitly demand variation in *τ* to be slow as compared to the rapid variations seen in R*I*, we can approximately determine *τ* (*t*) from R*I* and R*G* in a self-consistent manner. First we calculate the ratio ![Graphic][34] and fit a low order polynomial *q*(*t*) to it, such that high-frequency components are cut off. Then we set *τ*1(*t*) = *τ**q*(*t*), where *τ* is the constant *τ* we used initially calculating R*I* from eq. (6). This yields a new ![Graphic][35], which we use, together with *τ*1(*t*), to recalculate ![Graphic][36] from eq. (14). This is repeated until the result of the polynomial fit has become stable. ![Graphic][37] is then a reproduction index which is self-consistently matched to the one obtained from R*G* and should, just as the latter, provide high accuracy (in particular close to unity) while keeping all the fine details we found to be present in R*I*, which we initially calculated from eq. (6). The result we obtained using fourth-order polynomials for *q*(*t*) is presented in Fig. 2. After five iterations, the result for *q*(*t*), did not change anymore. As one can clearly see, ![Graphic][38] has been rather well matched to![Graphic][39], aside from some high-frequency variations. We have exploited the separation of time scales between *τ* and R*I*. One may admit faster variations to *τ* in order to achieve an even better match between ![Graphic][40] and ![Graphic][41], but this shall not concern us here. *q*(*t*) is indicated by the dotted curve in the figure, to be read off the left scale. The scale to the right shows the corresponding values of *τ*. Note that the pronounced peak around day 240 reaches a very high value, deviating from unity by about four times more strongly than R*G* (dotted curve in Fig. 1). This highlights the importance of our method to calculate R*I* for predictions of near-future epidemic spreading dynamics. Around day 240, using R*G* would have (or actually has) over-estimated the doubling time of infections by a factor of 4.5. ![FIG. 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/04/16/2022.04.11.22273599/F2.medium.gif) [FIG. 2:](http://medrxiv.org/content/early/2022/04/16/2022.04.11.22273599/F2) FIG. 2: Self-consistent matching of the low-frequency components of R*I* (*t*) (black) and ![Graphic][42] (grey). Matching has been achieved by setting *τ* = *q*(*t*)*τ*, where *τ* = 14 has been used for the initial calculation of R*I* (*t*), and *q*(*t*) is a fourth-order polynomial. It is indicated by the dotted curve. The corresponding values of *τ* (*t*) can be read off the scale to the right. ## V. REVEALING THE ORBITAL STRUCTURE OF EPIDEMIC WAVES Finally, it is instructive to elaborate on some additional aspects of data presentation and analysis. In Fig. 3a we plot the number of new infections during a period *τ* against the prevalence, i.e., the total number of currently ill individuals. The ordinate is calculated from the seven-day averaged incidence by multiplying with *τ/*7. Each of the small circles represents one day, with the symbol style representing the three epidemic waves (open, first wave. full grey, second wave. full black, third wave). The data are gathering into an elongated cloud along the first diagonal (dashed line). If we assume that about 3% of individuals infected with SARS-CoV2 need intensive care, we can estimate the maximum prevalence the society could bear. Since there are 16734 intensive care beds in Germany [19], we conclude that the displayed range of the abscissa represents the maximum “acceptable” range of prevalence (one half million infected individuals). ![FIG. 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/04/16/2022.04.11.22273599/F3.medium.gif) [FIG. 3:](http://medrxiv.org/content/early/2022/04/16/2022.04.11.22273599/F3) FIG. 3: Day-by-day epidemic trajectories (circles with polygons). First wave (open), *t* < 130. Second wave (filled grey), *t* ∈ [130, 351]. Third wave (filled black), *t* > 351. (a) Total number of new infections within a time interval *τ* versus total number of acutely infected. Each wave appears as a clockwise orbit, with smaller sub-orbits. Data points tend to group along the first diagonal (dashed line). Solid curve: a sample simulation of an epidemic trajectory, showig the generic clockwise orbit structure, very similar to the data from the first wave. Both are asymptotic to the dash-dotted line with slope R = 3.3. The second (lower) asymptote corresponds to the terminal value of *R* = 0.14. The horizontal axis spans the maximum range of SARS-CoV2 infections acceptable to the German health system (about one half million). (b) Reproduction index versus seven-days incidence. The orbit structure is more clearly revealed. This presentation allows to assess in which phase of an epidemic wave the system currently is. All data stay just above the asymptotic (dotted curve) solution of eq. (1), which enters the vertical axis at the terminal value of *R*1 = 0.14. Various incidents mentions in table I are marked along the trajectory. Note that R is generally a forerunner to incidence, and a new orbit is marked as a sharp increase of R. Hence one may miss important developments when monitoring incidence alone. In order to analyse the internal structure of the data cloud, we compare with a numerically simulated sample trajectory of eq. (1), which is shown as the solid curve. It forms a lobe, starting off at the origin with a slope equal to R, proceeds clockwise (arrows) and re-enters the origin at an inferior slope of R1 = 0.14 (lower dash-dotted line) for *t→∞*. The initial slope, which is indicated by the upper dash-dotted line, follows eq. (4), since the ordinate and abscissa just represent *τ* |*∂**t**s*| and *i*, respectively. For the simulation, we have set R = 3.3 in order to match typical values assumed for SARS-CoV2. The size of the lobe corresponds to the number of people affected by the epidemic. As a consequence of the structure of the solutions to eq. (1), R1 is a function of R, and R1(3.3) = 0.14. In fact, the data representing the first wave (open circles close to the origin) exhibit just the same lobe shape, initially following the dash-dotted line, and as time proceeds is traversed in clockwise direction. A closer look at the second (grey) and third (black) wave reveals that their trajectories tend to form clockwise orbits as well, with smaller sub-orbits, thus exposing additional fine structure of the infection dynamics. Positioning data within the course of an orbit may be very useful to assess the stage of the epidemic spreading. This orbit structure is revealed even more clearly in a different presentation, when we plot R*I* versus incidence, as shown in Fig. 3b. What we can clearly see, for instance, is that the situation was already on a relaxing path when the “lockdown light” came into effect on day 245. When the hard December lockdown was announced on day 275, the orbit was in fact almost finished and back towards the origin. Things would thus have probably eased off by themselves, if no major mistakes would have been made. After Mardi Gras (day 347), incidence stayed calm, and R*G* remained featureless within noise level, as Fig. 1 (top) shows. Hence no countermeasures were taken. The trajectory in Fig. 3b, however, clearly shows that this was when a new orbit had formed. This very probably launched the third wave. Most importantly, we see from Fig. 3b that the widespread exclusive use of the incidence (abscissa) for assessing the state of the epidemic is void of any sound basis, as it does not reflect the orbital structure of epidemic waves. In particular, the use of threshold values for incidence in legislation on public mitigation interventions is clearly not appropriate. ## VI. CONCLUSIONS We have presented analysis tools for different aspects of epidemic mitigation and management interventions. First, we showed that from data of daily incidence (and prevalence derived therefrom) one can derive, via eq. (6), an accurate value for the reproduction index. R*I* (*t*) provides a very sensitive seismograph of the currrent state of the epidemic. Second, we have shown that a similar quantity, ![Graphic][43], can be derived via eq. (14) from the generation-based reproduction index, R*G*. It has the same general behaviour as R*I* (*t*), but has the particular merit of being very accurate whenever it is close to unity. This is of great interest for forecasting the epidemic development, which abruptly changes at R = 1. The fact that it lacks much of the fine details showing up in R*I* (*t*) makes it particularly suited for polynomial extrapolation towards unity. Third, we have shown that a much improved (corrected) form of R(*t*) can be obtained in a self-consistent manner, thereby also providing the so far only poorly known, slowly varying, *τ* (*t*). Finally, we have demonstrated that a presentation of standard epidemic data in the plane spanned by the reproduction index and the incidence exploits the internal orbit structure of epidemic waves in a way beneficial for assessing the current epidemic state of affairs. It should be a very useful tool for policy makers during dangerous epidemics, such as the recent outbreak of COVID-19. ## Data Availability All data are available under "daily reports of the Robert-Koch-Institute / Germany" [https://www.rki.de/EN/Home/homepage\_node.html;jsessionid=6B1BE151A64136D8C071E6DF5EF019DC.internet082](https://www.rki.de/EN/Home/homepage_node.html;jsessionid=6B1BE151A64136D8C071E6DF5EF019DC.internet082) ## VII. ACKNOWLEDGEMENTS The authors are grateful to Paul Leiderer, Christoph Klein, and Armin Lambacher for inspiring discussions and critical reading of the manuscript. This paper is dedicated to all those who died in the Covid19 pandemic, especially to those who passed away in loneliness. ## Footnotes * * Electronic address: robert.conradt{at}ds.mpg.de * Received April 11, 2022. * Revision received April 11, 2022. * Accepted April 16, 2022. * © 2022, 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].W. H. O. (WHO), Tech. Rep. (2020). 2. [2]. M. Enserink and K. Kupferschmidt, Science 367, 1414 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE1OiIzNjcvNjQ4NS8xNDE0LWEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNC8xNi8yMDIyLjA0LjExLjIyMjczNTk5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 3. [3]. H. M. Baniamin, R. Mizanur, and M. T. Hasan, Asia Pacific Journal of Public Administration 42, 153 (2020). 4. [4]. N. Haug, L. Geyrhofer, A. Londei, E. Dervic, A. Desvars-Larrive, V. Loreto, B. Pinior, S. Thurner, and P. Klimek, Nature Human Behaviour 4, 1303 (2020). 5. [5]. T.-C. Ng, H.-Y. Cheng, H.-H. Chang, C.-C. Liu, C.-C. Yang, S.-W. Jian, D.-P. Liu, T. Cohen, and H.-H. Lin, JAMA International Medicine 181, 913 (2021). 6. [6]. G. Lu, O. Razum, A. Jahn, Y. Zhang, B. Sutton, D. Sridhar, K. Aryoshi, L. v. Seidlein, and O. Muüller, Global Health Action 14, 1 (2021). 7. [7]. A. Cascio, M. Bosilkovski, A. J. Rodriguez-Morales, and G. Pappas, Clinical Microbiology and Infection 17, 336 (2021). 8. [8]. S. S. Morse, J. A. K. Mazet, M. Woolhouse, C. R. Parrish, D. Carroll, W. B. Karesh, C. Zambrana-Torrelio, W. I. Lipkin, and P. Daszak, The Lancet 380, 1956 (2012). 9. [9]. L.-F. Wang and G. Crameri, Rev. sci. tech. Off. int. Epiz. 33, 569 (2014). 10. [10]. J. Dehning, J. Zierenberg, F. P. Spitzner, M. Wibral, J. P. Neto, M. Wilczek, and V. Priesemann (2020), URL [https://doi.org/10.1101/2020.09.16.20187484](https://doi.org/10.1101/2020.09.16.20187484);. 11. [11]. H. W. Hethcote, SIAM Review 42, 599 (2000). [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%2F2022%2F04%2F16%2F2022.04.11.22273599.atom) 12. [12]. P. Godara, S. Herminghaus, and K. M. Heidemann, PLoS ONE 16, e0247445 (2021). 13. [13]. A. Cori, N. M. Ferguson, C. Fraser, and S. Cauchemez, Am. J. Epidemiol. 178, 1505 (2013), ISSN 00029262. N/Neuartiges{\_}Coronavirus/Situationsberichte/ Gesamt.html. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwt133&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24043437&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F16%2F2022.04.11.22273599.atom) 14. [14].Robert Koch Institut, Aktueller Lage-/Situationsbericht des RKI zu COVID-19 (2020), URL [https://www.rki.de/DE/Content/InfAZ/](https://www.rki.de/DE/Content/InfAZ/) 15. [15]. J. Dehning, F. P. Spitzner, M. C. Linden, S. B. Mohr, J. P. Neto, J. Zierenberg, M. Wibral, M. Wilczek, and V. Priesemann, Science 369 (2020), ISSN 0036-8075, URL [https://science.sciencemag.org/content/early/2020/05/14/science.abb9789](https://science.sciencemag.org/content/early/2020/05/14/science.abb9789). 16. [16].RKI, Tech. Rep., Robert Koch Institut, Berlin (2020). 17. [17]. M. an der Heiden and O. Hamouda, Epidemiol. Bull. pp. 10–16 (2020), ISSN 2569-5266. 18. [18].Robert Koch Institut, Nowcasting und R-Schätzung: Schätzung URL [https://www.rki.de/DE/Content/InfAZ/N/Neuartiges{\_}Coronavirus/Projekte{\_}RKI/Nowcasting.html](https://www.rki.de/DE/Content/InfAZ/N/Neuartiges{_}Coronavirus/Projekte{_}RKI/Nowcasting.html). 19. [19]. A. Stang, M. Stang, and K.H. Jöckel, Dtsch. Arztebl. Int. 117, 329 (2020), ISSN 18660452. [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/graphic-1.gif [4]: /embed/graphic-2.gif [5]: /embed/graphic-3.gif [6]: /embed/graphic-4.gif [7]: /embed/inline-graphic-3.gif [8]: /embed/inline-graphic-4.gif [9]: /embed/graphic-5.gif [10]: /embed/inline-graphic-5.gif [11]: /embed/inline-graphic-6.gif [12]: /embed/graphic-6.gif [13]: /embed/inline-graphic-7.gif [14]: /embed/graphic-7.gif [15]: /embed/graphic-8.gif [16]: /embed/graphic-9.gif [17]: /embed/inline-graphic-8.gif [18]: /embed/graphic-12.gif [19]: /embed/graphic-13.gif [20]: /embed/graphic-14.gif [21]: /embed/inline-graphic-9.gif [22]: /embed/inline-graphic-10.gif [23]: /embed/inline-graphic-11.gif [24]: /embed/graphic-15.gif [25]: /embed/graphic-16.gif [26]: /embed/graphic-17.gif [27]: /embed/inline-graphic-12.gif [28]: /embed/inline-graphic-13.gif [29]: /embed/inline-graphic-14.gif [30]: /embed/inline-graphic-15.gif [31]: /embed/inline-graphic-16.gif [32]: /embed/inline-graphic-17.gif [33]: /embed/inline-graphic-18.gif [34]: /embed/inline-graphic-19.gif [35]: /embed/inline-graphic-20.gif [36]: /embed/inline-graphic-21.gif [37]: /embed/inline-graphic-22.gif [38]: /embed/inline-graphic-23.gif [39]: /embed/inline-graphic-24.gif [40]: /embed/inline-graphic-25.gif [41]: /embed/inline-graphic-26.gif [42]: F2/embed/inline-graphic-27.gif [43]: /embed/inline-graphic-28.gif