qtl (version 1.01-9)

max.scanone: Maximum peak in genome scan

Description

Print the row of the output from scanone that corresponds to the maximum LOD, genome-wide.

Usage

## S3 method for class 'scanone':
max(\dots, chr, na.rm=TRUE)

Arguments

...

An object of the form output by the function scanone: a data.frame whose third column is the LOD score. Note that while multiple objects may be input, only the first is considered.

chr

Numeric or character vector indicating chromosomes on which to find the maximum LOD score.

na.rm

A logical indicating whether missing values should be removed.

Value

• An object of class summary.scanone, to be printed by print.summary.scanone. This is a data.frame with one row, corresponding to the maximum LOD peak either genome-wide or for the particular chromosome specified.

See Also

scanone, plot.scanone, summary.scanone

Examples

data(listeria)
listeria <- subset(listeria,chr=c(5,13))
listeria <- calc.genoprob(listeria, step=2.5)
out <- scanone(listeria, model="2part")
# Maximum peak for LOD(p,mu)
max(out)
max(out,chr=13)

# Maximum peak for LOD(p)
max(out[,-3])

# Maximum peak for LOD(mu)
max(out[,-(3:4)])