ape (version 0.2-1)

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)

Arguments

s
a numeric giving the number of segregating sites.
n
a numeric giving the length of the sequences (i.e. the number of sites under study).
variance
a logical indicating whether the variance of the estimated THETA should be returned (TRUE), thedefault being FALSE.

Value

  • a numeric vector of length one with the estimated theta (the default), or of length two if the standard error is returned (standard.error = 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