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), the default 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
(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.