Abstract
The analysis of COVID-19 statistics in different regions of the world with the intention of understanding the trend of progression of the pandemic is a task of paramount importance. Publicly available data includes cumulative number of cases, new cases each day, and the mortality. Extracting information from this data necessitates mathematical modeling. In this note a simple technique is adopted to determine the trend towards stabilized elimination of the infection. Results pertaining to several representative regions of the world are presented.
Introduction
A question of utmost importance to individual nations and the whole is how the COVID pandemic would progress and whether available statistics can be utilized to determine the future trend in terms of mathematical models. In this note, we present a simple technique of analyzing publicly available COVID data based on the cumulative number and new cases each day for several geopolitical regions. Analysis indicates a current general trend of cumulative COVID cases approaching an equilibrium – possibly a sign that the pandemic is slowly moving towards a halt, perhaps because of the effectiveness of control measures. In few regions, slow exponential growth seems to continue.
Method and Discussion
The publicly available data on progression of COVID in different geographical regions of the world are:
Cumulative number of cases N at the beginning of a given date (defined as at the instant of time t).
New cases dN reported at the end of that day (after a time interval dt = 1, therefore dN is effectively ). The rate of growth can also be computed from plots of N vs t. OriginLab data analysis software is used plot curves and perform correlation analysis.
The very initial growth of the epidemic in many regions has been exponential1-2, therefore we assume that it evolves according to the dynamical system equation, where the growth constant F(N) is a function of N.
From (1)
The pandemic stops brewing and reach an equilibrium when after attaining a certain cumulative number N = NE. The point N = NE is asymptotically stable3 at this point, provided . Clearly this condition satisfied if
Eqn. (1) also embodies the requirement that point that point N = 0 corresponds to an unstable equilibrium. To illustrate the point, consider the logistic model, where F(N) = a - bN (a and b positive constants), with the equilibrium point NE = a/b as the equilibrium is stable. The cumulative number of cases in some epidemics approximately agree with the logistic model4-7. Our analysis shows that COVID-19 cumulative numbers deviate from the logistic variation.
The plot of F(N) vs N is same as the plot of vs N, if F(N) turns out be monotonically decreasing and reaching zero, either cutting the N-axis or approaching N-axis asymptotically, the pandemic will practically reach an equilibrium. In the simple logistic model F(N) is a negative slope straight line, cutting N-axis at the point of equilibrium. Possible scenarios are depicted in the Fig. 1 below, which portray different modes of variation of F(N) with N.
In the figures (a) and (b) F(N) decreases continuously, and in (a) a stable equilibrium is reached at N = NE, whereas in (b) equilibrium is approached asymptomatically. In figures (c) and (d), the monotonic decrease of N occurs after reaching a peak. In dealing with real data, the behaviour of F(N) at points close to N = 0 is obviously erratic, but does interfere with the subsequent trend.
The plots of F(N) vs N for several representative geopolitical regions are presented below. All computations are based on publicly available data8.
In most the above plots the later trend in variation of F(N) with N is monotonic decrease, indicating that remains negative and approach N-axis towards a stable equilibrium. However, there is no signal yet that curves will cut the N-axis, abruptly halting the pandemic in the near future. In few plots the current trend approximates to a straight line nearly parallel to N-axis, indicating that slow exponential growth continue.
Fig. 3 shows the plot of F(N) vs. N for Sri Lanka. Although there is no matured COVID epidemic in Sri Lanka, the current trend is a monotonically decreasing F(N).
Relaxation of control measures, mutations of the virus or influences of the climate may alter the observed trends F(N) shown in Figs. 2 and 3.
A limitation of the model is assumption of the equation (1) where F(N) not an explicit function of t. Dynamical systems interpretations, apply only under the above assumption.
Conclusion
The analysis suggests that fitting COVID-19 cumulative cases N to an equation of the form gives insight into the progression pandemic in different regions. As is essentially the number of cases reported each day, the function F(N) can be easily plotted. Since this quantity can be readily ascertained from logarithmic plots of the cumulative number The variation of F(N) with N depicts the progression towards stabilization of the cumulative number. If F(N) approaches N-axis with negative slope, indication is that the equilibrium point NE within reach. This tendency is visible majority of the curves presented - most probably a signature of success of the containment measures. More information about evolution of the pandemic could follow on updating of the plots.
Data Availability
Public available data: All Data Retrieved from JHU Data Repository:(Last access 5/06/2020)
Conflicts of Interest
Authors declare that there are no conflicts of interest
Funding
This study is not funded by organization or an individual.
Disclaimer
Authors are not responsible for the authenticity of data used and making decisions based on this study.
Acknowledgement
Authors wish to thank Prof. D. Tantrigoda, University of Sri Jayewardhanapura, Sri Lanka for valuable discussions, constructive comments and reading the manuscript.