theta.Watt.dosage: Estimates theta_{Watterson} from dosage data

View source: R/div_dosage.R

theta.Watt.dosageR 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


jgx65/hierfstat documentation built on Nov. 21, 2024, 9:06 p.m.