rbind.scan1 {qtl2} | R Documentation |
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 class 'scan1'
rbind(...)
Arguments
... |
Genome scan objects of class |
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.
Value
An object of class '"scan1", like the inputs, but with the results for different sets of chromosomes combined.
See Also
Examples
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)
[Package qtl2 version 0.36 Index]