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

Description Usage Arguments Details Value See Also Examples

View source: R/cbind_scan1.R

Description

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

Usage

1
2
## S3 method for class 'scan1'
cbind(...)

Arguments

...

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

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.

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.

See Also

rbind.scan1(), scan1()

Examples

1
2
3
4
5
6
7
8
9
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)

qtl2 documentation built on Oct. 18, 2021, 1:06 a.m.