# The Weighted Gini-Simpson Index: Revitalizing an Old Index of Biodiversity

Authors: Radu Cornel Guiasu; Silviu Guiasu
Journal: International Journal of Ecology
 (2012)

---

## Abstract

The distribution of biodiversity at multiple sites of a region has been traditionally investigated through the additive partitioning of the regional biodiversity into the average within-site biodiversity and the biodiversity among sites. The standard additive partitioning of diversity requires the use of a measure of diversity, which is a concave function of the relative abundance of species, such as the Gini-Simpson index, for instance. Recently, it was noticed that the widely used Gini-Simpson index does not behave well when the number of species is very large. The objective of this paper is to show that the new weighted Gini-Simpson index preserves the qualities of the classic Gini-Simpson index and behaves very well when the number of species is large. The weights allow us to take into account the abundance of species, the phylogenetic distance between species, and the conservation values of species. This measure may also be generalized to pairs of species and, unlike Rao’s index, this measure proves to be a concave function of the joint distribution of the relative abundance of species, being suitable for use in the additive partitioning of biodiversity. The weighted Gini-Simpson index may be easily transformed for use in the multiplicative partitioning of biodiversity as well.


---

## Body

## 1. Introduction

Measuring biodiversity is a major, and much debated, topic in ecology and conservation biology. The simplest measure of biodiversity is the number of species from a given community, habitat, or site. Obviously, this ignores how many individuals each species has. The best known measures of biodiversity, that also take into account the relative abundance of species, are the Gini-Simpson indexGS and the Shannon entropy H. Both measures have been imported into biology from other fields Thus, Gini [1] introduced his formula in statistics, in 1912. Much later, after 37 years, Simpson [2] pleaded convincingly in favour of using Gini’s formula as a measure of biodiversity. Shannon was an engineer who introduced his discrete entropy in information theory [3], in 1948, as a measure of uncertainty, inspired by Boltzmann’s continuous entropy from classical statistical mechanics [4], defined half a century earlier. Shannon’s formula was adopted by biologists about 17 years later [5–8], as a measure of specific diversity. This import of mathematical formulas has continued. Rényi, a probabilist, introduced his own entropy [9], in order to unify several generalizations of the Shannon entropy. He was a pure mathematician without any interest in applications, but later, Hill [10] claimed that by taking the exponential of Rényi’s entropy we obtain a class of suitable measures of biodiversity, called Hill’s numbers, which were praised by Jost [11–13] as being the “true” measures of biodiversity. In 1982, Rao [14], a statistician, introduced the so-called quadratic entropy R, which in fact has nothing to do with the proper entropy and depends not only on the relative abundance of species but also on the phylogenetic distance between species. This function has also been quickly adopted by biologists as a measure of dissimilarity between the pairs of species. In the last 20 years, a lot of other measures of diversity have been proposed. According to Ricotta [15], there is currently a “jungle” of biological measures of diversity. However, as mentioned by S. Hoffmann and A. Hoffmann [16], there is no unique “true” measure of diversity.Starting with MacArthur [5], MacArthur and Wilson [7], and Whittaker [8], the distribution of biodiversity at multiple sites of a region has been traditionally investigated through the partitioning of the regional or total biodiversity, called γ-diversity, into the average within-site biodiversity, called α-diversity, and the between-site biodiversity or diversity turnover, called β-diversity. All these diversities, namely, α-diversity, β-diversity, and γ-diversity, should be nonnegative numbers. Unlike α-diversity and γ-diversity, there is no consensus about how to interpret and calculate β-diversity. According to Whittaker [8], who introduced the terminology, β-diversity is the ratio between γ-diversity and α-diversity. This is the multiplicative partitioning of diversity. According to MacArthur [5], Lande [17], and, more recently, Veech at al. [18], β-diversity is the difference between γ-diversity and α-diversity. This is the additive partitioning of diversity.Let us assume that in a certain region there aren species, m sites, and θk is the distribution of the relative abundance of species at site k. Let λk>0be an arbitrary parameter assigned to site k, such that λ1+⋯+λm=1. These parameters may be used to make adjustments for differences (in size, altitude, etc.) between the sites. If no adjustment is made, we take these parameters to be equal, that is, λk=1/m, for every k. If μ is a nonnegative measure of diversity, which assigns a nonnegative number to each distribution of the relative abundance of the n species, then the corresponding γ-diversity is γ=μ(∑kλkθk) and the α-diversity is α=∑kλkμ(θk). The β-diversity is taken to be β=γ-α, in the additive partitioning of diversity, and β=γ/α, in the multiplicative partitioning of diversity. In general, a measure of biodiversity ought to be nonnegative, in which case the corresponding α-diversity and γ-diversity calculated by using such a measure are also nonnegative, as they should be. From a systemic point of view, the β-diversity shows to what extent the total, or regional diversity differs from the average diversity of the communities/habitats/sites taken together, as a system, reflecting the dissimilarity, or differentiation between communities/habitats/sites of the region with respect to the individual species. If the measure of biodiversity is a concave function of the distribution of the relative abundance of species θk, then the corresponding β-diversity is β≥0, in the additive partitioning of diversity, and β≥1, in the multiplicative partitioning of diversity, for arbitrary parameters λk>0,λ1+⋯+λm=1. If a measure of diversity μ is not a concave function of the distribution of the relative abundance of species θk, then the corresponding β-diversity could be negative, in the additive partitioning of biodiversity, or less than 1, in the multiplicative partitioning of biodiversity, for some parameters satisfying λk>0,λ1+⋯+λm=1, which is absurd. As discussed by Jost [11, 12], if a measure of diversity μ is not a concave function of the distribution of the relative abundance of species θk, we can still attempt the partitioning of biodiversity into α-, β-, and γ-diversity if a new kind of α-diversity may be introduced. This new type of α-diversity would be based on a different way of averaging the diversities of the individual communities/habitats/sites instead of the simple, golden mean value α=∑kλkμ(θk) from statistics, which works so well for the concave measures of diversity. However, finding such an unorthodox, nonstandard α-diversity when the measure of biodiversity μ is not concave is not easy. It is also difficult to find a mathematical interpretation for such a new kind of α-diversity. In spite of the passage of time, the most popular measures of biodiversity are still GS, H, and R. Both GS and H are concave functions of the distribution of the relative abundance of species and therefore can be used for doing the additive partitioning of biodiversity. The first two Hill’s numbers are mathematical transformations of GS and H, namely, Hl1=1/(1-GS) and Hl2=exp(H), and are used in the multiplicative partitioning of biodiversity. Recently, however, it was noted [13, 19] that both Shannon’s entropy and the Gini-Simpson index do not behave well when the number of species is very large. On the other hand, when a distance between species, such as the phylogenetic distance for instance, is also taken into account along with the relative abundance of species, Rao’s index [14] is a widely used measure of dissimilarity. But, unfortunately, Rao’s index is not a concave function of the distribution of the relative abundance of species, for an arbitrary distance matrix between species. Consequently, it proves to be suitable for use in the standard additive partitioning of diversity only in some special cases, but not in general. The objective of this paper is to show that the weighted Gini-Simpson quadratic index, a generalization of the classic Gini-Simpson index of biodiversity, offers a solution to both of the drawbacks just mentioned. Unlike Shannon’s entropy and the classic Gini-Simpson index, this new weighted measure of biodiversity behaves very well even if the number of species is very large. The weights allow us to measure biodiversity when a distance between species and/or conservation values of the species are taken into account, along with the abundance of species. When the phylogenetic distance between species is taken as the weight, the corresponding weighted Gini-Simpson index, unlike Rao’s index, is a concave function of the distribution of the relative abundance of the pairs of species, being suitable for use in the additive partitioning of biodiversity. A simple algebraic transformation makes the weighted Gini-Simpson index suitable for use in the multiplicative partitioning of biodiversity as well.In Methodology, the weighted Gini-Simpson quadratic index is defined both for individual species and for pairs of species. This new measure of biodiversity is used for calculating the average within-site biodiversity (α-diversity), the intersite biodiversity (β-diversity), and the regional or total biodiversity (γ-diversity). It is also shown that the weighted Gini-Simpson quadratic index may be easily modified, by a simple algebraic transformation, to get a measure of biodiversity suitable for use in the multiplicative partitioning of biodiversity as well. In Section 3, a numerical example is presented, which illustrates how the mathematical formalism should be applied from a practical standpoint.


## 2. Methodology

### 2.1. The Weighted Measure of Diversity with Respect to Individual Species

Let us assume that there aren species in a certain community/habitat/site and let pi be the relative abundance of species i (the number of individuals of species i divided by the total number of individuals in that community/habitat/site). We have diversity if species i is present at that location but other species are found there as well. The probability that the species i is present and there are other species present as well is pi(1-pi). If we take all possible values of pi from the unit interval [0,1] into account, the wave function pi(1-pi) corresponding to the species i is a nonnegative, symmetric, bell-shaped, concave function, reaching its maximum value 1/4 at pi=1/2. If we sum up these wave functions, for all n species, we obtain the classic Gini-Simpson index GS(θ) corresponding to the given distribution of the relative abundance of the species θ=(p1,…,pn). Since 1949, this has been considered to be a very good measure of biodiversity. In order to generalize it, we may assign an amplitude wi≥0 to the wave function of the species i, and the resulting new wave function wipi(1-pi) continues to be a nonnegative, symmetric, bell-shaped, concave function of pi, but this time its maximum value is wi/4. Summing up these wave functions for all the species, we get the weighted Gini-Simpson index:(1)GSw(θ)=∑iwipi(1-pi),
which depends both on the distribution of the relative abundance of species θ and on the nonnegative weights w=(w1,…,wn). The concavity of GSw(θ) was proven in [20, 21]. The weight wi could be anything which contributes to the increase in the diversity induced by the species i. However, the weights may not depend on the relative abundance of species. If wi=n, for each i, then (1) becomes the so-called Rich-Gini-Simpson index GSn(θ),introduced in [22], which is essentially dependent on the species richness of the respective community/habitat/site. If there are some conservation values assigned to the species v=(v1,…,vn), which are positive numbers on a certain scale of values, and the weights are wi=nvi, the corresponding weighted Gini-Simpson index is denoted by GSn,v(θ). Obviously, if wi=1, for each species i, then (1) is the classic Gini-Simpson index GS(θ). An upper bound for GSw(θ), which depends only on the maximum weight and the number of species, is(2)0≤GSw(θ)≤(max⁡i⁡wi)∑ipi(1-pi)≤(max⁡i⁡wi)(1-1n).
Denoting by B1 the bound from the right-hand side of the inequality (2), the relative weighted Gini-Simpson index for individual species is 0≤GSw(θ)/B1≤1. In Appendix A, another bound of GSw(θ) is given, denoted by B2, which depends on all the weights assigned to the species. If wi=1, for each species i, we get maxθGS(θ)=1-1/n, and this maximum biodiversity is obtained when all species have the same relative abundance pi=1/n. The fact that the maximum value of GS(θ) is almost insensitive to the increase of the number of species, tending very slowly to 1 when n increases, allowed Jost [12, 13] and Jost et al. [19] to give some examples showing that the Gini-Simpson index does not behave well when the number of species n is very large. R. C. Guiasu and S. Guiasu [22] showed, however, that the Rich-Gini-Simpson index GSn(θ) has no such problem. Indeed, if we take wi=n, for each species i, we get from (2), maxθGSn(θ)=n-1, whose value sensibly increases when the number of species increases, which makes inapplicable the criticism of the classic Gini-Simpson index.


### 2.2. The Additive Partitioning of Biodiversity with Respect to the Individual Species

Let us assume that in a certain region there aren species and m sites. In what follows, the subscripts i and j refer to species (i,j=1,…,n) and the subscripts k and r refer to sites, (k,r=1,…,m). Let θk=(p1,k,…,pn,k) be the vector whose components are the relative abundances of the individual species at site k, such that pi,k≥0,(i=1,…,n),∑ipi,k=1, for each k=1,…,m. Let w=(w1,…,wn) be nonnegative weights assigned to the species. In dealing with species diversity, a good measure of the differentiation, or dissimilarity, among the sites in a certain region has to be nonnegative and equal to zero if and only if there is no such difference. We assign a parameter λk to each site k, such that(3)λk≥0,(k=1,…,m),∑kλk=1.
These parameters may be used to make adjustments for differences (in size, altitude, etc.) between the sites, as shown in [23]. If no adjustment is made and we focus only on the species abundance, we take these parameters to be equal, that is, λk=1/m,for every (k=1,…,m). As GSw(θ) is a concave function of the distribution of the relative abundance θ, it may be used in the additive partitioning of biodiversity. The corresponding γ-diversity, reflecting the total or regional biodiversity, the α-diversity, interpreted as the within-site diversity or the average diversity of the sites, and the β-diversity, as a measure of between-site diversity, are given by(4)γ=GSw(∑kλkθk),α=∑kλkGSw(θk),β=γ-α.Theβ-diversity may be interpreted as a measure of dissimilarity or differentiation between the sites of the respective region with respect to the individual species. As shown in Appendix C, taking into account (4), the β-diversity has the expression(5)β=∑iwi∑k<rλkλr(pi,k-pi,r)2<2∑iwi∑k<rλkλr.
The upper bound in (5) is quite loose. As α is nonnegative and cannot exceed γ, a better upper bound for β is the value B2 from (A.1), corresponding to n and the given weights w. Thus, the relative between-site diversity, or the dissimilarity between sites is 0≤β/B2≤1, with respect to the individual species. From (5), we can see that if the species have the same abundance in each site, which means that pi,k=pi, for each site k and each species i, the β-diversity is equal to zero, reflecting the fact that in such a case there is no dissimilarity between the sites.


### 2.3. The Multiplicative Partitioning of Biodiversity with Respect to Individual Species

Dealing with the multiplicative partitioning of diversity, Whittaker [8] suggested the use of the exponential of the Shannon entropy as a measure of biodiversity. The weighted Gini-Simpson index GSw(θ), given by (1), which can be used in the additive partitioning of diversity induced by individual species, may also be transformed into the measure of biodiversity (R. C. Guiasu and S. Guiasu [21]):(6)1[∑iwipi-GSw(θ)]=(∑iwipi2)-1,
which can be used in the multiplicative partitioning of diversity induced by the individual species. This measure of biodiversity may be viewed as being the weighted version of the classic Hill number of first degree from [10]. The corresponding multiplicative γ-diversity and α-diversity are(7)γ=[∑iwi(∑kλkpi,k)2]-1,α=[∑kλk∑iwipi,k2]-1.
Due to the convexity of the function ∑iwipi2, as a function of the distribution of the relative abundance of species, the γ-diversity cannot be smaller than the α-diversity, as it should be, and, consequently, the multiplicative β-diversity satisfies the inequality: β=γ/α≥1. In the additive partitioning of diversity, the γ-diversity, α-diversity, and β-diversity are entities of the same kind and may be expressed in the same units. In the multiplicative partitioning of diversity, the β-diversity is simply a ratio between the total, regional biodiversity γ and the average within-site biodiversity α, a numerical indicator showing to what extent the regional biodiversity, as a whole, exceeds the average biodiversities of the sites of the respective region. Obviously, if the sites have the same species and the same abundance of these species, which means that pi,k=pi, for each species i, then β=1.


### 2.4. The Weighted Measure of Diversity with Respect to the Pairs of Species

LetD=[dij] be an n×n matrix whose entries are the distances between the pairs of n species, such that dij≥0,dii=0,(i,j=1,…,n). This could be the matrix of the phylogenetic distance between species, for instance. When Rao introduced his quadratic index, improperly called quadratic entropy [14]:(8)RD=∑i,jdijpipj,
he had to focus on the pairs of species instead of the individual species, using the distance between the pairs of distinct species, along with their relative abundance, in order to measure the dissimilarity between species. Rao’s indicator is very simple and may be easily interpreted as the average dissimilarity between two individuals belonging to two different species when the phylogenetic distance is taken into account. There have been numerous attempts, such as [24], for instance, at using Rao’s index RD in the additive partitioning of diversity. Unfortunately, RD is not a concave function of the distribution of the relative abundance of species θ=(p1,…,pn) for an arbitrary distance matrix D. Thus, RD can be applied to the additive partitioning of biodiversity only for some special kinds of such matrices D (as mentioned in [24], for instance), but not in general. On the other hand, there is no generally accepted proposal of a new kind of nonstandard α-divergence which could be defined for such a measure which is not a concave function of the distribution of the relative abundance of species. This paper shows that the generalization of the weighted Gini-Simpson index to the pairs of species provides a concave measure of diversity which could indeed be used both for the additive partitioning and the multiplicative partitioning of diversity when the phylogenetic distance between species is taken into account. Therefore, this provides a suitable replacement of Rao’s index in the partitioning of biodiversity.LetΘ=[πij] be an n×n matrix where πij is a joint probability of the pair of species (i,j), in this order. As πji is the probability of the pair (j,i), in this order, the probability of the subset of species {i,j} is πij+πji. We haveπij≥0,(i,j=1,…,n);∑i,jπij=1. Let W=[wij] be an n×n matrix whose entries are arbitrary nonnegative weights assigned to the pairs of species. However, these weights may not depend on the joint distribution [πij]. The weighted Gini-Simpson quadratic index of the pairs of species is(9)GSW(Θ)=∑i,jwijπij(1-πij),
If wij=1, for all pairs of species, then (9) becomes the generalization of the classic Gini-Simpson index to the pairs of species and is denoted by GS(Θ). As a function of the joint distribution Θ, the weighted Gini-Simpson index GSW(Θ) is nonnegative and concave, as shown in [20, 21]. An upper bound for GSW(Θ) which depends only on the maximum weight and the number of species is(10)0≤GSW(Θ)≤(max⁡i,j⁡wij)∑i,jπij(1-πij)≤(max⁡i,j⁡wij)(1-1n2).
Denoting by B3 the bound from the right-hand side of the inequality (10), the relative weighted Gini-Simpson index for pairs of species is 0≤GSW(Θ)/B3≤1. In Appendix B, another bound for GSW(Θ), denoted by B5, is given, which depends on all the weights assigned to the species. If wij=1, for each pair of species (i,j), we get maxΘGS(Θ)=1-1/n2, and this maximum biodiversity is obtained when all species have the same relative abundance πij=1/n2. The special cases of interest are the following ones.(a)
If the species are independent, which meansπij=pipj and the weights are wij=dij, then the weighted Gini-Simpson index (9) is denoted by
(11)GSD(Θ)=∑i,jdijpipj(1-pipj),
which generalizes Rao’s index RD given by (8).(b)
If there are the positive numbersv=(v1,…,vn), representing conservation values of the individual species, the species are independent, which means πij=pipj, and the weights are
(12)wij=n(n-1)212(vi+vj)dij,
where n(n-1)/2 is the number of distinct pairs of species (i,j), such that i<j, or the number of the pairs of species {i,j}, and (1/2)(vi+vj) is the average value of the pair of species (i,j), then the corresponding weighted Gini-Simpson index (9) is denoted by:
(13)GSn,v,D(Θ)=n(n-1)2∑i<j(vi+vj)dijpipj(1-pipj),
which takes into account all the information available, namely, the species richness n, the relative abundance θ of species, the matrix D of the distance between species, and the conservation values v of the species. As the measure given in (13) is a nonnegative concave function of the distribution of the relative abundance of the pairs of species Θ, for an arbitrary distance matrix D, and also depends explicitly on the species richness, the distance between species, and the conservation values of the species, all these are sufficient reasons to suggest that it could more than adequately replace the use of Rao’s index (8).


### 2.5. The Additive Partitioning of Biodiversity with Respect to the Pairs of Species

Let us assume that in a certain region there aren species and m sites. Again, in what follows, the subscripts i and j refer to species (i,j=1,…,n) and the subscripts k and r refer to sites, (k,r=1,…,m). Let Θk=[πij,k] be an arbitrary joint probability distribution of the pairs of species within site k, where πij,k is the probability of the pair of species (i,j), in this order, within site k, such that πij,k≥0,∑i,jπij,k=1. Let W=[wij] be the matrix whose entries are nonnegative weights assigned to the pairs of species. We assign a parameter λk to each site k, satisfying (3). As GSW(Θ), given by (9), is a concave function of the joint distribution Θ assigned to the pairs of species, it may be used in the additive partitioning of biodiversity. The corresponding γ-diversity, reflecting the total or regional biodiversity, the α-diversity, interpreted as the within-site diversity or the average diversity of the sites, and the β-diversity, as a measure of between-site diversity, with respect to the pairs of species, are given by(14)γ=GSW(∑kλkΘk),α=∑kλkGSW(Θk),β=γ-α.Theβ-diversity may be interpreted as a measure of dissimilarity or differentiation between the sites of the respective region with respect to the pairs of species. As shown in [20, 21], taking into account (14), the β-diversity has the following expression:(15)β=∑i,jwij∑k<rλkλr(πij,k-πij,r)2<∑i,jwij∑k<rλkλr.
The upper bound in (15) is quite loose. As α is nonnegative and cannot exceed γ, a better upper bound for β is the value B5 from (B.4) of Appendix B, corresponding to n and the given weights W. Thus, the relative between-site diversity, or the relative dissimilarity between sites, is 0≤β/B5≤1, with respect to the pairs of species.Letθk=(p1,k,…,pn,k) be the vector whose components are the relative abundances of the individual species at site k.If the species are independent, πij,k=pi,kpj,k. Let also D=[dij] be the matrix of the distances between species. Then, for the weights (12), the corresponding β-diversity (15) is(16)β=n(n-1)2∑i<j(vi+vj)dij∑k<rλkλr(pi,kpj,k-pi,rpj,r)2,
measuring the between-site diversity with respect to the distinct pairs of species when the species richness, the distance between species, and the conservation values are all taken into account along with the relative abundance of the species. If the species have the same abundance in each site, which means that pi,k=pi, for each site k and each species i, the corresponding β-diversity is equal to zero, reflecting the fact that in such a case there is no dissimilarity between sites. If we divide (16) by the upper bound B5 from the inequality (B.4) of Appendix B, taking the weights (12), we obtain the relative β-diversity, which has the advantage of always being a number between 0 and 1.At the same time, let us notice that Rao’s indexRD given by (8) is a linear function of the joint distribution of the relative abundance of the pairs of species Θ=[pipj], and, consequently, the corresponding β-diversity induced by the pairs of species is equal to zero, for an arbitrary distance matrix D and an arbitrary distribution of the relative abundance of individual species. Therefore, Rao’s index is not suitable for use in the standard additive partitioning of diversity induced by pairs of species when an arbitrary dissimilarity distance between species is taken into account. Attempts were made [24, 25] to either find particular cases of distance matrices or define new kinds of nonstandard α-divergence for which the partitioning of diversity can still be performed.


### 2.6. The Multiplicative Partitioning of Biodiversity with Respect to the Pairs of Species

The weighted Gini-Simpson quadratic indexGSW given by (9), which can be used in the additive partitioning of diversity induced by pairs of species, may be transformed into the measure of diversity [21]:(17)1[∑i,jwijπij-GSw(Θ)]=(∑i,jwijπij2)-1,
which can be used in the multiplicative partitioning of diversity induced by the pairs of species. This measure of biodiversity may be viewed as being the weighted version for pairs of species of the classic Hill number of first degree from [10]. Using the notations from the previous Section 2.5, the corresponding multiplicative γ-diversity, α-diversity, and β-diversity are(18)γ=[∑i,jwij(∑kλkπij,k)2]-1,α=[∑kλk∑i,jwijπij,k2]-1,β=γα.
Due to the convexity of the function ∑i,jwijπij2, as a function of the joint distribution Θ=[πij], the γ-diversity cannot be smaller than the α-diversity, as it should be, and, consequently, the multiplicative β-diversity satisfies the inequality: β=γ/α≥1. In the additive partitioning of diversity, the γ-diversity, α-diversity, and β-diversity are entities of the same kind and may be expressed in the same units. In the multiplicative partitioning of diversity, the β-diversity is simply a ratio between the total, regional biodiversity γ and the average within-site biodiversity α, a numerical indicator showing to what extent the regional biodiversity, as a whole, exceeds the average biodiversity of the sites of the respective region.Letθk=(p1,k,…,pn,k) be the vector whose components are the relative abundances of the individual species at site k.If the species are independent, πij,k=pi,kpj,k. Let also D=[dij] be the matrix of the distances between species. Then, for the weights (12), the corresponding β-diversity from (18) is(19)β=[∑kλk∑i<j(vi+vj)dijpi,k2pj,k2][∑i<j(vi+vj)dij(∑kλkpi,kpj,k)2]≥1,
measuring the ratio between the regional biodiversity and the average biodiversity of the sites with respect to the pairs of species. Obviously, if the sites have the same species and the same abundance of these species, which means that pi,k=pi, for each species i, then β=1.


### 2.7. The Weighted Shannon Entropy

The weighted Shannon entropy was introduced in [26]. If we have n species such that the distribution of the relative abundance of these species is θ=(p1,…,pn) and the nonnegative weights assigned to the species are w=(w1,…,wn), then the weighted entropy is the nonnegative, concave function Hw(θ)=-∑iwipilnpi. Similarly, if W=[wij] is a matrix of nonnegative weights and Θ=[πij] a joint probability distribution assigned to the pairs of species, the joint weighted entropy is HW(Θ)=-∑i,jwijπijlnπij. It is possible, in principle, to remake the analysis from Sections 2.1–2.6 using the weighted Shannon entropies Hw(θ) and HW(Θ) instead of the weighted Gini-Simpson indices GSw(θ) and GSW(Θ). However, the Shannon entropy is actually a measure of uncertainty and we cannot justify its use as a measure of diversity, as we did for the Gini-Simpson index at the beginning of Section 2.1. Also, since the Shannon entropy is a logarithmic function, it is much more difficult to obtain simple analytical formulas for its maximum values subject to given constraints. The weighted Gini-Simpson index is a simpler and more effective tool in measuring biodiversity.


## 2.1. The Weighted Measure of Diversity with Respect to Individual Species

Let us assume that there aren species in a certain community/habitat/site and let pi be the relative abundance of species i (the number of individuals of species i divided by the total number of individuals in that community/habitat/site). We have diversity if species i is present at that location but other species are found there as well. The probability that the species i is present and there are other species present as well is pi(1-pi). If we take all possible values of pi from the unit interval [0,1] into account, the wave function pi(1-pi) corresponding to the species i is a nonnegative, symmetric, bell-shaped, concave function, reaching its maximum value 1/4 at pi=1/2. If we sum up these wave functions, for all n species, we obtain the classic Gini-Simpson index GS(θ) corresponding to the given distribution of the relative abundance of the species θ=(p1,…,pn). Since 1949, this has been considered to be a very good measure of biodiversity. In order to generalize it, we may assign an amplitude wi≥0 to the wave function of the species i, and the resulting new wave function wipi(1-pi) continues to be a nonnegative, symmetric, bell-shaped, concave function of pi, but this time its maximum value is wi/4. Summing up these wave functions for all the species, we get the weighted Gini-Simpson index:(1)GSw(θ)=∑iwipi(1-pi),
which depends both on the distribution of the relative abundance of species θ and on the nonnegative weights w=(w1,…,wn). The concavity of GSw(θ) was proven in [20, 21]. The weight wi could be anything which contributes to the increase in the diversity induced by the species i. However, the weights may not depend on the relative abundance of species. If wi=n, for each i, then (1) becomes the so-called Rich-Gini-Simpson index GSn(θ),introduced in [22], which is essentially dependent on the species richness of the respective community/habitat/site. If there are some conservation values assigned to the species v=(v1,…,vn), which are positive numbers on a certain scale of values, and the weights are wi=nvi, the corresponding weighted Gini-Simpson index is denoted by GSn,v(θ). Obviously, if wi=1, for each species i, then (1) is the classic Gini-Simpson index GS(θ). An upper bound for GSw(θ), which depends only on the maximum weight and the number of species, is(2)0≤GSw(θ)≤(max⁡i⁡wi)∑ipi(1-pi)≤(max⁡i⁡wi)(1-1n).
Denoting by B1 the bound from the right-hand side of the inequality (2), the relative weighted Gini-Simpson index for individual species is 0≤GSw(θ)/B1≤1. In Appendix A, another bound of GSw(θ) is given, denoted by B2, which depends on all the weights assigned to the species. If wi=1, for each species i, we get maxθGS(θ)=1-1/n, and this maximum biodiversity is obtained when all species have the same relative abundance pi=1/n. The fact that the maximum value of GS(θ) is almost insensitive to the increase of the number of species, tending very slowly to 1 when n increases, allowed Jost [12, 13] and Jost et al. [19] to give some examples showing that the Gini-Simpson index does not behave well when the number of species n is very large. R. C. Guiasu and S. Guiasu [22] showed, however, that the Rich-Gini-Simpson index GSn(θ) has no such problem. Indeed, if we take wi=n, for each species i, we get from (2), maxθGSn(θ)=n-1, whose value sensibly increases when the number of species increases, which makes inapplicable the criticism of the classic Gini-Simpson index.


## 2.2. The Additive Partitioning of Biodiversity with Respect to the Individual Species

Let us assume that in a certain region there aren species and m sites. In what follows, the subscripts i and j refer to species (i,j=1,…,n) and the subscripts k and r refer to sites, (k,r=1,…,m). Let θk=(p1,k,…,pn,k) be the vector whose components are the relative abundances of the individual species at site k, such that pi,k≥0,(i=1,…,n),∑ipi,k=1, for each k=1,…,m. Let w=(w1,…,wn) be nonnegative weights assigned to the species. In dealing with species diversity, a good measure of the differentiation, or dissimilarity, among the sites in a certain region has to be nonnegative and equal to zero if and only if there is no such difference. We assign a parameter λk to each site k, such that(3)λk≥0,(k=1,…,m),∑kλk=1.
These parameters may be used to make adjustments for differences (in size, altitude, etc.) between the sites, as shown in [23]. If no adjustment is made and we focus only on the species abundance, we take these parameters to be equal, that is, λk=1/m,for every (k=1,…,m). As GSw(θ) is a concave function of the distribution of the relative abundance θ, it may be used in the additive partitioning of biodiversity. The corresponding γ-diversity, reflecting the total or regional biodiversity, the α-diversity, interpreted as the within-site diversity or the average diversity of the sites, and the β-diversity, as a measure of between-site diversity, are given by(4)γ=GSw(∑kλkθk),α=∑kλkGSw(θk),β=γ-α.Theβ-diversity may be interpreted as a measure of dissimilarity or differentiation between the sites of the respective region with respect to the individual species. As shown in Appendix C, taking into account (4), the β-diversity has the expression(5)β=∑iwi∑k<rλkλr(pi,k-pi,r)2<2∑iwi∑k<rλkλr.
The upper bound in (5) is quite loose. As α is nonnegative and cannot exceed γ, a better upper bound for β is the value B2 from (A.1), corresponding to n and the given weights w. Thus, the relative between-site diversity, or the dissimilarity between sites is 0≤β/B2≤1, with respect to the individual species. From (5), we can see that if the species have the same abundance in each site, which means that pi,k=pi, for each site k and each species i, the β-diversity is equal to zero, reflecting the fact that in such a case there is no dissimilarity between the sites.


## 2.3. The Multiplicative Partitioning of Biodiversity with Respect to Individual Species

Dealing with the multiplicative partitioning of diversity, Whittaker [8] suggested the use of the exponential of the Shannon entropy as a measure of biodiversity. The weighted Gini-Simpson index GSw(θ), given by (1), which can be used in the additive partitioning of diversity induced by individual species, may also be transformed into the measure of biodiversity (R. C. Guiasu and S. Guiasu [21]):(6)1[∑iwipi-GSw(θ)]=(∑iwipi2)-1,
which can be used in the multiplicative partitioning of diversity induced by the individual species. This measure of biodiversity may be viewed as being the weighted version of the classic Hill number of first degree from [10]. The corresponding multiplicative γ-diversity and α-diversity are(7)γ=[∑iwi(∑kλkpi,k)2]-1,α=[∑kλk∑iwipi,k2]-1.
Due to the convexity of the function ∑iwipi2, as a function of the distribution of the relative abundance of species, the γ-diversity cannot be smaller than the α-diversity, as it should be, and, consequently, the multiplicative β-diversity satisfies the inequality: β=γ/α≥1. In the additive partitioning of diversity, the γ-diversity, α-diversity, and β-diversity are entities of the same kind and may be expressed in the same units. In the multiplicative partitioning of diversity, the β-diversity is simply a ratio between the total, regional biodiversity γ and the average within-site biodiversity α, a numerical indicator showing to what extent the regional biodiversity, as a whole, exceeds the average biodiversities of the sites of the respective region. Obviously, if the sites have the same species and the same abundance of these species, which means that pi,k=pi, for each species i, then β=1.


## 2.4. The Weighted Measure of Diversity with Respect to the Pairs of Species

LetD=[dij] be an n×n matrix whose entries are the distances between the pairs of n species, such that dij≥0,dii=0,(i,j=1,…,n). This could be the matrix of the phylogenetic distance between species, for instance. When Rao introduced his quadratic index, improperly called quadratic entropy [14]:(8)RD=∑i,jdijpipj,
he had to focus on the pairs of species instead of the individual species, using the distance between the pairs of distinct species, along with their relative abundance, in order to measure the dissimilarity between species. Rao’s indicator is very simple and may be easily interpreted as the average dissimilarity between two individuals belonging to two different species when the phylogenetic distance is taken into account. There have been numerous attempts, such as [24], for instance, at using Rao’s index RD in the additive partitioning of diversity. Unfortunately, RD is not a concave function of the distribution of the relative abundance of species θ=(p1,…,pn) for an arbitrary distance matrix D. Thus, RD can be applied to the additive partitioning of biodiversity only for some special kinds of such matrices D (as mentioned in [24], for instance), but not in general. On the other hand, there is no generally accepted proposal of a new kind of nonstandard α-divergence which could be defined for such a measure which is not a concave function of the distribution of the relative abundance of species. This paper shows that the generalization of the weighted Gini-Simpson index to the pairs of species provides a concave measure of diversity which could indeed be used both for the additive partitioning and the multiplicative partitioning of diversity when the phylogenetic distance between species is taken into account. Therefore, this provides a suitable replacement of Rao’s index in the partitioning of biodiversity.LetΘ=[πij] be an n×n matrix where πij is a joint probability of the pair of species (i,j), in this order. As πji is the probability of the pair (j,i), in this order, the probability of the subset of species {i,j} is πij+πji. We haveπij≥0,(i,j=1,…,n);∑i,jπij=1. Let W=[wij] be an n×n matrix whose entries are arbitrary nonnegative weights assigned to the pairs of species. However, these weights may not depend on the joint distribution [πij]. The weighted Gini-Simpson quadratic index of the pairs of species is(9)GSW(Θ)=∑i,jwijπij(1-πij),
If wij=1, for all pairs of species, then (9) becomes the generalization of the classic Gini-Simpson index to the pairs of species and is denoted by GS(Θ). As a function of the joint distribution Θ, the weighted Gini-Simpson index GSW(Θ) is nonnegative and concave, as shown in [20, 21]. An upper bound for GSW(Θ) which depends only on the maximum weight and the number of species is(10)0≤GSW(Θ)≤(max⁡i,j⁡wij)∑i,jπij(1-πij)≤(max⁡i,j⁡wij)(1-1n2).
Denoting by B3 the bound from the right-hand side of the inequality (10), the relative weighted Gini-Simpson index for pairs of species is 0≤GSW(Θ)/B3≤1. In Appendix B, another bound for GSW(Θ), denoted by B5, is given, which depends on all the weights assigned to the species. If wij=1, for each pair of species (i,j), we get maxΘGS(Θ)=1-1/n2, and this maximum biodiversity is obtained when all species have the same relative abundance πij=1/n2. The special cases of interest are the following ones.(a)
If the species are independent, which meansπij=pipj and the weights are wij=dij, then the weighted Gini-Simpson index (9) is denoted by
(11)GSD(Θ)=∑i,jdijpipj(1-pipj),
which generalizes Rao’s index RD given by (8).(b)
If there are the positive numbersv=(v1,…,vn), representing conservation values of the individual species, the species are independent, which means πij=pipj, and the weights are
(12)wij=n(n-1)212(vi+vj)dij,
where n(n-1)/2 is the number of distinct pairs of species (i,j), such that i<j, or the number of the pairs of species {i,j}, and (1/2)(vi+vj) is the average value of the pair of species (i,j), then the corresponding weighted Gini-Simpson index (9) is denoted by:
(13)GSn,v,D(Θ)=n(n-1)2∑i<j(vi+vj)dijpipj(1-pipj),
which takes into account all the information available, namely, the species richness n, the relative abundance θ of species, the matrix D of the distance between species, and the conservation values v of the species. As the measure given in (13) is a nonnegative concave function of the distribution of the relative abundance of the pairs of species Θ, for an arbitrary distance matrix D, and also depends explicitly on the species richness, the distance between species, and the conservation values of the species, all these are sufficient reasons to suggest that it could more than adequately replace the use of Rao’s index (8).


## 2.5. The Additive Partitioning of Biodiversity with Respect to the Pairs of Species

Let us assume that in a certain region there aren species and m sites. Again, in what follows, the subscripts i and j refer to species (i,j=1,…,n) and the subscripts k and r refer to sites, (k,r=1,…,m). Let Θk=[πij,k] be an arbitrary joint probability distribution of the pairs of species within site k, where πij,k is the probability of the pair of species (i,j), in this order, within site k, such that πij,k≥0,∑i,jπij,k=1. Let W=[wij] be the matrix whose entries are nonnegative weights assigned to the pairs of species. We assign a parameter λk to each site k, satisfying (3). As GSW(Θ), given by (9), is a concave function of the joint distribution Θ assigned to the pairs of species, it may be used in the additive partitioning of biodiversity. The corresponding γ-diversity, reflecting the total or regional biodiversity, the α-diversity, interpreted as the within-site diversity or the average diversity of the sites, and the β-diversity, as a measure of between-site diversity, with respect to the pairs of species, are given by(14)γ=GSW(∑kλkΘk),α=∑kλkGSW(Θk),β=γ-α.Theβ-diversity may be interpreted as a measure of dissimilarity or differentiation between the sites of the respective region with respect to the pairs of species. As shown in [20, 21], taking into account (14), the β-diversity has the following expression:(15)β=∑i,jwij∑k<rλkλr(πij,k-πij,r)2<∑i,jwij∑k<rλkλr.
The upper bound in (15) is quite loose. As α is nonnegative and cannot exceed γ, a better upper bound for β is the value B5 from (B.4) of Appendix B, corresponding to n and the given weights W. Thus, the relative between-site diversity, or the relative dissimilarity between sites, is 0≤β/B5≤1, with respect to the pairs of species.Letθk=(p1,k,…,pn,k) be the vector whose components are the relative abundances of the individual species at site k.If the species are independent, πij,k=pi,kpj,k. Let also D=[dij] be the matrix of the distances between species. Then, for the weights (12), the corresponding β-diversity (15) is(16)β=n(n-1)2∑i<j(vi+vj)dij∑k<rλkλr(pi,kpj,k-pi,rpj,r)2,
measuring the between-site diversity with respect to the distinct pairs of species when the species richness, the distance between species, and the conservation values are all taken into account along with the relative abundance of the species. If the species have the same abundance in each site, which means that pi,k=pi, for each site k and each species i, the corresponding β-diversity is equal to zero, reflecting the fact that in such a case there is no dissimilarity between sites. If we divide (16) by the upper bound B5 from the inequality (B.4) of Appendix B, taking the weights (12), we obtain the relative β-diversity, which has the advantage of always being a number between 0 and 1.At the same time, let us notice that Rao’s indexRD given by (8) is a linear function of the joint distribution of the relative abundance of the pairs of species Θ=[pipj], and, consequently, the corresponding β-diversity induced by the pairs of species is equal to zero, for an arbitrary distance matrix D and an arbitrary distribution of the relative abundance of individual species. Therefore, Rao’s index is not suitable for use in the standard additive partitioning of diversity induced by pairs of species when an arbitrary dissimilarity distance between species is taken into account. Attempts were made [24, 25] to either find particular cases of distance matrices or define new kinds of nonstandard α-divergence for which the partitioning of diversity can still be performed.


## 2.6. The Multiplicative Partitioning of Biodiversity with Respect to the Pairs of Species

The weighted Gini-Simpson quadratic indexGSW given by (9), which can be used in the additive partitioning of diversity induced by pairs of species, may be transformed into the measure of diversity [21]:(17)1[∑i,jwijπij-GSw(Θ)]=(∑i,jwijπij2)-1,
which can be used in the multiplicative partitioning of diversity induced by the pairs of species. This measure of biodiversity may be viewed as being the weighted version for pairs of species of the classic Hill number of first degree from [10]. Using the notations from the previous Section 2.5, the corresponding multiplicative γ-diversity, α-diversity, and β-diversity are(18)γ=[∑i,jwij(∑kλkπij,k)2]-1,α=[∑kλk∑i,jwijπij,k2]-1,β=γα.
Due to the convexity of the function ∑i,jwijπij2, as a function of the joint distribution Θ=[πij], the γ-diversity cannot be smaller than the α-diversity, as it should be, and, consequently, the multiplicative β-diversity satisfies the inequality: β=γ/α≥1. In the additive partitioning of diversity, the γ-diversity, α-diversity, and β-diversity are entities of the same kind and may be expressed in the same units. In the multiplicative partitioning of diversity, the β-diversity is simply a ratio between the total, regional biodiversity γ and the average within-site biodiversity α, a numerical indicator showing to what extent the regional biodiversity, as a whole, exceeds the average biodiversity of the sites of the respective region.Letθk=(p1,k,…,pn,k) be the vector whose components are the relative abundances of the individual species at site k.If the species are independent, πij,k=pi,kpj,k. Let also D=[dij] be the matrix of the distances between species. Then, for the weights (12), the corresponding β-diversity from (18) is(19)β=[∑kλk∑i<j(vi+vj)dijpi,k2pj,k2][∑i<j(vi+vj)dij(∑kλkpi,kpj,k)2]≥1,
measuring the ratio between the regional biodiversity and the average biodiversity of the sites with respect to the pairs of species. Obviously, if the sites have the same species and the same abundance of these species, which means that pi,k=pi, for each species i, then β=1.


## 2.7. The Weighted Shannon Entropy

The weighted Shannon entropy was introduced in [26]. If we have n species such that the distribution of the relative abundance of these species is θ=(p1,…,pn) and the nonnegative weights assigned to the species are w=(w1,…,wn), then the weighted entropy is the nonnegative, concave function Hw(θ)=-∑iwipilnpi. Similarly, if W=[wij] is a matrix of nonnegative weights and Θ=[πij] a joint probability distribution assigned to the pairs of species, the joint weighted entropy is HW(Θ)=-∑i,jwijπijlnπij. It is possible, in principle, to remake the analysis from Sections 2.1–2.6 using the weighted Shannon entropies Hw(θ) and HW(Θ) instead of the weighted Gini-Simpson indices GSw(θ) and GSW(Θ). However, the Shannon entropy is actually a measure of uncertainty and we cannot justify its use as a measure of diversity, as we did for the Gini-Simpson index at the beginning of Section 2.1. Also, since the Shannon entropy is a logarithmic function, it is much more difficult to obtain simple analytical formulas for its maximum values subject to given constraints. The weighted Gini-Simpson index is a simpler and more effective tool in measuring biodiversity.


## 3. Discussion

It seems to be much easier to discuss the significance of the concepts introduced in Section2 by showing a representative numerical example. Let us assume that in a certain region there are three sites (m=3) and three species (n=3). If Aik denotes the absolute abundance (number of individuals) of species i within site k, let us assume that(20)A11=2,A21=24,A31=14;A12=32,A22=4,A32=14;A13=24,A23=36,A33=20.
The corresponding relative abundance is (21)p1,1=0.05,p2,1=0.60,p3,1=0.35;p1,2=0.64,p2,2=0.08,p3,2=0.28;p1,3=0.30,p2,3=0.45,p3,3=0.25.
Thus, in this example, θ1=(0.05,0.60,0.35), θ2=(0.64,0.08,0.28), and θ3=(0.30,0.45,0.25).


### 3.1. Biodiversity with Respect to the Individual Species

Using the Rich-Gini-Simpson indexGSn(θ), given by (1) with the weights wi=n, to calculate the amount of diversity with respect to the individual species, in each site, we obtain GS3(θ1)=1.5450, GS3(θ2)=1.5168, and GS3(θ3)=1.9350. The maximum biodiversity in this case would be n-1=2. We can see that the first two sites have almost the same biodiversity, both a little smaller than the biodiversity of the third site which is close to the maximum value, when only the richness and the abundance of species are taken into account.Let us assume now that the three species have the following conservation values:v1=6, v2=3, and v3=3. These conservation values v=(6,3,3) contribute to the diversity of the three sites. Taking the weights wi=nvi, we have w1=18,w2=9,w3=9. Therefore, w=(18,9,9). Using the weighted Gini-Simpson index GSw(θ) given by (1), we obtain the following values of the biodiversity of each site: GSw(θ1)=5.0625, GSw(θ2)=6.624, and GSw(θ3)=7.695. When the species have these conservation values, the biodiversity of the second and third sites are closer and higher than the biodiversity of the first site. But in order to have a better understanding of these numbers, we have to compare them with the bounds B1 and B2 from the inequalities (2) and (A.1), respectively. For the weights w=(18,9,9), the loose upper bound B1 for GSw, which takes into account only the number of species n=3 and the maximum weight maxiwi=18, has the value 12. For the much better upper bound B2 for GSw from (A.1), mentioned in Appendix A, which takes into account the number of species n=3 and all the weights w=(18,9,9), we get the value 8.1. Therefore, we can see that the bound B2 is obviously better than B1. With respect to B2, the second and third sites have 81.78% and 95% of the maximum biodiversity for the given weights, whereas the first site has only 62.5%. If we do not discriminate among sites with respect to size, altitude, or any other factor, then the parameters assigned to the three sites are λ1=λ2=λ3=1/3. In such a case, we have(22)∑kλkθk=13(0.05,0.60,0.35)+13(0.64,0.08,0.28)+13(0.30,0.45,0.25)=(0.3300,0.3767,0.2933)=(q1,q2,q3).
According to (4), the γ-diversity and α-diversity, with respect to the single species, are(23)γ=GSw(∑kλkθk)=∑iwiqi(1-qi)=7.9584,α=∑kλkGSw(θk)=13(5.0625+6.624+7.695)=6.4605.
Thus, in the additive partitioning of diversity, the β-diversity is β=γ-α=1.4979. For the weights w=(18,9,9) and n=3, according to the formula (A.1) from Appendix A, the maximum value of GSw is B2=8.1. Therefore, the biodiversity γ of the entire region is 98.25% of the maximum and the average within-site biodiversity α is 79.76%. The value of the between-site diversity β shows the average differentiation between sites corresponding to a difference of 18.49% between the values of γ and α. We note that for identical sites, the value of β would be equal to zero, as could be seen from (5). The advantage of the use of the additive partitioning of biodiversity is that the values of α,β, and γ are expressed on the same scale of values.Doing the multiplicative partitioning of biodiversity forλi=1/3,(i=1,2,3), and w=nv=(18,9,9), from (7) we get γ=0.2493 and α=0.1815. Consequently, β=γ/α=1.3736.


### 3.2. Biodiversity with Respect to the Pairs of Species

Let us assume that we have the matrix of the phylogenetic distances between the three speciesD=[dij], where d12=3,d13=2,andd23=2. If we assume that within each site the species are supposed to be independent from the point of view of their relative abundance, then the relative abundance of the pair of species (i,j), in this order, is the product of the relative abundance of the corresponding individual species, namely, pi,kpj,k, within every site k. Therefore, the matrices Θk=[pi,kpj,k] are:(24)Θ1=[0.00250.03000.01750.03000.36000.21000.01750.21000.1225],Θ2=[0.40960.05120.17920.05120.00640.02240.17920.02240.0784],Θ3=[0.09000.13500.07500.13500.20250.11250.07500.11250.0625].
If we do not discriminate among sites with respect to size, altitude, or any other factor, then the parameters assigned to the three sites are λ1=λ2=λ3=1/3. In such a case, we have(25)λ1Θ1+λ2Θ2+λ3Θ3=[0.16740.07210.09060.07210.18960.11500.09060.11500.0878].
Let us use Rao’s index (8) for doing the additive partitioning of diversity with respect to the pairs of species. Successively, we obtain RD(Θ1)=1.0900,RD(Θ2)=1.1136, and RD(Θ3)=1.5600. The corresponding α-diversity is α=λ1RD(Θ1)+λ2RD(Θ2)+λ3RD(Θ3)=1.255, and the γ-diversity is γ=RD(λ1Θ1+λ2Θ2+λ3Θ3)=1.255. Consequently, the β-diversity is β=γ-α=0, which is not surprising because Rao’s index is a linear function of the joint distribution of the pairs of species.If we use the weighted Gini-Simpson index (11) with the weights wij=dij, we obtain(26)GSD(Θ1)=0.9070,GSD(Θ2)=0.9674,GSD(Θ3)=1.3775,
and the corresponding α-diversity is
(27)α=λ1GSD(Θ1)+λ2GSD(Θ2)+λ3GSD(Θ3)=1.0840,
the γ-diversity is γ=GSD(λ1Θ1+λ2Θ2+λ3Θ3)=1.1381, and the β-diversity is β=γ-α=0.0541. Calculating the upper bound B5 of GSW given in the inequality (B.4) from Appendix B, for the weights wij=dij, which means W=D, we obtain maxGSD=2.75. Compared to this maximum value, GSD(Θ1) represents 32.98%; GSD(Θ2)=35.18%; GSD(Θ3)=50.09%; γ=41.39%; α=39.42%; β=1.97%.We take now into account the number of speciesn=3, the parameters assigned to the sites λ1=λ2=λ3=1/3, the phylogenetic distances between species d12=3,d13=2,d23=2, and the conservation values of the species v1=6, v2=3, v3=3. The computation of the weighted Gini-Simpson index given by (13), with the weights wij=(n(n-1)/2)(1/2)(vi+vj)dij, gives(28)GSn,v,D(Θ1)=9.2580,GSn,v,D(Θ2)=12.6659,GSn,v,D(Θ3)=16.7994,
and the corresponding α-diversity is
(29)α=λ1GSn,v,D(Θ1)+λ2GSn,v,D(Θ2)+λ3GSn,v,D(Θ3)=12.9078,
while the γ-diversity is γ=GSn,v,D(λ1Θ1+λ2Θ2+λ3Θ3)=13.5321, which gives the β-diversity: β=γ-α=0.6243. Calculating the upper bound B5 of GSn,v,D given in the inequality (B.4) from Appendix B, for the weights wij=(n(n-1)/2)(1/2)(vi+vj)dij, we obtain maxGSn,v,D=34.2237. Compared to this maximum value, GSn,v,D(Θ1) represents 27.05%; GSn,v,D(Θ2)=37.01%; GSn,v,D(Θ3)=49.09%; γ=39.54%; α=37.72%; β=1.82%.Doing the multiplicative partitioning of biodiversity forλi=1/3,(i=1,2,3), and wij=dij, from (18) and (19), we get γ=8.56 and α=5.86. Consequently, β=γ/α=1.46. Doing the multiplicative partitioning of biodiversity for the site parameters λi=1/3,(i=1,2,3), and the weights wij=(n(n-1)/2)(1/2)(vi+vj)dij, from (18), we get γ=0.75 and α=0.51. Consequently, β=γ/α=1.47.


## 3.1. Biodiversity with Respect to the Individual Species

Using the Rich-Gini-Simpson indexGSn(θ), given by (1) with the weights wi=n, to calculate the amount of diversity with respect to the individual species, in each site, we obtain GS3(θ1)=1.5450, GS3(θ2)=1.5168, and GS3(θ3)=1.9350. The maximum biodiversity in this case would be n-1=2. We can see that the first two sites have almost the same biodiversity, both a little smaller than the biodiversity of the third site which is close to the maximum value, when only the richness and the abundance of species are taken into account.Let us assume now that the three species have the following conservation values:v1=6, v2=3, and v3=3. These conservation values v=(6,3,3) contribute to the diversity of the three sites. Taking the weights wi=nvi, we have w1=18,w2=9,w3=9. Therefore, w=(18,9,9). Using the weighted Gini-Simpson index GSw(θ) given by (1), we obtain the following values of the biodiversity of each site: GSw(θ1)=5.0625, GSw(θ2)=6.624, and GSw(θ3)=7.695. When the species have these conservation values, the biodiversity of the second and third sites are closer and higher than the biodiversity of the first site. But in order to have a better understanding of these numbers, we have to compare them with the bounds B1 and B2 from the inequalities (2) and (A.1), respectively. For the weights w=(18,9,9), the loose upper bound B1 for GSw, which takes into account only the number of species n=3 and the maximum weight maxiwi=18, has the value 12. For the much better upper bound B2 for GSw from (A.1), mentioned in Appendix A, which takes into account the number of species n=3 and all the weights w=(18,9,9), we get the value 8.1. Therefore, we can see that the bound B2 is obviously better than B1. With respect to B2, the second and third sites have 81.78% and 95% of the maximum biodiversity for the given weights, whereas the first site has only 62.5%. If we do not discriminate among sites with respect to size, altitude, or any other factor, then the parameters assigned to the three sites are λ1=λ2=λ3=1/3. In such a case, we have(22)∑kλkθk=13(0.05,0.60,0.35)+13(0.64,0.08,0.28)+13(0.30,0.45,0.25)=(0.3300,0.3767,0.2933)=(q1,q2,q3).
According to (4), the γ-diversity and α-diversity, with respect to the single species, are(23)γ=GSw(∑kλkθk)=∑iwiqi(1-qi)=7.9584,α=∑kλkGSw(θk)=13(5.0625+6.624+7.695)=6.4605.
Thus, in the additive partitioning of diversity, the β-diversity is β=γ-α=1.4979. For the weights w=(18,9,9) and n=3, according to the formula (A.1) from Appendix A, the maximum value of GSw is B2=8.1. Therefore, the biodiversity γ of the entire region is 98.25% of the maximum and the average within-site biodiversity α is 79.76%. The value of the between-site diversity β shows the average differentiation between sites corresponding to a difference of 18.49% between the values of γ and α. We note that for identical sites, the value of β would be equal to zero, as could be seen from (5). The advantage of the use of the additive partitioning of biodiversity is that the values of α,β, and γ are expressed on the same scale of values.Doing the multiplicative partitioning of biodiversity forλi=1/3,(i=1,2,3), and w=nv=(18,9,9), from (7) we get γ=0.2493 and α=0.1815. Consequently, β=γ/α=1.3736.


## 3.2. Biodiversity with Respect to the Pairs of Species

Let us assume that we have the matrix of the phylogenetic distances between the three speciesD=[dij], where d12=3,d13=2,andd23=2. If we assume that within each site the species are supposed to be independent from the point of view of their relative abundance, then the relative abundance of the pair of species (i,j), in this order, is the product of the relative abundance of the corresponding individual species, namely, pi,kpj,k, within every site k. Therefore, the matrices Θk=[pi,kpj,k] are:(24)Θ1=[0.00250.03000.01750.03000.36000.21000.01750.21000.1225],Θ2=[0.40960.05120.17920.05120.00640.02240.17920.02240.0784],Θ3=[0.09000.13500.07500.13500.20250.11250.07500.11250.0625].
If we do not discriminate among sites with respect to size, altitude, or any other factor, then the parameters assigned to the three sites are λ1=λ2=λ3=1/3. In such a case, we have(25)λ1Θ1+λ2Θ2+λ3Θ3=[0.16740.07210.09060.07210.18960.11500.09060.11500.0878].
Let us use Rao’s index (8) for doing the additive partitioning of diversity with respect to the pairs of species. Successively, we obtain RD(Θ1)=1.0900,RD(Θ2)=1.1136, and RD(Θ3)=1.5600. The corresponding α-diversity is α=λ1RD(Θ1)+λ2RD(Θ2)+λ3RD(Θ3)=1.255, and the γ-diversity is γ=RD(λ1Θ1+λ2Θ2+λ3Θ3)=1.255. Consequently, the β-diversity is β=γ-α=0, which is not surprising because Rao’s index is a linear function of the joint distribution of the pairs of species.If we use the weighted Gini-Simpson index (11) with the weights wij=dij, we obtain(26)GSD(Θ1)=0.9070,GSD(Θ2)=0.9674,GSD(Θ3)=1.3775,
and the corresponding α-diversity is
(27)α=λ1GSD(Θ1)+λ2GSD(Θ2)+λ3GSD(Θ3)=1.0840,
the γ-diversity is γ=GSD(λ1Θ1+λ2Θ2+λ3Θ3)=1.1381, and the β-diversity is β=γ-α=0.0541. Calculating the upper bound B5 of GSW given in the inequality (B.4) from Appendix B, for the weights wij=dij, which means W=D, we obtain maxGSD=2.75. Compared to this maximum value, GSD(Θ1) represents 32.98%; GSD(Θ2)=35.18%; GSD(Θ3)=50.09%; γ=41.39%; α=39.42%; β=1.97%.We take now into account the number of speciesn=3, the parameters assigned to the sites λ1=λ2=λ3=1/3, the phylogenetic distances between species d12=3,d13=2,d23=2, and the conservation values of the species v1=6, v2=3, v3=3. The computation of the weighted Gini-Simpson index given by (13), with the weights wij=(n(n-1)/2)(1/2)(vi+vj)dij, gives(28)GSn,v,D(Θ1)=9.2580,GSn,v,D(Θ2)=12.6659,GSn,v,D(Θ3)=16.7994,
and the corresponding α-diversity is
(29)α=λ1GSn,v,D(Θ1)+λ2GSn,v,D(Θ2)+λ3GSn,v,D(Θ3)=12.9078,
while the γ-diversity is γ=GSn,v,D(λ1Θ1+λ2Θ2+λ3Θ3)=13.5321, which gives the β-diversity: β=γ-α=0.6243. Calculating the upper bound B5 of GSn,v,D given in the inequality (B.4) from Appendix B, for the weights wij=(n(n-1)/2)(1/2)(vi+vj)dij, we obtain maxGSn,v,D=34.2237. Compared to this maximum value, GSn,v,D(Θ1) represents 27.05%; GSn,v,D(Θ2)=37.01%; GSn,v,D(Θ3)=49.09%; γ=39.54%; α=37.72%; β=1.82%.Doing the multiplicative partitioning of biodiversity forλi=1/3,(i=1,2,3), and wij=dij, from (18) and (19), we get γ=8.56 and α=5.86. Consequently, β=γ/α=1.46. Doing the multiplicative partitioning of biodiversity for the site parameters λi=1/3,(i=1,2,3), and the weights wij=(n(n-1)/2)(1/2)(vi+vj)dij, from (18), we get γ=0.75 and α=0.51. Consequently, β=γ/α=1.47.


## 4. Conclusion

Using a measure of biodiversity, as a mathematical tool, the distribution of biodiversity at multiple sites of a region has been traditionally investigated through the partitioning of the regional biodiversity, calledγ-diversity, into the average within-site biodiversity, or α-diversity, and the biodiversity among sites, or β-diversity. According to Whittaker [8], who introduced the terminology, β-diversity is the ratio between γ-diversity and α-diversity. This is the multiplicative partitioning of diversity. According to MacArthur [5], MacArthur and Wilson [7], and Lande [17], β-diversity is the difference between γ-diversity and α-diversity. This is the additive partitioning of diversity. All these diversities, namely, α-diversity, β-diversity, and γ-diversity, should be nonnegative numbers. In general, a measure of biodiversity ought to be nonnegative, in which case the corresponding α-diversity and γ-diversity, calculated by using such a measure, are nonnegative as well, as they should be. But the corresponding β-diversity is also nonnegative, in the additive partitioning of the biodiversity, or larger than 1, in the multiplicative partitioning of biodiversity, if the measure of biodiversity used is a concave function of the distribution of the relative abundance of species.The best known measures of biodiversity are Shannon’s entropy and the Gini-Simpson index. Both of them measure the biodiversity taking into account only the relative abundance of species. The widely used Rao’s index measures the dissimilarity between species taking into account not only the relative abundance of species but also a distance between species, such as the phylogenetic distance, for instance. Both Shannon’s entropy and the classic Gini-Simpson index satisfy the mathematical properties (nonnegativity and concavity) that allow them to be successfully used in the additive partitioning of biodiversity. Unfortunately, as was pointed out recently [12, 13], these two measures do not give good results when the number of species is very large. On the other hand, Rao’s index of dissimilarity is not a concave function of the relative abundance of species for arbitrary distances between species and, consequently, can be used in the additive partitioning of biodiversity only for some particular distance matrices, but not in general. The main objective of this paper is to show that the weighted Gini-Simpson quadratic index GSD given by (11), which is a generalization of the classic Gini-Simpson index GS to the pairs of species, is a suitable measure for use in the standard additive partitioning of biodiversity because, unlike the commonly used Rao’s index of dissimilarity R, it is a concave function of the relative abundance of the pairs of species. Unlike the classic Gini-Simpson index GS, the weighted Gini-Simpson quadratic index GSn,D behaves very well when the number of species is very large. The index GSn,D may be generalized to get the diversity measure GSn.v,D, given by (13), which takes into account not only the number of species, the relative abundance of the pairs of species, and the matrix D of the distances between species, but also a vector v of values assigned to the individual species, such as some conservation values for instance. The algebraic transformations (6) and (17) of the weighted Gini-Simpson quadratic indices GSw, for single species, and GSW, for pairs of species, given by (1) and (9), respectively, provide two measures of biodiversity which are suitable for use in the multiplicative partitioning of biodiversity. A detailed numerical example shows how the formulas should be implemented in applications.From a practical point of view, the new weighted Gini-Simpson measure of biodiversityGSn,v,D, which is a positive concave function of the relative abundance of the pairs of species, which essentially depends both on the matrix D of the distances between species and on the conservation values v of the species, is proposed as a suitable and improved replacement for the well-known Rao’s index in the partitioning of biodiversity.



---

*Source: 478728-2012-02-16.xml*