data(hyper) # Hyper dataset
hyper <- hyper[1:5]
hyperfilled <- fill.geno(hyper)
# Automatic cofactors every third marker
cofactors <- mqmsetcofactors(hyperfilled,3)
result <- mqmscan(hyperfilled,cofactors) # Backward model selection
mqmgetmodel(result)
#Manual cofactors at markers 3,6,9,12,40 and 60
cofactors <- mqmsetcofactors(hyperfilled,cofactors=c(3,6,9,12,40,60))
result <- mqmscan(hyperfilled,cofactors) # Backward model selection
mqmgetmodel(result)
Run the code above in your browser using DataLab