Learn R Programming

qtl (version 1.38-4)

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 class 'scantwo':
subset(x, chr, lodcolumn, \dots)

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
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.

Value

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

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