qtl2 (version 0.24)

cbind.scan1: Join genome scan results for different phenotypes.

Description

Join multiple scan1() results for different phenotypes; must have the same map.

Usage

# S3 method for scan1
cbind(...)

Arguments

...

Genome scan objects of class "scan1", as produced by scan1(). Must have the same map.

Value

An object of class `"scan1", like the inputs, but with the lod score columns from the inputs combined as multiple columns in a single object.

Details

If components addcovar(), Xcovar, intcovar, weights do not match between objects, we omit this information.

If hsq present but has differing numbers of rows, we omit this information.

See Also

rbind.scan1(), scan1()

Examples

# NOT RUN {
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)
phe1 <- grav2$pheno[,1,drop=FALSE]
phe2 <- grav2$pheno[,2,drop=FALSE]

out1 <- scan1(probs, phe1) # phenotype 1
out2 <- scan1(probs, phe2) # phenotype 2
out <- cbind(out1, out2)

# }