Estimation of Daily Reproduction rates in COVID-19 Outbreak =========================================================== * Jacques Demongeot * Kayode Oshinubi * Hervé Seligmann * Florence Thuderoz ## Abstract **(1) Background** The estimation of daily reproduction rates throughout the infectivity period is rarely considered and only their sum *R**o* is calculated to quantify the level of virulence of an infectious agent; **(2) Methods** We give the equation of the discrete dynamics of epidemic growth and we obtain an estimation of the daily reproduction rates, by using a technique of deconvolution of the series of observed new cases of Covid-19; **(3) Results** We give both simulation results as well as estimations for several countries for the Covid-19 outbreak; **(4) Conclusions** We discuss the role of the noise on the precision of the estimation and we open on perspectives of forecasting methods to predict the distribution of daily reproduction rates along the infectivity period. Keywords * daily reproduction rate * Covid-19 outbreak * discrete epidemic growth equation * discrete deconvolution ## 1 Introduction After the Severe Acute Respiratory Syndrome outbreak by coronavirus SARS CoV in 2002 and the Middle East Respiratory Syndrome outbreak by coronavirus MERS CoV in 2012, the Covid-19 disease by coronavirus SARS CoV-2 is the third coronavirus outbreak in the past two decades. Human coronaviruses including 229E, OC43, NL63, and HKU1 are a group of viruses that cause a significant percentage of all common colds in humans. The SARS CoV-2 can be transmitted from person to person by respiratory droplets and through contact and fomites. Therefore, the severity of disease symptoms such as cough, sputum, and their viral load are most often important factors in the virus’s ability to spread, and these factors can change rapidly in just a few days during the period of infectivity of an individual. This ability to spread is quantified by the classic basic reproduction number *R**o* (also called basic reproduction ratio or rate), an epidemiology parameter which describes the transmissibility of an infectious agent and is equal to the mean number of susceptible individuals, an infected individual can contaminate during his infectivity period. *R**o* is not a biological constant for a pathogen: in fact, *R**o* is affected by numerous exogenous factors like geoclimate, demographic and socio-behavioural factors that govern pathogen transmission [1,2,3,4,5,6]. Because of these environmental factors, *R**o* might change seasonally, but these factors variations are not significant if a very short time in considered. *R**o* depends also on endogenous factors like the viral load of the infected individuals during their infectivity period; the variations of this viral load [8,9,10,11,12,13,14] is often neglected in the theoretical and applied studies on the Covid-19 outbreak, in which the authors estimate a unique reproduction number *R**o* linked to the Malthusian growth parameter of the exponential phase of the epidemic during which *R**o* is greater than 1 (Figure 1) and they rarely consider the distribution of partial daily reproduction numbers at day *j* of the infectivity period, denoted ![Graphic][1] [1]. When this distribution is considered, it is more convenient i) to define for each age class the distribution of the marginal daily reproduction numbers, ii) to estimate in each case its entropy and simulate the dynamics either using a Markovian model like that defined in Delbrück’s approach [2] or an ODE SIR model. In the Markovian case, *R**o* has to be replaced by the evolutionary entropy defined by L. Demetrius as the Kolmogorov-Sinaï entropy of the corresponding random process, which has properties concerning the stability of its Markovian invariant measure, analogue to the properties of the Malthusian parameter for the stability of the ODE’s steady state [3,4]. We compare these two modelling approaches in Section 2, then present in Section 3 as illustration results of estimation of daily reproduction rates in different countries at different periods. In Section 4, we discuss our method in the framework of the classical population dynamics (Leslie model) and eventually, in Section 5, we conclude and open the approach to the multi-age model. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/04/2020.12.30.20249010/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/F1) Figure 1. Spread of an epidemic disease from a first infected “patient zero”(in red) located on the centre of his influence sphere made of the successive generations of infected, for the same value of the reproduction number *R**o* = 3, with a deterministic dynamics (left) and a stochastic one (right) with standard deviation *σ* of the uniform distribution on an interval centred on *R**o* and with a random time interval *i* between infected generations (after [7]). ## 2 Materials and Methods ### 2.1 Relationship between Markovian and ODE SIR approach #### 2.1.1 Obtaining SIR equation from a discrete mechanism Let consider the classical reproduction rate *R**o*. If the model is deterministic, if *X**j* denotes the number of new cases at day *j*, and the contagious period is made of *r* consecutive days, with *R**k* the marginal reproduction rate at day *k* of the contagious period, we have: ![Formula][2] It is easy to show that, if *X**o* = 1 and *r* = 5 (inside the estimated interval of the duration of the maximal contagion period for the COVID-19 [5,6,7,8,9,10,11,12,13,14,15,16,17], we obtain: ![Formula][3] If *R*2 and *R*3 are dominant and equal to *R/*2, then *X*5 behaves as 2*R*2*R*3 = *R*2/2, which shows the difficulty to estimate *R**o*, which is the mean value of the ![Graphic][4] on the contagious period. The length of this period can be estimated from the ARIMA series of the stationary random variables ![Graphic][5], equal to the ![Graphic][6] without their trend: one can take the interval on which the auto-correlation function remains more than a certain threshold, e.g., 0.1 [15]. More generally if *R*1 = *a, R*2 = *b* and *R*3 = *c*, we get: ![Formula][7] If *R*1 and *R*2 equal respectively *a* and *b*, and if *a* = *b* = *R/*2, *c* = 0, then *X*5 behaves like: ![Formula][8] If *R* = 1, {*X**j*}*j*=1,*∞* is the Fibonacci sequence, and more generally, for *R >* 1, the generalised Fibonacci sequence. Let suppose now that *b* = 0 and *a* depends on the day *j* such as *a**j* = *νC*(*j*), where *C*(*j*) represents the number of possible susceptible individuals, which can be recruitable by one infectious individual at day *j*. If cumulated infected individuals (supposed to be all infectious) at day *j* are denoted by *I**j*, we have: ![Formula][9] Suppose that the first infected is recruited at the centre of his sphere of influence and that the secondary infected individuals remain in this sphere, by widening the radius on day *j*, therefore the susceptible individuals *C*(*j*), where each infectious on day *j* − 1 can recruit are located on a decreasing part of the sphere of influence of the first infected 0 (Figure 2). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/04/2020.12.30.20249010/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/F2) Figure 2. Spread of an epidemic disease from a first infected (located at his influence circle centre) progressively infecting all the neighbours in some regions (rectangles) on successive spatial layers. The function *C*(*j*) decreases due to the bulk on the successive spheres and we can consider the following functional form *C*(*j*) = *S*(*j*)/(*c* + *S*(*j*)), where *S*(*j*) is the number of susceptible individuals at day *j*. Then, we can write the following equation taking into account the mortality: ![Formula][10] The corresponding continuous equation is close to the SIR equation, if *c* is greater before *S*: ![Formula][11] #### 2.1.2 Second obtention of the SIR equation from a discrete mechanism Another way to derive the SIR equation from a probabilistic approach is to start from the microscopic equation of molecular shocks by Delbrück [2] which corresponds to a classical birth-and-death process: if at least one event (contact *ν*, birth *f*, death *μ* or recovering *ρ*) occurs in (*t, t* + *dt*), we have, if births compensate deaths, leaving constant the total size *N* of the population: ![Formula][12] Hence, we have, if *P**k*(*t*) denotes *P* (*S*(*t*) = *k, I*(*t*) = *N* − *k*): ![Formula][13] and we obtain: ![Formula][14] Then, by multiplying by *s**k* and summing over *k*, we obtain the characteristic function of the random variable *S*, which is proven to be a Poisson random variable if the coefficients *ν, f, μ* and *ρ* are sufficiently small. If births do not compensate deaths, we have: ![Formula][15] If *S* and *I* are independent and if the coefficients *ν, f, μ* and *ρ* are sufficiently small, they are Poisson random variables, whose expectations *E*(*S*) and *E*(*I*) verify: ![Formula][16] which leads to the SIR equations for the variables *S, I* and *R* considered as being deterministic: ![Formula][17] ### 2.2 Distribution of the rate of transmission along the infectiousness period of an individual Figure 3 gives the average transmission rate *R**o* calculated on the 25th of October 2020 just before the second lockdown (re-confinement) in France [16], but because at this day the second wave of the epidemic is still in its exponential phase, it is more convenient i) to consider the distribution of the marginal daily reproduction numbers, and ii) to calculate its entropy and simulate the epidemic dynamics using a Markovian model [2]. If *R**o* denotes the average transmission rate (or mean reproduction number) among the studied population, we can estimate the distribution *V* (whose coefficients are *V**j* = *R**j**/R**o*) of the daily reproduction numbers *R**j* along the infectiousness period, by remarking that the number *X**j* of new infectious cases at day *j*, equal to *X**j* = *I*(*j*)–*I*(*j* − 1), where *I*(*j*) is the number of infectious at day *j*, verifies the discrete convolution equation: ![Formula][18] where *r* is the duration of the contagion period, estimated by 1/(*k* +*μ*), where *k* is the recovering rate and *μ* the death rate in SIR equations: ![Formula][19] where *S* and *I* are respectively the size of the susceptible and infectious populations. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/04/2020.12.30.20249010/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/F3) Figure 3. Top: estimation of the average transmission rate *R**o* for the 20th and the 25th October 2020 [18]. Bottom left: daily new cases in France between February 15 and Oct 27 [6]. Bottom right: V-shape of the evolution of the daily ![Graphic][20] along the infectious 3-day period of an individual. If *r* and *S* can be considered as constant during the first exponential phase of the pandemic, we can also assume that the distribution *V* is constant and then, *V* can be estimated by solving the linear system: ![Formula][21] which can be written as *X* = *MR*, hence giving: ![Formula][22] and the equation (13) can be solved numerically, if the pandemic is observed during a time greater than 1/(*k* + *μ*). Then, the entropy of *V* = *R/R**o* can be calculated, as the Kolmogorov-Sinaï entropy of the Markovian Delbrück scheme [2] ruling the ![Graphic][23] and giving new parameters for characterising pandemic dynamics, namely for quantifying its robustness and stability [3,4]. ### 2.3 The biphasic pattern of the virulence curve in coronaviruses Mostly, the clinical course of patients with seasonal influenza shows a biphasic occurrence of symptoms with two distinct peaks. Patients have a classic influenza disease followed by an improvement period and a recurrence of the symptoms [1,2,3]. The influenza RNA virus shedding (the time during which a person might be infectious to another person) increases sharply one-half to one day after infection, peaks on day 2 and persists for an average total duration of 4.5 days, between 3 and 6 days, which explains that we will choose in the following as infectivity duration these extreme values, i.e., either 3 or 6 days depending on the positivity of the estimated daily reproduction rates. It is common to see this biphasic influenza clinically: after incubation of one day, there is a high fever, then a drop in temperature before rising, hence the term “*V*” fever. The other symptoms like coughing often also have this improvement on the second day of the flu attack: after a first feverish rise (39 − 40°*C*), the temperature drops to 38°*C* on 2nd day, then rises before disappearing on the 5th day, the fever being accompanied by respiratory signs (coughing, sneezing, clear rhinorrhea, etc.). By looking the shape of virulence curves observed in coronaviruses patients [5,6,7,8,9,10,11,12], we often see this biphasic pattern. ## 3 Results ### 3.1 Distribution of the daily reproduction numbers R*j*’s along the infectiousness period of an individual. A theoretical deterministic example #### Deterministic Toy Model Let’s suppose that *a priori* ![Graphic][24] are equal to: ![Formula][25] then the inverse of the matrix *M* is given by: ![Formula][26] and the deconvolution gives the *a posteriori* ![Graphic][27]: ![Formula][28] thanks to the following calculation: ![Formula][29] and we obtain for the *a posteriori* distribution of the daily reproduction rates the exact replica of the *a priori* distribution. ### 3.2 Distribution of the daily reproduction numbers R*j*’s. A simulated stochastic example Let’s consider a stochastic version of the deterministic toy model given in Section 3.1., by introducing an increasing noise on the ![Graphic][30], e.g., a uniform distribution on the three following intervals:[2−*a*, 2+*a*], [1−*a/*2, 1+*a/*2], [2− *a*, 2 + *a*], with increasing values of *a*, from 0.1 to 1, in order to see when the deconvolution would give negative *a posteriori* ![Graphic][31], with conservation of the average of their sum *R**o*, if you repeat the random choice of the values of the ![Graphic][32] at each generation? We will give the results of the random simulations for increasing values of *a* from *a* = 0.1 to *a* = 1. 1. **For a= 0.1**, let’s choose the *a priori* distribution of the daily reproduction numbers *R*1 in the interval [1.9, 2.1], *R*2 in [0.95, 1.05] and *R*3 in [1.9, 2.1] as *R*1 = 2.1, *R*2 = 0.95, *R*3 = 2.1. Then, transition matrix *M*1 is equal to: ![Formula][33] and we have: ![Formula][34] From *X*6 = 113.491, *X*5 = 41.7391, *X*4 = 15.351, *a posteriori* ![Graphic][35] can be calculated: *R*1 = 2.1, *R*2 = 0.95, *R*3 = 2.1 The next *a priori* ![Graphic][36] are chosen as: *R*1 = 2, *R*2 = 0.95, *R*3 = 1.9 and we have: ![Formula][37] Then, we get the matrices *M*2 and ![Graphic][38] ![Formula][39] The *a posteriori* ![Graphic][40] equal: *R*1 = 2.0279, *R*2 = 7.6158, *R*3 = −16.426 The following *a priori* ![Graphic][41] are: *R*1 = 2, *R*2 = 1.05, *R*3 = 1.9 and we have: ![Formula][42] Then, we get the matrices *M*3 and ![Graphic][43] ![Formula][44] The *a posteriori* ![Graphic][45] equal: *R*1 = 2.486, *R*2 = −2.33, *R*3 = 7.38769 The next *a priori* ![Graphic][46] are: *R*1 = 1.9, *R*2 = 1.05, *R*3 = 1.9 and we have: ![Formula][47] Then, we get the matrices *M*4 and ![Graphic][48] ![Formula][49] The *a posteriori* ![Graphic][50] equal: *R*1 = 2.69, *R*2 = −16.72, *R*3 = 43.809 2. **For a=1**, let’s choose the *a priori R*1 in [1.0, 3.0], *R*2 in [0.5, 1.5] and *R*3 in [1, 3] e.g., *R*1 = 1, *R*2 = 1.355, *R*3 = 1.1. Then, transition matrix is equal to: ![Formula][51] and its inverse is given by: ![Formula][52] The new cases are: *X*6 = 18.209, *X*5 = 9.101, *X*4 = 4.81, *X*3 = 2.355, *X*2 = 1, *X*1 = 1, and by deconvolution, we get *a posteriori* ![Graphic][53] equal: *R*1 = 1, *R*2 = 1.355, *R*3 = 1.1, i.e., the exact *a priori* distribution. Let’s consider now a new *a priori R*1 = 1, *R*2 = 1, *R*3 = 1. That gives a new matrix *M*2, with new *X*7 and *X*8 calculated from the new *a priori* ![Graphic][54], by using the former values of *X*6, &, *X*2 : ![Formula][55] Hence, we get: ![Formula][56] Then, by inverting *M*2, we can calculate the *a posteriori* ![Graphic][57]: ![Formula][58] and *a posteriori* ![Graphic][59] equal: *R*1 = 2.90, *R*2 = 5.4888, *R*3 = −14.696 We continue the process by calculating *X*9 and *X*10 using new *a priori R*1 = 3, *R*2 = 0.5, *R*3 = 2.9 : ![Formula][60] Hence, we get: ![Formula][61] and *a posteriori* ![Graphic][62] equal: *R*1 = 3.66898, *R*2 = −33.857, *R*3 = 61.32 Let’s choose now a new *a priori R*1 = 2.6, *R*2 = 0.7, *R*3 = 2.6 Then, we have: ![Formula][63] and we get the matrices *M*4 and![Graphic][64] : ![Formula][65] Then, *a posteriori* ![Graphic][66] equal: *R*1 = 2.99859, *R*2 = −1.785, *R*3 = 7.139 More precise simulation results are given in Table 1, which summarises computations made for random choices of *a priori* distributions, for *a* = 0.1 and *a* = 1. These simulations show a great sensitivity to the noise, but a qualitative conservation of the down biphasic (D-B) shape of their distribution along the infectivity period of individuals. View this table: [Table 1.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/T1) Table 1. Simulation results obtained for extreme noises *a* = 0.1 and *a* = 1, showing great variations of the deconvoluted *a posteriori* distribution of the daily reproduction numbers ![Graphic][67] and a qualitative conservation of the down biphasic(D-B) shape of their distribution along the infectivity period. Note: INV means Inverted ### 3.3 Distribution of the daily reproduction numbers R*j*’s. The real example of France The Figure 3 gives the average transmission rate *R**o* calculated the 25th of October 2020 just before the second lockdown (re-confinement) in France [18]. Because the second wave of the epidemic is still in its exponential phase, it is more convenient i) to consider the distribution of the marginal daily reproduction numbers, and ii) to calculate its entropy and simulate the epidemic dynamics using a Markovian model [2]. By using the daily new infectious cases given by [6], we can calculate *M* −1 for the period from 20th to 25th October 2020, by choosing 3 days for the duration of the infectiousness period and the following raw data for the new infected cases (Figure 3): Oct 25 : 52010, 45422, 42032, 41622, 26676, 20468 : Oct 20 We have: ![Formula][68] Hence, we can deduce the daily ![Graphic][69] i.e., the vector (*R*1, *R*2, *R*3) : ![Formula][70] The average transmission rate is equal to *R**o* ≈ 1, 174, value close to that calculated directly (Figure 3), giving *V* = (0.7, 0.085, 0.215), with a maximal daily reproduction rate the first day of the infectiousness period. The entropy *H* of *V* is equal to: *H* = −∑*k*=1,*r* *V**k**Log*(*V**k*) = 0.25 + 0.21 + 0.33 = 0.79. ### 3.4 Calculations of the R*j*’s for different countries #### 3.4.1 Chile By using the daily new infectious cases given by [6], we can calculate *M* −1 for the period from 1st to 12th November 2020, by choosing 6 days for the duration of the infectiousness period and the following 7-day moving average data for the new infected cases (Figure 4): ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/04/2020.12.30.20249010/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/F4) Figure 4. Top: estimation of the average transmission rate *R**o* for the 1st and the 12th November 2020 [18]. Bottom left: Daily new cases in Chile between November 1 and November 12 [6]. Bottom right: V-shape of the evolution of the daily ![Graphic][71] along the infectious 6-day period of an individual. Nov 12 : 1408, 1394, 1387, 1384, 1385, 1389, 1405, 1362, 1359, 1382, 1370, 1400 : Nov1 We have: ![Formula][72] Hence, after deconvolution we have: ![Formula][73] The average transmission rate is equal to *R**o* ≈ 1, 011, value close to that calculated directly, with a maximal daily reproduction rate the last day of the infectiousness period. Because of the negativity of *R*1, we cannot derive the distribution *V* and therefore calculate its entropy. The quasi-endemic situation in Chile since the end of August, which corresponds to the increase of temperature and drought at this period of the year [15], gives a constancy of new cases of infected and a periodicity of their occurrence corresponding to the length of the infectious period of about 6 days, analogue to the cyclic phenomenon observed in simulated stochastic data of Section 3.2. with a similar M-shaped distribution of the ![Graphic][74]. #### 3.4.2 Russia By using the daily new infectious cases given by [6], we can calculate *M* −1 for the period from September 30 to October 5, 2020, by choosing 3 days for the duration of the infectiousness period and the following 7-day moving average data for the new infected cases (Figure 5): ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/04/2020.12.30.20249010/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/F5) Figure 5. Top: estimation of the average transmission rate *R**o* for the 30th September and the 5th October 2020 [18]. Bottom left: Daily new cases in Russia between September 30 and October 5 [6]. Bottom right: V-shape of the evolution of the daily ![Graphic][75] along the infectious 3-day period of an individual. Oct 05 : 9473, 9081, 8704, 8371, 8056, 7721 : Sept 30 We have: ![Formula][76] and ![Formula][77] where: ![Formula][78] The average transmission rate is equal to *R**o* ≈ 1, 073, value close to that calculated directly, with a maximal daily reproduction rate the first day of the infectiousness period. Because of the negativity of *R*2, we cannot derive the distribution *V* and therefore calculate its entropy. The period studied corresponds to a local slow increase of new infectious cases at the start of the second wave in Russia, which seems to have a succession of slightly inclined 4-day plateaus followed by a step. #### 3.4.3 Nigeria By using the daily new infectious cases given by [6], we can calculate *M* −1 for the period from November 5 to November 10 2020, by choosing 3 days for the duration of the infectiousness period and the following raw data for the new infected cases (Figure 6) : ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/04/2020.12.30.20249010/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/F6) Figure 6. Top: estimation of the average transmission rate *R**o* for the 5th and the 10th November 2020 [18]. Bottom left: Daily new cases in Nigeria between November 5 and November 10 [6]. Bottom right: increasing evolution of the daily ![Graphic][79] along the infectious 3-day period of an individual. Nov 10 : 166, 164, 161, 133, 149, 141 : Nov5 We have: ![Formula][80] After deconvolution, we get: ![Formula][81] The average transmission rate is equal to *R**o* ≈ 1.129, value close to that calculated directly, with a maximal daily reproduction rate the last day of the infectiousness period. *V* = (0.143, 0.342, 0.515) and the entropy *H* of *V* is equal to : *H* = −∑*k*=1,*r* *V**k**Log*(*V**k*) = 0.29 + 0.37 + 0.34 = 1. ### 3.5 Weekly patterns in daily infection cases Daily new infectious cases are highly affected by weekdays, such that case numbers are lowest at the start of the week, and increase afterwards. This pattern is observed at the global level, as well as at the level of almost each single country or USA state. Hence, in order to estimate biologically meaningful infection rates, clean of weekly patterns due to administrative constraints, analyses have to be restricted to specific periods shorter than a week, or at the rare occasions when patterns escape the administrative constraints on weekly patterns. This weekly phenomenon occurs during exponential increase as well as decrease phases of the pandemic, and during stagnation periods in numbers of daily cases. In addition, the daily new infection case record is discontinuous for many countries/regions, which frequently publish on Monday or Tuesday a cumulative count for that day and the weekend days. For example, Sweden typically publishes over one week only four numbers, the one on Tuesday cumulating cases for Saturday, Sunday and the two first weekdays. Discontinuity in records further limit the availability of data enabling detailed analyses of daily infection rates and can be considered as extreme weekday effects on new case records due to various administrative constraints. We calculated Pearson correlation coefficients *r* between a running window of daily new case numbers of 20 consecutive days, and a running window of identical duration with different intervals between the two running windows. There is typically a peak with a time-lag of seven days between the two running windows. This could reflect a biological phenomenon of seven infection days. However, examinations of frequency distributions of time-lags for r maxima reveals, besides the median at time-lag 7 days, local maxima for multiples of 7 (14, 21, 28, 35, etc). About 50 percent of all local maxima in *r* involve time-lags that are multiples of seven (seven included). This excludes a biological causation. We tried to control for week days using two methods, and combinations thereof. For the first method, we calculated z-scores for each weekday, considering the mean number of cases for each weekday, and subtracted that mean from the observed number for a day. This delta was then divided by the standard deviation for the number of cases for that day. This was done for each weekday. The second method implies data smoothing using a running window of 5 consecutive days, where the mean number of cases calculated across the five days is subtracted from the observed number of cases for the third day. Hence data for a given day are compared to a mean including two previous, and two later days. We constructed two further datasets, one where the second data-smoothing method is applied to the z-scores from the first method, and the other, where the z-scores from the first method are applied to the data after data-smoothing from the second method (Figures 10 and 11). ![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/04/2020.12.30.20249010/F7.medium.gif) [Figure 7.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/F7) Figure 7. Daily new confirmed cases for the whole world as a function of days since February 24 until August 23.+ indicate Sundays, x indicate Mondays. ![Figure 8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/04/2020.12.30.20249010/F8.medium.gif) [Figure 8.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/F8) Figure 8. Z-transformed scores of daily new confirmed cases for the whole world, from Figure 7, as a function of days since February 24, until August 23.+ indicate Sundays, x indicate Mondays. Z-transformations are specific to each weekday. ![Figure 9.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/04/2020.12.30.20249010/F9.medium.gif) [Figure 9.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/F9) Figure 9. Smoothed daily new confirmed cases for the whole world, from Figure 7, as a function of days since February 24, until August 23. + indicate Sundays, x indicate Mondays. For each specific day *j*, the mean number of new confirmed cases calculated for days *j* − 1, *j* − 2, *j* + 1 and *j* + 2 is subtracted from the number for day *j*. ![Figure 10.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/04/2020.12.30.20249010/F10.medium.gif) [Figure 10.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/F10) Figure 10. Smoothed daily new confirmed cases for the whole world applied to *Z* scores from Figure 8, as a function of days since February 24, until August 23. + indicate Sundays, x indicate Mondays. For each specific day *j*, the mean number of new confirmed cases calculated for days *j* − 1, *j* − 2, *j* + 1 and *j* + 2 is subtracted from the number for day *j*. ![Figure 11.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/04/2020.12.30.20249010/F11.medium.gif) [Figure 11.](http://medrxiv.org/content/early/2021/01/04/2020.12.30.20249010/F11) Figure 11. Z-transformed scores of daily new confirmed cases for the whole world, smoothed data from Figure 9, as a function of days since February 24, until August 23.+ indicate Sundays, x indicate Mondays. Z-transformations are specific to each weekday. These four datasets transformed according to different methods and combinations thereof designed to control for weekday were analysed using the running window method. Despite attempts at controlling for weekday effects, the median time-lag was always seven days across all four transformed datasets, and local maxima in time-lag distributions were multiples of seven. Despite data transformations, about 50 percent of all local maxima were time-lags that are multiples of seven, seven included. Visual inspection of plots of these transformed data versus time for daily new infection cases from the whole world systematically show systematic local biases in daily new infection cases (after transformations) on Sundays and Mondays, for all four transformed datasets, with Sundays and/or Mondays as local minima and/or local maxima, along which method or combination thereof. Hence, the methods we used failed to neutralise the weekly patterns in daily new cases due to administrative constraints. This issue highly limits the data available for detailed analyses of daily new cases aimed at estimating biologically relevant estimates of infection rates at the level of short temporal scales, but it reinforces the hypotheses done in previous sections on the existence of an heterogeneity of the virulence along the infectiousness period, independent of administrative constraints observed in collecting new daily cases. ## 4 Discussion Rhodes and Demetrius have pointed out the interest of the distribution *V* of the daily reproduction numbers [5]. In particular, they found that this distribution was generally not uniform, which we have confirmed here by showing many cases with the biphasic form of the virulence observed in respiratory viruses, such as those of influenza. The entropy of the distribution *V* makes it possible to evaluate the intensity of this biphasic character. By taking into account the mortality due to the Covid-19, the discrete dynamics of new cases can be considered as a Leslie dynamics governed by the matrix equation **X***j* = L**X***j*−1, where **X***j* is the vector of the new cases living at day *j* and *L* is the Leslie matrix given by: ![Formula][82] and ![Formula][83] where *b**j* = 1 − *μ**j* ≤ 1, ∀*i* = 1, &, *r*, is the survival recovered probability between days *j* and *j* + 1. The dynamical *X**k* stability for *L*2 distance to the stationary infection age pyramid **P**= *lim**j***X***j**/∑* *k*=*j,j*−*r*+1 *X**k* is related to |*λ* − *λ*,|, the modulus of the difference between the dominant and sub-dominant eigenvalues of *L*, namely *λ* = *e**R* and *λ*, (where *R* is the Malthusian growth rate), where **P** is the left eigenvector of *L* corresponding to *λ*. The dynamical stability for the distance (or symmetrized divergence) of Kullback-Leibler to **P** considered as stationary distribution is related to the population entropy *H* [4,19,20,21,22,23,24], which is defined if *I**j* = ∏ *i*=1,*j*−1 *b**i* and *p**j* = *I**j**R**j**/λ**j*, as follows: ![Formula][84] The mathematical characterisation by population entropy defined in (18) of the stochastic stability of the dynamics described by equation (17) has its origin in the theory of large deviations [25,26,27]. This notion of stability pertains to the rate at which the system returns to its steady state conditions after a random exogenous and/or endogenous perturbation and it could be useful to quantify the variations of the distribution of the daily reproduction rates observed for many countries. ## 5 Conclusions and Perspectives The article has used a discrete approach for describing the dynamics of the Covid-19 outbreak. In the ODE SIR case [28,29,30], marginal ![Graphic][85] can be estimated in the same way, but by using continuous deconvolution equation. An improvement in this framework could be the introduction of the demographic notion of age class, because the daily reproduction rates distribution is depending on the state of the immune defences both of the infected and susceptible individuals, and on the mode of transmission, this mode being associated in social settings recording age and infection schedule with more secondary cases for young and adult than for elderly at home [31,32]. Then, for each age class *k* a reproduction number *R**ok* can be calculated by summing the marginal ![Graphic][86], which correspond to this age class *k*, and eventually a global *R**o* can be estimated from the ![Graphic][87], weighted by the proportion of their age class in the whole population. Both global *R**o* and dedicated to age classes ![Graphic][88] could serve to follow the effects of public health measures as lockdowns and quarantines. ## Data Availability All data comes from public data bases. ## 6 Author Contributions Conceptualisation, J.D.; methodology, J.D., H.S.; K.O. and F.T. have performed calculations and Figures; all authors have equally participated to the other steps of the article elaboration. ## 8 Conflicts of Interest The authors declare no conflict of interest ## 7 Acknowledgements The authors hereby give their thanks to the International Excellence Fellowship of Karlsruhe Institute of Technology (KIT) funded within the Framework of the University of Excellence Concept “The Research University in the Helmholtz Association I Living the Change”. ## Footnotes * Jacques.Demongeot{at}univ-grenoble-alpes.fr, Kayode.Oshinubi{at}univ-grenoble-alpes.fr, florence.thuderoz{at}gmail.com, varanuseremius{at}gmail.com * Received December 30, 2020. * Revision received December 30, 2020. * Accepted January 4, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. [1].Ouassou, I.; Hobbad, L.; Alahiane, M.; Gaudart, J.; Iggui, S.; Rachdi, M.; Demongeot, J. The ARIMA model to analyse incidence pattern and estimate short-term forecasts for retro-predicting the first wave of the COVID-19 outbreak. Math. Population Studies (under review). 2. [2].Delbrück, M. Statistical fluctuations in autocatalytic reactions. J. Chem. Physics 1940, 8, 120–124. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1063/1.1750549&link_type=DOI) 3. [3].Demetrius, L. Boltzmann; Darwin and the directionality theory. Physics Reports 2013,530, 1–86. 4. [4].Demongeot, J.; Demetrius, L. Complexity and Stability in Biological Systems. Int. J. Bifurcation and Chaos 2015, 25, 40013. 5. [5].Rhodes, C.J.; Demetrius, L. Evolutionary Entropy Determines Invasion Success in Emergent Epidemics. PLoS ONE 2010, 5, e12951. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20886082&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F04%2F2020.12.30.20249010.atom) 6. [6].Worldometers. Available online: [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/) (accessed on 2nd November 2020) 2020. 7. [7].Lacoude, P. Available online: [https://www.contrepoints.org/2020/07/22/376624-covid-19-le-debut-de-la-fin-1](https://www.contrepoints.org/2020/07/22/376624-covid-19-le-debut-de-la-fin-1) (accessed on 22nd November 2020) 2020. 8. [8].Carrat, F.; Vergu, E.; Ferguson, N.M.; Lematre, M.; Cauchemez, S.; Leach, S; Valleron, A.J. Lines of Infection and Disease in Human In-fluenza: A Review of Volunteer. Am. J. Epidemiology 2008, 167, 775–785. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwm375&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18230677&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F04%2F2020.12.30.20249010.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000254469500002&link_type=ISI) 9. [9].Pan, Y.; Zhang, D.; Yang, P.; Poon, L.L.M.; Wang, Q. Viral load of SARS-CoV-2 in clinical samples. The Lancet Infectious Diseases 2020, 20, 411–412. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30113-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32105638&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F04%2F2020.12.30.20249010.atom) 10. [10].Wölfel, R.; Corman, V.M.; Guggemos, W.; Seilmaier, M.; Zange, S.; Muller, M.A.; Niemeyer, D.; Jones, T.C.; Vollma, P.; Rothe, C.; Hoelscher, M.; Bleicker, T.; Brunink, S.; Schneider, J.; Ehmann, R.; Zwirglmaier, K.; Drosten, C.; Wendtner, C. Virological Assessment of hospitalized patients with COVID-2019. Nature 2020, 581, 465–469. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2196-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F04%2F2020.12.30.20249010.atom) 11. [11].Liu, W.D.; Chang, S. Y.; Wang, J.T.; Tsai, M.J.; Hung, C. C; Hsu, C.L.; Chang, S.C. Prolonged Virus shedding even after Seroconversion in a Patient with COVID-19. J. Infection 2020, 81, 318–356. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jinf.2020.04.025&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F04%2F2020.12.30.20249010.atom) 12. [12].Ferretti, L.; Wymant, C.; Kendall, M.; Shao, L.; Nurtay, A.; Abeler-Dorner, L.; Parker, M. D. Bonsall, D.; Fraser, C. Quantifying SARS-CoV-2 transmission suggets epidemic control with digital contact tracing. Science 2020, 368, eabb6936. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNjgvNjQ5MS9lYWJiNjkzNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAxLzA0LzIwMjAuMTIuMzAuMjAyNDkwMTAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 13. [13].Cheng, H.Y.; Jian, S.W.; Liu, D.P.; Ng, T.C.; Huang, W.T.; Lin, H.H. Contact Tracing Assessment of COVID-19 Transmission Dynamics in Taiwan and Risk at Different Exposure Periods before and after Symptom Onset. JAMA Intern. Med. 2020, 180, 1156–1163. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F04%2F2020.12.30.20249010.atom) 14. [14].He, X.; Lau, E.H.Y.; Wu, P.; Deng, X.; Wang, J.; Hao, X.; Lau, Y.C.; Wong, J.Y.; Guan, Y.; Tan, X.; Mo, X.; Chen, Y.; Liao, B.; Chen, W.; Hu, F.; Zhang, Q.; Zhong, M.; Wu, Y.; Zhao, L.; Zhang, F.; Cowling, B.J.; Li, F.; Leung, G.M. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nat. Med. 2020, 26, 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%2F2021%2F01%2F04%2F2020.12.30.20249010.atom) 15. [15].Demongeot, J.; Flet-Berliac, Y.; and Seligmann, H. Temperature decreases spread parameters of the new COVID-19 cases dynamics. Biology (Basel) 2020, 9,94. 16. [16].Seligmann, H.; Iggui, S.; Rachdi, M.; Vuillerme, N.; Demongeot, J. Inverted covariate effects for mutated 2nd vs 1st wave Covid-19: high temperature spread biased for young. Biology (Basel) 2020, 9, 226. 17. [17].Seligmann, H.; Vuillerme, N.; Demongeot, J. Summer COVID-19 third wave: faster high altitude spread suggests high UV adaptation. MedRxiv 2020, [https://doi.org/10.1101/2020.08.17.20176628](https://doi.org/10.1101/2020.08.17.20176628). 18. [18].Renkulab. Available online: [https://renkulab.shinyapps.ioCOVID-19-Epidemic-Forecastingw](https://renkulab.shinyapps.ioCOVID-19-Epidemic-Forecastingw) (accessed on 22nd November 2020) 2020. 19. [19].Demetrius, L. Demographic parameters and natural selection. Proc. Natl Acad. Sci. USA 1974, 71, 4645–4647. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiNzEvMTIvNDY0NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAxLzA0LzIwMjAuMTIuMzAuMjAyNDkwMTAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 20. [20].Demetrius, L. Statistical mechanics and population biology. Journal of Statistical Physics 1983, 30, 709–750. 21. [21].Demongeot, J.; Demetrius, L. La dérive démographique et la sélection naturelle : Etude empirique de la France (1850-1965). Population 1989, 2, 231–248. 22. [22].Demongeot, J. Biological boundaries and biological age. Acta Biotheoretica 2009, 57, 397–419. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19907923&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F04%2F2020.12.30.20249010.atom) 23. [23].Gaudart, J.; Ghassani, M.; Mintsa, J.; Rachdi, M.; Waku, J.; Demongeot, J. Demography and Diffusion in epidemics: Malaria and Black Death spread. Acta Biotheoretica 2010, 58, 277–305. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10441-010-9103-z&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20706773&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F04%2F2020.12.30.20249010.atom) 24. [24].Demongeot, J.; Hansen, O.; Hessami, H.; Jannot, A.S.; Mintsa, J.; Rachdi, M.; Taramasco, C. Random modelling of contagious diseases. Acta Biotheoretica 2013, 61, 141–172. 25. [25].Wentzell, A.D.; Freidlin, M.I. On small random perturbations of dynamical systems. Russ. Math. Surv. 1970, 25, 1–55. 26. [26].Donsker, M.D.; Varadhan, S.R.S. Asymptotic evaluation of certain Markov process expectations for large time. I. Comm. Pure Appl. Math. 1975, 28, 1–47. 27. [27].Freidlin, M.I.; Wentzell, A.D. Random Perturbations of Dynamical Systems. Springer, New York 1984. 28. [28].Scarpino, S.V.; Petri, G. On the predictability of infectious disease outbreaks. Nature Commun 2019, 10, 898. 29. [29].Liu, Z.; Magal, P.; Seydi, O.; Webb, G. Understanding Unreported Cases in the covid-19 Epidemic Outbreak in Wuhan, China, and Importance of Major Public Health Interventions. Biology (Basel) 2020, 9, 50. 30. [30].Demongeot, J.; Griette, Q.; Magal, P. Computations of the transmission rates in SI epidemic model applied to COVID-19 data in mainland China. Royal Society Open Science 2020, 7, 201878. 31. [31].Adam, D.C.; Wu, P.; Wong, J.Y.; Lau, E.H.Y; Tsang, T.K.; Cauchemez, S.; Leung, G.M.; Cowling, B.J. Clustering and superspreading potential of SARS-CoV-2 infections in Hong Kong. Nat Med 2020, 26, 1714–1719. 32. [32].Nishiura, H.; Lintona, N.M.; Akhmetzhanov, A.R. Serial interval of novel coronavirus (COVID-19) infections. International Journal of Infectious Diseases 2020, 93, 284–286. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.02.060&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32145466&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F04%2F2020.12.30.20249010.atom) [1]: /embed/inline-graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/inline-graphic-2.gif [5]: /embed/inline-graphic-3.gif [6]: /embed/inline-graphic-4.gif [7]: /embed/graphic-4.gif [8]: /embed/graphic-5.gif [9]: /embed/graphic-6.gif [10]: /embed/graphic-8.gif [11]: /embed/graphic-9.gif [12]: /embed/graphic-10.gif [13]: /embed/graphic-11.gif [14]: /embed/graphic-12.gif [15]: /embed/graphic-13.gif [16]: /embed/graphic-14.gif [17]: /embed/graphic-15.gif [18]: /embed/graphic-17.gif [19]: /embed/graphic-18.gif [20]: F3/embed/inline-graphic-5.gif [21]: /embed/graphic-19.gif [22]: /embed/graphic-20.gif [23]: /embed/inline-graphic-6.gif [24]: /embed/inline-graphic-7.gif [25]: /embed/graphic-21.gif [26]: /embed/graphic-22.gif [27]: /embed/inline-graphic-8.gif [28]: /embed/graphic-23.gif [29]: /embed/graphic-24.gif [30]: /embed/inline-graphic-9.gif [31]: /embed/inline-graphic-10.gif [32]: /embed/inline-graphic-11.gif [33]: /embed/graphic-25.gif [34]: /embed/graphic-26.gif [35]: /embed/inline-graphic-12.gif [36]: /embed/inline-graphic-13.gif [37]: /embed/graphic-27.gif [38]: /embed/inline-graphic-14.gif [39]: /embed/graphic-28.gif [40]: /embed/inline-graphic-15.gif [41]: /embed/inline-graphic-16.gif [42]: /embed/graphic-29.gif [43]: /embed/inline-graphic-17.gif [44]: /embed/graphic-30.gif [45]: /embed/inline-graphic-18.gif [46]: /embed/inline-graphic-19.gif [47]: /embed/graphic-31.gif [48]: /embed/inline-graphic-20.gif [49]: /embed/graphic-32.gif [50]: /embed/inline-graphic-21.gif [51]: /embed/graphic-33.gif [52]: /embed/graphic-34.gif [53]: /embed/inline-graphic-22.gif [54]: /embed/inline-graphic-23.gif [55]: /embed/graphic-35.gif [56]: /embed/graphic-36.gif [57]: /embed/inline-graphic-24.gif [58]: /embed/graphic-37.gif [59]: /embed/inline-graphic-25.gif [60]: /embed/graphic-38.gif [61]: /embed/graphic-39.gif [62]: /embed/inline-graphic-26.gif [63]: /embed/graphic-40.gif [64]: /embed/inline-graphic-27.gif [65]: /embed/graphic-41.gif [66]: /embed/inline-graphic-28.gif [67]: T1/embed/inline-graphic-29.gif [68]: /embed/graphic-43.gif [69]: /embed/inline-graphic-30.gif [70]: /embed/graphic-44.gif [71]: F4/embed/inline-graphic-31.gif [72]: /embed/graphic-46.gif [73]: /embed/graphic-47.gif [74]: /embed/inline-graphic-32.gif [75]: F5/embed/inline-graphic-33.gif [76]: /embed/graphic-49.gif [77]: /embed/graphic-50.gif [78]: /embed/graphic-51.gif [79]: F6/embed/inline-graphic-34.gif [80]: /embed/graphic-53.gif [81]: /embed/graphic-54.gif [82]: /embed/graphic-60.gif [83]: /embed/graphic-61.gif [84]: /embed/graphic-62.gif [85]: /embed/inline-graphic-35.gif [86]: /embed/inline-graphic-36.gif [87]: /embed/inline-graphic-37.gif [88]: /embed/inline-graphic-38.gif