mutRate: Mutation Rate

In aumath-advancedr2019/simDNA: Simulate and analyze DNA

Description

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

Usage

mutRate(SFS)

Arguments

SFS	vector with the site frequency spectrum.

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.

Value

List containing the following components:

Watterson	Watterson's estimator.
pairwDiff	The pairwise difference estimator.

References

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

Examples

mutRate(c(2,1,0,0,0,1,0))

aumath-advancedr2019/simDNA documentation built on Nov. 27, 2019, 11 a.m.