iron <- read_cross2(system.file("extdata", "iron.zip", package="qtl2"))
iron <- iron[,"8"]
gmap <- insert_pseudomarkers(iron$gmap, step=1)
pr <- calc_genoprob(iron, gmap, error_prob=0.002)
pmar <- find_marker(gmap, 8, 40)
pr_8_40 <- pull_genoprobpos(pr, pmar)
pr_8_40_alt <- pull_genoprobpos(pr, gmap, 8, 40)
Run the code above in your browser using DataLab