data(hyper)
hyper <- subset(hyper, chr=1:4)
hyperfilled <- fill.geno(hyper)
cofactors <- mqmsetcofactors(hyperfilled,4)
result <- mqmscan(hyperfilled,cofactors)
mqmgetmodel(result)
plot(mqmgetmodel(result))
Run the code above in your browser using DataLab