This function is a follow-up of
getPeaks() by merging closely located ambiguous peaks together.
findRobustPeaks(fit, thr = 1.5)
a density fit returned by
threshold to separate adjacent peaks. See details.
If two peaks are closely located, one way to tell if they are actually one peak or otherwise is by comparing the peak height with peak width. By assuming both peaks are normally distributed, this is equivalent to compare the peak height with the height of the local minima (h) between them.
thr is the minimum ratio of the lower peak height over h. If this ratio is larger than
thr, the two peaks are considered separated.
a numeric matrix in the same format as returned by
1 2 3 4 5 6 7