Learn R Programming

qtl (version 1.70)

subset.scantwo: Subsetting the results of a 2-d genome scan

Description

Pull out a specified set of chromosomes and/or LOD columns from scantwo output.

Usage

# S3 method for scantwo
subset(x, chr, lodcolumn, ...)

Value

The input scantwo object, but with only the specified subset of the data.

Arguments

x

An object of class scantwo, output from scantwo.

chr

Optional vector specifying which chromosomes to keep. This should be a vector of character strings referring to chromosomes by name; numeric values are converted to strings. Refer to chromosomes with a preceding - to have all chromosomes but those considered. A logical (TRUE/FALSE) vector may also be used.

lodcolumn

A vector specifying which LOD columns to keep (or, if negative), omit. These should be between 1 and the number of LOD columns in the input x.

...

Ignored at this point.

Author

Karl W Broman, broman@wisc.edu

See Also

summary.scantwo, scantwo

Examples

Run this code
data(fake.bc)
fake.bc <- subset(fake.bc, chr=16:19)
fake.bc <- calc.genoprob(fake.bc)
out <- scantwo(fake.bc, method="hk", pheno.col=1:2)

summary(subset(out, chr=18:19))

Run the code above in your browser using DataLab