lod.peaks {paramlink}

LOD score peaks

Description

Identify LOD score peaks

Usage

lod.peaks(x, threshold, width = 1)

Arguments

x

a linkres object

threshold

a single numeric

width

a positive integer

Details

The function first transforms x to a data frame (using as.data.frame.linkres with sort=T. A peak is defined a run of at least width consecutive markers with LOD score above or equal to threshold. If possible, one flanking marker is included on each side of the peak.

Value

A list of data frames.

See Also

linkres, lod, merlin,

Examples

## minimal example
x = linkdat(toyped, model=1)
res = lod(x)
peak1 = lod.peaks(res, threshold=0)
peak2 = lod.peaks(res, threshold=0, width=2)
peak3 = lod.peaks(res, threshold=1)
stopifnot(length(peak1)==1, nrow(peak1[[1]])==1, length(peak2)==0, length(peak3)==0)