View source: R/get.largest.mode.R
get.last.mode | R Documentation |
This function grabs the last peak or shoulder.
get.last.mode(x, y, rare.prop = 1/5, verbose = TRUE, remove.tail = TRUE, ...)
x |
values of midpoints from |
y |
values of counts from |
rare.prop |
Proportion used to remove rare intensity tail |
verbose |
print diagnostic output |
remove.tail |
Remove rare intensity tail |
... |
arguments to be passed to |
Value of x that is the last peak
data(t1.voi.hist) system.time({ y = t1.voi.hist$counts x = t1.voi.hist$mids x = x[!is.na(y)]; y = y[!is.na(y)] # 20 used for speed of example nawm_peak = get.last.mode(x, y, k=20) plot(t1.voi.hist, border="red") abline(v=nawm_peak) })
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.