theta.s {ape} | R Documentation |
This function computes the population parameter THETA using the
number of segregating sites s
in a sample of n
DNA sequences.
theta.s(s, n, variance = FALSE)
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 . |
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
).
The number of segregating sites needs to be computed beforehand.
Emmanuel Paradis paradis@isem.univ-montp2.fr
Watterson, G. (1975) On the number of segragating sites in genetical models without recombination. Theoretical Population Biology, 7, 256276.
Tajima, F. (1989) Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics, 123, 585595.