hierfstat: theta.Watt.dosage – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

theta.Watt.dosage

Estimates θ_{Watterson} from dosage data

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

Estimation and Tests of Hierarchical F-Statistics

v0.5-10

GPL (>= 2)

Authors

Jerome Goudet [aut, cre], Thibaut Jombart [aut], Zhian N. Kamvar [ctb], Eric Archer [ctb], Olivier Hardy [ctb]

Initial release

2021-11-16