data(multitrait)
multitrait <- subset(multitrait, chr=1,
ind=!apply(multitrait$pheno,
1, function(a) any(is.na(a))))
data(locations)
multiloc <- addloctocross(multitrait,locations)
multiloc <- calc.genoprob(multiloc)
results <- scanall(multiloc, method="hk")
mqmplot.cistrans(results, multiloc, 5, FALSE, TRUE)
Run the code above in your browser using DataLab