genepeak: Estimating modal distribution

In Fang0828/SCMarker: What the package does (short line)

Description

Estimating the number of peak corresponding to non-zero value from probability density distribution

Usage

genepeak(x, width = 1)

Arguments

x	A vector
width	The numeric width which is used to estimate probability density distribution function

See Also

peak

Examples

x=c(rnorm(50,4,0.5),rep(0,length=20))
genepeak(x)
genepeak(x,width=2)

## The function is currently defined as
function (x, width = 1) 
{
    if (length(x[x != 0]) > 1) {
        pp = density(x[x != 0], width = width)
        Hpeak = peak(pp)
        return(dim(Hpeak)[1])
    }
    else {
        return(0)
    }
  }

Fang0828/SCMarker documentation built on May 13, 2019, 12:51 p.m.