Calculate heterozygosities
Calculate heterozygosites, by individual or by marker
calc_het(probs, by = c("individual", "marker"), omit_x = TRUE, cores = 1)probs |
List of arrays of genotype probabilities, as
calculated by |
by |
Whether to summarize by individual or marker |
omit_x |
If TRUE, omit the X chromosome. |
cores |
Number of CPU cores to use, for parallel calculations.
(If |
The result is a vector of estimated heterozygosities
iron <- read_cross2(system.file("extdata", "iron.zip", package="qtl2"))
p <- calc_genoprob(iron, err=0.002)
# heterozygosities by individual
het_ind <- calc_het(p)
# heterozygosities by marker
het_mar <- calc_het(p, "marker")Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.