get.h {rflowcyt} | R Documentation

## Estimate the critical bandwidth for specific number of modes

### Description

get.h finds the critical bandwidth for specific number of modes. That is, it finds the smallest bandwidth for which "m" modes are present for a kernel density estimator.

### Usage

    get.h(x, m = 1, prec = 0.001, hmin = 0, hmax = 1)

### Arguments

`x` | the data vector in which to find the critical bandwidth
`m` | the number of modes for the critical bandwidth
`prec` | the precision for the resulting bandwidth
`hmin` | the minimum value to start searching for the critical bandwidth, h
`hmax` | the maximum value to start searching for the critical bandwidth, h

### Details

get.h uses the Gaussian kernel to estimate the density of a data vector given by x. The bandwidth determines the spread of each data point. Thus a larger bandwidth leads to a smoother density estaimate. get.h finds the smallest bandwidth in which "m" modes are still present.

### Value

`h` | the critical bandwidth, rescaled for the standardized x-values for direct comparison

### Author(s)

Kevin Rader

### References

B.W. Silverman (1981),Using Kernel Density Estimates to Investigate Multimodatlity. J.R. Statist. Soc. B,43,1,97-99.

### See Also

`get.p`, `emp.f`, `get.num.modes`

### Examples

     set.seed(12345)
     x <- c(rnorm(20,0),rnorm(20,3))
     get.h(x)

[Package rflowcyt version 1.8.0 Index]