data(fake.bc)
fake.bc <- calc.genoprob(fake.bc)
operm1 <- scantwo(fake.bc, pheno.col=1, method="hk", n.perm=50)
operm2 <- scantwo(fake.bc, pheno.col=2, method="hk", n.perm=50)
operm <- cbind(operm1, operm2)
Run the code above in your browser using DataLab