Misplaced Pages

Qualitative variation

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
Statistical dispersion in nominal distributions
This article has multiple issues. Please help improve it or discuss these issues on the talk page. (Learn how and when to remove these messages)
This article needs additional citations for verification. Please help improve this article by adding citations to reliable sources. Unsourced material may be challenged and removed.
Find sources: "Qualitative variation" – news · newspapers · books · scholar · JSTOR (April 2024) (Learn how and when to remove this message)
This article may contain an excessive amount of intricate detail that may interest only a particular audience. Please help by spinning off or relocating any relevant information, and removing excessive detail that may be against Misplaced Pages's inclusion policy. (April 2024) (Learn how and when to remove this message)
This article needs attention from an expert in Statistics. Please add a reason or a talk parameter to this template to explain the issue with the article. WikiProject Statistics may be able to help recruit an expert. (April 2024)
(Learn how and when to remove this message)

An index of qualitative variation (IQV) is a measure of statistical dispersion in nominal distributions. Examples include the variation ratio or the information entropy.

Properties

There are several types of indices used for the analysis of nominal data. Several are standard statistics that are used elsewhere - range, standard deviation, variance, mean deviation, coefficient of variation, median absolute deviation, interquartile range and quartile deviation.

In addition to these several statistics have been developed with nominal data in mind. A number have been summarized and devised by Wilcox (Wilcox 1967), (Wilcox 1973), who requires the following standardization properties to be satisfied:

  • Variation varies between 0 and 1.
  • Variation is 0 if and only if all cases belong to a single category.
  • Variation is 1 if and only if cases are evenly divided across all categories.

In particular, the value of these standardized indices does not depend on the number of categories or number of samples.

For any index, the closer to uniform the distribution, the larger the variance, and the larger the differences in frequencies across categories, the smaller the variance.

Indices of qualitative variation are then analogous to information entropy, which is minimized when all cases belong to a single category and maximized in a uniform distribution. Indeed, information entropy can be used as an index of qualitative variation.

One characterization of a particular index of qualitative variation (IQV) is as a ratio of observed differences to maximum differences.

Wilcox's indexes

Wilcox gives a number of formulae for various indices of QV (Wilcox 1973), the first, which he designates DM for "Deviation from the Mode", is a standardized form of the variation ratio, and is analogous to variance as deviation from the mean.

ModVR

The formula for the variation around the mode (ModVR) is derived as follows:

M = i = 1 K ( f m f i ) {\displaystyle M=\sum _{i=1}^{K}(f_{m}-f_{i})}

where fm is the modal frequency, K is the number of categories and fi is the frequency of the i group.

This can be simplified to

M = K f m N {\displaystyle M=Kf_{m}-N}

where N is the total size of the sample.

Freeman's index (or variation ratio) is

v = 1 f m N {\displaystyle v=1-{\frac {f_{m}}{N}}}

This is related to M as follows:

( f m N ) 1 K N K ( K 1 ) N = M N ( K 1 ) {\displaystyle {\frac {({\frac {f_{m}}{N}})-{\frac {1}{K}}}{{\frac {N}{K}}{\frac {(K-1)}{N}}}}={\frac {M}{N(K-1)}}}

The ModVR is defined as

ModVR = 1 K f m N N ( K 1 ) = K ( N f m ) N ( K 1 ) = K v K 1 {\displaystyle \operatorname {ModVR} =1-{\frac {Kf_{m}-N}{N(K-1)}}={\frac {K(N-f_{m})}{N(K-1)}}={\frac {Kv}{K-1}}}

where v is Freeman's index.

Low values of ModVR correspond to small amount of variation and high values to larger amounts of variation.

When K is large, ModVR is approximately equal to Freeman's index v.

RanVR

This is based on the range around the mode. It is defined to be

RanVR = 1 f m f l f m = f l f m {\displaystyle \operatorname {RanVR} =1-{\frac {f_{m}-f_{l}}{f_{m}}}={\frac {f_{l}}{f_{m}}}}

where fm is the modal frequency and fl is the lowest frequency.

AvDev

This is an analog of the mean deviation. It is defined as the arithmetic mean of the absolute differences of each value from the mean.

AvDev = 1 1 2 N K K 1 i = 1 K | f i N K | {\displaystyle \operatorname {AvDev} =1-{\frac {1}{2N}}{\frac {K}{K-1}}\sum _{i=1}^{K}\left|f_{i}-{\frac {N}{K}}\right|}

MNDif

This is an analog of the mean difference - the average of the differences of all the possible pairs of variate values, taken regardless of sign. The mean difference differs from the mean and standard deviation because it is dependent on the spread of the variate values among themselves and not on the deviations from some central value.

MNDif = 1 1 N ( K 1 ) i = 1 K 1 j = i + 1 K | f i f j | {\displaystyle \operatorname {MNDif} =1-{\frac {1}{N(K-1)}}\sum _{i=1}^{K-1}\sum _{j=i+1}^{K}|f_{i}-f_{j}|}

where fi and fj are the i and j frequencies respectively.

The MNDif is the Gini coefficient applied to qualitative data.

VarNC

This is an analog of the variance.

VarNC = 1 1 N 2 K K 1 ( f i N K ) 2 {\displaystyle \operatorname {VarNC} =1-{\frac {1}{N^{2}}}{\frac {K}{K-1}}\sum \left(f_{i}-{\frac {N}{K}}\right)^{2}}

It is the same index as Mueller and Schussler's Index of Qualitative Variation and Gibbs' M2 index.

It is distributed as a chi square variable with K – 1 degrees of freedom.

StDev

Wilson has suggested two versions of this statistic.

The first is based on AvDev.

StDev 1 = 1 i = 1 K ( f i N K ) 2 ( N N K ) 2 + ( K 1 ) ( N K ) 2 {\displaystyle \operatorname {StDev} _{1}=1-{\sqrt {\frac {\sum _{i=1}^{K}\left(f_{i}-{\frac {N}{K}}\right)^{2}}{\left(N-{\frac {N}{K}}\right)^{2}+(K-1)\left({\frac {N}{K}}\right)^{2}}}}}

The second is based on MNDif

StDev 2 = 1 i = 1 K 1 j = i + 1 K ( f i f j ) 2 N 2 ( K 1 ) {\displaystyle \operatorname {StDev} _{2}=1-{\sqrt {\frac {\sum _{i=1}^{K-1}\sum _{j=i+1}^{K}(f_{i}-f_{j})^{2}}{N^{2}(K-1)}}}}

HRel

This index was originally developed by Claude Shannon for use in specifying the properties of communication channels.

HRel = p i log 2 p i log 2 K {\displaystyle \operatorname {HRel} ={\frac {-\sum p_{i}\log _{2}p_{i}}{\log _{2}K}}}

where pi = fi / N.

This is equivalent to information entropy divided by the log 2 ( K ) {\displaystyle \log _{2}(K)} and is useful for comparing relative variation between frequency tables of multiple sizes.

B index

Wilcox adapted a proposal of Kaiser based on the geometric mean and created the B' index. The B index is defined as

B = 1 1 [ i = 1 k f i K N k ] 2 {\displaystyle B=1-{\sqrt {1-\left{\prod _{i=1}^{k}{\frac {f_{i}K}{N}}}}\,\right]^{2}}}}

R packages

Several of these indices have been implemented in the R language.

Gibb's indices and related formulae

Gibbs & Poston Jr (1975) proposed six indexes.

M1

The unstandardized index (M1) (Gibbs & Poston Jr 1975, p. 471) is

M 1 = 1 i = 1 K p i 2 {\displaystyle M1=1-\sum _{i=1}^{K}p_{i}^{2}}

where K is the number of categories and p i = f i / N {\displaystyle p_{i}=f_{i}/N} is the proportion of observations that fall in a given category i.

M1 can be interpreted as one minus the likelihood that a random pair of samples will belong to the same category, so this formula for IQV is a standardized likelihood of a random pair falling in the same category. This index has also referred to as the index of differentiation, the index of sustenance differentiation and the geographical differentiation index depending on the context it has been used in.

M2

A second index is the M2 (Gibbs & Poston Jr 1975, p. 472) is:

M 2 = K K 1 ( 1 i = 1 K p i 2 ) {\displaystyle M2={\frac {K}{K-1}}\left(1-\sum _{i=1}^{K}p_{i}^{2}\right)}

where K is the number of categories and p i = f i / N {\displaystyle p_{i}=f_{i}/N} is the proportion of observations that fall in a given category i. The factor of K K 1 {\displaystyle {\frac {K}{K-1}}} is for standardization.

M1 and M2 can be interpreted in terms of variance of a multinomial distribution (Swanson 1976) (there called an "expanded binomial model"). M1 is the variance of the multinomial distribution and M2 is the ratio of the variance of the multinomial distribution to the variance of a binomial distribution.

M4

The M4 index is

M 4 = i = 1 K | X i m | 2 i = 1 K X i {\displaystyle M4={\frac {\sum _{i=1}^{K}|X_{i}-m|}{2\sum _{i=1}^{K}X_{i}}}}

where m is the mean.

M6

The formula for M6 is

M 6 = K [ 1 i = 1 K | X i m | 2 N ] {\displaystyle M6=K\left}

· where K is the number of categories, Xi is the number of data points in the i category, N is the total number of data points, || is the absolute value (modulus) and

m = i = 1 K X i N {\displaystyle m={\frac {\sum _{i=1}^{K}X_{i}}{N}}}

This formula can be simplified

M 6 = K [ 1 i = 1 K | p i 1 N | 2 ] {\displaystyle M6=K\left}

where pi is the proportion of the sample in the i category.

In practice M1 and M6 tend to be highly correlated which militates against their combined use.

Related indices

The sum

i = 1 K p i 2 {\displaystyle \sum _{i=1}^{K}p_{i}^{2}}

has also found application. This is known as the Simpson index in ecology and as the Herfindahl index or the Herfindahl-Hirschman index (HHI) in economics. A variant of this is known as the Hunter–Gaston index in microbiology

In linguistics and cryptanalysis this sum is known as the repeat rate. The incidence of coincidence (IC) is an unbiased estimator of this statistic

IC = f i ( f i 1 ) n ( n 1 ) {\displaystyle \operatorname {IC} =\sum {\frac {f_{i}(f_{i}-1)}{n(n-1)}}}

where fi is the count of the i grapheme in the text and n is the total number of graphemes in the text.

M1

The M1 statistic defined above has been proposed several times in a number of different settings under a variety of names. These include Gini's index of mutability, Simpson's measure of diversity, Bachi's index of linguistic homogeneity, Mueller and Schuessler's index of qualitative variation, Gibbs and Martin's index of industry diversification, Lieberson's index. and Blau's index in sociology, psychology and management studies. The formulation of all these indices are identical.

Simpson's D is defined as

D = 1 i = 1 K n i ( n i 1 ) n ( n 1 ) {\displaystyle D=1-\sum _{i=1}^{K}{\frac {n_{i}(n_{i}-1)}{n(n-1)}}}

where n is the total sample size and ni is the number of items in the i category.

For large n we have

u 1 i = 1 K p i 2 {\displaystyle u\sim 1-\sum _{i=1}^{K}p_{i}^{2}}

Another statistic that has been proposed is the coefficient of unalikeability which ranges between 0 and 1.

u = c ( x , y ) n 2 n {\displaystyle u={\frac {c(x,y)}{n^{2}-n}}}

where n is the sample size and c(x,y) = 1 if x and y are unalike and 0 otherwise.

For large n we have

u 1 i = 1 K p i 2 {\displaystyle u\sim 1-\sum _{i=1}^{K}p_{i}^{2}}

where K is the number of categories.

Another related statistic is the quadratic entropy

H 2 = 2 ( 1 i = 1 K p i 2 ) {\displaystyle H^{2}=2\left(1-\sum _{i=1}^{K}p_{i}^{2}\right)}

which is itself related to the Gini index.

M2

Greenberg's monolingual non weighted index of linguistic diversity is the M2 statistic defined above.

M7

Another index – the M7 – was created based on the M4 index of Gibbs & Poston Jr (1975)

M 7 = i = 1 K j = 1 L | R i R | 2 R i {\displaystyle M7={\frac {\sum _{i=1}^{K}\sum _{j=1}^{L}|R_{i}-R|}{2\sum R_{i}}}}

where

R i j = O i j E i j = O i j n i p j {\displaystyle R_{ij}={\frac {O_{ij}}{E_{ij}}}={\frac {O_{ij}}{n_{i}p_{j}}}}

and

R = i = 1 K j = 1 L R i j i = 1 K n i {\displaystyle R={\frac {\sum _{i=1}^{K}\sum _{j=1}^{L}R_{ij}}{\sum _{i=1}^{K}n_{i}}}}

where K is the number of categories, L is the number of subtypes, Oij and Eij are the number observed and expected respectively of subtype j in the i category, ni is the number in the i category and pj is the proportion of subtype j in the complete sample.

Note: This index was designed to measure women's participation in the work place: the two subtypes it was developed for were male and female.

Other single sample indices

These indices are summary statistics of the variation within the sample.

Berger–Parker index

The Berger–Parker index, named after Wolfgang H. Berger and Frances Lawrence Parker, equals the maximum p i {\displaystyle p_{i}} value in the dataset, i.e. the proportional abundance of the most abundant type. This corresponds to the weighted generalized mean of the p i {\displaystyle p_{i}} values when q approaches infinity, and hence equals the inverse of true diversity of order infinity (1/D).

Brillouin index of diversity

This index is strictly applicable only to entire populations rather than to finite samples. It is defined as

I B = log ( N ! ) i = 1 K ( log ( n i ! ) ) N {\displaystyle I_{B}={\frac {\log(N!)-\sum _{i=1}^{K}(\log(n_{i}!))}{N}}}

where N is total number of individuals in the population, ni is the number of individuals in the i category and N! is the factorial of N. Brillouin's index of evenness is defined as

E B = I B / I B ( max ) {\displaystyle E_{B}=I_{B}/I_{B(\max )}}

where IB(max) is the maximum value of IB.

Hill's diversity numbers

Hill suggested a family of diversity numbers

N a = 1 [ i = 1 K p i a ] a 1 {\displaystyle N_{a}={\frac {1}{\left^{a-1}}}}

For given values of a, several of the other indices can be computed

  • a = 0: Na = species richness
  • a = 1: Na = Shannon's index
  • a = 2: Na = 1/Simpson's index (without the small sample correction)
  • a = 3: Na = 1/Berger–Parker index

Hill also suggested a family of evenness measures

E a , b = N a N b {\displaystyle E_{a,b}={\frac {N_{a}}{N_{b}}}}

where a > b.

Hill's E4 is

E 4 = N 2 N 1 {\displaystyle E_{4}={\frac {N_{2}}{N_{1}}}}

Hill's E5 is

E 5 = N 2 1 N 1 1 {\displaystyle E_{5}={\frac {N_{2}-1}{N_{1}-1}}}

Margalef's index

I Marg = S 1 log e N {\displaystyle I_{\text{Marg}}={\frac {S-1}{\log _{e}N}}}

where S is the number of data types in the sample and N is the total size of the sample.

Menhinick's index

I M e n = S N {\displaystyle I_{\mathrm {Men} }={\frac {S}{\sqrt {N}}}}

where S is the number of data types in the sample and N is the total size of the sample.

In linguistics this index is the identical with the Kuraszkiewicz index (Guiard index) where S is the number of distinct words (types) and N is the total number of words (tokens) in the text being examined. This index can be derived as a special case of the Generalised Torquist function.

Q statistic

This is a statistic invented by Kempton and Taylor. and involves the quartiles of the sample. It is defined as

Q = 1 2 ( n R 1 + n R 2 ) + j = R 1 + 1 R 2 1 n j log ( R 2 / R 1 ) {\displaystyle Q={\frac {{\frac {1}{2}}(n_{R1}+n_{R2})+\sum _{j=R_{1}+1}^{R_{2}-1}n_{j}}{\log(R_{2}/R_{1})}}}

where R1 and R2 are the 25% and 75% quartiles respectively on the cumulative species curve, nj is the number of species in the jth category, nRi is the number of species in the class where Ri falls (i = 1 or 2).

Shannon–Wiener index

This is taken from information theory

H = log e N 1 N n i p i log ( p i ) {\displaystyle H=\log _{e}N-{\frac {1}{N}}\sum n_{i}p_{i}\log(p_{i})}

where N is the total number in the sample and pi is the proportion in the i category.

In ecology where this index is commonly used, H usually lies between 1.5 and 3.5 and only rarely exceeds 4.0.

An approximate formula for the standard deviation (SD) of H is

SD ( H ) = 1 N [ p i [ log e ( p i ) ] 2 H 2 ] {\displaystyle \operatorname {SD} (H)={\frac {1}{N}}\left^{2}-H^{2}\right]}

where pi is the proportion made up by the i category and N is the total in the sample.

A more accurate approximate value of the variance of H(var(H)) is given by

var ( H ) = p i [ log ( p i ) ] 2 [ p i log ( p i ) ] 2 N + K 1 2 N 2 + 1 + p i 2 p i 1 log ( p i ) + p i 1 p i log ( p i ) 6 N 3 {\displaystyle \operatorname {var} (H)={\frac {\sum p_{i}^{2}-\left^{2}}{N}}+{\frac {K-1}{2N^{2}}}+{\frac {-1+\sum p_{i}^{2}-\sum p_{i}^{-1}\log(p_{i})+\sum p_{i}^{-1}\sum p_{i}\log(p_{i})}{6N^{3}}}}

where N is the sample size and K is the number of categories.

A related index is the Pielou J defined as

J = H log e ( S ) {\displaystyle J={\frac {H}{\log _{e}(S)}}}

One difficulty with this index is that S is unknown for a finite sample. In practice S is usually set to the maximum present in any category in the sample.

Rényi entropy

The Rényi entropy is a generalization of the Shannon entropy to other values of q than unity. It can be expressed:

q H = 1 1 q ln ( i = 1 K p i q ) {\displaystyle {}^{q}H={\frac {1}{1-q}}\;\ln \left(\sum _{i=1}^{K}p_{i}^{q}\right)}

which equals

q H = ln ( 1 i = 1 K p i p i q 1 q 1 ) = ln ( q D ) {\displaystyle {}^{q}H=\ln \left({1 \over {\sqrt{\sum _{i=1}^{K}p_{i}p_{i}^{q-1}}}}\right)=\ln({}^{q}\!D)}

This means that taking the logarithm of true diversity based on any value of q gives the Rényi entropy corresponding to the same value of q.

The value of q D {\displaystyle {}^{q}\!D} is also known as the Hill number.

McIntosh's D and E

McIntosh proposed measure of diversity:

I = i = 1 K n i 2 {\displaystyle I={\sqrt {\sum _{i=1}^{K}n_{i}^{2}}}}

where ni is the number in the i category and K is the number of categories.

He also proposed several normalized versions of this index. First is D:

D = N I N N {\displaystyle D={\frac {N-I}{N-{\sqrt {N}}}}}

where N is the total sample size.

This index has the advantage of expressing the observed diversity as a proportion of the absolute maximum diversity at a given N.

Another proposed normalization is E — ratio of observed diversity to maximum possible diversity of a given N and K (i.e., if all species are equal in number of individuals):

E = N I N N K {\displaystyle E={\frac {N-I}{N-{\frac {N}{K}}}}}

Fisher's alpha

This was the first index to be derived for diversity.

K = α ln ( 1 + N α ) {\displaystyle K=\alpha \ln(1+{\frac {N}{\alpha }})}

where K is the number of categories and N is the number of data points in the sample. Fisher's α has to be estimated numerically from the data.

The expected number of individuals in the r category where the categories have been placed in increasing size is

E ( n r ) = α X r r {\displaystyle \operatorname {E} (n_{r})=\alpha {\frac {X^{r}}{r}}}

where X is an empirical parameter lying between 0 and 1. While X is best estimated numerically an approximate value can be obtained by solving the following two equations

N = α X 1 X {\displaystyle N={\frac {\alpha X}{1-X}}}
K = α ln ( 1 X ) {\displaystyle K=-\alpha \ln(1-X)}

where K is the number of categories and N is the total sample size.

The variance of α is approximately

var ( α ) = α ln ( X ) ( 1 X ) {\displaystyle \operatorname {var} (\alpha )={\frac {\alpha }{\ln(X)(1-X)}}}

Strong's index

This index (Dw) is the distance between the Lorenz curve of species distribution and the 45 degree line. It is closely related to the Gini coefficient.

In symbols it is

D w = m a x [ c i K i N ] {\displaystyle D_{w}=max}

where max() is the maximum value taken over the N data points, K is the number of categories (or species) in the data set and ci is the cumulative total up and including the ith category.

Simpson's E

This is related to Simpson's D and is defined as

E = 1 / D K {\displaystyle E={\frac {1/D}{K}}}

where D is Simpson's D and K is the number of categories in the sample.

Smith & Wilson's indices

Smith and Wilson suggested a number of indices based on Simpson's D.

E 1 = 1 D 1 1 K {\displaystyle E_{1}={\frac {1-D}{1-{\frac {1}{K}}}}}
E 2 = log e ( D ) log e ( K ) {\displaystyle E_{2}={\frac {\log _{e}(D)}{\log _{e}(K)}}}

where D is Simpson's D and K is the number of categories.

Heip's index

E = e H 1 K 1 {\displaystyle E={\frac {e^{H}-1}{K-1}}}

where H is the Shannon entropy and K is the number of categories.

This index is closely related to Sheldon's index which is

E = e H K {\displaystyle E={\frac {e^{H}}{K}}}

where H is the Shannon entropy and K is the number of categories.

Camargo's index

This index was created by Camargo in 1993.

E = 1 i = 1 K j = i + 1 K p i p j K {\displaystyle E=1-\sum _{i=1}^{K}\sum _{j=i+1}^{K}{\frac {p_{i}-p_{j}}{K}}}

where K is the number of categories and pi is the proportion in the i category.

Smith and Wilson's B

This index was proposed by Smith and Wilson in 1996.

B = 1 2 π arctan ( θ ) {\displaystyle B=1-{\frac {2}{\pi }}\arctan(\theta )}

where θ is the slope of the log(abundance)-rank curve.

Nee, Harvey, and Cotgreave's index

This is the slope of the log(abundance)-rank curve.

Bulla's E

There are two versions of this index - one for continuous distributions (Ec) and the other for discrete (Ed).

E c = O 1 K 1 1 K {\displaystyle E_{c}={\frac {O-{\frac {1}{K}}}{1-{\frac {1}{K}}}}}
E d = O 1 K K 1 N 1 1 K K 1 N {\displaystyle E_{d}={\frac {O-{\frac {1}{K}}-{\frac {K-1}{N}}}{1-{\frac {1}{K}}-{\frac {K-1}{N}}}}}

where

O = 1 1 2 | p i 1 K | {\displaystyle O=1-{\frac {1}{2}}\left|p_{i}-{\frac {1}{K}}\right|}

is the Schoener–Czekanoski index, K is the number of categories and N is the sample size.

Horn's information theory index

This index (Rik) is based on Shannon's entropy. It is defined as

R i k = H max H o b s H max H min {\displaystyle R_{ik}={\frac {H_{\max }-H_{\mathrm {obs} }}{H_{\max }-H_{\min }}}}

where

X = x i j {\displaystyle X=\sum x_{ij}}
X = x k j {\displaystyle X=\sum x_{kj}}
H ( X ) = x i j X log X x i j {\displaystyle H(X)=\sum {\frac {x_{ij}}{X}}\log {\frac {X}{x_{ij}}}}
H ( Y ) = x k j Y log Y x k j {\displaystyle H(Y)=\sum {\frac {x_{kj}}{Y}}\log {\frac {Y}{x_{kj}}}}
H min = X X + Y H ( X ) + Y X + Y H ( Y ) {\displaystyle H_{\min }={\frac {X}{X+Y}}H(X)+{\frac {Y}{X+Y}}H(Y)}
H max = ( x i j X + Y log X + Y x i j + x k j X + Y log X + Y x k j ) {\displaystyle H_{\max }=\sum \left({\frac {x_{ij}}{X+Y}}\log {\frac {X+Y}{x_{ij}}}+{\frac {x_{kj}}{X+Y}}\log {\frac {X+Y}{x_{kj}}}\right)}
H o b s = x i j + x k j X + Y log X + Y x i j + x k j {\displaystyle H_{\mathrm {obs} }=\sum {\frac {x_{ij}+x_{kj}}{X+Y}}\log {\frac {X+Y}{x_{ij}+x_{kj}}}}

In these equations xij and xkj are the number of times the j data type appears in the i or k sample respectively.

Rarefaction index

In a rarefied sample a random subsample n in chosen from the total N items. In this sample some groups may be necessarily absent from this subsample. Let X n {\displaystyle X_{n}} be the number of groups still present in the subsample of n items. X n {\displaystyle X_{n}} is less than K the number of categories whenever at least one group is missing from this subsample.

The rarefaction curve, f n {\displaystyle f_{n}} is defined as:

f n = E [ X n ] = K ( N n ) 1 i = 1 K ( N N i n ) {\displaystyle f_{n}=\operatorname {E} =K-{\binom {N}{n}}^{-1}\sum _{i=1}^{K}{\binom {N-N_{i}}{n}}}

Note that 0 ≤ f(n) ≤ K.

Furthermore,

f ( 0 ) = 0 ,   f ( 1 ) = 1 ,   f ( N ) = K . {\displaystyle f(0)=0,\ f(1)=1,\ f(N)=K.}

Despite being defined at discrete values of n, these curves are most frequently displayed as continuous functions.

This index is discussed further in Rarefaction (ecology).

Caswell's V

This is a z type statistic based on Shannon's entropy.

V = H E ( H ) SD ( H ) {\displaystyle V={\frac {H-\operatorname {E} (H)}{\operatorname {SD} (H)}}}

where H is the Shannon entropy, E(H) is the expected Shannon entropy for a neutral model of distribution and SD(H) is the standard deviation of the entropy. The standard deviation is estimated from the formula derived by Pielou

S D ( H ) = 1 N [ p i [ log e ( p i ) ] 2 H 2 ] {\displaystyle SD(H)={\frac {1}{N}}\left^{2}-H^{2}\right]}

where pi is the proportion made up by the i category and N is the total in the sample.

Lloyd & Ghelardi's index

This is

I L G = K K {\displaystyle I_{LG}={\frac {K}{K'}}}

where K is the number of categories and K' is the number of categories according to MacArthur's broken stick model yielding the observed diversity.

Average taxonomic distinctness index

This index is used to compare the relationship between hosts and their parasites. It incorporates information about the phylogenetic relationship amongst the host species.

S T D = 2 i < j ω i j s ( s 1 ) {\displaystyle S_{TD}=2{\frac {\sum \sum _{i<j}\omega _{ij}}{s(s-1)}}}

where s is the number of host species used by a parasite and ωij is the taxonomic distinctness between host species i and j.

Index of qualitative variation

Several indices with this name have been proposed.

One of these is

I Q V = K ( 100 2 i = 1 K p i 2 ) 100 2 ( K 1 ) = K K 1 ( 1 i = 1 K ( p i / 100 ) 2 ) {\displaystyle IQV={\frac {K(100^{2}-\sum _{i=1}^{K}p_{i}^{2})}{100^{2}(K-1)}}={\frac {K}{K-1}}(1-\sum _{i=1}^{K}(p_{i}/100)^{2})}

where K is the number of categories and pi is the proportion of the sample that lies in the i category.

Theil's H

This index is also known as the multigroup entropy index or the information theory index. It was proposed by Theil in 1972. The index is a weighted average of the samples entropy.

Let

E a = i = 1 a p i l o g ( p i ) {\displaystyle E_{a}=\sum _{i=1}^{a}p_{i}log(p_{i})}

and

H = i = 1 r n i ( E E i ) N E {\displaystyle H=\sum _{i=1}^{r}{\frac {n_{i}(E-E_{i})}{NE}}}

where pi is the proportion of type i in the a sample, r is the total number of samples, ni is the size of the i sample, N is the size of the population from which the samples were obtained and E is the entropy of the population.

Indices for comparison of two or more data types within a single sample

Several of these indexes have been developed to document the degree to which different data types of interest may coexist within a geographic area.

Index of dissimilarity

Let A and B be two types of data item. Then the index of dissimilarity is

D = 1 2 i = 1 K | A i A B i B | {\displaystyle D={\frac {1}{2}}\sum _{i=1}^{K}\left|{\frac {A_{i}}{A}}-{\frac {B_{i}}{B}}\right|}

where

A = i = 1 K A i {\displaystyle A=\sum _{i=1}^{K}A_{i}}
B = i = 1 K B i {\displaystyle B=\sum _{i=1}^{K}B_{i}}

Ai is the number of data type A at sample site i, Bi is the number of data type B at sample site i, K is the number of sites sampled and || is the absolute value.

This index is probably better known as the index of dissimilarity (D). It is closely related to the Gini index.

This index is biased as its expectation under a uniform distribution is > 0.

A modification of this index has been proposed by Gorard and Taylor. Their index (GT) is

G T = D ( 1 A A + B ) {\displaystyle GT=D\left(1-{\frac {A}{A+B}}\right)}

Index of segregation

The index of segregation (IS) is

S I = 1 2 i = 1 K | A i A t i A i T A | {\displaystyle SI={\frac {1}{2}}\sum _{i=1}^{K}\left|{\frac {A_{i}}{A}}-{\frac {t_{i}-A_{i}}{T-A}}\right|}

where

A = i = 1 K A i {\displaystyle A=\sum _{i=1}^{K}A_{i}}
T = i = 1 K t i {\displaystyle T=\sum _{i=1}^{K}t_{i}}

and K is the number of units, Ai and ti is the number of data type A in unit i and the total number of all data types in unit i.

Hutchen's square root index

This index (H) is defined as

H = 1 i = 1 K j = 1 i p i p j {\displaystyle H=1-\sum _{i=1}^{K}\sum _{j=1}^{i}{\sqrt {p_{i}p_{j}}}}

where pi is the proportion of the sample composed of the i variate.

Lieberson's isolation index

See also: isolation index

This index ( Lxy ) was invented by Lieberson in 1981.

L x y = 1 N i = 1 K X i Y i X t o t {\displaystyle L_{xy}={\frac {1}{N}}\sum _{i=1}^{K}{\frac {X_{i}Y_{i}}{X_{\mathrm {tot} }}}}

where Xi and Yi are the variables of interest at the i site, K is the number of sites examined and Xtot is the total number of variate of type X in the study.

Bell's index

This index is defined as

I R = p x x p x 1 p x {\displaystyle I_{R}={\frac {p_{xx}-p_{x}}{1-p_{x}}}}

where px is the proportion of the sample made up of variates of type X and

p x x = i = 1 K x i p i N x {\displaystyle p_{xx}={\frac {\sum _{i=1}^{K}x_{i}p_{i}}{N_{x}}}}

where Nx is the total number of variates of type X in the study, K is the number of samples in the study and xi and pi are the number of variates and the proportion of variates of type X respectively in the i sample.

Index of isolation

The index of isolation is

I I = i = 1 K A i A A i t i {\displaystyle II=\sum _{i=1}^{K}{\frac {A_{i}}{A}}{\frac {A_{i}}{t_{i}}}}

where K is the number of units in the study, Ai and ti is the number of units of type A and the number of all units in ith sample.

A modified index of isolation has also been proposed

M I I = I I A T 1 A T {\displaystyle MII={\frac {II-{\frac {A}{T}}}{1-{\frac {A}{T}}}}}

The MII lies between 0 and 1.

Gorard's index of segregation

This index (GS) is defined as

G S = 1 2 i = 1 K | A i A t i T | {\displaystyle GS={\frac {1}{2}}\sum _{i=1}^{K}\left|{\frac {A_{i}}{A}}-{\frac {t_{i}}{T}}\right|}

where

A = i = 1 K A i {\displaystyle A=\sum _{i=1}^{K}A_{i}}
T = i = 1 K t i {\displaystyle T=\sum _{i=1}^{K}t_{i}}

and Ai and ti are the number of data items of type A and the total number of items in the i sample.

Index of exposure

This index is defined as

I E = i = 1 K A i A B i t i {\displaystyle IE=\sum _{i=1}^{K}{\frac {A_{i}}{A}}{\frac {B_{i}}{t_{i}}}}

where

A = i = 1 K A i {\displaystyle A=\sum _{i=1}^{K}A_{i}}

and Ai and Bi are the number of types A and B in the i category and ti is the total number of data points in the i category.

Ochiai index

This is a binary form of the cosine index. It is used to compare presence/absence data of two data types (here A and B). It is defined as

O = a ( a + b ) ( a + c ) {\displaystyle O={\frac {a}{\sqrt {(a+b)(a+c)}}}}

where a is the number of sample units where both A and B are found, b is number of sample units where A but not B occurs and c is the number of sample units where type B is present but not type A.

Kulczyński's coefficient

This coefficient was invented by Stanisław Kulczyński in 1927 and is an index of association between two types (here A and B). It varies in value between 0 and 1. It is defined as

K = a 2 ( 1 a + b + 1 a + c ) {\displaystyle K={\frac {a}{2}}\left({\frac {1}{a+b}}+{\frac {1}{a+c}}\right)}

where a is the number of sample units where type A and type B are present, b is the number of sample units where type A but not type B is present and c is the number of sample units where type B is present but not type A.

Yule's Q

This index was invented by Yule in 1900. It concerns the association of two different types (here A and B). It is defined as

Q = a d b c a d + b c {\displaystyle Q={\frac {ad-bc}{ad+bc}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. Q varies in value between -1 and +1. In the ordinal case Q is known as the Goodman-Kruskal γ.

Because the denominator potentially may be zero, Leinhert and Sporer have recommended adding +1 to a, b, c and d.

Yule's Y

This index is defined as

Y = a d b c a d + b c {\displaystyle Y={\frac {{\sqrt {ad}}-{\sqrt {bc}}}{{\sqrt {ad}}+{\sqrt {bc}}}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present.

Baroni–Urbani–Buser coefficient

This index was invented by Baroni-Urbani and Buser in 1976. It varies between 0 and 1 in value. It is defined as

B U B = a d + a a d + a + b + c = a d + a N + a d d = 1 N ( a d ) N + a d d {\displaystyle BUB={\frac {{\sqrt {ad}}+a}{{\sqrt {ad}}+a+b+c}}={\frac {{\sqrt {ad}}+a}{N+{\sqrt {ad}}-d}}=1-{\frac {N-(a-d)}{N+{\sqrt {ad}}-d}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size.

When d = 0, this index is identical to the Jaccard index.

Hamman coefficient

This coefficient is defined as

H = ( a + d ) ( b + c ) a + b + c + d = ( a + d ) ( b + c ) N {\displaystyle H={\frac {(a+d)-(b+c)}{a+b+c+d}}={\frac {(a+d)-(b+c)}{N}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size.

Rogers–Tanimoto coefficient

This coefficient is defined as

R T = a + d a + 2 ( b + c ) + d = a + d N + b + c {\displaystyle RT={\frac {a+d}{a+2(b+c)+d}}={\frac {a+d}{N+b+c}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size

Sokal–Sneath coefficient

This coefficient is defined as

S S = 2 ( a + d ) 2 ( a + d ) + b + c = 2 ( a + d ) N + a + d {\displaystyle SS={\frac {2(a+d)}{2(a+d)+b+c}}={\frac {2(a+d)}{N+a+d}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size.

Sokal's binary distance

This coefficient is defined as

S B D = b + c a + b + c + d = b + c N {\displaystyle SBD={\sqrt {\frac {b+c}{a+b+c+d}}}={\sqrt {\frac {b+c}{N}}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size.

Russel–Rao coefficient

This coefficient is defined as

R R = a a + b + c + d = a N {\displaystyle RR={\frac {a}{a+b+c+d}}={\frac {a}{N}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size.

Phi coefficient

This coefficient is defined as

φ = a d b c ( a + b ) ( a + c ) ( b + c ) ( c + d ) {\displaystyle \varphi ={\frac {ad-bc}{\sqrt {(a+b)(a+c)(b+c)(c+d)}}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present.

Soergel's coefficient

This coefficient is defined as

S = b + c b + c + d = b + c N a {\displaystyle S={\frac {b+c}{b+c+d}}={\frac {b+c}{N-a}}}

where b is the number of samples where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size.

Simpson's coefficient

This coefficient is defined as

S = a a + min ( b , c ) {\displaystyle S={\frac {a}{a+\min(b,c)}}}

where b is the number of samples where type A is present but not type B, c is the number of samples where type B is present but not type A.

Dennis' coefficient

This coefficient is defined as

D = a d b c ( a + b + c + d ) ( a + b ) ( a + c ) = a d b c N ( a + b ) ( a + c ) {\displaystyle D={\frac {ad-bc}{\sqrt {(a+b+c+d)(a+b)(a+c)}}}={\frac {ad-bc}{\sqrt {N(a+b)(a+c)}}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size.

Forbes' coefficient

This coefficient was proposed by Stephen Alfred Forbes in 1907. It is defined as

F = a N ( a + b ) ( a + c ) {\displaystyle F={\frac {aN}{(a+b)(a+c)}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size (N = a + b + c + d).

A modification of this coefficient which does not require the knowledge of d has been proposed by Alroy

F A = a ( n + n ) a ( n + n ) + 3 2 b c = 1 3 b c 2 a ( n + n ) + 3 b c {\displaystyle F_{A}={\frac {a(n+{\sqrt {n}})}{a(n+{\sqrt {n}})+{\frac {3}{2}}bc}}=1-{\frac {3bc}{2a(n+{\sqrt {n}})+3bc}}}

Where n = a + b + c.

Simple match coefficient

This coefficient is defined as

S M = a + d a + b + c + d = a + d N {\displaystyle SM={\frac {a+d}{a+b+c+d}}={\frac {a+d}{N}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size.

Fossum's coefficient

This coefficient is defined as

F = ( a + b + c + d ) ( a 0.5 ) 2 ( a + b ) ( a + c ) = N ( a 0.5 ) 2 ( a + b ) ( a + c ) {\displaystyle F={\frac {(a+b+c+d)(a-0.5)^{2}}{(a+b)(a+c)}}={\frac {N(a-0.5)^{2}}{(a+b)(a+c)}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size.

Stile's coefficient

This coefficient is defined as

S = log [ n ( | a d b c | n 2 ) 2 ( a + b ) ( a + c ) ( b + d ) ( c + d ) ] {\displaystyle S=\log \left}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A, d is the sample count where neither type A nor type B are present, n equals a + b + c + d and || is the modulus (absolute value) of the difference.

Michael's coefficient

This coefficient is defined as

M = 4 ( a d b c ) ( a + d ) 2 + ( b + c ) 2 {\displaystyle M={\frac {4(ad-bc)}{(a+d)^{2}+(b+c)^{2}}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present.

Peirce's coefficient

In 1884 Charles Peirce suggested the following coefficient

P = a b + b c a b + 2 b c + c d {\displaystyle P={\frac {ab+bc}{ab+2bc+cd}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present.

Hawkin–Dotson coefficient

In 1975 Hawkin and Dotson proposed the following coefficient

H D = 1 2 ( a a + b + c + d b + c + d ) = 1 2 ( a N d + d N a ) {\displaystyle HD={\frac {1}{2}}\left({\frac {a}{a+b+c}}+{\frac {d}{b+c+d}}\right)={\frac {1}{2}}\left({\frac {a}{N-d}}+{\frac {d}{N-a}}\right)}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size.

Benini coefficient

In 1901 Benini proposed the following coefficient

B = a ( a + b ) ( a + c ) a + min ( b , c ) ( a + b ) ( a + c ) {\displaystyle B={\frac {a-(a+b)(a+c)}{a+\min(b,c)-(a+b)(a+c)}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B and c is the number of samples where type B is present but not type A. Min(b, c) is the minimum of b and c.

Gilbert coefficient

Gilbert proposed the following coefficient

G = a ( a + b ) ( a + c ) a + b + c ( a + b ) ( a + c ) = a ( a + b ) ( a + c ) N ( a + b ) ( a + c ) d {\displaystyle G={\frac {a-(a+b)(a+c)}{a+b+c-(a+b)(a+c)}}={\frac {a-(a+b)(a+c)}{N-(a+b)(a+c)-d}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the sample count where neither type A nor type B are present. N is the sample size.

Gini index

The Gini index is

G = a ( a + b ) ( a + c ) ( 1 ( a + b ) 2 ) ( 1 ( a + c ) 2 ) {\displaystyle G={\frac {a-(a+b)(a+c)}{\sqrt {(1-(a+b)^{2})(1-(a+c)^{2})}}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B and c is the number of samples where type B is present but not type A.

Modified Gini index

The modified Gini index is

G M = a ( a + b ) ( a + c ) 1 | b c | 2 ( a + b ) ( a + c ) {\displaystyle G_{M}={\frac {a-(a+b)(a+c)}{1-{\frac {|b-c|}{2}}-(a+b)(a+c)}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B and c is the number of samples where type B is present but not type A.

Kuhn's index

Kuhn proposed the following coefficient in 1965

I = 2 ( a d b c ) K ( 2 a + b + c ) = 2 ( a d b c ) K ( N + a d ) {\displaystyle I={\frac {2(ad-bc)}{K(2a+b+c)}}={\frac {2(ad-bc)}{K(N+a-d)}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B and c is the number of samples where type B is present but not type A. K is a normalizing parameter. N is the sample size.

This index is also known as the coefficient of arithmetic means.

Eyraud index

Eyraud proposed the following coefficient in 1936

I = a ( a + b ) ( a + c ) ( a + c ) ( a + d ) ( b + d ) ( c + d ) {\displaystyle I={\frac {a-(a+b)(a+c)}{(a+c)(a+d)(b+d)(c+d)}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the number of samples where both A and B are not present.

Soergel distance

This is defined as

SD = b + c b + c + d = b + c N a {\displaystyle \operatorname {SD} ={\frac {b+c}{b+c+d}}={\frac {b+c}{N-a}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the number of samples where both A and B are not present. N is the sample size.

Tanimoto index

This is defined as

T I = 1 a b + c + d = 1 a N a = N 2 a N a {\displaystyle TI=1-{\frac {a}{b+c+d}}=1-{\frac {a}{N-a}}={\frac {N-2a}{N-a}}}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A and d is the number of samples where both A and B are not present. N is the sample size.

Piatetsky–Shapiro's index

This is defined as

P S I = a b c {\displaystyle PSI=a-bc}

where a is the number of samples where types A and B are both present, b is where type A is present but not type B, c is the number of samples where type B is present but not type A.

Indices for comparison between two or more samples

Czekanowski's quantitative index

This is also known as the Bray–Curtis index, Schoener's index, least common percentage index, index of affinity or proportional similarity. It is related to the Sørensen similarity index.

C Z I = min ( x i , x j ) ( x i + x j ) {\displaystyle CZI={\frac {\sum \min(x_{i},x_{j})}{\sum (x_{i}+x_{j})}}}

where xi and xj are the number of species in sites i and j respectively and the minimum is taken over the number of species in common between the two sites.

Canberra metric

The Canberra distance is a weighted version of the L1 metric. It was introduced by introduced in 1966 and refined in 1967 by G. N. Lance and W. T. Williams. It is used to define a distance between two vectors – here two sites with K categories within each site.

The Canberra distance d between vectors p and q in a K-dimensional real vector space is

d ( p , q ) = i = 1 n | p i q i | | p i | + | q i | {\displaystyle d(\mathbf {p} ,\mathbf {q} )=\sum _{i=1}^{n}{\frac {|p_{i}-q_{i}|}{|p_{i}|+|q_{i}|}}}

where pi and qi are the values of the i category of the two vectors.

Sorensen's coefficient of community

This is used to measure similarities between communities.

C C = 2 c s 1 + s 2 {\displaystyle CC={\frac {2c}{s_{1}+s_{2}}}}

where s1 and s2 are the number of species in community 1 and 2 respectively and c is the number of species common to both areas.

Jaccard's index

This is a measure of the similarity between two samples:

J = A A + B + C {\displaystyle J={\frac {A}{A+B+C}}}

where A is the number of data points shared between the two samples and B and C are the data points found only in the first and second samples respectively.

This index was invented in 1902 by the Swiss botanist Paul Jaccard.

Under a random distribution the expected value of J is

J = 1 A ( 1 A + B + C ) {\displaystyle J={\frac {1}{A}}\left({\frac {1}{A+B+C}}\right)}

The standard error of this index with the assumption of a random distribution is

S E ( J ) = A ( B + C ) N ( A + B + C ) 3 {\displaystyle SE(J)={\sqrt {\frac {A(B+C)}{N(A+B+C)^{3}}}}}

where N is the total size of the sample.

Dice's index

This is a measure of the similarity between two samples:

D = 2 A 2 A + B + C {\displaystyle D={\frac {2A}{2A+B+C}}}

where A is the number of data points shared between the two samples and B and C are the data points found only in the first and second samples respectively.

See also: Sørensen–Dice coefficient

Match coefficient

This is a measure of the similarity between two samples:

M = N B C N = 1 B + C N {\displaystyle M={\frac {N-B-C}{N}}=1-{\frac {B+C}{N}}}

where N is the number of data points in the two samples and B and C are the data points found only in the first and second samples respectively.

Morisita's index

Masaaki Morisita's index of dispersion ( Im ) is the scaled probability that two points chosen at random from the whole population are in the same sample. Higher values indicate a more clumped distribution.

I m = x ( x 1 ) n m ( m 1 ) {\displaystyle I_{m}={\frac {\sum x(x-1)}{nm(m-1)}}}

An alternative formulation is

I m = n x 2 x ( x ) 2 x {\displaystyle I_{m}=n{\frac {\sum x^{2}-\sum x}{\left(\sum x\right)^{2}-\sum x}}}

where n is the total sample size, m is the sample mean and x are the individual values with the sum taken over the whole sample. It is also equal to

I m = n   I M C n m 1 {\displaystyle I_{m}={\frac {n\ IMC}{nm-1}}}

where IMC is Lloyd's index of crowding.

This index is relatively independent of the population density but is affected by the sample size.

Morisita showed that the statistic

I m ( x 1 ) + n x {\displaystyle I_{m}\left(\sum x-1\right)+n-\sum x}

is distributed as a chi-squared variable with n − 1 degrees of freedom.

An alternative significance test for this index has been developed for large samples.

z = I m 1 2 / n m 2 {\displaystyle z={\frac {I_{m}-1}{2/nm^{2}}}}

where m is the overall sample mean, n is the number of sample units and z is the normal distribution abscissa. Significance is tested by comparing the value of z against the values of the normal distribution.

Morisita's overlap index

Morisita's overlap index is used to compare overlap among samples. The index is based on the assumption that increasing the size of the samples will increase the diversity because it will include different habitats

C D = 2 i = 1 S x i y i ( D x + D y ) X Y {\displaystyle C_{D}={\frac {2\sum _{i=1}^{S}x_{i}y_{i}}{(D_{x}+D_{y})XY}}}
xi is the number of times species i is represented in the total X from one sample.
yi is the number of times species i is represented in the total Y from another sample.
Dx and Dy are the Simpson's index values for the x and y samples respectively.
S is the number of unique species

CD = 0 if the two samples do not overlap in terms of species, and CD = 1 if the species occur in the same proportions in both samples.

Horn's introduced a modification of the index

C H = 2 i = 1 S x i y i ( i = 1 S x i 2 X 2 + i = 1 S y i 2 Y 2 ) X Y {\displaystyle C_{H}={\frac {2\sum _{i=1}^{S}x_{i}y_{i}}{\left({\sum _{i=1}^{S}x_{i}^{2} \over X^{2}}+{\sum _{i=1}^{S}y_{i}^{2} \over Y^{2}}\right)XY}}}

Standardised Morisita's index

Smith-Gill developed a statistic based on Morisita's index which is independent of both sample size and population density and bounded by −1 and +1. This statistic is calculated as follows

First determine Morisita's index ( Id ) in the usual fashion. Then let k be the number of units the population was sampled from. Calculate the two critical values

M u = χ 0.975 2 k + x x 1 {\displaystyle M_{u}={\frac {\chi _{0.975}^{2}-k+\sum x}{\sum x-1}}}
M c = χ 0.025 2 k + x x 1 {\displaystyle M_{c}={\frac {\chi _{0.025}^{2}-k+\sum x}{\sum x-1}}}

where χ is the chi square value for n − 1 degrees of freedom at the 97.5% and 2.5% levels of confidence.

The standardised index ( Ip ) is then calculated from one of the formulae below

When IdMc > 1

I p = 0.5 + 0.5 ( I d M c k M c ) {\displaystyle I_{p}=0.5+0.5\left({\frac {I_{d}-M_{c}}{k-M_{c}}}\right)}

When Mc > Id ≥ 1

I p = 0.5 ( I d 1 M u 1 ) {\displaystyle I_{p}=0.5\left({\frac {I_{d}-1}{M_{u}-1}}\right)}

When 1 > IdMu

I p = 0.5 ( I d 1 M u 1 ) {\displaystyle I_{p}=-0.5\left({\frac {I_{d}-1}{M_{u}-1}}\right)}

When 1 > Mu > Id

I p = 0.5 + 0.5 ( I d M u M u ) {\displaystyle I_{p}=-0.5+0.5\left({\frac {I_{d}-M_{u}}{M_{u}}}\right)}

Ip ranges between +1 and −1 with 95% confidence intervals of ±0.5. Ip has the value of 0 if the pattern is random; if the pattern is uniform, Ip < 0 and if the pattern shows aggregation, Ip > 0.

Peet's evenness indices

These indices are a measure of evenness between samples.

E 1 = I I min I max I min {\displaystyle E_{1}={\frac {I-I_{\min }}{I_{\max }-I_{\min }}}}
E 2 = I I max {\displaystyle E_{2}={\frac {I}{I_{\max }}}}

where I is an index of diversity, Imax and Imin are the maximum and minimum values of I between the samples being compared.

Loevinger's coefficient

Loevinger has suggested a coefficient H defined as follows:

H = p max ( 1 p min ) p min ( 1 p max ) {\displaystyle H={\sqrt {\frac {p_{\max }(1-p_{\min })}{p_{\min }(1-p_{\max })}}}}

where pmax and pmin are the maximum and minimum proportions in the sample.

Tversky index

The Tversky index is an asymmetric measure that lies between 0 and 1.

For samples A and B the Tversky index (S) is

S = | A B | | A B | + α | A B | + β | B A | {\displaystyle S={\frac {|A\cap B|}{|A\cap B|+\alpha |A-B|+\beta |B-A|}}}

The values of α and β are arbitrary. Setting both α and β to 0.5 gives Dice's coefficient. Setting both to 1 gives Tanimoto's coefficient.

A symmetrical variant of this index has also been proposed.

S 1 = | A B | | A B | + β ( α a + ( 1 α ) b ) {\displaystyle S_{1}={\frac {|A\cap B|}{|A\cap B|+\beta \left(\alpha a+(1-\alpha )b\right)}}}

where

a = min ( | X Y | , | Y X | ) {\displaystyle a=\min \left(|X-Y|,|Y-X|\right)}
b = max ( | X Y | , | Y X | ) {\displaystyle b=\max \left(|X-Y|,|Y-X|\right)}

Several similar indices have been proposed.

Monostori et al. proposed the SymmetricSimilarity index

S S ( A , B ) = | d ( A ) d ( B ) | | d ( A ) + d ( B ) | {\displaystyle SS(A,B)={\frac {|d(A)\cap d(B)|}{|d(A)+d(B)|}}}

where d(X) is some measure of derived from X.

Bernstein and Zobel have proposed the S2 and S3 indexes

S 2 = | d ( A ) d ( B ) | min ( | d ( A ) | , | d ( B ) ) | {\displaystyle S2={\frac {|d(A)\cap d(B)|}{\min(|d(A)|,|d(B))|}}}
S 3 = 2 | d ( A ) d ( B ) | | d ( A ) + d ( B ) | {\displaystyle S3={\frac {2|d(A)\cap d(B)|}{|d(A)+d(B)|}}}

S3 is simply twice the SymmetricSimilarity index. Both are related to Dice's coefficient

Metrics used

A number of metrics (distances between samples) have been proposed.

Euclidean distance

While this is usually used in quantitative work it may also be used in qualitative work. This is defined as

d j k = i = 1 N ( x i j x i k ) 2 {\displaystyle d_{jk}={\sqrt {\sum _{i=1}^{N}(x_{ij}-x_{ik})^{2}}}}

where djk is the distance between xij and xik.

Gower's distance

This is defined as

G D = Σ i = 1 n w i d i Σ i = 1 n w i {\displaystyle GD={\frac {\Sigma _{i=1}^{n}w_{i}d_{i}}{\Sigma _{i=1}^{n}w_{i}}}}

where di is the distance between the i samples and wi is the weighing give to the i distance.

Manhattan distance

While this is more commonly used in quantitative work it may also be used in qualitative work. This is defined as

d j k = i = 1 N | x i j x i k | {\displaystyle d_{jk}=\sum _{i=1}^{N}|x_{ij}-x_{ik}|}

where djk is the distance between xij and xik and || is the absolute value of the difference between xij and xik.

A modified version of the Manhattan distance can be used to find a zero (root) of a polynomial of any degree using Lill's method.

Prevosti's distance

This is related to the Manhattan distance. It was described by Prevosti et al. and was used to compare differences between chromosomes. Let P and Q be two collections of r finite probability distributions. Let these distributions have values that are divided into k categories. Then the distance DPQ is

D P Q = 1 r j = 1 r i = 1 k | p j i q j i | {\displaystyle D_{PQ}={\frac {1}{r}}\sum _{j=1}^{r}\sum _{i=1}^{k}|p_{ji}-q_{ji}|}

where r is the number of discrete probability distributions in each population, kj is the number of categories in distributions Pj and Qj and pji (respectively qji) is the theoretical probability of category i in distribution Pj (Qj) in population P(Q).

Its statistical properties were examined by Sanchez et al. who recommended a bootstrap procedure to estimate confidence intervals when testing for differences between samples.

Other metrics

Let

A = x i j {\displaystyle A=\sum x_{ij}}
B = x i k {\displaystyle B=\sum x_{ik}}
J = min ( x i j , x j k ) {\displaystyle J=\sum \min(x_{ij},x_{jk})}

where min(x,y) is the lesser value of the pair x and y.

Then

d j k = A + B 2 J {\displaystyle d_{jk}=A+B-2J}

is the Manhattan distance,

d j k = A + B 2 J A + B {\displaystyle d_{jk}={\frac {A+B-2J}{A+B}}}

is the Bray−Curtis distance,

d j k = A + B 2 J A + B J {\displaystyle d_{jk}={\frac {A+B-2J}{A+B-J}}}

is the Jaccard (or Ruzicka) distance and

d j k = 1 1 2 ( J A + J B ) {\displaystyle d_{jk}=1-{\frac {1}{2}}\left({\frac {J}{A}}+{\frac {J}{B}}\right)}

is the Kulczynski distance.

Similarities between texts

HaCohen-Kerner et al. have proposed a variety of metrics for comparing two or more texts.

Ordinal data

If the categories are at least ordinal then a number of other indices may be computed.

Leik's D

Leik's measure of dispersion (D) is one such index. Let there be K categories and let pi be fi/N where fi is the number in the i category and let the categories be arranged in ascending order. Let

c a = i = 1 a p i {\displaystyle c_{a}=\sum _{i=1}^{a}p_{i}}

where aK. Let da = ca if ca ≤ 0.5 and 1 − ca ≤ 0.5 otherwise. Then

D = 2 a = 1 K d a K 1 {\displaystyle D=2\sum _{a=1}^{K}{\frac {d_{a}}{K-1}}}

Normalised Herfindahl measure

This is the square of the coefficient of variation divided by N − 1 where N is the sample size.

H = 1 N 1 s 2 m 2 {\displaystyle H={\frac {1}{N-1}}{\frac {s^{2}}{m^{2}}}}

where m is the mean and s is the standard deviation.

Potential-for-conflict Index

The potential-for-conflict Index (PCI) describes the ratio of scoring on either side of a rating scale's centre point. This index requires at least ordinal data. This ratio is often displayed as a bubble graph.

The PCI uses an ordinal scale with an odd number of rating points (−n to +n) centred at 0. It is calculated as follows

P C I = X t Z [ 1 | i = 1 r + X + X t i = 1 r X X t | ] {\displaystyle PCI={\frac {X_{t}}{Z}}\left}

where Z = 2n, |·| is the absolute value (modulus), r+ is the number of responses in the positive side of the scale, r is the number of responses in the negative side of the scale, X+ are the responses on the positive side of the scale, X are the responses on the negative side of the scale and

X t = i = 1 r + | X + | + i = 1 r | X | {\displaystyle X_{t}=\sum _{i=1}^{r_{+}}|X_{+}|+\sum _{i=1}^{r_{-}}|X_{-}|}

Theoretical difficulties are known to exist with the PCI. The PCI can be computed only for scales with a neutral center point and an equal number of response options on either side of it. Also a uniform distribution of responses does not always yield the midpoint of the PCI statistic but rather varies with the number of possible responses or values in the scale. For example, five-, seven- and nine-point scales with a uniform distribution of responses give PCIs of 0.60, 0.57 and 0.50 respectively.

The first of these problems is relatively minor as most ordinal scales with an even number of response can be extended (or reduced) by a single value to give an odd number of possible responses. Scale can usually be recentred if this is required. The second problem is more difficult to resolve and may limit the PCI's applicability.

The PCI has been extended

P C I 2 = i = 1 K j = 1 i k i k j d i j δ {\displaystyle PCI_{2}={\frac {\sum _{i=1}^{K}\sum _{j=1}^{i}k_{i}k_{j}d_{ij}}{\delta }}}

where K is the number of categories, ki is the number in the i category, dij is the distance between the i and i categories, and δ is the maximum distance on the scale multiplied by the number of times it can occur in the sample. For a sample with an even number of data points

δ = N 2 2 d max {\displaystyle \delta ={\frac {N^{2}}{2}}d_{\max }}

and for a sample with an odd number of data points

δ = N 2 1 2 d max {\displaystyle \delta ={\frac {N^{2}-1}{2}}d_{\max }}

where N is the number of data points in the sample and dmax is the maximum distance between points on the scale.

Vaske et al. suggest a number of possible distance measures for use with this index.

D 1 : d i j = | r i r j | 1 {\displaystyle D_{1}:d_{ij}=|r_{i}-r_{j}|-1}

if the signs (+ or −) of ri and rj differ. If the signs are the same dij = 0.

D 2 : d i j = | r i r j | {\displaystyle D_{2}:d_{ij}=|r_{i}-r_{j}|}
D 3 : d i j = | r i r j | p {\displaystyle D_{3}:d_{ij}=|r_{i}-r_{j}|^{p}}

where p is an arbitrary real number > 0.

D p i j : d i j = [ | r i r j | ( m 1 ) ] p {\displaystyle Dp_{ij}:d_{ij}=^{p}}

if sign(ri ) ≠ sign(ri ) and p is a real number > 0. If the signs are the same then dij = 0. m is D1, D2 or D3.

The difference between D1 and D2 is that the first does not include neutrals in the distance while the latter does. For example, respondents scoring −2 and +1 would have a distance of 2 under D1 and 3 under D2.

The use of a power (p) in the distances allows for the rescaling of extreme responses. These differences can be highlighted with p > 1 or diminished with p < 1.

In simulations with a variates drawn from a uniform distribution the PCI2 has a symmetric unimodal distribution. The tails of its distribution are larger than those of a normal distribution.

Vaske et al. suggest the use of a t test to compare the values of the PCI between samples if the PCIs are approximately normally distributed.

van der Eijk's A

This measure is a weighted average of the degree of agreement the frequency distribution. A ranges from −1 (perfect bimodality) to +1 (perfect unimodality). It is defined as

A = U ( 1 S 1 K 1 ) {\displaystyle A=U\left(1-{\frac {S-1}{K-1}}\right)}

where U is the unimodality of the distribution, S the number of categories that have nonzero frequencies and K the total number of categories.

The value of U is 1 if the distribution has any of the three following characteristics:

  • all responses are in a single category
  • the responses are evenly distributed among all the categories
  • the responses are evenly distributed among two or more contiguous categories, with the other categories with zero responses

With distributions other than these the data must be divided into 'layers'. Within a layer the responses are either equal or zero. The categories do not have to be contiguous. A value for A for each layer (Ai) is calculated and a weighted average for the distribution is determined. The weights (wi) for each layer are the number of responses in that layer. In symbols

A o v e r a l l = w i A i {\displaystyle A_{\mathrm {overall} }=\sum w_{i}A_{i}}

A uniform distribution has A = 0: when all the responses fall into one category A = +1.

One theoretical problem with this index is that it assumes that the intervals are equally spaced. This may limit its applicability.

Related statistics

Birthday problem

If there are n units in the sample and they are randomly distributed into k categories (nk), this can be considered a variant of the birthday problem. The probability (p) of all the categories having only one unit is

p = i = 1 n ( 1 i k ) {\displaystyle p=\prod _{i=1}^{n}\left(1-{\frac {i}{k}}\right)}

If c is large and n is small compared with k then to a good approximation

p = exp ( n 2 2 k ) {\displaystyle p=\exp \left({\frac {-n^{2}}{2k}}\right)}

This approximation follows from the exact formula as follows:

log e ( 1 i k ) i k {\displaystyle \log _{e}\left(1-{\frac {i}{k}}\right)\approx -{\frac {i}{k}}}
Sample size estimates

For p = 0.5 and p = 0.05 respectively the following estimates of n may be useful

n = 1.2 k {\displaystyle n=1.2{\sqrt {k}}}
n = 2.448 k 2.5 k {\displaystyle n=2.448{\sqrt {k}}\approx 2.5{\sqrt {k}}}

This analysis can be extended to multiple categories. For p = 0.5 and p 0.05 we have respectively

n = 1.2 1 i = 1 k 1 c i {\displaystyle n=1.2{\sqrt {\frac {1}{\sum _{i=1}^{k}{\frac {1}{c_{i}}}}}}}
n 2.5 1 i = 1 k 1 c i {\displaystyle n\approx 2.5{\sqrt {\frac {1}{\sum _{i=1}^{k}{\frac {1}{c_{i}}}}}}}

where ci is the size of the i category. This analysis assumes that the categories are independent.

If the data is ordered in some fashion then for at least one event occurring in two categories lying within j categories of each other than a probability of 0.5 or 0.05 requires a sample size (n) respectively of

n = 1.2 k 2 j + 1 {\displaystyle n=1.2{\sqrt {\frac {k}{2j+1}}}}
n 2.5 k 2 j + 1 {\displaystyle n\approx 2.5{\sqrt {\frac {k}{2j+1}}}}

where k is the number of categories.

Birthday-death day problem

Whether or not there is a relation between birthdays and death days has been investigated with the statistic

log 10 ( 1 + 2 d 365 ) , {\displaystyle -\log _{10}\left({\frac {1+2d}{365}}\right),}

where d is the number of days in the year between the birthday and the death day.

Rand index

The Rand index is used to test whether two or more classification systems agree on a data set.

Given a set of n {\displaystyle n} elements S = { o 1 , , o n } {\displaystyle S=\{o_{1},\ldots ,o_{n}\}} and two partitions of S {\displaystyle S} to compare, X = { X 1 , , X r } {\displaystyle X=\{X_{1},\ldots ,X_{r}\}} , a partition of S into r subsets, and Y = { Y 1 , , Y s } {\displaystyle Y=\{Y_{1},\ldots ,Y_{s}\}} , a partition of S into s subsets, define the following:

  • a {\displaystyle a} , the number of pairs of elements in S {\displaystyle S} that are in the same subset in X {\displaystyle X} and in the same subset in Y {\displaystyle Y}
  • b {\displaystyle b} , the number of pairs of elements in S {\displaystyle S} that are in different subsets in X {\displaystyle X} and in different subsets in Y {\displaystyle Y}
  • c {\displaystyle c} , the number of pairs of elements in S {\displaystyle S} that are in the same subset in X {\displaystyle X} and in different subsets in Y {\displaystyle Y}
  • d {\displaystyle d} , the number of pairs of elements in S {\displaystyle S} that are in different subsets in X {\displaystyle X} and in the same subset in Y {\displaystyle Y}

The Rand index - R {\displaystyle R} - is defined as

R = a + b a + b + c + d = a + b ( n 2 ) {\displaystyle R={\frac {a+b}{a+b+c+d}}={\frac {a+b}{n \choose 2}}}

Intuitively, a + b {\displaystyle a+b} can be considered as the number of agreements between X {\displaystyle X} and Y {\displaystyle Y} and c + d {\displaystyle c+d} as the number of disagreements between X {\displaystyle X} and Y {\displaystyle Y} .

Adjusted Rand index

The adjusted Rand index is the corrected-for-chance version of the Rand index. Though the Rand Index may only yield a value between 0 and +1, the adjusted Rand index can yield negative values if the index is less than the expected index.

The contingency table

Given a set S {\displaystyle S} of n {\displaystyle n} elements, and two groupings or partitions (e.g. clusterings) of these points, namely X = { X 1 , X 2 , , X r } {\displaystyle X=\{X_{1},X_{2},\ldots ,X_{r}\}} and Y = { Y 1 , Y 2 , , Y s } {\displaystyle Y=\{Y_{1},Y_{2},\ldots ,Y_{s}\}} , the overlap between X {\displaystyle X} and Y {\displaystyle Y} can be summarized in a contingency table [ n i j ] {\displaystyle \left} where each entry n i j {\displaystyle n_{ij}} denotes the number of objects in common between X i {\displaystyle X_{i}} and Y j {\displaystyle Y_{j}}  : n i j = | X i Y j | {\displaystyle n_{ij}=|X_{i}\cap Y_{j}|} .

X\Y Y 1 {\displaystyle Y_{1}} Y 2 {\displaystyle Y_{2}} {\displaystyle \ldots } Y s {\displaystyle Y_{s}} Sums
X 1 {\displaystyle X_{1}} n 11 {\displaystyle n_{11}} n 12 {\displaystyle n_{12}} {\displaystyle \ldots } n 1 s {\displaystyle n_{1s}} a 1 {\displaystyle a_{1}}
X 2 {\displaystyle X_{2}} n 21 {\displaystyle n_{21}} n 22 {\displaystyle n_{22}} {\displaystyle \ldots } n 2 s {\displaystyle n_{2s}} a 2 {\displaystyle a_{2}}
{\displaystyle \vdots } {\displaystyle \vdots } {\displaystyle \vdots } {\displaystyle \ddots } {\displaystyle \vdots } {\displaystyle \vdots }
X r {\displaystyle X_{r}} n r 1 {\displaystyle n_{r1}} n r 2 {\displaystyle n_{r2}} {\displaystyle \ldots } n r s {\displaystyle n_{rs}} a r {\displaystyle a_{r}}
Sums b 1 {\displaystyle b_{1}} b 2 {\displaystyle b_{2}} {\displaystyle \ldots } b s {\displaystyle b_{s}}

Definition

The adjusted form of the Rand Index, the Adjusted Rand Index, is

AdjustedIndex = Index ExpectedIndex MaxIndex ExpectedIndex , {\displaystyle {\text{AdjustedIndex}}={\frac {{\text{Index}}-{\text{ExpectedIndex}}}{{\text{MaxIndex}}-{\text{ExpectedIndex}}}},}

more specifically

ARI = i j ( n i j 2 ) [ i ( a i 2 ) j ( b j 2 ) ] / ( n 2 ) 1 2 [ i ( a i 2 ) + j ( b j 2 ) ] [ i ( a i 2 ) j ( b j 2 ) ] / ( n 2 ) {\displaystyle {\text{ARI}}={\frac {\sum _{ij}{\binom {n_{ij}}{2}}-\left.\left\right/{\binom {n}{2}}}{{\frac {1}{2}}\left-\left.\left\right/{\binom {n}{2}}}}}

where n i j , a i , b j {\displaystyle n_{ij},a_{i},b_{j}} are values from the contingency table.

Since the denominator is the total number of pairs, the Rand index represents the frequency of occurrence of agreements over the total pairs, or the probability that X {\displaystyle X} and Y {\displaystyle Y} will agree on a randomly chosen pair.

Evaluation of indices

Different indices give different values of variation, and may be used for different purposes: several are used and critiqued in the sociology literature especially.

If one wishes to simply make ordinal comparisons between samples (is one sample more or less varied than another), the choice of IQV is relatively less important, as they will often give the same ordering.

Where the data is ordinal a method that may be of use in comparing samples is ORDANOVA.

In some cases it is useful to not standardize an index to run from 0 to 1, regardless of number of categories or samples (Wilcox 1973, pp. 338), but one generally so standardizes it.

See also

Notes

  1. This can only happen if the number of cases is a multiple of the number of categories.
  2. Freemen LC (1965) Elementary applied statistics. New York: John Wiley and Sons pp. 40–43
  3. Kendal MC, Stuart A (1958) The advanced theory of statistics. Hafner Publishing Company p. 46
  4. Mueller JE, Schuessler KP (1961) Statistical reasoning in sociology. Boston: Houghton Mifflin Company. pp. 177–179
  5. Wilcox (1967), p. .
  6. Kaiser HF (1968) "A measure of the population quality of legislative apportionment." The American Political Science Review 62 (1) 208
  7. Joel Gombin (August 18, 2015). "qualvar: Initial release (Version v0.1)". Zenodo. doi:10.5281/zenodo.28341.
  8. Gibbs & Poston Jr (1975).
  9. Lieberson (1969), p. 851.
  10. IQV at xycoon
  11. Hunter, PR; Gaston, MA (1988). "Numerical index of the discriminatory ability of typing systems: an application of Simpson's index of diversity". J Clin Microbiol. 26 (11): 2465–2466. doi:10.1128/jcm.26.11.2465-2466.1988. PMC 266921. PMID 3069867.
  12. Friedman WF (1925) The incidence of coincidence and its applications in cryptanalysis. Technical Paper. Office of the Chief Signal Officer. United States Government Printing Office.
  13. Gini CW (1912) Variability and mutability, contribution to the study of statistical distributions and relations. Studi Economico-Giuricici della R. Universita de Cagliari
  14. Simpson, EH (1949). "Measurement of diversity". Nature. 163 (4148): 688. Bibcode:1949Natur.163..688S. doi:10.1038/163688a0.
  15. Bachi R (1956) A statistical analysis of the revival of Hebrew in Israel. In: Bachi R (ed) Scripta Hierosolymitana, Vol III, Jerusalem: Magnus press pp 179–247
  16. Mueller JH, Schuessler KF (1961) Statistical reasoning in sociology. Boston: Houghton Mifflin
  17. Gibbs, JP; Martin, WT (1962). "Urbanization, technology and division of labor: International patterns". American Sociological Review. 27 (5): 667–677. doi:10.2307/2089624. JSTOR 2089624.
  18. Lieberson (1969), p. .
  19. Blau P (1977) Inequality and Heterogeneity. Free Press, New York
  20. Perry M, Kader G (2005) Variation as unalikeability. Teaching Stats 27 (2) 58–60
  21. Greenberg, JH (1956). "The measurement of linguistic diversity". Language. 32 (1): 109–115. doi:10.2307/410659. JSTOR 410659.
  22. Lautard EH (1978) PhD thesis.
  23. Berger, WH; Parker, FL (1970). "Diversity of planktonic Foramenifera in deep sea sediments". Science. 168 (3937): 1345–1347. Bibcode:1970Sci...168.1345B. doi:10.1126/science.168.3937.1345. PMID 17731043. S2CID 29553922.
  24. ^ Hill, M O (1973). "Diversity and evenness: a unifying notation and its consequences". Ecology. 54 (2): 427–431. Bibcode:1973Ecol...54..427H. doi:10.2307/1934352. JSTOR 1934352.
  25. Margalef R (1958) Temporal succession and spatial heterogeneity in phytoplankton. In: Perspectives in marine biology. Buzzati-Traverso (ed) Univ Calif Press, Berkeley pp 323–347
  26. Menhinick, EF (1964). "A comparison of some species-individuals diversity indices applied to samples of field insects". Ecology. 45 (4): 859–861. Bibcode:1964Ecol...45..859M. doi:10.2307/1934933. JSTOR 1934933.
  27. Kuraszkiewicz W (1951) Nakladen Wroclawskiego Towarzystwa Naukowego
  28. Guiraud P (1954) Les caractères statistiques du vocabulaire. Presses Universitaires de France, Paris
  29. Panas E (2001) The Generalized Torquist: Specification and estimation of a new vocabulary-text size function. J Quant Ling 8(3) 233–252
  30. Kempton, RA; Taylor, LR (1976). "Models and statistics for species diversity". Nature. 262 (5571): 818–820. Bibcode:1976Natur.262..818K. doi:10.1038/262818a0. PMID 958461. S2CID 4168222.
  31. Hutcheson K (1970) A test for comparing diversities based on the Shannon formula. J Theo Biol 29: 151–154
  32. McIntosh RP (1967). An Index of Diversity and the Relation of Certain Concepts to Diversity. Ecology, 48(3), 392–404
  33. Fisher RA, Corbet A, Williams CB (1943) The relation between the number of species and the number of individuals in a random sample of an animal population. Animal Ecol 12: 42–58
  34. Anscombe (1950) Sampling theory of the negative binomial and logarithmic series distributions. Biometrika 37: 358–382
  35. Strong, WL (2002). "Assessing species abundance uneveness within and between plant communities" (PDF). Community Ecology. 3 (2): 237–246. doi:10.1556/comec.3.2002.2.9.
  36. Camargo JA (1993) Must dominance increase with the number of subordinate species in competitive interactions? J. Theor Biol 161 537–542
  37. Smith, Wilson (1996)
  38. Bulla, L (1994). "An index of evenness and its associated diversity measure". Oikos. 70 (1): 167–171. Bibcode:1994Oikos..70..167B. doi:10.2307/3545713. JSTOR 3545713.
  39. Horn, HS (1966). "Measurement of 'overlap' in comparative ecological studies". Am Nat. 100 (914): 419–423. doi:10.1086/282436. S2CID 84469180.
  40. Siegel, Andrew F (2006) "Rarefaction curves." Encyclopedia of Statistical Sciences 10.1002/0471667196.ess2195.pub2.
  41. Caswell H (1976) Community structure: a neutral model analysis. Ecol Monogr 46: 327–354
  42. Poulin, R; Mouillot, D (2003). "Parasite specialization from a phylogenetic perspective: a new index of host specificity". Parasitology. 126 (5): 473–480. CiteSeerX 10.1.1.574.7432. doi:10.1017/s0031182003002993. PMID 12793652. S2CID 9440341.
  43. Theil H (1972) Statistical decomposition analysis. Amsterdam: North-Holland Publishing Company>
  44. Duncan OD, Duncan B (1955) A methodological analysis of segregation indexes. Am Sociol Review, 20: 210–217
  45. Gorard S, Taylor C (2002b) What is segregation? A comparison of measures in terms of 'strong' and 'weak' compositional invariance. Sociology, 36(4), 875–895
  46. Massey, DS; Denton, NA (1988). "The dimensions of residential segregation". Social Forces. 67 (2): 281–315. doi:10.1093/sf/67.2.281.
  47. Hutchens RM (2004) One measure of segregation. International Economic Review 45: 555–578
  48. Lieberson S (1981). "An asymmetrical approach to segregation". In Peach C, Robinson V, Smith S (eds.). Ethnic segregation in cities. London: Croom Helm. pp. 61–82.
  49. Bell, W (1954). "A probability model for the measurement of ecological segregation". Social Forces. 32 (4): 357–364. doi:10.2307/2574118. JSTOR 2574118.
  50. Ochiai A (1957) Zoogeographic studies on the soleoid fishes found in Japan and its neighbouring regions. Bull Jpn Soc Sci Fish 22: 526–530
  51. Kulczynski S (1927) Die Pflanzenassoziationen der Pieninen. Bulletin International de l'Académie Polonaise des Sciences et des Lettres, Classe des Sciences
  52. Yule GU (1900) On the association of attributes in statistics. Philos Trans Roy Soc
  53. Lienert GA and Sporer SL (1982) Interkorrelationen seltner Symptome mittels Nullfeldkorrigierter YuleKoeffizienten. Psychologische Beitrage 24: 411–418
  54. Baroni-Urbani, C; Buser, MW (1976). "similarity of binary Data". Systematic Biology. 25 (3): 251–259. doi:10.2307/2412493. JSTOR 2412493.
  55. Forbes SA (1907) On the local distribution of certain Illinois fishes: an essay in statistical ecology. Bulletin of the Illinois State Laboratory of Natural History 7:272–303
  56. Alroy J (2015) A new twist on a very old binary similarity coefficient. Ecology 96 (2) 575-586
  57. Carl R. Hausman and Douglas R. Anderson (2012). Conversations on Peirce: Reals and Ideals. Fordham University Press. p. 221. ISBN 9780823234677.
  58. Lance, G. N.; Williams, W. T. (1966). "Computer programs for hierarchical polythetic classification ("similarity analysis")". Computer Journal. 9 (1): 60–64. doi:10.1093/comjnl/9.1.60.
  59. Lance, G. N.; Williams, W. T. (1967). "Mixed-data classificatory programs I.) Agglomerative Systems". Australian Computer Journal: 15–20.
  60. Jaccard P (1902) Lois de distribution florale. Bulletin de la Socíeté Vaudoise des Sciences Naturelles 38:67-130
  61. Archer AW and Maples CG (1989) Response of selected binomial coefficients to varying degrees of matrix sparseness and to matrices with known data interrelationships. Mathematical Geology 21: 741–753
  62. ^ Morisita M (1959) Measuring the dispersion and the analysis of distribution patterns. Memoirs of the Faculty of Science, Kyushu University Series E. Biol 2:215–235
  63. Lloyd M (1967) Mean crowding. J Anim Ecol 36: 1–30
  64. Pedigo LP & Buntin GD (1994) Handbook of sampling methods for arthropods in agriculture. CRC Boca Raton FL
  65. Morisita M (1959) Measuring of the dispersion and analysis of distribution patterns. Memoirs of the Faculty of Science, Kyushu University, Series E Biology. 2: 215–235
  66. Horn, HS (1966). "Measurement of "Overlap" in comparative ecological studies". The American Naturalist. 100 (914): 419–424. doi:10.1086/282436. S2CID 84469180.
  67. Smith-Gill SJ (1975). "Cytophysiological basis of disruptive pigmentary patterns in the leopard frog Rana pipiens. II. Wild type and mutant cell specific patterns". J Morphol. 146 (1): 35–54. doi:10.1002/jmor.1051460103. PMID 1080207. S2CID 23780609.
  68. Peet (1974) The measurements of species diversity. Annu Rev Ecol Syst 5: 285–307
  69. Tversky, Amos (1977). "Features of Similarity" (PDF). Psychological Review. 84 (4): 327–352. doi:10.1037/0033-295x.84.4.327.
  70. Jimenez S, Becerra C, Gelbukh A SOFTCARDINALITY-CORE: Improving text overlap with distributional measures for semantic textual similarity. Second Joint Conference on Lexical and Computational Semantics (*SEM), Volume 1: Proceedings of the main conference and the shared task: semantic textual similarity, p194-201. June 7–8, 2013, Atlanta, Georgia, USA
  71. Monostori K, Finkel R, Zaslavsky A, Hodasz G and Patke M (2002) Comparison of overlap detection techniques. In: Proceedings of the 2002 International Conference on Computational Science. Lecture Notes in Computer Science 2329: 51-60
  72. Bernstein Y and Zobel J (2004) A scalable system for identifying co-derivative documents. In: Proceedings of 11th International Conference on String Processing and Information Retrieval (SPIRE) 3246: 55-67
  73. Prevosti, A; Ribo, G; Serra, L; Aguade, M; Balanya, J; Monclus, M; Mestres, F (1988). "Colonization of America by Drosophila subobscura: experiment in natural populations that supports the adaptive role of chromosomal inversion polymorphism". Proc Natl Acad Sci USA. 85 (15): 5597–5600. Bibcode:1988PNAS...85.5597P. doi:10.1073/pnas.85.15.5597. PMC 281806. PMID 16593967.
  74. Sanchez, A; Ocana, J; Utzetb, F; Serrac, L (2003). "Comparison of Prevosti genetic distances". Journal of Statistical Planning and Inference. 109 (1–2): 43–65. doi:10.1016/s0378-3758(02)00297-5.
  75. HaCohen-Kerner Y, Tayeb A and Ben-Dror N (2010) Detection of simple plagiarism in computer science papers. In: Proceedings of the 23rd International Conference on Computational Linguistics pp 421-429
  76. Leik R (1966) A measure of ordinal consensus. Pacific sociological review 9 (2): 85–90
  77. Manfredo M, Vaske, JJ, Teel TL (2003) The potential for conflict index: A graphic approach tp practical significance of human dimensions research. Human Dimensions of Wildlife 8: 219–228
  78. ^ Vaske JJ, Beaman J, Barreto H, Shelby LB (2010) An extension and further validation of the potential for conflict index. Leisure Sciences 32: 240–254
  79. Van der Eijk C (2001) Measuring agreement in ordered rating scales. Quality and quantity 35(3): 325–341
  80. Von Mises R (1939) Uber Aufteilungs-und Besetzungs-Wahrcheinlichkeiten. Revue de la Facultd des Sciences de de I'Universite d'lstanbul NS 4: 145−163
  81. Sevast'yanov BA (1972) Poisson limit law for a scheme of sums of dependent random variables. (trans. S. M. Rudolfer) Theory of probability and its applications, 17: 695−699
  82. Hoaglin DC, Mosteller, F and Tukey, JW (1985) Exploring data tables, trends, and shapes, New York: John Wiley
  83. ^ W. M. Rand (1971). "Objective criteria for the evaluation of clustering methods". Journal of the American Statistical Association. 66 (336): 846–850. arXiv:1704.01036. doi:10.2307/2284239. JSTOR 2284239.
  84. Lawrence Hubert and Phipps Arabie (1985). "Comparing partitions". Journal of Classification. 2 (1): 193–218. doi:10.1007/BF01908075. S2CID 189915041.
  85. Nguyen Xuan Vinh, Julien Epps and James Bailey (2009). "Information Theoretic Measures for Clustering Comparison: Is a Correction for Chance Necessary?" (PDF). ICML '09: Proceedings of the 26th Annual International Conference on Machine Learning. ACM. pp. 1073–1080. Archived from the original (PDF) on 25 March 2012.PDF.
  86. Wagner, Silke; Wagner, Dorothea (12 January 2007). "Comparing Clusterings - An Overview" (PDF). Retrieved 14 February 2018.

References

  • Lieberson, Stanley (December 1969), "Measuring Population Diversity", American Sociological Review, 34 (6): 850–862, doi:10.2307/2095977, JSTOR 2095977
  • Swanson, David A. (September 1976), "A Sampling Distribution and Significance Test for Differences in Qualitative Variation", Social Forces, 55 (1): 182–184, doi:10.2307/2577102, JSTOR 2577102
  • Wilcox, Allen R. (June 1973). "Indices of Qualitative Variation and Political Measurement". The Western Political Quarterly. 26 (2): 325–343. doi:10.2307/446831. JSTOR 446831.
Categories: