max.scanone {qtl} | R Documentation |
Print the row of the output from scanone
that
corresponds to the maximum LOD, genome-wide.
max.scanone(..., chr, na.rm=TRUE)
... |
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. |
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.
Karl W Broman, kbroman@jhsph.edu
scanone
, plot.scanone
,
summary.scanone
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)])