max.scanone {qtl}R Documentation

Maximum peak in genome scan

Description

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

Usage

max.scanone(..., 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.

Author(s)

Karl W Broman, kbroman@jhsph.edu

See Also

scanone, plot.scanone, summary.scanone

Examples

data(listeria)

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)])

[Package qtl version 0.98-57 Index]