theta.Watt.dosage {hierfstat} | R Documentation |
Estimates \theta_{Watterson}
from dosage data
Description
Estimates \theta_{Watterson}=S/a
, where S
is the number of segregating sites
in a set of sequences and a=1/\sum_i^{n-1} i
.
Usage
theta.Watt.dosage(dos,L=NULL)
Arguments
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 |
Value
if L=NULL
(default), returns \theta_{Watterson}
, else return \theta_{Watterson}/L
[Package hierfstat version 0.5-11 Index]