getHeterozygosityBySample: Compute heterozygosity by sample

Description Usage Arguments Details Value

View source: R/seqarray_process_popgen.R

Description

Compute heterozygosity by sample

Usage

1

Arguments

gdsfile

a SeqVarGDSClass object

Details

Follows the method outlined in Manske et al. 2012. Briefly, for each sample, the reference and and alternate allele frequencies are computed as the proportion of reads covering each allele. Then heterozygosity at a SNP is 1 - (raf^2 + aaf^2)

Value

a numeric matrix of size l by n where l is the number of samples and n is the number of SNPs.


bahlolab/moimix documentation built on April 16, 2020, 8:40 a.m.