Next Article in Journal
Personalized Image Classification by Semantic Embedding and Active Learning
Next Article in Special Issue
Selection of Embedding Dimension and Delay Time in Phase Space Reconstruction via Symbolic Dynamics
Previous Article in Journal
Thermodynamic Evaluation and Sensitivity Analysis of a Novel Compressed Air Energy Storage System Incorporated with a Coal-Fired Power Plant
Previous Article in Special Issue
Classification of Actigraphy Records from Bipolar Disorder Patients Using Slope Entropy: A Feasibility Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy Ratio and Entropy Concentration Coefficient, with Application to the COVID-19 Pandemic

Institute of Mathematics, University of Greifswald, 17487 Greifswald, Germany
Submission received: 7 October 2020 / Revised: 13 November 2020 / Accepted: 14 November 2020 / Published: 18 November 2020
(This article belongs to the Special Issue Information theory and Symbolic Analysis: Theory and Applications)

Abstract

:
In order to study the spread of an epidemic over a region as a function of time, we introduce an entropy ratio U describing the uniformity of infections over various states and their districts, and an entropy concentration coefficient C = 1 U . The latter is a multiplicative version of the Kullback-Leibler distance, with values between 0 and 1. For product measures and self-similar phenomena, it does not depend on the measurement level. Hence, C is an alternative to Gini’s concentration coefficient for measures with variation on different levels. Simple examples concern population density and gross domestic product. Application to time series patterns is indicated with a Markov chain. For the Covid-19 pandemic, entropy ratios indicate a homogeneous distribution of infections and the potential of local action when compared to measures for a whole region.
MSC:
62B10; 94A17; 62P10; 28A80

1. Introduction

In any statistical data analysis, a first step is to determine two parameters that characterize center (mean, median) and variation of the data. For the Covid-19 pandemic, the mean in any region is the case incidence, given by weekly averages of case numbers per capita. How can we measure variation of infection intensity? In abstract terms, we have a basic probability measure Q, which expresses population size, and another probability measure P, which counts infections. We want to know whether P is uniformly distributed with respect to Q , or whether it is concentrated on certain hotspots.
Similar questions are frequently studied in economy, where P represents market share, wealth, or income, while the basic measure Q again expresses population size. The Gini concentration coefficient G is the standard parameter of variation [1]. It has a clear, interpretable scale, with values that are between zero and one. Zero describes totally uniform distribution, and one is approached when there are many actors, while all income or production is concentrated in the hands of one monopolist. Gini’s coefficient was used for spatial patterns in medicine [2], and it is applied here to Covid-19 case numbers. However, it does not really fit the data structure and the fractal character of the infection measure P .
Shannon entropy [3] and the relative entropy introduced by Kullback and Leibler [4] have also been used as measures of uniform distribution and concentration. Their values grow with the size of the data, which makes interpretation difficult. In 1971, their use was rejected by Hart [5], since, for the lognormal distribution, a standard model of income distribution, they lead to the ordinary variance. Fractals were not known at this time. However, already in the 1950s, Rényi [6,7] developed entropy concepts that later provided adequate description of fractal measures [8].
Here, we introduce the entropy ratio U which is related to Rényi’s entropy dimension of a measure on a metric space. While Rényi’s theoretical concept is defined as a limit, our parameter U is simply a descriptive statistics of the measures P , Q . It may vary with the size of the partition on which data of P and Q are given. Our entropy concentration coefficient C = 1 U has the same scale and interpretation as the Gini coefficient. It can be considered as a version of the Kullback–Leibler distance on multiplicative scale.
Entropy concepts indeed seem to be appropriate to describe the spread of an epidemic since there is a strong analogy with Boltzmann’s classical thermodynamical paradigm. A closed system, say molecules in a box with fixed temperature, will approach in the course of time a final state of maximum entropy, which is characterized by uniform distribution of molecules in state space. Under certain conditions, the distance between the actual distribution and the limit distribution is a relative entropy that can be interpreted as a measure of free energy or available work [9].
In the beginning of an epidemic wave, the virus itself will resist this paradigm by creating a fractal structure of hotspots. But if we let the virus do its job, spots will overlap and their number will increase so fast that uniform distribution is approached soon [10]. Of course, reality is not a closed system, and people do their best to avoid this scenario with various containment and lockdown measures. In the language of physics, external forces drive our system away from equilibrium. In this connection, the interpretation of relative entropy as available work seems to be useful for Covid-19. Local and federal authorities in many countries argue on the right type of containment. Measures for a whole country are easier to implement, while local actions can adapt better to a specific situation. The entropy concentration coefficient as a distance to uniform distribution may help to settle such arguments. For small concentration, general measures are justified. When concentration rises, local action becomes more preferable.
We shall now define our concepts, prove their basic properties and compare with Gini’s coefficient. Although Covid-19 was the motivation for our study, it turns out that entropy ratio and entropy concentration coefficient are appropriate parameters for many other phenomena. For this reason, we consider, in Section 3, two small applications to population density and gross domestic product. Section 4 and Section 5 study product spaces and show how entropy ratios are connected with entropy rate of random processes and entropy dimension of measures in Euclidean space. The entropy ratio does not change when the measures P , Q are replaced by their products P k , Q k for some k . We introduce Covid-19 data in Section 6 and study the pandemic in the USA and Germany in Section 7 and Section 8, while using data from Johns Hopkins University and Robert Koch Institut, respectively. The paper is completed with conclusions and Appendix A on the Gini coefficient.

2. Basic Concepts

Independence of intensity and concentration. We take a finite basic space Ω = { 1 , 2 , , n } . In our application, the n points are the subregions of a region, for example, states of the USA. On this space, two measures P 0 and Q 0 are given. In the example, they count the number of infections and inhabitants in each subregion. P 0 ( Ω ) and Q 0 ( Ω ) denote the total number of infections and inhabitants in the region. They are important parameters, and P 0 ( Ω ) / Q 0 ( Ω ) is the incidence, or intensity, of the disease in the region.
The parameters of uniformity and concentration should be formally independent of the intensity of the observed phenomenon. For this reason, we divide the measures P 0 , Q 0 by their total mass P 0 ( Ω ) , Q 0 ( Ω ) , respectively. We obtain two probability measures P , Q . In other words, we obtain numbers
p 1 , , p n a n d q 1 , , q n w i t h i = 1 n p i = i = 1 n q i = 1 .
When we have twice or ten times the number of cases in the same proportions, we obtain the same p i and q i . Intensity has been eliminated. In the sequel, Q will be the measure that is considered as `uniform distribution’. For this reason, we require q i > 0 for i = 1 , , n . Our standard reference method, the Gini coefficient, is described in the Appendix A of the paper.
Entropy and relative entropy. The Shannon entropy [3] of the distribution P is
H = H ( P ) = i = 1 n p i log p i .
Because 0 log 0 = 0 , all terms of H are non-negative. Thus, H 0 , and H = 0 only if p i = 1 for one i . The relative entropy of P with respect to Q was defined in [4].
D = D ( P , Q ) = i = 1 n p i log p i q i = i = 1 n p i log q i H .
Jensen’s inequality for the logarithm implies that D is always positive, except for the case P = Q where it is zero, cf. [11] (Lemma 2.3.1). For this reason, D is considered to be a kind of distance from P to Q , and is called Kullback-Leibler distance, or Kullback-Leibler divergence, or information gain [7]. Note that D is finite since we assumed q i > 0 . If we would interchange P and Q , we would get D = as soon as p i = 0 for some i (that is, no infections in a subregion).
On the one hand, D is a measure of non-uniformity or concentration of P . It measures how far the proportions of the p i deviate from those of the q i . On the other hand, D does not possess a nice bounded scale: when p i = 1 and q i = 1 / K for some i , , then D = log K . As a rule, D admits larger values when the number n of items increases.
The entropy ratio. Instead of D , we recommend to take the entropy ratio
U = U ( P , Q ) = i = 1 n p i log p i i = 1 n p i log q i = H ( P ) i = 1 n p i log q i = H H + D
as a measure of uniformity, and the entropy concentration coefficient
C = 1 U = D H + D
as a measure of non-uniformity or concentration. It turns out that C has the same properties as Gini’s concentration coefficient.
Proposition 1.
(Basic properties of entropy ratio and entropy concentration coefficient)
(i) 
0 U 1 and 0 C 1
(ii) 
U = 1 and C = 0 if and only if P = Q .
(iii) 
U = 0 and C = 1 if and only if there is an i with p i = 1 .
(iv) 
The values of U and C do not depend on the base of the logarithm.
Proof. 
(i). Because log q < 0 for 0 < q < 1 and q i > 0 for all i , the denominator of U is always positive. Thus U 0 . And U 1 follows directly from D 0 (see above). (ii). U = 1 is equivalent to D = 0 , which is equivalent to P = Q . (iii). U = 0 means H = 0 , which happens only if some p i = 1 . (lv). When we first work with natural logarithm and then change to base b , we have to divide both denominator and numerator by log b .

3. Examples

The case of two items. Let n = 2 . Setting p 1 = p and q 1 = q , a small calculation gives G = | p q | . For q = 1 2 , we obtain C = 1 p log 2 p ( 1 p ) log 2 ( 1 p ) . Figure 1 shows G and C as functions of p while q is fixed at 1 2 and 0.8 .
For p = q , we have P = Q , and both coefficients indicate concentration zero. Deviations from this minimum point will be punished stronger by Gini’s coefficient, which is what mathematicians call an L 1 -parameter, cf. [1]. Our coefficient C is an L 2 -parameter with a quadratic minimum, which assumes small values in the vicinity of the minimum point.
We see that Gini’s coefficient is not perfect in recognizing monopolists. It approaches the value 1 only if p i = 1 and also q i 0 . Our coefficient C will always be 1 for p i = 1 , , as shown in Proposition 1.
Population density in Germany. For the 16 federal states and 412 districts of Germany we take the area as q i and the population size as p i . For both states and districts, we show D , G and C in Table 1. We want to understand to which extend population is concentrated on small areas.
When compared to the world, Germany has only one large town: Berlin, with 3.6 milliion inhabitants, more than 4000 per km 2 . It is surrounded by five fairly rural states of East Germany, with a population density of between 70 and 220 people per km 2 .
All three parameters reflect this situation. Deleting Berlin from the statistics must decrease the concentration, and must do so drastically for East Germany. The decrease on the district level is smaller than on state level, since districts of Berlin are well comparable with many urban districts in Germany.
For D , all values of districts are much larger than corresponding values of states. In other words, D is biased by sample size n . Moreover, a value of 0.971 does not say anything, since the scale of D is unbounded. Accordingly, D is rather useless as a concentration measure. The bias of D is known, cf. [11] (Lemma 2.3.3), and this paper started as an attempt to standardize D , cf. Proposition 1 (v).
However, G is also biased by sample size, with values around 0.5 for districts and around 0.3 for states. According to G , the concentration of districts without Berlin is larger than concentration of states with Berlin. This is simply not true.
The values of C in Germany are nearly the same on the state and district level. In the East with Berlin, concentration of districts is clearly smaller than concentration of states, because of the size of Berlin as a state. Without Berlin, the eastern states are very similar, while the districts show a similar structure as all over Germany. On the whole, C describes the real situation much better than G . However, in agreement with Figure 1, the values of C are generally smaller than those of G, since we are fairly near to the minimum point P = Q .
Gross domestic product. For all countries with more than one million inhabitants, we took the 2020 population projection from the UN population division [12] and the 2019 gross domestic product (GDP) according to the World Bank database [13] (for a few countries, data of a previous year had to be taken). We study the concentration of the GDP among countries, and among collections of countries defined by the UN as sustainable development goal regions. Table 2 shows the Gini and entropy concentration coefficients determined from these data.
This is a typical application of the Gini coefficient. Nevertheless, G is biased by sample size. The values for the world are larger than for the continents, values for regions are smaller than corresponding values for countries. We think that C gives a better picture by showing that inequality in Asia is as large as inequality all over the world, while Africa, Europe, and America without US and Canada show a great degree of homogeneity. When the region US+Canada is included, inequality in America becomes very large, due to the immense economic power of the USA.
The values of C are again generally smaller than those of G . Hence, C may not be welcome as a measure of economic inequality. It takes small deviations from the equality point less seriously than Gini’s G , raising suspicion to downplay social differences.
In the Appendix A, we sketch a proof indicating how G depends on sample size. Actually, G was not designed for large lists of data of varying length. The standard applications either involve a small number of enterprises that share a market, or an income distribution in classes like 0–1000 $, 1000–4000 $ etc. which is expected to approximate a smooth density function. The geographical distribution of wealth, as measured in regions and subregions, is much more fragmented.

4. Product Measures

This note studies measures P , Q which are given by their values on nested partitions. Unfortunately, our data provide only two partitions—states and districts. However, in stochastic processes and time series analysis, we have an unbounded sequence of nested partitions, indexed by words of length 1,2,... over some alphabet. This setting seems to bear a large potential for the application of entropy ratios. We now introduce product measures and show a basic property of the entropy ratio.
If P is defined on Ω = { 1 , , n } , and the probability measure P is defined on Ω = { 1 , . . . , n } , then the product measure M = P × P is defined on Ω × Ω by m i j = p i p j for i in Ω and j in Ω . The k-th power P k is a probability measure on the set Ω k of all words w = i 1 i k of length k over the alphabet Ω . The measure P k assigns, to the word w, the probability p w k = p i 1 p i 2 · · p i k .
Product measures are at the core of entropy theory. When Shannon [3] characterized entropy axiomatically, the following property played a key role.
H ( P × P ) = H ( P ) + H ( P ) which implies H ( P k ) = k · H ( P ) .
Kullback and Leibler [4] (Theorem 3.1) proved the product formula for relative entropy, and Rényi [7] used it as first axiom for characterizing relative entropy:
D ( P × P , Q × Q ) = D ( P , Q ) + D ( P , Q ) which implies D ( P k , Q k ) = k · D ( P , Q ) .
Proposition 2.
(Product invariance of entropy ratio)
U ( P k , Q k ) = U ( P , Q ) a n d C ( P k , Q k ) = C ( P , Q )
for all probability measures P , Q on Ω and each positive integer k .
Proof. 
The formula directly follows from Proposition 1 (v) and (6),(7), since k cancels out. For our finite spaces Ω there is a direct proof by calculation with sums. However, the derivation from (6) and (7) works for arbitrary probability measures. □
It is custom to also consider infinite product measures P , Q on the space X = Ω of all sequences x = x 1 x 2 x k over the alphabet Ω [3,11]. The basic sets in X are the cylinder sets
C w = C y 1 . . . y k = { x 1 x 2 . . . X | x j = y j f o r j = 1 , . . , k } .
for words w = y 1 y k in Ω k . For each level k , we can consider the partition P k of X into the cylinder sets of all words of length k . These partitions form a nested sequence: each P k refines P k 1 . Proposition 2 implies
Proposition 3.
(Entropy ratio on infinite products) The measures P , Q on the infinite product space X have the same entropy ratio with respect to all partitions P k , k = 1 , 2 , . . .
It should be noted that, even in Ω 2 , the entropy ratio can change if we use other partitions. It is not allowed to split partition sets into small pieces where p i / q i is large, and leaves the cylinder sets big where p i / q i is small. This is known from everyday news: if there are ten reports from districts with more infections than average, and only one report from a district below average, we shall overestimate the danger.
There are various ways to realize infinite product measures in the plane. THe construction of self-similar measures by iterated function systems is a standard method which includes the decimal and the binary number system and can generate fancy fractals [14,15]. The above proposition says that under certain self-similarity assumptions, the entropy ratio of two measures P , Q in the plane will be the same for a sequence of nested partitions.
We do not pretend that a division into states and districts, together with the Covid-19 infection distribution, will fulfil any self-similarity conditions. However, it is natural to expect that the entropy ratios for two nested partitions coincide – in situations where the measure is fractal and the subdivision of the sets of the first partition has nothing to do with the measures P , Q . This is a kind of null hypothesis which says that P and Q show the same amount of statistical variation on different levels. In reality, there are deviations from the hypothesis, as discussed in Section 6, Section 7 and Section 8.

5. Entropy Ratios for Processes and Fractal Measures

To show connections with theory, in this section we consider random processes as in Shannon’s work [3]. Let Ω = { 0 , 1 } and q 0 = q 1 = 1 2 . Subsequently, X = Ω is the space of 01-sequences, and Q ˜ = Q models sequences of flips of a fair coin. This is the standard uniform distribution: each cylinder set of level k has measure 2 k . The denominator D + H = i = 1 n p ˜ i log q ˜ i in the entropy ratio then equals log 2 k = k · log 2 for the partition P k of X into cylinder sets of level k , for any probability measure P ˜ on X . In particular, when P ˜ = P is a product measure itself, then, for the partition P k , we have U ( P ˜ , Q ˜ ) = H ( P k ) / ( k log 2 ) , which is H ( P ) with base 2 logarithms, a well-known fact. The limit of 1 k H ( P ˜ , P k ) for k is called entropy rate [3,11]. For finite alphabets, we have
Proposition 4.
(Entropy rate from entropy ratios) Let Q be the equidistribution on a finite alphabet Ω and P ˜ a probability measure on X = Ω . Subsequently, the entropy ratios U ( P ˜ , Q ) on P k are the entropies per letter in words of length k . They converge to the entropy rate of P ˜ .
Incidentally, in the first version of this paper, we suggested to take the entropy H ( Q ) as denominator in order to standardize the relative entropy D . When all partition sets have the same Q-measure, then the quantity D / H ( Q ) agrees with U . Unfortunately, this does not hold in applications: adminstrative districts differ a lot in population size. Therefore, H ( Q ) was replaced by D + H ( P ) .
Using the correspondence between 01-sequences x = x 1 x 2 and binary numbers 0 . x 1 x 2 , the space X = { 0 , 1 } with its uniform distribution Q ˜ can be represented as the unit interval [ 0 , 1 ] with the ordinary length measure. The cylinder sets correspond to intervals [ m 2 k , ( m + 1 ) 2 k ] . The product measure P ˜ corresponds to a fractal measure π with dense support in [ 0 , 1 ] . This correspondence is one-to-one, except for a countable set, which has measure zero. We shall see now that the entropy rate of P ˜ coincides with the Hausdorff dimension of the measure π .
The entropy dimension of a finite measure μ on R d was defined by Rényi [6] as
dim ent μ = lim k 1 k d log 2 · A P k μ ( A ) log μ ( A ) ,
where P k denotes the partition of R d into d-dimensional dyadic cubes of side length 2 k . Here, we have d = 1 . Obviously, Rényi had in mind the connection between dyadic cubes and 01-words. In 1982, Young proved that the entropy dimension and Hausdorff dimension of a measure coincide if the measure is homogeneous in the sense that pointwise dimension is the same almost everywhere [8] (Section 4). Without going into further detail (see [15] for Hausdorff and pointwise dimension), we just note that on the unit cube, where the volume measure λ is a probability, the term behind the limit is the entropy ratio U ( μ , λ ) . Thus
Proposition 5.
(Entropy ratios approximate fractal dimension)
Let μ be a probability measure on [ 0 , 1 ] d which fulfils the homogeneity condition of [8]. Let λ denote the volume measure and P k the partition of the unit cube into dyadic cubes of side length 2 k . Subsequently, the entropy ratios U k ( μ , λ ) on these partitions converge to the Hausdorff dimension of μ .
It should be added that a quotient of logarithms, like U, is the typical form of a fractal dimension [15]. Great effort was invested in proving the existence of entropy rates and dimension limits [11,15]. Our point here is that, in real-world situations, limits are hard to determine, but the first terms of the sequence have a descriptive power and can be calculated.
Young’s homogeneity condition is fulfilled for the product measures and for the Markov chains that we consider now. We assume that P ˜ is generated by a Markov chain, which does not permit the word 11. The transition probabilities are p 11 = 0 , p 10 = 1 and p 01 = t , p 00 = 1 t . The stationary distribution is given as p 0 = 1 1 + t and p 1 = t 1 + t , On the interval, the corresponding measure π lives on those binary numbers that have no two consecutive digits 1. They form a Cantor set A = A 0 A 1 [ 0 , 2 3 ] . Here, A 0 , A 1 denote the subsets of A that consist of binary numbers with first digit 0 and 1, respectively. They fulfil the equations
A 0 = f 0 ( A 0 ) f 0 ( A 1 ) and A 1 = f 1 ( A 0 ) , thus A 0 = f 0 ( A 0 ) f 0 f 1 ( A 0 )
where f 0 ( x ) = x / 2 and f 1 ( x ) = ( x + 1 ) / 2 assign to 0 . x 1 x 2 the values 0 . 0 x 1 x 2 and 0 . 0 x 1 x 2 , respectively.
The similarity maps f 0 and f 0 f 1 have factors 1 2 and 1 4 , respectively. Accordingly, the Hausdorff dimension α of A 0 can be determined from the equation ( 1 2 ) α + ( 1 4 ) α = 1 [14,15]. We get α = log 2 τ 0.694 where τ 1.618 is the golden mean. For every t between 0 and 1, π is a measure on A, which has dimension smaller than or equal to α .
The measure π has maximal dimension when t = 1 / τ , or t 2 = 1 t . In this case π is the normalized Hausdorff measure on A . This is the uniform distribution on A , so it seems clear that it has the smallest distance to the uniform distribution on [ 0 , 1 ] , which is, minimal concentration 0.306 . In this case, all of the cylinder sets have measure t k for some k , and entropy ratios can be expressed with Fibonacci numbers F k as U k = ( F k + 1 + F k k + 1 k t ) t k log 2 τ . For any Markov chain, we can easily determine entropy ratios by matrix calculation, however. Table 3 shows how the concentration values C k = 1 U k approach their limit. The first step away from the uniform distribution is decisive, and the convergence slows down after the second step.
It seems interesting to determine similar distances between empirical distributions and model distributions, taking, for instance, time series of order patterns derived from high-frequency biomedical measurements, as in [16].

6. The Global Covid-19 Pandemic

We study Covid data for two reasons. On the one hand, even though they are given only on two levels, they give a good testbed, since they are determined every day. We can determine time series of concentration coefficients, and their graphs give much more information than the tables in Section 3. On the other hand, we provide a little tool for the investigation of previous waves of infection which perhaps can help to prevent the next one. We refrain from evaluating countries or political measures. That would require much more care, including study of social and economic aspects.
The numbers of positively tested cases and deaths are taken as seven-day sums per 100,000 inhabitants of the region. The data come from the database of Johns Hopkins University [17]. Weekly averages smooth out random variations of the daily numbers as well as systematic variations caused by weekends where administrations have no office hours. Sums are taken instead of averages, since seven days is the typical period for which an infected person remains infectious [18]. Thus, the sum over one week is a measure of the infection potential or Covid activity in a region. Such sums have become the basis of government decisions in Germany [19], and a main criterion for defining hotspots of the pandemic in the USA [20]. The seven-day sums are updated every day, which results in a daily time series for each country and province. See [17,21] for more details.
Figure 2 shows a global view, summarizing the values of 55 large countries, including 23 from Europe, three from North, and six from South America. The number of cases on the whole world is continuously increasing, with slow growth in August and September. The worldwide number of deaths remained almost constant at about six Covid-related deaths per one million people and week. In Europe, there was a wave of infections in March and April, and, now, there is a much larger second wave of cases. The first wave came with large numbers of deaths which lasted only for few weeks. The death rate of the second wave is now twice as large as in America. The number of cases in America increased from March to July, declined in August and is rising again. The number of deaths in America remained on a level four times as large as the world’s average for a long time, higher in July and August, and lower in October. The two rectangular peaks come from irregularities in the data, when a statistical office reported one thousand more deaths on a single day which they did not count before. We did not smooth the data. Because we take seven-day sums, irregularities of a day remain in the time series for 13 days. As a matter of fact, errors, like negative daily numbers of deaths, appear quite frequently in Europe (for instance, France, Spain, Germany) and in the USA, but not in Brazil or Iran.
Time series of concentration coefficients were calculated for cases only. For deaths, they look similar, but more distorted due to the irregularities in the data. Gini’s G and our C show a similar time course. The pandemic becomes increasingly spread over the world, but the concentration increased since October because of the second wave in Europe. In Europe, raise of concentration from April to June indicated desrease of cases in various countries, while the second wave now hits all of Europe, causing uniform distribution of infections. In America, the situation has been constant for several months, but since October the concentration rises, probably due to an improvement in Brazil and Peru and increase of cases in Argentina and the US.
At second glance, we recognize the sample size bias of G . The worldwide concentration is large because of Asian countries, like China, Vietnam, and Thailand, which have almost no cases. However, it is not as large as indicated by G . The difference between Gini concentration levels of America and Europe is too large and seems to be caused by the fact that we had 23 countries from Europe and only nine from America. The American values have less variation, when compared with the world and Europe. Therefore, the initial concentration of the pandemic in the north-east of the US in March/April is underestimated. For entropy concentration, the variation is larger and the initial concentration in New York is more correct. For the world, the initial maximum concentration (without counting Wuhan) was rather in the beginning of March, as shown by C , than in the middle of March, as indicated by G . Altogether, entropy concentration gives a more correct description than Gini’s coefficient. However, the values of C are considerably smaller than those of G .

7. States and Counties in the USA

The most comprehensive public dataset of the pandemic, maintained by the Johns Hopkins University, contains case and death numbers of 3143 counties of the 51 federal states of the USA. For Figure 3, we have chosen four states that represent the time development. New York suffered a lot in the beginning of the pandemic and has now shown the smallest numbers of cases for some months. Pennsylvania had a similar development, less dramatic in the beginning, and with more recent infections. In California and Texas, the peak of the pandemic came in summer. Curiously, all four, and a number iof other states, had approximately the same number of cases in the beginning of June, and the same number of deaths in the end of June. This is reflected by the concentration coefficients.
In March, the pandemic first concentrated near Seattle and then in New York and the northeast of the US, which results in concentration values that are near one. For this reason, C and G are shown by logarithmic scale on [ 0.022 , 1 ] and [ 0.22 , 1 ] , respectively. With this modification, it is surprising how similar C and G are as functions of time. Once again, it becomes apparent how much G is pushed by sample size. The picture of C with values of the same magnitude for states and counties seems to be correct. The concentration values for deaths in counties look almost constant since June. They are greater than corresponding values of states and greater than values for cases. Different conditions of cities and countryside could be a reason. Among others, large intensive care units are usually in cities. The two peaks in the concentration functions for deaths come from data irregularities in New York and New Jersey, respectively. Because these were spread over various counties in one state, they disturb the function of states more than the function of counties.

8. States and Districts in Germany

We use the case-related public file RKI COVID19 [22] of Robert Koch Institut, the health authority in charge of epidemic. These data allow for determining accurate case fatality rates, the distribution of daily cases in 16 federal states and 412 districts, and their joint distribution with six age classes.
Figure 4 shows the cases and deaths on a logarithmic scale. Germany is one of the European countries with mild impact of the first Covid-19 wave. The case fatality rate was about 5% in April and then decreased tremendously, below 0.5%. In June, the cases reached a minimum. Just at this time, there was a big outbreak in the main meat enterprise, which caused the big peak in concentration functions. Since July, cases have been continuously rising, and deaths followed since September.
The entropy concentration coefficients for states and districts do almost coincide since July, and they show a continuous approach to uniform distribution, with exception of the very last week. However, in spring, the concentration of districts was higher than concentration of states. This is true, because many districts with small population density recovered very fast from the first wave, so that for several weeks there were between 100 and 200 districts each day without cases. On state level, the differences were smaller. The outbreak in the meat factory appeared in the federal state with the largest population, but in a district of average size that explains why the peak is so much bigger for districts than for states.
Gini coefficients show similar features as the entropy concentration coefficients C , but their bias on sample size is obvious. In the search for concentration measures we also calculated standard deviations of case numbers or probabilities p i with respect to the weights q i They gave no reasonable result, no matter whether we took ordinary or logarithmic scale. The only variant worth mentioning is given by coefficients of variation of the case numbers (for the p i the parameter is the same). That is, the standard deviation divided by mean. Figure 4 shows that they admit no natural scale and have a bias. Actually they do not show much change in time, except for the meat factory outbreak. The parameter C clearly performs the best, and our null hypothesis of equal concentration for the two partitions is confirmed by the remarkable coincidence of the two time series.
Because Covid deaths are the most frequent among old people, it is very important to study the age dependence of cases, as shown on the left of Figure 5. There are four age classes: children below 15, the working population in the age of 15 up to 59, the elderly with 60 to 79 years, who have case fatality rate of about 2%, and the very vulnerable group of age 80 and above, who now have a case fatality of about 12% (in spring, fatality rates were much higher since even among old people, only cases with serious symptoms were tested). The last group, represented by nearly six-million inhabitants of Germany, needs the most protection.
Until August, the seven-day incidence of people above 80 was very low, only 2 per 100,000. In September, the number raised to eight and, on November 10, it is 138. Simple calculation shows that this is more or less the number of deaths in this group, which we can expect all over Germany two weeks later. Approximately 60% of all Covid deaths come from this age class. So it is possible to predict deaths in Germany two weeks in advance. It should be noted that the rise of infections of old people is not due to their careless behavior or lack of protection. More likely, the number of infections have risen in all age classes, but younger people do not easily get tested anymore.
For our study, it is interesting to consider the concentration functions C in different age classes. We use the districts that give more information. It can be seen that cases in the group of working age 15–59 are most uniformly spread, certainly due to their mobility. The outbreak in the meat factory is almost restricted to this group. The infections in the group of children and the group of elderly above 80 were more concentrated in certain districts. However, with the progress of the second wave, they have also become uniformly spread. While local action was perhaps still good in September, it is quite reasonable that a partial lockdown for all Germany took effect on November 2. There is an upwards trend in the most recent concentration coefficients, which indicates that, like in March, districts with a small population density are already improving. However, incidences are still growing.

9. Conclusions

We intended to describe the spread of the Covid-19 pandemic, given by daily numbers of confirmed cases and deaths in various regions. The question is how uniform the distribution of infections is with respect to the distribution of population. The Kullback–Leibler distance can be applied, but it is biased by sample size and has an unbounded scale. The well-established Gini coefficient can be used and it has a unit scale, but it is strongly influenced by sample size. It was not designed for data sampled in many regions and districts. When countries are compared by Gini coefficients, with respect to the inequality of income say, then a smooth income distribution is calculated for every country. The Covid data structure of cases and deaths reflects the spatial distribution of infections, which is more fragmented. Many variables can be sampled in this way: air and water pollution, precipitation, wealth, crime, etc.
An entropy concentration coefficient C was introduced. It has the same scale and basic properties as the Gini coefficient and it seems to be much less dependent on sample size. The new coefficient is a multiplicative version of the Kullback–Leibler divergence and can be interpreted as distance of a probability measure P from a uniform distribution. When the uniform distribution is volume in Euclidean space, and P is homogeneous, then U approximates the Hausdorff dimension of P . When the uniform distribution is the product measure of an equidistribution on a finite alphabet, then U = 1 C is connected with the entropy rate. In the setting of random processes and empirical distributions of long time series, we anticipate applications of C .
We analysed data of the global pandemic, of counties in the US and of districts in Germany. The concentration coefficient C as a function of time becomes small when infections become widespread. Its growth can indicate improvement, sometimes earlier than the decline of cases. Large values of C also appear in the initial stage of an epidemic, and for big local outbreaks. In Germany, age groups—in particular, the vulnerable group of elderly above 80—develop differently with respect to concentration. A referee suggested that cross-correlation between countries and age classes could be further studied by forecasting and backcasting analysis, while using detrended cross-correlation coefficients, as in [23]. All in all, the entropy concentration coefficient is an interesting new parameter.

Funding

This research received no external funding.

Acknowledgments

The author is grateful to all three referees for their constructive criticism of the first version of the paper. He thanks Mariano Matilla for the invitation to join the special issue on Information theory and Symbolic Analysis, and Bernd Pompe for a number of helpful discussions.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. The Gini Coefficient

Definition of G . Given the p i , q i in (1), we renumber the n items in such a way that p 1 / q 1 < p 2 / q 2 < . . . < p n / q n . For an epidemic, this means ordering subregions by increasing incidence. To save notation, we assume the renumbering has already been done. Now, take the cumulative probabilities x k = i = 1 k q i and y k = i = 1 k p i for k = 1 , . . . , n . Let ( x 0 , y 0 ) = ( 0 , 0 ) and note that ( x n , y n ) = ( 1 , 1 ) . The polygonal line which connects the successive points ( x k , y k ) , k = 0 , 1 , . . . , n in the x y -plane is called Lorenz curve. Due to the assumption of increasing quotients, the Lorenz curve runs below the line y = x . The Gini coefficient G is twice the area between Lorenz curve and line y = x . For details, see any statistics textbook for economists or [1].
Now, we show that the Gini coefficient increases when the partition determining the p i , q i is refined. It is enough to prove this for a basic step of refinement, where one partition set is subdivided into two. It would be interesting to modify the geometric proof so that under stochastic model assumptions the increase can be quantified. The area described below can easily be calculated.
Proposition A1.
(Gini coefficient increases under refinement)
Let G be the Gini coefficient for probabilities p i , q i with p i / q i < p i + 1 / q i + 1 for i = 1 , . . . , n . Now suppose that one pair of probabilties is splitted: p k = p k + p k and q k = q k + q k . Then the new Gini coefficient G is larger than G , except for the case p k / q k = p k / q k where G = G .
Proof. 
We can assume p k / q k p k / q k . For the Lorenz curve connecting the points ( x i , y i ) defined above, we add one more vertex ( x , y ) with x = i < k q i + q k and y = i < k p i + p k . When p k / q k = p k / q k then ( x , y ) is on the Lorenz curve and G = G . Otherwise, ( x , y ) is below the Lorenz curve, and connecting the points with increasing x , the area between the new Lorenz curve and the line y = x is enlarged by the triangle with vertices A = ( x k 1 , y k 1 ) , B = ( x , y ) , and C = ( x k , y k ) .
However, in that case, we have to reorder the points, which will lead to even more enlargement. To save notation, we only consider one step of reordering. We assume that p k 1 / q k 1 > p k / q k so that k and k 1 must interchange order. In that case the point A of the Lorenz curve is replaced by A = ( x k 2 , y k 2 ) + A B , and the parallelogram with vertices A , B , A , ( x k 2 , y k 2 ) will be added to the area between Lorenz curve and line y = x .
When k has to be arranged still further down, then more of these parallelograms will be added. The same holds for ordering k above k . Those parallelograms are constructed with vector C B instead of A B . The more p k / q k and p k / q k differ, the more parallelograms have to be added in order to obtain G from G .

References

  1. Wikipedia. Gini Coefficient. Available online: https://en.wikipedia.org/wiki/Gini_coefficient (accessed on 6 October 2020).
  2. Brown, M. Using Gini-Style indices to evaluate the spatial patterns of health practitioners. Soc. Sci. Med. 1994, 38, 1243–1256. [Google Scholar] [CrossRef]
  3. Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 379–423, 623–656. [Google Scholar] [CrossRef] [Green Version]
  4. Kullback, S.; Leibler, R.A. On information and sufficiency. Ann. Math. Stat. 1951, 22, 79–86. [Google Scholar] [CrossRef]
  5. Hart, P.E. Entropy and other measures of concentration. J. R. Stat. Soc. Ser. 1971, 134, 73–85. [Google Scholar] [CrossRef]
  6. Rényi, A. Dimension, entropy and information. In Transactions of the second Prague Conference on Information Theory, Statistical Decision Functions, Random Processes; Czech Academy of Sciences: Prague, Czech Republic, 1957; pp. 545–556. [Google Scholar]
  7. Rényi, A. Wahrscheinlichkeitsrechnung mit einem Anhang über Informationstheorie; VEB Deutscher Verlag der Wissenschaften: Berlin, Germany, 1962; English translation: Probability Theory, North Holland, The Netherlands, 1970. [Google Scholar]
  8. Young, L.S. Dimension, entropy and Lyapunov exponents. Ergod. Theory Dyn. Syst. 1982, 2, 109–124. [Google Scholar] [CrossRef] [Green Version]
  9. Müller, I.; Weiss, W. Entropy and Energy; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  10. Ferguson, N.M.; Laydon, D.; Nedjati-Gilani, G.; Imai, N.; Ainslie, K.; Baguelin, M.; Bhatia, S.; Boonyasiri, A.; Cucunuba, P.Z.; Cuomo-Dannenburg, G.; et al. Impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand. Imp. Coll. Lond. 2020. [Google Scholar] [CrossRef]
  11. Gray, R.M. Entropy and Information Theory; Springer: Berlin/Heidelberg, Germany, 1990. [Google Scholar]
  12. United Nations, Department of Economic and Social Affairs, Population Division. World Population Prospects, Standard Projections. Available online: https://population.un.org/wpp/Download/Standard/Population/ (accessed on 30 October 2020).
  13. The World Bank. World Bank Open Data, Gross Domestic Product. Available online: https://data.worldbank.org/indicator/NY.GDP.MKTP.CD (accessed on 30 October 2020).
  14. Barnsley, M.F. Fractals Everywhere, 2nd ed.; Academic Press: Cambridge, MA, USA, 1993. [Google Scholar]
  15. Falconer, K.J. Fractal Geometry: Mathematical Foundations and Applications, 3rd ed.; J. Wiley & Sons: Hoboken, NJ, USA, 2014. [Google Scholar]
  16. Bandt, C. A new kind of permutation entropy used to classify sleep stages from invisible EEG microstructure. Entropy 2017, 19, 197. [Google Scholar] [CrossRef] [Green Version]
  17. Johns Hopkins University. Coronavirus COVID-19 Global Cases by the Center for Systems Science and Engineering. Available online: https://gisanddata.maps.arcgis.com/apps/opsdashboard/index.html;https://github.com/CSSEGISandData/COVID-19/tree/master/csse_covid_19_data/csse_covid_19_time_series (accessed on 10 November 2020).
  18. World Health Organization. Criteria for Releasing COVID-19 Patients from Isolation. Available online: https://www.who.int/news-room/commentaries/detail/criteria-for-releasing-covid-19-patients-from-isolation (accessed on 6 November 2020).
  19. Robert Koch Institut. Daily Situation Reports. Available online: https://www.rki.de/EN/Home/homepage_node.html (accessed on 10 November 2020).
  20. Oster, A.M.; Kang, G.J.; Cha, A.E.; Beresovsky, V.; Rose, C.E.; Rainisch, G.; Porter, L.; Valverde, E.E.; Peterson, E.B.; Driscoll, A.K.; et al. Trends in Number and Distribution of COVID-19 Hotspot Counties in the United States, March 8 to July 15, 2020. In Morbidity and Mortality Weekly Reports; Centers for Disease Control and Prevention: Atlanta, GA, USA, 21 August 2020; p. 69. [Google Scholar]
  21. Roser, M.; Ritchie, H.; Ortiz-Ospina, E.; Hasell, J. Coronavirus Pandemic (Covid-19). Available online: https://ourworldindata.org/coronavirus (accessed on 10 November 2020).
  22. Robert Koch Institut. RKI COVID19 Database. Available online: https://npgeo-corona-npgeo-de.hub.arcgis.com/datasets/dd4580c810204019a7b8eb3e0b329dd6_0 (accessed on 10 November 2020).
  23. Contreras-Reyes, J.E.; Idrovo-Aguirre, B.J. Backcasting and forecasting time series using detrended cross-correlation analysis. Phys. A 2020, 560, 125109. [Google Scholar] [CrossRef]
Figure 1. Gini coefficient and entropy concentration coefficient as functions of p = p 1 for fixed values of q . Left: q = 0.5 . Right: q = 0.8 .
Figure 1. Gini coefficient and entropy concentration coefficient as functions of p = p 1 for fixed values of q . Left: q = 0.5 . Right: q = 0.8 .
Entropy 22 01315 g001
Figure 2. Time series of the pandemic. The upper row describes the intensity of the pandemic, with cases on the left and deaths on the right. The lower row specifies the non-uniformity of the distribution of cases, with Gini coefficients on the left and entropy concentration coefficients on the right. Data from Johns Hopkins University [17], updated 10 November 2020.
Figure 2. Time series of the pandemic. The upper row describes the intensity of the pandemic, with cases on the left and deaths on the right. The lower row specifies the non-uniformity of the distribution of cases, with Gini coefficients on the left and entropy concentration coefficients on the right. Data from Johns Hopkins University [17], updated 10 November 2020.
Entropy 22 01315 g002
Figure 3. Upper row: cases and deaths for four representative states of the USA. Middle row: entropy concentration coefficients for the partitions into 51 states and 3143 counties on logarithmic scale. Lower row: Gini coefficients are very similar, but larger and not adapted to sample size. Data from Johns Hopkins University [17], 10 November 2020.
Figure 3. Upper row: cases and deaths for four representative states of the USA. Middle row: entropy concentration coefficients for the partitions into 51 states and 3143 counties on logarithmic scale. Lower row: Gini coefficients are very similar, but larger and not adapted to sample size. Data from Johns Hopkins University [17], 10 November 2020.
Entropy 22 01315 g003
Figure 4. Upper left: cases and deaths in Germany. Upper right: entropy concentration coefficients for the 16 federal states and 412 districts. Lower row: Gini coefficient and coefficient of variation are biased by sample size. Data from Robert Koch Institut [22], updated 10 Nov 2020.
Figure 4. Upper left: cases and deaths in Germany. Upper right: entropy concentration coefficients for the 16 federal states and 412 districts. Lower row: Gini coefficient and coefficient of variation are biased by sample size. Data from Robert Koch Institut [22], updated 10 Nov 2020.
Entropy 22 01315 g004
Figure 5. Left: case incidence with respect to age. Right: entropy concentration coefficients of cases for age groups. These figures were updated in proof on 16 November since the most recent concentration coefficients already indicate improvement in part of the country while case incidences are still rising.
Figure 5. Left: case incidence with respect to age. Right: entropy concentration coefficients of cases for age groups. These figures were updated in proof on 16 November since the most recent concentration coefficients already indicate improvement in part of the country while case incidences are still rising.
Entropy 22 01315 g005
Table 1. Concentration of population in German states and districts.
Table 1. Concentration of population in German states and districts.
ParameterGermanywithout BerlinEast with BerlinEast without Berlin
D for states0.2600.1840.6180.089
D for districts0.6500.5800.9710.469
G for states0.3360.3090.3900.227
G for districts0.5280.5090.5390.414
C for states0.0980.0730.2620.054
C for districts0.1010.0920.1830.101
Table 2. Concentration of gross domestic product in countries and regions.
Table 2. Concentration of gross domestic product in countries and regions.
ParameterWorldAsiaAfricaEuropeAmericawithout US, Ca
number of regions/countries19/1595/475/494/354/253/23
G for regions0.5750.3790.2910.2820.4410.012
G for countries0.6210.4960.4120.3400.4710.190
C for regions0.2180.2100.0920.1020.3820.0004
C for countries0.1820.1770.0880.0670.2690.034
Table 3. Concentration of Markov Fibonacci measure on words of length k.
Table 3. Concentration of Markov Fibonacci measure on words of length k.
k123456
C k 0.04060.23820.25880.26710.27390.27830.30576
C k C k 1 0.19770.02060.00830.00680.0044
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bandt, C. Entropy Ratio and Entropy Concentration Coefficient, with Application to the COVID-19 Pandemic. Entropy 2020, 22, 1315. https://0-doi-org.brum.beds.ac.uk/10.3390/e22111315

AMA Style

Bandt C. Entropy Ratio and Entropy Concentration Coefficient, with Application to the COVID-19 Pandemic. Entropy. 2020; 22(11):1315. https://0-doi-org.brum.beds.ac.uk/10.3390/e22111315

Chicago/Turabian Style

Bandt, Christoph. 2020. "Entropy Ratio and Entropy Concentration Coefficient, with Application to the COVID-19 Pandemic" Entropy 22, no. 11: 1315. https://0-doi-org.brum.beds.ac.uk/10.3390/e22111315

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop