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

In rqtl/qtl2scan: Genome Scans for QTL Experiments

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 as produced by scan1(). Must have the same lod score columns.

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

A single genome scan object with the results for different sets of chromosomes combined.

Examples

library(qtl2geno)
grav2 <- read_cross2(system.file("extdata", "grav2.zip", package="qtl2geno"))
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)

rqtl/qtl2scan documentation built on May 28, 2019, 2:36 a.m.