theta.Watt.dosage | R Documentation |
Estimates θ_{Watterson}=S/a, where S is the number of segregating sites in a set of sequences and a=1/∑_i^{n-1} i.
theta.Watt.dosage(dos,L=NULL)
dos |
a ni X nl dosage matrix containing the number of derived/alternate alleles each individual carries at each SNP |
L |
the length of the sequence |
if L=NULL
(default), returns θ_{Watterson}, else return θ_{Watterson}/L
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.