#Simulated F2 Population
f2qtl <- c(3,15,1,0) # QTL at chromosome 3
data(map10) # Mouse genetic map
map10 <- map10[1:3]
f2cross <- sim.cross(map10,f2qtl,n=100,type="f2") # Simulate a F2 Cross
f2cross <- mqmaugment(f2cross)
f2result <- mqmscan(f2cross) # Do a MQM scan of the genome
mqmplot.singletrait(f2result) # Use our fancy plotting routine
Run the code above in your browser using DataLab