| pHetErrOnStateCount {diemr} | R Documentation |
Hybrid index, heterozygosity, error rate
Description
Using genotype allele counts, calculates the hybrid index, heterozygosity and error rate in a single individual.
Usage
pHetErrOnStateCount(sCount)
Arguments
sCount |
a numeric vector of length 4 with allele counts for missing data, homozygots for allele 1, heterozygots, and homozygots for allele 2. |
Details
Allele counts are genomic state counts multiplied by ploidy. As different compartments might have different ploidies (e.g. autosomal markers, sex chromosomes, mitochondrial markers), allele counts should be calculated per compartment and then summarised to obtain the correct genomic allele counts. When all individuals in each compartmenst have the same ploidy, state counts do not need to be corrected.
Value
Returns a named numeric vector with three values: p - hybrid index, Het - heterozygosity, Err - error rate.
Examples
pHetErrOnStateCount(sCount = c(2, 4, 2, 6))
# p Het Err
# 0.5833333 0.1666667 0.1428571