ape (version 1.8-4)

theta.s: Population Parameter THETA using Segregating Sites in DNA Sequences

Description

This function computes the population parameter THETA using the number of segregating sites s in a sample of n DNA sequences.

Usage

theta.s(s, n, variance = FALSE)

Value

a numeric vector of length one with the estimated theta (the default), or of length two if the standard error is returned (variance = TRUE).

References

Watterson, G. (1975) On the number of segragating sites in genetical models without recombination. Theoretical Population Biology, 7, 256--276.

Tajima, F. (1989) Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics, 123, 585--595.

See Also

theta.h, theta.k, seg.sites, nuc.div

Examples

data(woodmouse)
y <- seg.sites(woodmouse)
s <- length(y)
n <- length(woodmouse)
theta.s(s, n)
theta.s(s, n, variance = TRUE)