bkfe(x, drv, bandwidth, gridsize=401, range.x=range(x), binned=F, truncate=F)
x
| vector of observations from the distribution whose density is to be estimated. Missing values are not allowed. |
drv
| order of derivative in the density functional. Must be a non-negative even integer. |
bandwidth
| the kernel bandwidth smoothing parameter. |
gridsize
| the number of equally-spaced points over which binning is performed. |
range.x
|
vector containing the minimum and maximum values of x
at which to compute the estimate.
The default is the minimum and maximum data values.
|
binned
|
logical flag: if TRUE , then x and y are taken to be grid counts
rather than raw data.
|
truncate
|
logical flag: if TRUE , data with x values outside the
range specified by range.x are ignored.
|
drv
is the integral of the
product of the density and its drv
th derivative.
The kernel estimates
of such quantities are computed using a binned implementation,
and the kernel is the standard normal density.Wand, M. P. and Jones, M. C. (1995). Kernel Smoothing. Chapman and Hall, London.
data(geyser) x <- geyser$duration est <- bkfe(x,drv=4,bandwidth=0.3)