paramlink (version 1.1-6)

lod.peaks: LOD score peaks

Description

Identify LOD score peaks

Usage

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

Value

A list of data frames.

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.

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)