lod.peaks | R Documentation |
Identify LOD score peaks
lod.peaks(x, threshold, width = 1)
x |
a |
threshold |
a single numeric |
width |
a positive integer |
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.
A list of data frames.
linkres
, lod
, merlin
,
## 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.