if (FALSE) data(hyper)
hyper <- calc.genoprob(hyper, step=2.5)
out <- scanone(hyper)
out.cim <- cim(hyper, n.marcovar=3)
plot(out, out.cim, chr=c(1,4,6,15), col=c("blue", "red"))
add.cim.covar(out.cim, chr=c(1,4,6,15))
Run the code above in your browser using DataLab