Mutation Rate

The function calculates two different estimates of the mutation rate: Watterson's estimator and the pairwise difference estimator.

mutRate(SFS)

Arguments

SFS	vector with the site frequency spectrum.

Value

List containing the following components:

Watterson

Watterson's estimator.

pairwDiff

The pairwise difference estimator.

Details

See SFS for explanation of the site frequency spectrum.

If one runs mutRate with a vector that contains anything other than natural numbers, an error will occur.

For details about Watterson's estimator and the pairwise difference estimator, consult the analyzeDNA vignette.

References

Wakeley J. (2009) Coalescent Theory: An Introduction. Colorado: Roberts and Company Publishers.

Examples

mutRate(c(2,1,0,0,0,1,0))
#> $Watterson
#> [1] 1.5427
#> 
#> $pairwDiff
#> [1] 1.357143
#>