data(fake.f2)
fake.f2 <- calc.genoprob(fake.f2)
fake.f2 <- subset(fake.f2, chr=18:19)
out.hk <- scanone(fake.f2, method="hk")
out.np <- scanone(fake.f2, model="np")
out <- c(out.hk, out.np, labels=c("hk","np"))
plot(out, lod=1:2, col=c("blue", "red"))
Run the code above in your browser using DataLab