Join genome scan results for different phenotypes.
Join multiple scan1() results for different phenotypes;
must have the same map.
## S3 method for class 'scan1' cbind(...)
... |
Genome scan objects of class |
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.
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.
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)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.