# QTL at chromosome 3
f2qtl <- c(3,15,1,0)
# Mouse genetic map
data(map10)
map10 <- map10[1:3]
# Simulate a F2 Cross
f2cross <- sim.cross(map10,f2qtl,n=100,type="f2")
if (FALSE) # Bootstrap MQM mapping on the f2cross
f2result <- mqmpermutation(f2cross,scanfunction=mqmscan)
f2result <- mqmpermutation(f2cross,scanfunction=mqmscan, n.perm=2, batchsize=2)
# Create a permutation object
f2perm <- mqmprocesspermutation(f2result)
# What LOD score is considered significant?
summary(f2perm)
Run the code above in your browser using DataLab