LOD score peaks
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.
Magnus Dehli Vigeland
## 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)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.