SCMarker: What the package does (short line)

Description Usage Arguments Value See Also Examples

Estimating peak number of probabiliy density distribution

peak(pp)

`pp`	The output of density function

The value corresponding to each peak

genepeak, ModalFilter

x=c(rnorm(50,4,0.5),rep(0,length=20))
pp=density(x)
peak(pp)
## The function is currently defined as
function (pp) 
{
    y = pp$y
    x = pp$x
    Hpeak = c()
    for (i in 2:(length(y) - 1)) {
        if (y[i] >= y[i - 1] & y[i] > y[i + 1]) {
            Hpeak = rbind(Hpeak, c(x[i], y[i]))
        }
    }
    if (length(Hpeak) == 0) {
        Hpeak = rbind(Hpeak, c(0, 0))
    }
    Hpeak = as.data.frame(Hpeak)
    names(Hpeak) = c("H", "Density")
    return(Hpeak)
  }