paramlink (version 1.1-5)

lod.peaks: 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

Value

A list of data frames.

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

# NOT RUN {
## 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)

# }