data(fake.bc)
fake.bc <- subset(fake.bc, chr=18:19)
fake.bc <- calc.genoprob(fake.bc, step=0)
operm <- scantwo(fake.bc, method="hk", pheno.col=1:2, n.perm=5)
operm2 <- subset(operm, lodcolumn=2)
# alternatively
operm2alt <- operm[,2]
Run the code above in your browser using DataLab