# Simulated F2 Population
# 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")
f2cross <- calc.genoprob(f2cross)
if (FALSE) # Permutations to obtain significance threshold
f2result <- mqmpermutation(f2cross, n.perm=1000, method="permutation")
f2result <- mqmpermutation(f2cross, n.perm=2, batchsize=2, method="permutation")
# Plot results
mqmplot.permutations(f2result)
Run the code above in your browser using DataLab