R/modes.R

modes <-
function (y) {
    test <- try(d <- density(y, bw = "nrd", adjust = 3, n = 1000), silent = TRUE)
    if (!inherits(test, "try-error")) d$x[which.max(d$y)] else NA
}
drizopoulos/JMbayes documentation built on Feb. 2, 2021, 12:34 a.m.