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