Taylor's law (also known as Taylor’s power law) is an empirical law in ecology that relates the variance of the number of individuals of a species per unit area of habitat to the corresponding mean by a power law relationship. It is named after the ecologist who first proposed it in 1961, Lionel Roy Taylor (1924–2007). Taylor's original name for this relationship was the law of the mean.
- 1 Definition
- 2 History
- 3 The Tweedie Hypothesis
- 4 Mathematical formulation
- 5 Applications
- 6 Related analyses
- 6.1 Barlett-Iawo model
- 6.2 Nachman model
- 6.3 Kono-Sugino equation
- 6.4 Hughes-Madden equation
- 6.5 Negative binomial distribution model
- 6.6 Tests for a common dispersion parameter
- 6.7 Charlier coefficient
- 6.8 Cole's index of dispersion
- 6.9 Lloyd's indexes
- 6.10 Patchiness regression test
- 6.11 Morisita’s index of dispersion
- 6.12 Standardised Morisita’s index
- 6.13 Southwood's index of spatial aggregation
- 6.14 Fisher's index of dispersion
- 6.15 Index of cluster size
- 6.16 Green’s index
- 6.17 Binary dispersal index
- 6.18 Ma's population aggregation critical density
- 7 Related statistics
- 8 See also
- 9 References
This law was originally defined for ecological systems, specifically to assess the spatial clustering of organisms. For a population count Y with mean µ and variance var(Y), Taylor’s law is written,
where a and b are both positive constants. Taylor proposed this relationship in 1961, suggesting that the exponent b be considered a species specific index of aggregation. This power law has subsequently been confirmed for many hundreds of species.
Taylor’s law has also been applied to assess the time dependent changes of population distributions. Related variance to mean power laws have also been demonstrated in several non-ecological systems:
- cancer metastasis
- the numbers of houses built over the Tonami plain in Japan.
- measles epidemiology
- HIV epidemiology,
- the geographic clustering of childhood leukemia
- blood flow heterogeneity
- the genomic distributions of single-nucleotide polymorphisms (SNPs)
- gene structures
- in number theory with sequential values of the Mertens function and also with the distribution of prime numbers  
- from the eigenvalue deviations of Gaussian orthogonal and unitary ensembles of random matrix theory
Fisher in 1921 proposed the equation
Smith in 1938 while studying crop yields proposed a relationship similar to Taylor's. This relationship was
where Vx is the variance of yield for plots of x units, V1 is the variance of yield per unit area and x is the size of plots. The slope (b) is the index of heterogeneity. The value of b in this relationship lies between 0 and 1. Where the yield are highly correlated b tends to 0; when they are uncorrelated b tends to 1.
L.R. Taylor (1924–2007) was an English entomologist who worked on the Rothamsted Insect Survey for pest control. His 1961 paper used data from 24 papers published between 1936 and 1960. These papers considered a variety of biological settings: virus lesions, macro-zooplankton, worms and symphylids in soil, insects in soil, on plants and in the air, mites on leaves, ticks on sheep and fish in the sea. In these papers the b value lay between 1 and 3. Taylor proposed the power law as a general feature of the spatial distribution of these species. He also proposed a mechanistic hypothesis to explain this law. Among the papers cited were those of Bliss and Yates and Finney.
Initial attempts to explain the spatial distribution of animals had been based on approaches like Bartlett’s stochastic population models and the negative binomial distribution that could result from birth-death processes. Taylor’s novel explanation was based the assumption of a balanced migratory and congregatory behavior of animals. His hypothesis was initially qualitative, but as it evolved it became semi-quantitative and was supported by simulations. In proposing that animal behavior was the principal mechanism behind the clustering of organisms, Taylor though appeared to have ignored his own report of clustering seen with tobacco necrosis virus plaques.
Following Taylor’s initial publications several alternative hypotheses for the power law were advanced. Hanski proposed a random walk model, modulated by the presumed multiplicative effect of reproduction. Hanski’s model predicted that the power law exponent would be constrained to range closely about the value of 2, which seemed inconsistent with many reported values.
Anderson et al formulated a simple stochastic birth, death, immigration and emigration model that yielded a quadratic variance function. As a response to this model Taylor argued that such a Markov process would predict that the power law exponent would vary considerably between replicate observations, and that such variability had not been observed.
About this time concerns were, however, raised regarding the statistical variability with measurements of the power law exponent, and the possibility that observations of a power law might reflect more mathematical artefact than a mechanistic process. Taylor et al responded with an additional publication of extensive observations which he claimed refuted Downing’s concerns.
In addition, Thórarinsson published a detailed critique of the animal behavioral model, noting that Taylor had modified his model several times in response to concerns raised, and that some of these modifications were inconsistent with earlier versions. Thórarinsson also claimed that Taylor confounded animal numbers with density and that Taylor had incorrectly interpreted simulations that had been constructed to demonstrate his models as validation.
Kemp reviewed a number of discrete stochastic models based on the negative binomial, Neyman type A, and Polya-Aeppli distributions that with suitable adjustment of parameters could produce a variance to mean power law. Kemp, however, did not explain the parameterizations of his models in mechanistic terms. Other relatively abstract models for Taylor’s law followed.
A number of additional statistical concerns were raised regarding Taylor’s law, based on the difficulty with real data in distinguishing between Taylor’s law and other variance to mean functions, as well the inaccuracy of standard regression methods.
Reports also began to accumulate where Taylor’s law had been applied to time series data. Perry showed how simulations based on chaos theory could yield Taylor’s law, and Kilpatrick & Ives provided simulations which showed how interactions between different species might lead to Taylor's law.
Other reports appeared where Taylor’s law had been applied to the spatial distribution of plants and bacterial populations As with the observations of Tobacco necrosis virus mentioned earlier, these observations were not consistent with Taylor’s animal behavioral model.
Earlier it was mentioned that variance to mean power function had been applied to non-ecological systems, under the rubric of Taylor’s law. To provide a more general explanation for the range of manifestations of the power law a hypothesis was proposed based on the Tweedie distributions, a family of probabilistic models that express an inherent power function relationship between the variance and the mean. Details regarding this hypothesis will be provided in the next section.
A further alternative explanation for Taylor's law was proposed by Cohen et al, derived from the Lewontin Cohen growth model. This model was successfully used to describe the spatial and temporal variability of forest populations.
Another paper by Cohen and Xu that random sampling in blocks where the underling distribution is skewed with the first four moments finite gives rise to Taylor's law. Approximate formulae for the parameters and their variances were also derived. These estimates were tested again data from the Black Rock Forest and found to be in reasonable agreement.
In the physics literature Taylor’s law has been referred to as ﬂuctuation scaling. Eisler et al, in a further attempt to find a general explanation for fluctuation scaling, proposed a process they called impact inhomogeneity in which frequent events are associated with larger impacts. In appendix B of the Eisler article, however, the authors noted that the equations for impact inhomogeneity yielded the same mathematical relationships as found with the Tweedie distributions.
Another group of physicists, Fronczak and Fronczak, derived Taylor’s power law for fluctuation scaling from principles of equilibrium and non-equilibrium statistical physics. Their derivation was based on assumptions of physical quantities like free energy and an external field that caused the clustering of biological organisms. Direct experimental demonstration of these postulated physical quantities in relationship to animal or plant aggregation has yet to be achieved, though. Shortly thereafter, an analysis of Fronczak and Fronczak’s model was presented that showed their equations directly lead to the Tweedie distributions, a finding that suggested that Fronczak and Fronczak had possibly provided a maximum entropy derivation of these distributions.
Taylor's law has been shown to hold for prime numbers not exceeding a given real number. This result has been shown to hold for the first 11 million primes. If the Hardy-Littlewood twin primes conjecture is true then this law also holds for twin primes.
Naming of law
The law itself is named after the ecologist Lionel Roy Taylor (1924–2007). The name 'Taylor's law' was coined by Southwood in 1966. Taylor's original name for this relationship was the law of the mean.
The Tweedie Hypothesis
About the time that Taylor was substantiating his ecological observations, MCK Tweedie, a British statistician and medical physicist, was investigating a family of probabilistic models that are now known as the Tweedie distributions. As mentioned above, these distributions are all characterized by a variance to mean power law mathematically identical to Taylor’s law.
The Tweedie distribution most applicable to ecological observations is the compound Poisson-gamma distribution, which represents the sum of N independent and identically distributed random variables with a gamma distribution where N is a random variable distributed in accordance with a Poisson distribution. In the additive form its cumulant generating function (CGF) is:
where κb(θ) is the cumulant function,
the Tweedie exponent
s is the generating function variable, and θ and λ are the canonical and index parameters, respectively.
These last two parameters are analogous to the scale and shape parameters used in probability theory. The cumulants of this distribution can be determined by successive differentiations of the CGF and then substituting s=0 into the resultant equations. The first and second cumulants are the mean and variance, respectively, and thus the compound Poisson-gamma CGF yields Taylor’s law with the proportionality constant
The compound Poisson-gamma cumulative distribution function has been verified for limited ecological data through the comparison of the theoretical distribution function with the empirical distribution function. A number of other systems, demonstrating variance to mean power laws related to Taylor’s law, have been similarly tested for the compound Poisson-gamma distribution.
The main justification for the Tweedie hypothesis rests with the mathematical convergence properties of the Tweedie distributions. The Tweedie convergence theorem requires the Tweedie distributions to act as foci of convergence for a wide range of statistical processes. As a consequence of this convergence theorem, processes based on the sum of multiple independent small jumps will tend to express Taylor’s law and obey a Tweedie distribution. A limit theorem for independent and identically distributed variables, as with the Tweedie convergence theorem, might then be considered as being fundamental relative to the ad hoc population models, or models proposed on the basis of simulation or approximation.
This hypothesis remains controversial; more conventional population dynamic approaches seem preferred amongst ecologists, despite the fact that the Tweedie compound Poisson distribution can be directly applied to population dynamic mechanisms.
One difficulty with the Tweedie hypothesis is that the value of b does not range between 0 and 1. Values of b < 1 are rare but have been reported.
In logarithmic form
Taylor's law is scale invariant. If the unit of measurement is changed by a constant factor c, the exponent (b) remains unchanged.
To see this let y = cx. Then
Taylor's law expressed in the original variable (x) is
and in the rescaled variable (y) it is
It has been shown that Taylor's law is the only relationship between the mean and variance that is scale invariant.
Extensions and refinements
A refinement in the estimation of the slope b has been proposed by Rayner.
where r is the Pearson moment correlation coefficient between log(s2) and log m, f is the ratio of sample variances in log(s2) and log m and φ is the ratio of the errors in log(s2) and log m.
Ordinary least squares regression assumes that φ = ∞. This tends to underestimate the value of b because the estimates of both log(s2) and log m are subject to error.
An extension of Taylor's law has been proposed by Ferris et al when multiple samples are taken
where s2 and m are the variance and mean respectively, b, c and d are constants and n is the number of samples taken. To date, this proposed extension has not been verified to be as applicable as the original version of Taylor's law.
An extension to this law for small samples has been proposed by Hanski. For small samples the Poisson variation (P) - the variation that can be ascribed to sampling variation - may be significant. Let S be the total variance and let V be the biological (real) variance. Then
Assuming the validity of Taylor's law, we have
Because in the Poisson distribution the mean equals the variance, we have
This gives us
This closely resembles Barlett's original suggestion.
Slope values (b) significantly > 1 indicate clumping of the organisms.
For populations that are experiencing constant per capita environmental variability, the regression of log( variance ) versus log( mean abundance ) should have a line with b = 2.
Most populations that have been studied have b < 2 (usually 1.5–1.6) but values of 2 have been reported. Occasionally cases with b > 2 have been reported. b values below 1 are uncommon but have also been reported ( b = 0.93 ).
The origin of the slope (b) in this regression remains unclear. Two hypotheses have been proposed to explain it. One suggests that b arises from the species behavior and is a constant for that species. The alternative suggests that it is dependent on the sampled population. Despite the considerable number of studies carried out on this law (over 1000), this question remains open.
It is known that both a and b are subject to change due to age-specific dispersal, mortality and sample unit size.
Extension to cluster sampling of binary data
A form of Taylor's law applicable to binary data in clusters (e.q., quadrats) has been proposed. In a binomial distribution, the theoretical variance is
where (varbin) is the binomial variance, n is the sample size per cluster, and p is the proportion of individuals with a trait (such as disease), an estimate of the probability of an individual having that trait.
One difficulty with binary data is that the mean and variance, in general, have a particular relationship: as the mean proportion of individuals infected increases above 0.5, the variance deceases.
It is now known that the observed variance (varobs) changes as a power function of (varbin).
Hughes and Madden noted that if the distribution is Poisson, the mean and variance are equal. As this is clearly not the case in many observed proportion samples, they instead assumed a binomial distribution. They replaced the mean in Taylor's law with the binomial variance and then compared this theoretical variance with the observed variance. For binomial data, they showed that varobs = varbin with overdispersion, varobs > varbin.
In symbols, Hughes and Madden's modification to Tyalor's law was
In logarithmic form this relationship is
This latter version is known as the binary power law.
A key step in the derivation of the binary power law by Hughes and Madden was the observation made by Patil and Stiteler that the variance-to-mean ratio used for assessing over-dispersion of unbounded counts in a single sample is actually the ratio of two variances: the observed variance and the theoretical variance for a random distribution. For unbounded counts, the random distribution is the Poisson. Thus, the Taylor power law for a collection of samples can be considered as a relationship between the observed variance and the Poisson variance.
More broadly, Madden and Hughes considered the power law as the relationship between two variances, the observed variance and the theoretical variance for a random distribution. With binary data, the random distribution is the binomial (not the Poisson). Thus the Taylor power law and the binary power law are two special cases of a general power-law relationships for heterogeneity.
When both a and b are equal to 1, then a small-scale random spatial pattern is suggested and is best described by the binomial distribution. When b = 1 and a > 1, there is over-dispersion (small scale aggregation). When b is > 1, the degree of aggregation varies with p. Turechek et al have showed that the binary power law describes numerous data sets in plant pathology. In general, b is greater than 1 and less than 2.
The fit of this law has been tested by simulations. These results suggest that rather than a single regression line for the data set, a segmental regression may be a better model for genuinely random distributions. However, this segmentation only occurs for very short-range dispersal distances and large quadrat sizes. The break in the line occurs only at p very close to 0.
An extension to this law has been proposed. The original form of this law is symmetrical but it can be extended to an asymmetrical form. Using simulations the symmetrical form fits the data when there is positive correlation of disease status of neighbors. Where there is a negative correlation between the likelihood of neighbours being infected, the asymmetrical version is a better fit to the data.
Because of the ubiquitous occurrence of Taylor's law in biology it has found a variety of uses some of which are listed here.
Recommendations as to use
It has been recommended based on simulation studies in applications testing the validity of Taylor's law to a data sample that:
(1) the total number of organisms studied be > 15
(2) the minimum number of groups of organisms studied be > 5
(3) the density of the organisms should vary by at least 2 orders of magnitude within the sample
Randomly distributed populations
It is common assumed (at least initially) that a population is randomly distributed in the environment. If a population is randomly distributed then the mean ( m ) and variance ( s2 ) of the population are equal and the proportion of samples that contain at least one individual ( p ) is
When a species with a clumped pattern is compared with one that is randomly distributed with equal overall densities, p will be less for the species having the clumped distribution pattern. Conversely when comparing a uniformly and a randomly distributed species but at equal overall densities, p will be greater for the randomly distributed population. This can be graphically tested by plotting p against m.
Wilson and Room developed a binomial model that incorporates Taylor's law. The basic relationship is
where the log is taken to the base e.
Incorporating Taylor's law this relationship becomes
Dispersion parameter estimator
The common dispersion parameter (k) of the negative binomial distribution is
where m is the sample mean and s2 is the variance. If 1 / k is > 0 the population is considered to be aggregated; 1 / k = 0 ( s2 = m ) the population is considered to be randomly (Poisson) distributed and if 1 / k is < 0 the population is considered to be uniformly distributed. No comment on the distribution can be made if k = 0.
Wilson and Room assuming that Taylor's law applied to the population gave an alternative estimator for k:
where a and b are the constants from Taylor's law.
derived an estimator for the probability of a sample containing x individuals per sampling unit. Jones's formula is
where P( x ) is the probability of finding x individuals per sampling unit, k is estimated from the Wilon and Room equation and m is the sample mean. The probability of finding zero individuals P( 0 ) is estimated with the negative binomial distribution
Jones also gives confidence intervals for these probabilities.
where CI is the confidence interval, t is the critical value taken from the t distribution and N is the total sample size.
Katz family of distributions
Katz proposed a family of distributions (the Katz family) with 2 parameters ( w1, w2 ). This family of distributions includes the Bernoulli, Geometric, Pascal and Poisson distributions as special cases. The mean and variance of a Katz distribution are
where m is the mean and s2 is the variance of the sample. The parameters can be estimated by the method of moments from which we have
For a Poisson distribution w2 = 0 and w1 = λ the parameter of the Possion distribution. This family of distributions is also sometimes known as the Panjer family of distributions.
The Katz family is related to the Sundt-Jewel family of distributions:
If the population obeys a Katz distribution then the coefficients of Taylor's law are
Katz also introduced a statistical test
where Jn is the test statistic, s2 is the variance of the sample, m is the mean of the sample and n is the sample size. Jn is asymptotically normally distributed with a zero mean and unit variance. If the sample is Poisson distributed Jn = 0; values of Jn < 0 and > 0 indicate under and over dispersion respectively. Overdispersion is often caused by latent heterogeneity - the presence of multiple sub populations within the population the sample is drawn from.
This statistic is related to the Neyman-Scott statistic
which is known to be asymptotically normal and the conditional chi-squared statistic (Poisson dispersion test)
which is known to have an asymptotic chi squared distribution with n − 1 degrees of freedom when the population is Poisson distributed.
If the population obeys Taylor's law then
Time to extinction
If Taylor's law is assumed to apply it is possible to determine the mean time to local extinction. This model assumes a simple random walk in time and the absence of density dependent population regulation.
Let where Nt+1 and Nt are the population sizes at time t + 1 and t respectively and r is parameter equal to the annual increase (decrease in population). Then
where var( r ) is the variance of r.
Let K be a measure of the species abundance (organisms per unit area). Then
where TE is the mean time to local extinction.
The probability of extinction by time t is
Minimum population size required to avoid extinction
Given that H must be > 0 for the population to persist then rearranging we have
is the minimum size of population for the species to persist.
The assumption of a lognormal distribution appears to apply to about half of a sample of 544 species. suggesting that it is at least a plausible assumption.
Sampling size estimators
The degree of precision (D) is defined to be s / m where s is the standard deviation and m is the mean. The degree of precision is known as the coefficient of variation in other contexts. In ecology research it is recommended that D be in the range 10-25%. The desired degree of precision is important in estimating the required sample size where an investigator wishes to test if Taylor's law applies to the data. The required sample size has been estimated for a number of simple distributions but where the population distribution is not known or cannot be assumed more complex formulae may needed to determine the required sample size.
Where the population is Poisson distributed the sample size (n) needed is
If the population is distributed as a negative binomial distribution then the required sample size is
where k is the parameter of the negative binomial distribution.
A more general sample size estimator has also been proposed
where a and b are derived from Taylor's law.
An alternative has been proposed by Southwood
where n is the required sample size, a and b are the Taylor's law coefficients and D is the desired degree of precision.
where d is the ratio of half the desired confidence interval (CI) to the mean. In symbols
The second estimator is used in binomial (presence-absence) sampling. The desired sample size (n) is
where the dp is ratio of half the desired confidence interval to the proportion of sample units with individuals, p is proportion of samples containing individuals and q = 1 - p. In symbols
For binary (presence/absence) sampling, Schulthess et al modified Karandinos' equation
where N is the required sample size, p is the proportion of units containing the organisms of interest, t is the chosen level of significance and Dip is a parameter derived from Taylor's law.
Sequential analysis is a method of statistical analysis where the sample size is not fixed in advance. Instead samples are taken in accordance with a predefined stopping rule. Taylor's law has been used to derive a number of stopping rules.
A formula for fixed precision in serial sampling to test Taylor's law was derived by Green in 1970.
where T is the cumulative sample total, D is the level of precision, n is the sample size and a and b are obtained from Taylor's law.
As an aid to pest control Wilson et al developed a test that incorporated a threshold level where action should be taken. The required sample size is
where a and b are the Taylor coefficients, || is the absolute value, m is the sample mean, T is the threshold level and t is the critical level of the t distribution. The authors also provided a similar test for binomial (presence-absence) sampling
where p is the probability of finding a sample with pests present and q = 1 - p.
Green derived another sampling formula for sequential sampling based on Taylor's law
where D is the degree of precision, a and b are the Taylor's law coefficients, n is the sample size and T is the total number of individuals sampled.
Serra et al have proposed a stopping rule based on Taylor's law.
where a and b are the parameters from Taylor's law, D is the desired level of precision and Tn is the total sample size.
Serra et al also proposed a second stopping rule based on Iwoa's regression
where α and β are the parameters of the regression line, D is the desired level of precision and Tn is the total sample size.
The authors recommended that D be set at 0.1 for studies of population dynamics and D = 0.25 for pest control.
It is considered to be good practice to estimate at least one additional analysis of aggregation (other than Taylor's law) because the use of only a single index may be misleading. Although a number of other methods for detecting relationships between the variance and mean in biological samples have been proposed, to date none have achieved the popularity of Taylor's law. The most popular analysis used in conjunction with Taylor's law is probably Iowa's Patchiness regression test but all the methods listed here have been used in the literature.
where s is the variance in the ith sample and mi is the mean of the ith sample
When the population follows a negative binomial distribution, a = 1 and b = k (the exponent of the negative binomial distribution).
This alternative formulation has not been found to be as good a fit as Taylor's law in most studies.
Nachman proposed a relationship between the mean density and the proportion of samples with zero counts:
where p0 is the proportion of the sample with zero counts, m is the mean density, a is a scale parameter and b is a dispersion parameter. If a = b = 0 the distribution is random. This relationship is usually tested in its logarithmic form
Allsop used this relationship along with Taylor's law to derive an expression for the proportion of infested units in a sample
where D2 is the degree of precision desired, zα/2 is the upper α/2 of the normal distribution, a and b are the Taylor's law coefficients, c and d are the Nachman coefficients, n is the sample size and N is the number of infested units.
Binary sampling is not uncommonly used in ecology. In 1958 Kono and Sugino derived an equation that relates the proportion of samples without individuals to the mean density of the samples.
where p0 is the proportion of the sample with no individuals, m is the mean sample density, a and b are constants. Like Taylor's law this equation has been found to fit a variety of populations including ones that obey Taylor's law. Unlike the negative binomial distribution this model is independent of the mean density.
The derivation of this equation is straightforward. Let the proportion of empty units be p0 and assume that these are distributed exponentially. Then
Taking logs twice and re arranging, we obtain the equation above. This model is the same as that proposed by Nachman.
The advantage of this model is that it does not require counting the individuals but rather their presence or absence. Counting individuals may not be possible in many cases particularly where insects are the matter of study.
The equation was derived while examining the relationship between the proportion ( P ) of a series of rice hills infested and the mean severity of infestation ( m ). The model studied was
where a and b are empirical constants. Based on this model the constants a and b were derived and a table prepared relating the values of P and m
where var() is the variance of the sample unit means ( mi ) and m is the overall mean.
An alternative adjustment to the mean estimates is
where MSE is the mean square error of the regression.
This model may also be used to estimate stop lines for enumerative (sequential) sampling. The variance of the estimated means is
where MSE is the mean square error of the regression, α and β are the constant and slope of the regression respectively, sβ2 is the variance of the slope of the regression, N is the number of points in the regression, n is the number of sample units and p is the mean value of p0 in the regression. The parameters a and b are estimated from Taylor's law:
Hughes and Madden have proposed testing a similar relationship applicable to binary observations in cluster, where each cluster contains from 0 to n individuals.
where a, b and c are constants, varobs is the observed variance, and p is the proportion of individuals with a trait (such as disease), an estimate of the probability of an individual with a trait. In logarithmic form, this relationship is
In most cases, it is assumed that b = c, leading to a simple model
This relationship has been subjected to less extensive testing than Taylor's law. However, it has accurately described over 100 data sets, and there are no published examples reporting that it does not works.
A variant of this equation was proposed by Shiyomi et al. () who suggested testing the regression
where varobs is the variance, a and b are the constants of the regression, n here is the sample size (not sample per cluster) and p is the probability of a sample containing at least one individual.
Negative binomial distribution model
A negative binomial model has also been proposed. The dispersion parameter (k) using the method of moments is m2 / ( s2 - m ) and pi is the proportion of samples with counts > 0. The s2 used in the calculation of k are the values predicted by Taylor's law. pi is plotted against 1 - ( k ( k + m ) −1 )k and the fit of the data is visually inspected.
Perry and Taylor have proposed an alternative estimator of k based on Taylor's law.
where Ax is the total number of samples with more than x individuals, N is the total number of individuals, x is the number of individuals in a sample, m is the mean number of individuals per sample and k is the exponent. The value of k has to be estimated numerically.
Goodness of fit of this model can be tested in a number of ways including using the chi square test. As these may be biased by small samples an alternative is the U statistic - the difference between the variance expected under the negative binomial distribution and that of the sample. The expected variance of this distribution is m + m2 / k and
where s2 is the sample variance, m is the sample mean and k is the negative binomial parameter.
The variance of U is
where p = m / k, q = 1 + p, R = p / q and N is the total number of individuals in the sample. The expected value of U is 0. For large sample sizes U is distributed normally.
Note: The negative binomial is actually a family of distributions defined by the relation of the mean to the variance
where a and p are constants. When a = 0 this defines the Poisson distribution. With p = 1 and p = 2, the distribution is known as the NB1 and NB2 distribution respectively.
This model is a version of that proposed earlier by Barlett.
Tests for a common dispersion parameter
The dispersion parameter (k) is
where m is the sample mean and s2 is the variance. If k−1 is > 0 the population is considered to be aggregated; k−1 = 0 the population is considered to be random; and if k−1 is < 0 the population is considered to be uniformly distributed.
Southwood has recommended regressing k against the mean and a constant
where ki and mi are the dispersion parameter and the mean of the ith sample respectively to test for the existence of a common dispersion parameter (kc). A slope (b) value significantly > 0 indicates the dependence of k on the mean density.
An alternative method was proposed by Elliot who suggested plotting ( s2 - m ) against ( m2 - s2 / n ). kc is equal to 1/slope of this regression.
This coefficient (C) is defined as
If the population can be assumed to be distributed in a negative binomial fashion, then C = 100 (1/k)0.5 where k is the dispersion parameter of the distribution.
Cole's index of dispersion
This index (Ic) is defined as
The usual interpretation of this index is as follows: values of Ic < 1, 1, > 1 are taken to mean a uniform distribution, a random distribution or an aggregated distribution.
Because s2 = Σ x2 - (Σx)2, the index can also be written
If Taylor's law can be assumed to hold, then
Lloyd's index of mean crowding (IMC) is the average number of other points contained in the sample unit that contains a randomly chosen point.
where m is the sample mean and s2 is the variance.
Lloyd's index of patchiness (IP) is
It is a measure of pattern intensity that is unaffected by thinning (random removal of points). This index was also proposed by Pielou in 1988 and is sometimes known by this name also.
Because an estimate of the variance of IP is extremely difficult to estimate from the formula itself, LLyod suggested fitting a negative binomial distribution to the data. This method gives a parameter k
where SE(IP) is the standard error of the index of patchiness,var(k) is the variance of the parameter k and q is the number of quadrats sampled..
If the population obeys Taylor's law then
Patchiness regression test
yi here is Lloyd's index of mean crowding. Perform an ordinary least squares regression of mi against y.
In this regression the value of the slope (b) is an indicator of clumping: the slope = 1 if the data is Poisson-distributed. The constant (a) is the number of individuals that share a unit of habitat at infinitesimal density and may be < 0, 0 or > 0. These values represent regularity, randomness and aggregation of populations in spatial patterns respectively. A value of a < 1 is taken to mean that the basic unit of the distribution is a single individual.
Where the statistic s2 / m is not constant it has been recommended to use instead to regress Lloyd's index against am + bm2 where a and b are constants.
The sample size (n) for a given degree of precision (D) for this regression is given by
where a is the constant in this regression, b is the slope, m is the mean and t is the critical value of the t distribution.
Iawo has proposed a sequential sampling test based on this regression. The upper and lower limits of this test are based on critical densities mc where control of a pest requires action to be taken.
where Nu and Nl are the upper and lower bounds respectively, a is the constant from the regression, b is the slope and i is the number of samples.
Kuno has proposed an alternative sequential stopping test also based on this regression.
where Tn is the total sample size, D is the degree of precision, n is the number of samples units, a is the constant and b is the slope from the regression respectively.
Kuno's test is subject to the condition that n ≥ (b - 1) / D2
Parrella and Jones have proposed an alternative but related stop line
where a and b are the parameters from the regression, N is the maximum number of sampled units and n is the individual sample size.
Morisita’s index of dispersion
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.
An alternative formulation is
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
where IMC is Lloyd's index of crowding.
This index is relatively independent of the population density but is affected by the sample size. Values > 1 indicate clumping; values < 1 indicate a uniformity of distribution and a value of 1 indicates a random sample.
Morisita showed that the statistic
is distributed as a chi squared variable with n - 1 degrees of freedom.
A alternative significance test for this index has been developed for large samples.
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.
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
where χ2 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 Id ≥ Mc > 1
When Mc > Id ≥ 1
When 1 > Id ≥ Mu
When 1 > Mu > Id
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.
Southwood's index of spatial aggregation
Southwood's index of spatial aggregation (k) is defined as
where m is the mean of the sample and m* is Lloyd's index of crowding.
Fisher's index of dispersion
This index may be used to test for over dispersion of the population. It is recommended that in applications n > 5 and that the sample total divided by the number of samples is > 3. In symbols
where x is an individual sample value. The expectation of the index is equal to n and it is distributed as the chi-square distribution with n − 1 degrees of freedom when the population is Poisson distributed. It is equal to the scale parameter when the population obeys the gamma distribution.
It can be applied both to the overall population and to the individual areas sampled individually. The use of this test on the individual sample areas should also include the use of a Bonferroni correction factor.
If the population obeys Taylor's law then
Index of cluster size
The index of cluster size (ICS) was created by David and Moore. Under a random (Poisson) distribution ICS is expected to equal 0. Positive values indicate a clumped distribution; negative values indicate a uniform distribution.
where s2 is the variance and m is the mean.
If the population obeys Taylor's law
The ICS is also equal to Katz's test statistic divided by ( n / 2 )1/2 where n is the sample size. It is also related to Clapham's test statistic. It is also sometimes referred to as the clumping index.
Green’s index (GI) is a modification of the index of cluster size that is independent of n the number of sample units.
This index equals 0 if the distribution is random, 1 if it is maximally aggregated and -1 / ( nm - 1 ) if it is uniform.
The distribution of Green's index is not currently known so statistical tests have been difficult to devise for it.
If the population obeys Taylor's law
Binary dispersal index
Binary sampling (presence/absence) is frequently used where it is difficult to obtain accurate counts. The dispersal index (D) is used when the study population is divided into a series of equal samples ( number of units = N: number of units per sample = n: total population size = n x N ). The theoretical variance of a sample from a population with a binomial distribution is
where s2 is the variance, n is the number of units sampled and p is the mean proportion of sampling units with at least one individual present. The dispersal index (D) is defined as the ratio of observed variance to the expected variance. In symbols
where varobs is the observed variance and varbin is the expected variance. The expected variance is calculated with the overall mean of the population. Values of D > 1 are considered to suggest aggregation. D( n - 1 ) is distributed as the chi squared variable with n - 1 degrees of freedom where n is the number of units sampled.
An alternative test is the C test.
where D is the dispersal index, n is the number of units per sample and N is the number of samples. C is distributed normally. A statistically significant value of C indicates overdispersion of the population.
where T is the number of organisms per sample, p is the likelihood of the organism having the sought after property (diseased, pest free, etc), and xi is the number of organism in the ith unit with this property. T must be the same for all sampled units. In this case with n constant
where θ is the parameter of the distribution.
Ma's population aggregation critical density
Ma has proposed a parameter (m0) - the population aggregation critical density - to relate population density to Taylor's law.
A number of statistical tests are known that may be of use in applications.
de Oliveria's statistic
where t is the Poisson parameter, s2 is the variance, m is the mean and n is the sample size. The expected value of s2 - m is zero. This statistic is distributed normally.
If the Poisson parameter in this equation is estimated by putting t = m, after a little manipulation this statistic can be written
This is almost identical to Katz's statistic with ( n - 1 ) replacing n. Again OT is normally distributed with mean 0 and unit variance for large n.
de Oliveria actually suggested that the variance of s2 - m was ( 1 - 2t1/2 + 3t ) / n where t is the Poisson parameter. He suggested that t could be estimated by putting it equal to the mean (m) of the sample. Further investigation by Bohning showed that this estimate of the variance was incorrect. Bohning's correction is given in the equations above.
In 1936 Clapham proposed using the ratio of the variance to the mean as a test statistic (the relative variance). In symbols
For a Possion distribution this ratio equals 1. To test for deviations from this value he proposed testing its value against the chi square distribution with n degrees of freedom where n is the number of sample units. The distribution of this statistic was studied further by Blackman who noted that it was approximately normally distributed with a mean of 1 and a variance ( Vθ ) of
The derivation of the variance was re analysed by Bartlett who considered it to be
For large samples these two formulae are in approximate agreement. This test is related to the later Katz's Jn statistic.
If the population obeys Taylor's law then
A refinement on this test has also been published These authors noted that this test tends to detect overdispersion at higher scales even when this was not present in the data. They noted that the use of the multinomial distribution may be more appropriate than the use of a Poisson distribution for such data. The statistic θ is distributed
where N is the number of sample units, n is the total number of samples examined and xi are the individual data values.
The expectation and variance of θ are
For large N E(θ) is approximately 1 and
If the number of individuals sampled ( n ) is large this estimate of the variance is in agreement with those derived earlier. However, for smaller samples these latter estimates are more precise and should be used.
- Morisita's overlap index
- Natural exponential family
- Scaling pattern of occupancy
- Spatial ecology
- Watson's power law
- Density-mass allometry
- Variance-mass allometry
- Taylor, L. R. (1961). "Aggregation, variance and the mean". Nature. 189 (4766): 732–735. doi:10.1038/189732a0.
- Thomas R. E. Southwood (1966). Ecological methods, with particular reference to the study of insect populations. Methuen.
- Taylor, L. R.; Woiwod, I. P. (1980). "Temporal stability as a density-dependent species characteristic". Journal of Animal Ecology. 49 (1): 209–224. doi:10.2307/4285. JSTOR 4285.
- Taylor LR & Woiwod (1982) Comparative Synoptic Dynamics. I. Relationships Between Inter- and Intra-Specific Spatial and Temporal Variance/Mean Population Parameters. J Animal Ecol 51, 879–906
- Kendal WS & Frost P (1987) Experimental metastasis: a novel application of the variance-to-mean power function. J Natl Cancer Inst 79, 1113-1115
- Kendal WS (1995) A probabilistic model for the variance to mean power law in ecology. Ecological Modelling 80, 293-297
- Keeling M & Grenfell B (1999) Stochastic dynamics and a power law for measles variability. Phil Trans R Soc Lond B Biol Sci 354, 769-776
- Anderson RM & May RM (1989) Epidemiological parameters of HIV transmission. Nature 333, 514–519
- Philippe P (1999) The scale-invariant spatial clustering of leukemia in San Francisco. J Theor Biol 199, 371-381
- Bassingthwaighte JB (1989) Fractal nature of regional myocardial blood flow heterogeneity. Circ Res 65, 578-590
- Kendal WS (2001) A stochastic model for the self-similar heterogeneity of regional organ blood flow. Proc Natl Acad Sci U S A 98, 837-841
- Kendal WS (2003) An exponential dispersion model for the distribution of human single nucleotide polymorphisms. Mol Biol Evol 20 579-590
- Kendal WS (2004) A scale invariant clustering of genes on human chromosome 7. BMC Evol Biol 4(1), 3
- Kendal WS & Jørgensen B (2011) Taylor's power law and fluctuation scaling explained by a central-limit-like convergence. Phys. Rev. E 83,066115
- Kendal, W (2013). "Fluctuation scaling and 1/f noise: shared origins from the Tweedie family of statistical distributions". J Basic Appl Phys. 2 (2): 40–9.
- Kendal WS & Jørgensen B (2015) A scale invariant distribution of the prime numbers. Computation 2015,3,528-540
- Kendal WS & Jørgensen BR (2011) Tweedie convergence: a mathematical basis for Taylor's power law, 1/f noise and multifractality. Phys. Rev E 84, 066120
- Reynolds O (1879) On certain dimensional properties of matter in the gaseous state. Part I. Experimental researches on thermal transpiration of gases through porous plates and on the laws of transpiration and impulsion, including an experimental proof that gas is not a continuous plenum. Part II. On an extension of the dynamical theory of gas, which includes the stresses, tangential and normal, caused by a varying condition of gas, and affords an explanation of the phenomena of transpiration and impulsion. Philosophical Transactions of the Royal Society of London 170:727-845
- Pareto V (1897) Cours D'économie Politique. Volume 2. Lausanne: F. Rouge
- Fisher, RA (1918). "The correlation between relatives on the supposition of Mendelian inheritance". Transactions of the Royal Society of Edinburgh. 52 (2): 399–433. doi:10.1017/S0080456800012163.
- Fisher RA (1921) Studies in crop variation. I. An examination of the yield of dressed grain from Broadbalk. Journal of Agricultural Science 11:107-135
- Neyman J (1926) On the correlation of the mean and the variance in samples drawn from an "infinite" population. Biometrika 18:401-413
- Bartlett M (1936) Some notes on insecticide tests in the laboratory and in the field. J Roy Stat Soc Supplement 3: 185-194
- Smith HF (1938) An empirical law describing heterogeneity in the yield of agricultural crops. J Agric Sci 28: 1-23
- Bliss CI (1941) Statistical problems in estimating populations of Japanese beetle larve. J Econ Entomol 34, 221-232
- Fracker SB & Brischle HA (1944) Measuring the local distribution of Ribes. Ecology 25: 283-303
- Hayman BI & Lowe AD (1961) The transformation of counts of the cabbage aphid (Brevicovyne brassicae (L.)) NZ J Sci 4:271-278
- Taylor LR (1984) Anscombe’s hypothesis and the changing distributions of insect populations. Antenna 8, 62–67
- Taylor LR & Taylor RAJ (1977) Aggregation, migration and population mechanics. Nature 265, 415–421
- Hanski I (1980) Spatial patterns and movements in coprophagous beetles. Oikos 34, 293-310
- Anderson RD, Crawley GM & Hassell M (1982) Variability in the abundance of animal and plant species. Nature 296, 245–248
- Taylor LR, Taylor RAJ, Woiwod IP & Perry JN (1983) Behavioural dynamics. Nature 303, 801-804
- Downing JA (1986) Spatial heterogeneity: evolved behaviour or mathematical artefact? Nature 323, 255-257
- Taylor LR (1988) Specificity of the spatial power-law exponent in ecology and agriculture. Nature 332, 721-722
- Thórarinsson K (1986) Population density and movement: a critique of Δ-models. Oikos 189, 70–81
- Kemp AW (1987) Families of discrete distributions satisfying Taylor’s power law. Biometrics 43, 693-699
- Yamamura K (1990) Sampling scale dependence of Taylor’s power law. Oikos 59, 121-125
- Routledge RD & Swartz TB (1991) Taylor’s power law re-examined. Oikos 60, 107–112
- Tokeshi M (1995) On the mathematical basis of the variance–mean power relationship. Res Pop Ecol 37, 43–48
- Perry JN (1994) Chaotic dynamics can generate Taylor’s power law. Proc Roy Soc Lond Ser B Biol Sci 257, 221–226
- Kilpatrick AM & Ives AR (2003) Species interactions can explain Taylor’s power law for ecological time series. Nature 422, 65–68
- Clark S, Perry JJN & Marshall JP (1996) Estimating Taylor’s power law parameters for weeds and the effect of spatial scale. Weed Res 36, 405-417
- Ramsayer J, Fellous S, Cohen JE & Hochberg ME (2011) Taylor's Law holds in experimental bacterial populations but competition does not influence the slope. Biology Letters
- Jørgensen, Bent (1997). The theory of dispersion models. [Chapman & Hall]. ISBN 978-0412997112.
- Kendal WS (2002) Spatial aggregation of the Colorado potato beetle described by an exponential dispersion model. Ecol Model 151, 261–269
- Cohen, J E; Xu, m; Schuster, W S (2013). "Stochastic multiplicative population growth predicts and interprets Taylor's power law of fluctuation scaling". Proc R soc Lond B biol Sci. 280: 20122955.
- Cohen JE, Xu M (2015) Random sampling of skewed distributions implies Taylor's power law of fluctuation scaling. Proc Natl Acad Sci USA 112(25):7749-7754. doi: 10.1073/pnas.1503824112
- Eisler Z, Bartos I & Kertesz (2008) Fluctuation scaling in complex systems: Taylor’s law and beyond. Adv Phys 57, 89–142
- Fronczak A & Fronczak P (2010) Origins of Taylor’s power law for fluctuation scaling in complex systems. Phys Rev E 81, 066112
- Cohen JE (2016) Statistics of primes (and probably twin primes) satisfy Taylor's Law from ecology. The American Statistician 70 (4) 399-404
- Tweedie MCK (1984) An index which distinguishes between some important exponential families. In: Statistics: Applications and New Directions Proceedings of the Indian Statistical Institute Golden Jubilee International Conference pp 579-604 Eds: JK Ghosh & J Roy, Indian Statistical Institute, Calcutta
- Jørgensen B (1987)Exponential dispersion models. J Roy Stat Soc Ser B 49, 127–162}
- Jørgensen B, Marinez JR & Tsao M (1994) Asymptotic behaviour of the variance function. Scandinavian Journal of Statistics 21, 223–243
- Wilson LT & Room PM (1982) The relative efficiency and reliability of three methods for sampling arthropods in Australian cotton fields. J Aust Entomol Soc 21: 175-181
- Jørgensen B (1997) The theory of exponential dispersion models. Chapman & Hall. London
- Rayner JMV (1985) Linear relations in biomechanics: the statistics of scaling functions. J Zoo 206 (3) 415–439
- Ferris H, Mullens TA, Foord KE (1990) Stability and characteristics of spatial description parameters for nematode populations. J Nematol 22(4):427-439
- Hanski I(1982) On patterns of temporal and spatial variation in animal populations. Ann. zool. Fermici 19: 21—37
- Boag B, Hackett CA, Topham PB (1992) The use of Taylor's power law to describe the aggregated distribution of gastro-intestinal nematodes of sheep. Int J Parasitol 22(3):267–270
- Cohen J E, Xua M (2015) Random sampling of skewed distributions implies Taylor’s power law of fluctuation scaling.Proc. Natl. Acad. Sci. USA 2015 112 (25) 7749–7754
- Reply to Chen: Under specified assumptions, adequate random samples of skewed distributions obey Taylor's law (2015) Proc Natl Acad Sci USA 112 (25) E3157-E3158
- Random sampling of skewed distributions does not necessarily imply Taylor's law (2015) Proc Natl Acad Sci USA 112 (25) E3156
- Banerjee B (1976) Variance to mean ratio and the spatial distribution of animals. Experientia 32: 993-994
- Hughes G, and Madden LV (1992) Aggregation and incidence of disease. Plant Pathology 41(6):657–660
- Patil GP, and Stiteler WM (1974) Concepts of aggregation and their quantification: a critical review with some new results and applications. Researches on Population Ecology 15:238-254
- Turechek WW, Madden LV, Gent DH, Xu XM (2011) Comments regarding the binary power law for heterogeneity of disease incidence. Phytopathology 101(12):1396-1407
- Gosme M, Lucas P (2009) Disease spread across multiple scales in a spatial hierarchy: effect of host spatial structure and of inoculum quantity and distribution. Phytopathology 99(7):833-9. doi: 10.1094/PHYTO-99-7-0833
- Xu X-M, Madden LV (2013) The limits of the binary power law describing spatial variability for incidence data. Plant Pathology 63 (5) 973–982
- Clark SJ & Perry JN (1994) Small sample estimation for Taylor's power law. Environment Ecol Stats 1 (4) 287–302, doi:10.1007/BF00469426
- Wilson LT & Room PM (1983) Clumping patterns of fruit and arthropods in cotton with implications for binomial sampling. Environ Entomol 12:50-54
- Bliss CI and Fisher RA (1953) Fitting the negative binomial distribution to biological data (also includes note on the efficient fitting of the negative binomial). Biometrics 9: 177-200
- Jones VP (1991) Binomial sampling plans for tentiform leafminer (Lepidoptera: Gracillariidae) on apple in Utah. J Econ Entomol 84: 484-488
- Katz L (1965) United treatment of a broad class of discrete probability distributions. in Proceedings of the International Symposium on Discrete Distributions. Montreal
- Jewel W, Sundt B (1981) Improved approximations for the distribution of a heterogeneous risk portfolio. Bull Assoc Swiss Act 81: 221-240
- Foley P (1994) Predicting extinction times from environmental stochasticity and carrying capacity. Conserv Biol 8: 124-137
- Pertoldi C, Bach LA, Loeschcke V (2008) On the brink between extinction and persistence. Biol Direct 3: 47 doi: 10.1186/1745-6150-3-47
- Halley J, Inchausti P (2002) Lognormality in ecological time series. Oikos 99:518–530. doi: 10.1034/j.1600-0706.2002.11962.x
- Southwood TRE & Henderson PA (2000) Ecological methods. 3rd ed. Blackwood, Oxford
- Service MW (1971) Studies on sampling larval populations of the Anopheles gambiae complex. Bull World Health Organ 45: 169-180
- Southwood TRE (1978) Ecological methods. Chapman & Hall, London, England
- Karandinos MG (1976) Optimum sample size and comments on some published formulae. Bull Entomol Soc Am 22:417-421
- Ruesink WG (1980) Introduction to sampling theory, in Kogan M & Herzog DC (eds.) Sampling Methods in Soybean Entomology. Springer-Verlag New York, Inc, New York. pp 61-78
- Schulthess F, Bosque-Péreza NA, Gounoua S (1991) Sampling lepidopterous pests on maize in West Africa. Bull Entomol Res 81 (3) 297-301
- Bisseleua DHB, Yede, & Vida S (2011) Dispersion models and sampling of cacao mirid bug Sahlbergella singularis (Hemiptera: Miridae) on theobroma cacao in southern Cameroon. Environ Entomol 40(1): 111-119
- Wilson LT, Gonzalez D & Plant RE(1985) Predicting sampling frequency and economic status of spider mites on cotton. Proc. Beltwide Cotton Prod Res Conf, National Cotton Council of America, Memphis, TN pp 168-170
- Green RH (1970) On fixed precision level sequential sampling. Res Pop Ecol 12 (2) 249-251 doi:10.1007/BF02511568
- Serraa GV, La Porta NC, Avalos S, Mazzuferi V (2012) Fixed-precision sequential sampling plans for estimating alfalfa caterpillar, Colias lesbia, egg density in alfalfa, Medicago sativa, fields in Córdoba, Argentina. J Insect Sci 13: 41
- Myers JH (1978). Selecting a measure of dispersion. Environ Entomol 7: 619–621
- Iwao S (1968) A new regression method for analyzing the aggregation pattern of animal populations. Res Popul Ecol 10: 1–20
- Nachman G (1981) A mathematical model of the functional relationship between density and spatial distribution of a population. J Anim Ecol 50: 453-460
- Allsopp PG (1991) Binomial sequential sampling of adult Saccharicoccus sacchari on sugarcane. Entomologia Experimentalis et Applicata 60 (3) 213–218
- Kono T & Sugino T (1958) On the estimation of density of rice stem infested by the rice stem borer. Japanese J Appl Entomol and Zool 2: 184
- Binns MR & Bostonian NJ (1990) Robustness in empirically based binomial decision rules for integrated pest management. J Econ Entomol 83: 420-442
- Nachman G (1984) Estimates of mean population density and spatial distribution of Tetranychus urticae (Acarina: Tetranychidae) and Phytoseiulus persimilis (Acarina: Phytoseiidae) based upon the proportion of empty sampling units. J Appl Ecolog 21: 903-991
- Schaalje GB ,Butts RA, Lysyk TL (1991) Simulation studies of binomial sampling: a new variance estimator and density pre&ctor, with special reference to the Russian wheat aphid (Homoptera: Aphididae). J Econ Entomol 84: 140-147
- Shiyomi M, Egawa T, Yamamoto Y (1998) Negative hypergeometric series and Taylor's power law in occurrence of plant populations in semi-natural grassland in Japan. Proceedings of the 18th International Grassland Congress on grassland management. The Inner Mongolia Univ Press pp35-43 (1998)
- Wilson L T & Room PM (1983) Clumping patterns of fruit and arthropods in cotton, with implications for binomial sampling. Environ Entomol 12: 50-54
- Perry JN & Taylor LR(1986). Stability of real interacting populations in space and time: implications, alternatives and negative binomial. J Animal Ecol 55: 1053–1068
- Elliot JM (1977) Some methods for the statistical analysis of samples of benthic invertebrates. 2nd ed. Freshwater Biological Association, Cambridge, United Kingdom
- Cole LC (1946) A theory for analyzing contagiously distributed populations. Ecol 27: 329-341
- Lloyd M (1967) Mean crowding. J Anim Ecol 36: 1-30
- Iwao S & Kuno E (1968) Use of the regression of mean crowding on mean density for estimating sample size and the transformation of data for the analysis of variance. Res Pop Ecology, 10, 210–214
- Ifoulis AA, Savopoulou-Soultani M (2006) Developing optimum sample size and multistage sampling plans for Lobesia botrana (Lepidoptera: Tortricidae) larval infestation and injury in northern Greece. J Econ Entomol 99(5):1890–1898
- Ho CC (1993) Dispersion statistics and sample size estimates for Tetranychus kanzawai (Acari: Tetranychidae) on mulberry. Environ Entomol 22: 21-25
- Iwao S (1975) A new method of sequential sampling to classify populations relative to a critical density. Res Popul Ecol 16: 281-28
- Kuno E (1969) A new method of sequential sampling to obtain the population estimates with a fixed level of precision. Res. Pooul. Ecol 11: 127-136
- Parrella MP, Jones VP (1985) Yellow traps as monitoring tools for Liriomyza trifolii (Diptera: Agromyzidae) in chrysanthemum greenhouses. J Econ Entomol 78, 53-56
- Morisita M (1959) Measuring the dispersion and the analysis of distribution patterns. Memoires of the Faculty of Science, Kyushu University Series E. Biol 2:215-235
- Pedigo LP & Buntin GD (1994) Handbook of sampling methods for arthropods in agriculture. CRC Boca Raton FL
- 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, 35-54
- Elliot JM (1977) Statistical analysis of samples of benthic invertebrates. Freshwater Biological Association. Ambleside
- Fisher RA (1925) Statistical methods for research workers. Hafner, New York
- Hoel P (1943) On the indices of dispersion. Ann Math Statist 14: 155
- David FN & Moore PG (1954) Notes on contagious distributions in plant populations. Ann Bot of London 18:47-53
- Green RH (1966) Measurement of non-randomness in spatial distributions. Res Pop Ecol 8:1-7
- Gottwald TR, Bassanezi RB, Amorim L, Bergamin-Filho A (2007) Spatial pattern analysis of citrus canker-infected plantings in São Paulo, Brazil, and augmentation of infection elicited by the Asian leafminer. Phytopathology 97(6):674-683
- Hughes G & Madden LV (1993) Using the beta-binomial distribution to describe aggregated patterns of disease incidence. Phytopathology 83:759-763
- Fleiss JL (1981) Statistical methods for rates and proportions. 2nd ed. Wiley, Newy York, USA
- Ma ZS (1991) Further interpreted Taylor’s Power Law and population aggregation critical density. Trans Ecol Soc China (1991) 284–288
- de Oliveria T (1965) Some elementary tests for mixtures of discrete distributions, in Patil, GP ed., Classical and contagious discrete distributions. Calcutá, Calcutta Publishing Society pp379-384
- Bohning D (1994) A note on a test for Poisson overdispersion. Biometrica 81(2) 418-419
- Ping S (1995) Further study on the statistical test to detect spatial pattern. Biometrical J 37(2) 199–203
- Clapham AR (1936) Overdispersion in grassland communities and the use of statistical methods in plant ecology. J Ecol 14: 232
- Blackman GE (1942) Statistical and ecological studies on the distribution of species in plant communities. I. Dispersion as a factor in the study of changes in plant populations. Ann Bot N.s. vi: 351
- Greig-Smith P (1952) The use of random and contiguous quadrats in the study of the structure of plant communities. Ann. Bot. 16:293-316
- Gosset E, Louis B (1986) The binning analysis - Towards a better significance test. Astrophysics Space Sci 120 (2) 263-306