data(fake.bc)
fake.bc <- calc.genoprob(fake.bc)
fake.bc <- subset(fake.bc, chr=18:19)
out2a <- scantwo(fake.bc, method="hk")
out2b <- scantwo(fake.bc, pheno.col=2, method="hk")
out2 <- c(out2a, out2b)
Run the code above in your browser using DataLab