powered by
Creates a matrix of opposing homozygotes from the genotype matrix.
ohg(genotypeMatrix)
Returns a square matrix (sample X sample) with the pairwise counts of opposing homozygotes.
matrix genotype (Data should be numeric. Use 0, 1 and 2 respectively for AA, AB and BB. Use 9 for missing data)
matrix
Ferdosi, M. H., & Boerner, V. (2014). A fast method for evaluating opposing homozygosity in large SNP data sets. Livestock Science.
rpoh
genotype <- matrix(c( 2,1,0, 2,0,0, 0,0,2 ), byrow = TRUE, ncol = 3) ohg(genotype)
Run the code above in your browser using DataLab