The κ-statistics approach to epidemiology.

A great variety of complex physical, natural and artificial systems are governed by statistical distributions, which often follow a standard exponential function in the bulk, while their tail obeys the Pareto power law. The recently introduced [Formula: see text]-statistics framework predicts distribution functions with this feature. A growing number of applications in different fields of investigation are beginning to prove the relevance and effectiveness of [Formula: see text]-statistics in fitting empirical data. In this paper, we use [Formula: see text]-statistics to formulate a statistical approach for epidemiological analysis. We validate the theoretical results by fitting the derived [Formula: see text]-Weibull distributions with data from the plague pandemic of 1417 in Florence as well as data from the COVID-19 pandemic in China over the entire cycle that concludes in April 16, 2020. As further validation of the proposed approach we present a more systematic analysis of COVID-19 data from countries such as Germany, Italy, Spain and United Kingdom, obtaining very good agreement between theoretical predictions and empirical observations. For these countries we also study the entire first cycle of the pandemic which extends until the end of July 2020. The fact that both the data of the Florence plague and those of the Covid-19 pandemic are successfully described by the same theoretical model, even though the two events are caused by different diseases and they are separated by more than 600 years, is evidence that the [Formula: see text]-Weibull model has universal features.

Introduction

Many phenomena in a large variety of disciplines present apparent regularities described by well-known statistical distributions[1]. Some examples include the populations of cities, the frequencies of words in a text, the energy distribution in solids, the behaviour of financial markets, the statistical distribution of the kinetic energy in a gas, etc.[2,3]. Moreover, statistical approaches provide powerful tools for medical and epidemiological applications, since they allow predicting the behaviour of certain diseases[4,5]. Thus, such approaches are very useful for informing health policy and decision making, particularly regarding control and mitigation measures in response to the societal impacts of epidemics and pandemics.

While an epidemic is defined as “the occurrence in a community or region of cases of an illness clearly in excess of normal expectancy”[6,7], a pandemic is defined as “an epidemic occurring over a very wide area, crossing international boundaries, and usually affecting a large number of people”[6,7]. Pandemics are large-scale outbreaks of infectious diseases which cause a growth in mortality over a wide geographic area.

In addition to the obvious health consequences, both epidemics and pandemics also cause economic stress and social hardship. It has been emphasized that the probability of occurrence of pandemics is increasing due to the high inter-connectedness of the modern world, which is facilitated by the ease of travel and a continuous increase of urbanisation[6,8,9]. Consequently, the international community is increasing its efforts toward the mitigation of the impacts of pandemics[6]. Some recent examples of pandemics are the 2003 SARS (Severe Acute Respiratory Syndrome), the 2014 West Africa Ebola epidemic, and the present COVID-19 (COronaVIrus Disease) pandemic caused by the SARS-CoV-2 virus.

In this paper we formulate a statistical thermodynamics approach to epidemiology, which demonstrates the utility of -statistics for the analysis of epidemics and pandemics. First, we must consider certain facts which form the biomedical base of these phenomena. Below we summarise the principle causes of pandemics[6]:Recently, pandemics have also found their origin in zoonotic transmissions of pathogens from animals to humans[10].

Pathogens such as influenza viruses that are capable of efficient transmission among humans, with a high potential to cause global and severe pandemics. They are characterized by a relatively long and asymptomatic infectious period, during which it is not possible to detect the infected persons and their movements.

Pathogens such as Nipah virus which are capable of generating a moderate global threat. They are transmitted efficiently as a result of mutations and adaptation.

Pathogens such as Ebola that are capable of causing local epidemics, with a non-negligible risk of evolving into a global pandemic.

The risk of pandemic spread is conditioned by the following factors[11]:In order to better predict and mitigate the societal impact, the ability to quantify the morbidity and mortality associated with pandemics is very important. Consequently, statistical approaches based on available data that can lead to accurate models for predicting the behaviour of future pandemics are very useful. Historical data from past pandemics play a fundamental role, because they enable comparisons with theoretical models and can also be used to assess the performance of model-based forecasts. Historical records are often sparse and incomplete. Nonetheless, in all fields of research, scientists and engineers always search for “optimal” statistical distributions that can reliably predict the behaviour of different natural, engineered, and social systems[12-14].

pathogen factors such as genetic adaptation and mode of transmission;

human factors such as population density, susceptibility to infection, travel, migration, poverty, malnutrition, and caloric deficits;

factors related to public policy such as public health surveillance and measurements.

Since the beginning of the Covid-19 epidemic, various works in the scientific literature have presented statistical analyses of the epidemiological data. Most of these works focus on the number of infections[15-21] while the fatality curves have only been analyzed in a few studies. The analysis in the paper[16] focuses on the cumulative data (cumulative distribution function) of Covid-19, using the Richardson growth model which depends on three free parameters. The data related to Covid-19 deaths per day (i.e., the empirical probability density function) are analyzed in the paper[15] by adopting a four-parameter model which is based on non-extensive statistics. However, to our knowledge the development of a statistical model that admits explicit expressions for both the mortality cumulative distribution and the respective probability density function remains an open problem.

A second question concerns the mathematical features of suitable statistical models and in particular the behavior of the tails of the distribution functions. All of the existing models depend on three to five free parameters. However, a good fit of the data with a statistical model that depends on various free parameters is not sufficient proof of model validity. Two different model classes are mainly used to describe the Covid-19 data. The statistical distributions in the first class include the Weibull, Richardson, and SIR (Susceptible-Infected-Recovered) model functions. These depend on three to five free parameters and exhibit exponential or stretched exponential tails. The second class includes statistical distributions with two, three, or four free parameters that feature power-law tails. Models in this class include the Pareto, Zipf, Burr, Student and Tsallis distributions.

The present work aims to address the two problems that are described above, i.e., analytical tractability and tail behavior. We propose a statistical model for the analysis of the fatality curves that originated in physics but has also been successfully used in other fields such as econophysics, finance, and seismology. This model depends only on three free parameters but nonetheless provides excellent fits to the empirical curves obtained from the data. In addition, the model admits simple, closed-form analytic expressions for all the main statistical functions such as the probability density function, the cumulative distribution function, the survival function, the quantile function, the hazard function and the cumulative hazard function. Hence, this model overcomes the problem of analytical tractability and allows straightforward analysis of both daily and cumulative empirical data.

Regarding the second problem (tail behavior), the proposed -Weibull model is closely related to both classes mentioned above, but it does not strictly belong in either class. Indeed, the model’s distribution function is a one-parameter continuous deformation of the Weibull distribution, which belongs in the first model class. However, the tail of the deformed distribution follows Pareto’s power law as the distributions in the second class. Two of the three model parameters, namely and , correspond to the standard Weibull parameters. The third parameter, , is related to deformation of the tail of the distribution into a power law; is directly linked to the Pareto exponent n of the survival function through the simple relationship .

In order to test the validity of the -Weibull model we analyze various epidemics with different characteristics. First, we analyze the mortality data from the Florence plague that occurred in 1417[22-24]. These historical mortality data are of particular interest, because they were carefully recorded and allow testing the model with an epidemic process that evolved over several months. A second test of the -Weibull model involves the diffusion of Covid-19 in China[25,26]. The mortality curve derived from the Chinese Covid-19 data becomes flat after April 15, 2020. On the other hand, on April 17, 2020 China reported 1290 additional deaths which occurred during the entire cycle of the epidemic. However, no details were given regarding the temporal distribution of the fatalities. Hence, these data cannot be incorporated in our analysis. After April 15, 2020 only two cases of death have been reported to date (both in May 2020). However, these concerned infected individuals who had returned to China from abroad. For this reason, it is considered that the mortality cycle from Covid-19 in China concluded on April 16. The Chinese data analyzed herein are the “partial data” until April 16, 2020 which exclude the additional 1290 deaths that were a posteriori reported on April 17.

To further test the -Weibull model, we present a systematic analysis of COVID-19 data from other countries such as Germany, Italy, Spain and United Kingdom[25,26]. The analyzed data from these countries spans the time window from the beginning of the epidemic in February 2020 until the beginning of August 2020. This temporal range completely captures the entire cycle of the first wave of the epidemic. Florence’s data are completely different from those of Covid-19 in terms of the epidemiological cause, the duration, the fatality rate and the historical time. The fact that the -Weibull model can successfully capture the behavior of data so different as those of Florence and Covid-19, as well as data from other fields of science (economics, finance, seismology) indicates that it possesses universality features which make it suitable for a wide range of applications.

Methods

The fundamental difficulty faced by mathematical approaches to epidemiology is that all forecasts strongly depend on the model employed, the parameter estimates, as well as the choices of the initial conditions. There are both deterministic and stochastic epidemiological models[5,27]. The -deformed model that we propose herein has its roots in statistical mechanics and is therefore stochastic by construction.

The -deformed exponential function

Let us consider the function n(t) to represent the number density of deaths at any given time , where is the temporal domain of interest. In most cases of interest , where . The probability of death within a short time interval around t is given by , where f(t) is the probability density function (pdf). Then, the respective number density is given by , where N is the total number of deaths.

In statistical mechanics a general rate equation for f(t) is the first-order linear ordinary differential equation (ODE)where the function r(t) is the decay rate. The solution of the above ODE is the exponentialwith the standard normalisation condition which determines the constant c:The normalization also enforces a constraint on the number density, i.e.,

In the context of the exponential solution, the following three simple cases must be considered.The -exponential represents a very powerful tool which can be used to formulate a generalized statistical theory capable of treating systems described by distribution functions that exhibit power-law tails[46-48]. The mechanism generating the -exponential function is based on first principles from special relativity, and therefore the new function appears very promising for physical applications. Generalized statistical mechanics, based on the -exponential, preserves the main features of ordinary Boltzmann-Gibbs statistical mechanics which is based on the ordinary exponential through the Boltzmann factor. For this reason, it has attracted the interest of many researchers over the last two decades who have studied its foundations and mathematical aspects[49-58], the underlying thermodynamics[59,60], and specific applications of the theory to various fields. A non-exhaustive list of applications includes those in quantum statistics[61-64], in quantum entanglement[65,66], in plasma physics[67-73], in nuclear fission[74-77], in astrophysics[78-81], in quantum gravity[82-87], in geomechanics[88,89], in genomics[90,91], in complex networks[92-94], in economy[95-99], in finance[100-103], as well as in reliability analysis and seismology[104-106].

The Exponential Model: this model is fundamental in every branch of science; indeed, it allows us to describe a great variety of phenomena, from elasticity to electricity, from nuclear decay to thermal transient response, to name a few. The exponential model can be obtained for constant decay rate, i.e., which leads to the following exponential pdf: with .

Power-law Model (Pareto distribution): Zipf[28-31], in his studies on the size distribution of cities, incomes and word frequencies, pointed out the notion of regularity in the distribution of sizes. In a great variety of these cases the distributions follow a power law with an exponent close to , also known as Zipf’s Law[32-38]. Such power law distributions have been considered with increasing interest in the description of regular distributions. Distributions with exponents different from are known as Pareto’s law[3,36,39]. In general, distributions with power-law tails are known as heavy-tailed, fat-tailed or subexponential distributions, in juxtaposition to distributions whose tails decay exponentially[2].

Perline introduced the following classification of power laws[40]: Pareto obtained his Type I model, named the Pareto law by fitting the available data in his time on increasing social inequalities. He concluded that only economic growth can increase the income of the poor and decrease inequality[42]. Today, we know that his conclusions were partially right and partially wrong, because economic growth and equity are strictly related to the social relations of production, to the technological level and institutional structures of the economy, and to the composition, accumulation and distribution of human capital, as well as the quality and accessibility of the educational and financial structures in place[43]. Moreover, since Pareto disseminated his approach in 1897[44], the application of heavy-tailed distributions in economics has been developed[45].

A power law is strong if the distribution is a power law over the entire domain of definition;

A power law is weak if only part of the distribution is fitted by a power law;

A power law is false if only a highly truncated part of the distribution is approximated by a power law (in the scientific literature there are many examples of false power laws[41]).

The Pareto pdf f(t) is obtained from Eq. (1) by inserting the following time-dependent decay rate function which leads to a power-law solution for the pdf, i.e., in the domain , with :

-Exponential Model: This model, which is based on a fundamental approach derived from relativity[46], has proved useful in many applications. Experimental evidence suggests that probability density functions should resemble the exponential function for . However, for the Pareto pdf diverges. On the other hand, for high values of t many experimental results show a Pareto-like pdf with power law tails instead of exponential decay. Consequently, for it follows that while for it follows . So, the actual decay rate function r(t) should smoothly interpolate between these two regimes; a good proposal for r(t) has been introduced in the context of special relativity, where the function r(t) is given in terms of the Lorentz factor. We recall the expression of the Lorentz factor ; this expression involves the dimensionless momentum q where the parameter is the reciprocal of the dimensionless light speed c, i.e. . After posing or more explicitly it follows that for the decay rate r(t) approaches the exponential regime, i.e. . On the other hand, for it follows that r(t) approaches the decay rate of the Pareto model, i.e. .

The solution of the rate equation in this case yields the following pdf where the -deformed exponential function is given by with . It is important to note that in the limit and in the limit the function approaches the ordinary exponential , i.e. On the other hand the function for presents a power-law tail, i.e. Furthermore, the -exponential satisfies the following identity in analogy with the standard, non-deformed, exponential.

The -deformed statistical model

Given a pdf f(t) which represents the death rate, the cumulative distribution function (cdf) represents the probability of death between the initial time and the current time t. F(t), which is also known as the lifetime distribution, is given by means of the following integral

Conversely, the pdf f(t) is given by the derivative of F(t). In the following, we will assume without loss of generality that . The complement of F(t), i.e., is known as the survival function, and it represents the probability of survival at time t.

In most models of population dynamics the following time-dependent monomial is introducedThe expression (15) for T contains the real-valued parameters and . We can think of T(t) as the time measured by a nonlinear clock[107]. T(t) is regularly used in the definition of the survival function S(t) which becomes an implicit function of time based on the dependence of S on , i.e. .

The survival function of the Weibull model is then defined according toand it represents a stretched exponential function in time.

Popular models for empirical data that exhibit power-law tails include the Log-Logistic model[108], with survival function given bythe Burr type XII or Singh-Maddala model[109], with survival function given byand the Dagum model[110], with survival function given byVarious other postulates can be used for the analytical expression of the survival function. However, for physical applications it is extremely important to identify physical mechanisms or first principles which lead to such expressions for the survival function.

The -deformed statistical model can be viewed as a one-parameter generalization of the Weibull model, obtained by replacing the ordinary exponential in the definition (16) of the survival function by the -deformed exponential . Then, using the Weibull dependence for T given by Eq. (15), we obtain the following expression for the -deformed survival function or more explicitlyThe -deformed survival function reduces to the ordinary survival function of the Weibull model in the limit. The rate equation obeyed by is expressed in terms of T as a linear first-order ODEwith initial condition . The ODE (22) represents an interpolation between the rate equation of the exponential model and that of the Pareto model; moreover, the functional form in Eq. (22) is dictated by the first principles of special relativity.

The most important feature of the function is that it continuously interpolates between a power-law tail for large , i.e.,and exponential dependencefor .

The lifetime distribution function is given by the expressionwhile the pdf , defined , is given byThe rate equation for the survival function assumes the form of the first-order linear ODEwhere , and is the hazard function (hazard rate) defined throughHence, the hazard function assumes the expressionThe cumulative hazard function is defined by means of the integraland is linked with throughAfter taking into account that the -exponential function (10) can also be expressed in the formthe cumulative hazard function assumes the following explicit expressionand in the limit reduces to the standard Weibull cumulative hazard function . It is important to note that in the standard Weibull model the cumulative hazard function coincides with the function T(t).

Finally, the quantile function is defined as the inverse of the survival function follows . By expressing the quantile function in the form , one easily obtainswhere the -logarithm is the inverse function of , i.e. and is given byAfter observing that the function approaches the function in the limit, it follows that in the same limit reduces to the quantile function of the standard Weibull model.

Regarding the meaning of the parameters , and , we observe that the first two are the same as in the Weibull model and therefore have the same meaning. In particular, is the shape parameter while is linked to the scale parameter through the relationship . The deformation parameter , on the other hand, is linked to the Pareto exponent n of the survival function through the simple relationship and to the Pareto exponent p of the pdf via ; the former follows from the asymptotic dependence of the survival function given by Eq. (23), while the latter from the fact that .

The mode of the probability density function, i.e. the time where the pdf attains its maximum value can be easily obtained analytically as a function of the free parameters , and . More explicitly, after setting the first time derivative of the probability density function equal to zero we obtain the following biquadratic equationwith .

This equation can be solved to obtain as a function of and . It can also be solved to determine as a function of and , i.e.,This last expression is very useful in the analysis of empirical data.

Finally, the -deformed statistical model presented above, which represents a one-parameter continuous deformation of the Weibull model, has been successfully applied in econophysics for the analysis of personal income distribution[95] and in seismology[104]. We believe that the applicability of the model to such complex dynamical systems which are quite different from each other suggests that the model involves universal features.

Results

The purpose of this section is the validation of the -Weibull statistical model which is described in detail in the “Methods” section. The validation is performed by applying the model to various pandemic data. An intriguing first application regards the analysis of the number of deaths during the plague pandemic that ravaged the city of Florence in the XV century.

The Bubonic Plague or Black Death arrived in Europe in 1348 and in Italy in the spring of the same year[22,23]. It should be noted that the plague inspired Boccaccio to write his famous nouvelle Decameron only a couple of years after the end of the pandemic. Its consequence was the death of around 25-50% of Europe’s population by 1351. The pandemic is believed to have started in China and came to Europe via the trading routes though Asia and the Black Sea. In the year 1417, Florence had a population of 60,000 inhabitants and the Italian city’s Grain Office maintained a series of Books of Dead which recorded the number of human deaths caused by the bubonic plague. Florentines contracted the disease through contact with infected black rats and the fleas that they carried[23]. From a thermodynamic point of view, deaths which are due to the spread of an epidemic can be considered as an irreversible process inside an open system. Herein we show how the evolution of the plague in relation to the time and the number of deceased people can be explained using the -deformed statistics. In Table 1 the statistical data of the 1417 plague in Florence are summarised as a function of time[24].

Table 1

Temporal distribution of plague victims and population in Florence during the year 1417[24].

Month	Dead people	Population of Florence
May	600	59400
June	700	58700
July	2700	56000
August	5000	51000
September	2000	49000
October	600	48400
November	200	48200
December	100	48100

We consider the death occurrence probability as the ratio between the monthly number of deaths (given by the second column of the Table 1) and the cumulative number of deaths. We assume that the latter is equal to the total number of deaths that were recorded between May and December increased by an additional 100 deaths; we hypothesize that the latter occurred in the months following December, raising the total number of deaths to .

Figure 1 shows the quantile function (left) and the cumulative distribution function (right), versus the time t for the Florence epidemic. The theoretical quantile function of the -deformed model, defined in Eq. (34), is just the bisectrix of the plane represented by the continuous straight line. The continuous sigmoid curve represents the theoretical cumulative distribution function as given by Eq. (25). The empirical data related to the two functions (marked on the plots by dots) are deduced from the data in Table 1. The optimal fit parameters have the following values: , and . It is clear that the model describes the empirical data very well. Deviations of the empirical quantiles from the theoretical curves are likely due to detection errors since as it is evidenced in Table 1, the number of monthly deaths is approximated with an error margin of the order of 100 units. This statistical error is responsible for the dispersion of empirical quantiles around the bisectrix. Hence, the optimal parameter estimate clearly shows the difference (supported by the data) between the optimal -deformed model and the standard Weibull model corresponding to value .

Figure 1

Theoretical (continuous curve) and empirical (dots) plots of the quantile function (left) and the cumulative distribution function (right) versus time for the 1417 Florence plague epidemic. The theoretical curves are based on Eqs. (34) and (25), respectively.

As second application we analyze the data of Covid-19 mortality in China during the winter and spring of 2020[25,26]. The Chinese data on SARS-CoV-2 are important, as they represent the first statistical record for the entire evolutionary cycle of the COVID-19 pandemic in the world. It is important to note that the first and only cycle of Covid-19 ended in China with a total of 3,342 deaths around April 16, 2020. After April 17, 2020 while the Covid-19 cycle in China had already ended, the authorities reported 1290 additional deaths that occurred outside the hospitals during the entire cycle without disclosing their temporal distribution. For this reason the present analysis is applied to the 3,342 deaths which were reported by April 16, 2020.

Figure 2 reports the probability density function Eq. (26), the cumulative distribution function Eq. (25), the quantile function Eq. (34), the survival function Eq. (21), the hazard function Eq. (29) and the cumulative hazard function Eq. (33) versus the time t (measured in days) for the COVID-19 pandemic in China. The best fit parameters have the values: , and . The plots in this figure confirm the goodness of fit of the -deformed model. In particular, there is no systematic deviation of the empirical quantiles from the theoretical linear trend while the statistical nature of data leads to slight dispersion around the bisectrix. Furthermore the survival function versus time t is shown on a log-log plot in order to better explore differences between the empirical data and the theoretical predictions, especially in the tail of the distribution. The almost linear decay of the tail in this log-log plot is remarkable, showing clear signs that the Chinese Covid-19 data follows the Pareto power law in the distribution tail. The dispersion around the theoretical curves of the empirical data for both and reflects statistical fluctuations that are independent on the adopted theoretical model. These fluctuations are averaged out and become less evident in the cumulative functions , , and .

Figure 2

Theoretical (continuous curve) and empirical (dots) plots of the probability density function (top left), cumulative distribution function (top right), quantile function (middle left), survival function (middle right), hazard function (bottom left) and cumulative hazard function (bottom right) versus time for the Covid-19 mortality data related to China 2020 epidemic. The theoretical curves are based on Eqs. (26), (25), (34), (21), (29) and (33) respectively.

In the following we focus on the pandemic in Europe; in particular, we analyze the mortality Covid-19 data of Germany, Italy, Spain and United Kingdom, during the spring of 2020[25,26]. In these European countries the virus began to spread later than in China while the mortality exceeded the Chinese rate by more than 10 times.

Figure 3 shows the probability density function versus time related to Covid-19 mortality daily data for Germany, Italy, Spain and United Kingdom. The theoretical curve is given by Eq. (26). The parameter for the four countries assumes the values 3.83, 3.90, 4.00 and 4.15 respectively. The parameter for the four countries assumes the values 1.68, 1.67, 1.70 and 2.185 respectively. Finally the parameter for the four countries assumes the values , , and respectively. The comparison of the above optimal parameters for the four countries is performed below.

Figure 3

Theoretical (continuous curve) and empirical (dots) plots of the probability density function versus time for the Covid-19 mortality data related to Germany (top left), Italy (top right), Spain (bottom left) and United Kingdom (bottom right) 2020 pandemic. The theoretical curves are based on Eq. (26).

Figure 4 shows the cumulative distribution function versus time related to Covid-19 mortality data for Germany, Italy, Spain and United Kingdom. The theoretical curve is given by Eq. (25). The optimal values of the parameters , and for the four countries are the same as in Fig. 3. It is worth noting that the statistical fluctuations of the daily data which are evident in Fig. 3 are absorbed in the representation of the cumulative data shown in Fig. 4.

Figure 4

Theoretical (continuous curve) and empirical (dots) plots of the cumulative distribution function versus time for the Covid-19 mortality data related to Germany (top left), Italy (top right), Spain (bottom left) and United Kingdom (bottom right) 2020 pandemic. The theoretical curves are based on Eq. (25).

Figure 5 displays the survival function versus time on a log-log plot related to Covid-19 mortality data for Germany, Italy, Spain and United Kingdom. The continuous curve is the theoretical function defined in Eq. (21). The optimal values of the parameter , and for the four countries are the same as in Fig. 3. In this log-log plot, a remarkable agreement between the theoretical predictions and the empirical data clearly emerges. Furthermore, the almost linear decay of the tail of the survival function shown in this plot indicates a Pareto power-law tail for all the countries studied here, in agreement with the -deformed model.

Figure 5

Theoretical (continuous curve) and empirical (dots) plots of the survival function versus time for the Covid-19 mortality data related to Germany (top left), Italy (top right), Spain (bottom left) and United Kingdom (bottom right) 2020 pandemic. The theoretical curves are based on Eq. (21).

For the purpose of comparing the Covid-19 mortality data from China, Germany, Italy, Spain and United Kingdom we focus on the probability density function. First, we discuss the growth in mortality (i.e., the left tail) which is due to the initial phase of the virus spread. This is followed by a discussion of the final phase of the studied cycle, during which the decrease in the virus spread is described by the right tail of the curve.

The asymptotic behaviour of for (left tail) is given by where is the Weibull scale parameter. The shape parameter for China, Germany, Italy, Spain and United Kingdom has the values 3.83, 3.83, 3.90, 4.00 and 4.15 respectively. These estimates of are very close in value, with an average value equal to . It can therefore be concluded that for the five countries analyzed, the initial growing phase of the Covid-19 mortality presents quite similar dynamic behavior.

The asymptotic behaviour of for (right tail) is given by where is the pdf’s Pareto parameter (tail exponent). This parameter which governs the behavior of the right tail of the distribution takes the values 6.31, 3.28, 3.33, 3.35 and 2.83 for China, Germany, Italy, Spain and United Kingdom respectively. It can be observed that the dynamics of the final phase of the first cycle of the epidemic is essentially the same for the three continental European countries in the study. For these countries the Pareto parameter is approximately equal to . For the United Kingdom, on the other hand, the Pareto parameter is lower (equal to 2.83) indicating that the tail is thicker, i.e., that mortality is more persistent over time. On the contrary, the Paretto index is higher (equal to 6.31) for the Chinese data which indicates a faster decay of the mortality curve than the respective European curves. The differences in the Pareto parameters between the different countries are likely to reflect differences in the response strategies, in the compliance of the citizens to restrictive measures that can curb the spread of the virus, as well as in the readiness and quality of the health care system including emergency treatment units.

Discussion

Both the plague data from the pandemic of 1417 in Florence as well as the Covid-19 data of the 2020 pandemic from China, Germany, Italy, Spain and United Kingdom have been analyzed by means of the proposed -Weibull model, to obtain information about the spreading dynamics of these deadly disease outbreaks.

It is worth recalling that the Florence plague data are approximate and relatively sparse, and thus they cannot be used to reliably test a statistical model. However, the Florence data are of great historical importance as they represent the first quantitative and regularly collected record of a pandemic. The first two figures in this paper undoubtedly show that at least the cumulative Florence data are compatible with the proposed -deformed statistical model.

On the other hand, the Chinese data set on SARS-CoV-2 are important, as they represent the first statistical record for the entire evolutionary cycle of the COVID-19 pandemic in the world. The Chinese mortality data have been successfully used to perform a first validation test of the -deformed model.

It is important to note that the first and only (so far) cycle of Covid-19 ended in China with a total of 3,342 deaths. This is a relatively low mortality compared to other countries where the number of deaths was decidedly higher and mortality figures rose to a few dozen times that of China. A second noteworthy fact is that on April 17, 2020 while the Covid-19 cycle in China had already ended, the Chinese authorities reported 1290 additional deaths that occurred outside hospitals during the entire cycle without providing any information regarding their temporal distribution.

It is therefore evident that the data pertaining to the 3,342 deaths caused by Covid-19 in China are partial. For this reason, the validation of the -deformed model with complete data from other countries affected by Covid-19 is extremely important. Thus, we focused on Europe and analyzed the Covid-19 data of Italy and Spain, where the virus began to spread towards the beginning of March 2020, as well as the data from Germany and the United Kingdom, where the virus began to spread two weeks later. These data from the European countries and China differ with respect to the total number of deaths, which in the case of Germany is about 3 times the number of Chinese deaths while for Italy, Spain and United Kingdom the mortality is about 10 times that of China. Very good agreement has been obtained by fitting the -Weibull model to the data from these European countries, thereby providing further and more reliable validation of the theoretical model.

The -deformed, three-parameter model admits simple analytical forms for all the main statistical functions (probability density function, cumulative distribution function, survival function, quantile function, hazard function, cumulative hazard function) and can therefore be easily applied to Covid-19 data. Furthermore this model preserves all the universality properties already present in the original Weibull model. The results obtained by applying the -deformed model to mortality data from Covid-19 pertaining to the first cycle of the pandemic in China, Germany, Italy, Spain and United Kingdom, suggest that the -deformed model could also be used to model infection data from the first cycle of the Covid-19 pandemic as well as for the analysis of the second pandemic cycle which has already started in August 2020.