library(qtl2)
grav2 <- read_cross2(system.file("extdata", "grav2.zip", package="qtl2"))
map <- insert_pseudomarkers(grav2$gmap, step=1)
probs <- calc_genoprob(grav2, map, error_prob=0.002)
dir <- tempdir()
fprobs <- fst_genoprob(probs, "grav2", dir, overwrite=TRUE)
# subset probabilities
fprobs2 <- subset(fprobs, chr=1:2)
# use object to get the full probabilities back
fprobs5 <- fst_restore(fprobs2)
# clean up: remove all the files we created
unlink(fst_files(fprobs))
Run the code above in your browser using DataLab