Watterson estimator

In population genetics, the Watterson estimator is a method for estimating the population mutation rate, $heta\; =\; 4N\_emu$, where $N\_e$ is the effective population size and $mu$ is the per-generation mutation rate of the population of interest (harvtxt|Watterson|1975). The assumptions made are that there is a sample of "n" haploid individuals from the population of interest, that there are an infinitely many alleles possible, and that $n\; ll\; N\_e$.

The estimate of $heta$, often denoted as $\{hat\; heta\_w\}$, is

: $\{hat\; heta\_w\}\; =\; \{\; K\; over\; a\_n\; \},$

where "K" is the number of segregating sites in the sample and

: $a\_n\; =\; sum^\{n-1\}\_\{i=1\}\; \{1\; over\; i\}$

is the ("n" − 1)th harmonic number.

This estimate is based on coalescent theory. Watterson's estimator is commonly used for its simplicity. The estimator is unbiased and the variance of the estimator decreases with increasing sample size and/or recombination rate. However, the estimator can be biased by population structure. For example, $\{hat\; heta\_w\}$ is downwardly biased in an exponentially growing population.

See also

* Tajima estimator
* Fu estimator
* Coupon collector's problem
* Ewens sampling formula

References

*