Numerically statistical investigation of the partly super-exponential growth rate in the COVID-19 pandemic (throughout the world)

In [1] the multiplication process of the particles in a random medium σ with asymptotically exponential (via the time t) number n⁢(t,σ)≈C⁢(σ)⁢exp⁡(λ⁢(σ)⁢t) was considered. It was assumed that the distribution of random variable λ⁢(σ) is Gaussian, the limiting process n⁢(t,σ)→t→∞C⁢(σ)⁢exp⁡(λ⁢(σ)⁢t) is homogeneous over σ and the quantities C⁢(σ), exp⁡(λ⁢(σ)⁢t) are asymptotically weakly correlated. On this basis the following estimate of the function E⁢n⁢(t,σ)=n⁢(t) was obtained

where a=E⁢λ⁢(σ) is the mathematical expectation and d2=D⁢λ⁢(σ) is the variance of random λ. By the use of [2, integral formula (2.3.15), No. 11] the following expression was obtained:

Note that formula (1) may be the basis for numerical investigations of concrete problems. The law of increase, corresponding to formula (1), may be naturally called “super-exponential”. More total and practically convenient definition of super-exponentiality is related to the growth of the increase coefficientα⁢(t)=n⁢(t+Δ⁢t)n⁢(t) while t rises.

In real problems the relation -∞<λ1<λ⁢(σ)<λ2<+∞ holds. Therefore, the obtained asymptotics possibly approximates n⁢(t) only in the interval 0<T1≤t≤T2<+∞. So, when solving real problems, it is expedient to perform additional numerical investigations. Similar investigations were fulfilled in [1] for the special stochastic model of particles multiplication process.

The above results show that growth rate of mean number of the arbitrary nature particles (for example, microorganisms) may be super-exponential. This behavior corresponds to the rise of the increase coefficient α of the particle numbers {ni}, which are fixed at the equidistant moments {ti}, i.e., the ratio αi=ni+1ni growths. On the other hand if this rise is observed, then it is possible that the distribution of breeding grounds is related with it.

In particular, the novel coronavirus pandemic in the world behaved like this according to the WHO (World Health Organization) statistics [3] in the time interval from 9.03.2020 to 21.03.2020. The corresponding number of confirmed cases (via days) is approximated as in formula (1) with the error not exceeding 2 % as

for i=9,10,…,21. Some comparative results are given in Table 1.

Approximation (2) was obtained as follows. Let

and therefore

Solving these equations by the least-squares method (i.e., by the linear regression method) gives us coefficients in formula (2). Comparison of the results obtained by this formula with statistical data is given in Figure 1 in logarithmic scale. Comparatively with formula (2) the rise of WHO statistics weakens over time, in particular, due to the medical activities. When t>21, this growth can be exponentially approximated. This possibly happens because fluctuations of λ diminish and therefore the quantity D⁢λ decreases.

Figure 1

The number of confirmed cases of COVID-19 in the world (circles) and approximation by the formula (2) (solid curve)

It is remarkable that accordingly to expression (1) formula (2) may be obtained by averaging the random number n⁢(t,σ), for which λ has Gaussian distribution with parameters

So, if the estimates of the quantities E⁢λ, D⁢λ for t=t0 are known, then it is possible to extrapolate n⁢(t) by formula (2) significantly more accurate than those made by n⁢(t)=n⁢(t0)⁢exp⁡(E⁢λ⁢(t-t0)).