scan1() results for different phenotypes;
must have the same map.
Genome scan objects of class
weights do not match between objects, we
omit this information.
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.
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.