| 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]