## Not run:
# qtl = scanone(pheno = pheno, probs = probs, addcovar = addcovar, snps = anps)
# perms = scanone.perm(pheno = pheno, probs = probs, addcovar = addcovar, snps = anps)
# sig.thr = get.sig.thr(perms)
# plot(qtl)
# add.sig.thr(sig.thr, chrsum = cumsum(get.chr.lengths()))
# # Could also run:
# plot(qtl, sig.thr = sig.thr)
# ## End(Not run)
Run the code above in your browser using DataLab