cbind.feather_genoprob: Join genotype probabilities for different chromosomes

Description Usage Arguments Value Examples

Description

Join multiple genotype probability objects, as produced by feather_genoprob for different individuals.

Usage

1
2
## S3 method for class 'feather_genoprob'
cbind(..., fbase, fdir = NULL)

Arguments

...

Genotype probability objects as produced by feather_genoprob. Must have the same set of individuals.

fbase

Base of fileame for feather database. Needed if objects have different feather databases.

fdir

Directory for feather database.

Value

A single genotype probability object.

Examples

1
2
3
4
5
6
7
8
9
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,1:2], map, error_prob=0.002)
probsB <- calc_genoprob(grav2[1:5,3:4], map, error_prob=0.002)
dir <- tempdir()
fprobsA <- feather_genoprob(probsA, "exampleAc", dir)
fprobsB <- feather_genoprob(probsB, "exampleBc", dir)
fprobs <- cbind(fprobsA, fprobsB, fbase = "exampleABc")

byandell/qtl2feather documentation built on May 13, 2019, 9:30 a.m.