Description Usage Arguments Value See Also Examples
View source: R/rbind_fst_genoprob.R
Join multiple genotype probability objects, as produced by
fst_genoprob()
for different individuals.
1 2 |
... |
Genotype probability objects as produced by
|
fbase |
Base of fileame for fst database. Needed if objects have different fst databases. |
fdir |
Directory for fst database. |
overwrite |
If FALSE (the default), refuse to overwrite existing |
quiet |
If TRUE, don't show any messages. Passed to |
A single genotype probability object.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | library(qtl2)
grav2 <- read_cross2(system.file("extdata", "grav2.zip", package="qtl2"))
map <- insert_pseudomarkers(grav2$gmap, step=1)
probsA <- calc_genoprob(grav2[1:5,], map, error_prob=0.002)
probsB <- calc_genoprob(grav2[6:12,], map, error_prob=0.002)
dir <- tempdir()
fprobsA <- fst_genoprob(probsA, "exampleAr", dir, overwrite=TRUE)
fprobsB <- fst_genoprob(probsB, "exampleBr", dir, overwrite=TRUE)
# use rbind to combine probabilities for same chromosomes but different individuals
fprobs <- rbind(fprobsA, fprobsB, fbase = "exampleABr")
# clean up: remove all the files we created
unlink(fst_files(fprobsA))
unlink(fst_files(fprobsB))
unlink(fst_files(fprobs))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.