# NOT RUN {
data(ch3a)
# phenotype problem
pairs(ch3a$pheno)
ch3a$pheno[ch3a$pheno[,4]==0,] # individual 159
# individuals with similar genotypes
library(qtl)
cg <- comparegeno(ch3a)
hist(cg, breaks=200)
max(cg[cg < 1])
which(cg == max(cg[cg < 1]), arr.ind=TRUE)
# }
Run the code above in your browser using DataLab