qtl2 (version 0.24)

rbind.scan1: Join genome scan results for different chromosomes.

Description

Join multiple scan1() results for different chromosomes; must have the same set of lod score column.

Usage

# S3 method for scan1
rbind(...)

Arguments

...

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

Value

An object of class `"scan1", like the inputs, but with the results for different sets of chromosomes combined.

Details

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

If hsq present, we simply rbind() the contents.

See Also

cbind.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)
phe <- grav2$pheno[,1,drop=FALSE]

out1 <- scan1(probs[,1], phe) # chr 1
out2 <- scan1(probs[,5], phe) # chr 5
out <- rbind(out1, out2)

# }