Learn R Programming

hierfstat (version 0.5-11)

theta.Watt.dosage: 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\)