powered by
Creates a matrix of opposing homozygotes from the genotype matrix.
ohg(genotypeMatrix)
matrix genotype (Data should be numeric. Use 0, 1 and 2 respectively for AA, AB and BB. Use 9 for missing data)
matrix
Returns a square matrix (sample X sample) with the pairwise counts of opposing homozygotes.
rpoh
# NOT RUN { 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