h.ucv: Unbiased (Least-Squares) Cross-Validation for Bandwidth...

In kedd: Kernel Estimator and Bandwidth Selection for Density and Its Derivatives

Description

The (S3) generic function h.ucv computes the unbiased (least-squares) cross-validation bandwidth selector of r'th derivative of kernel density estimator one-dimensional.

Usage

h.ucv(x, ...)
## Default S3 method:
h.ucv(x, deriv.order = 0, lower = 0.1 * hos, upper = 2 * hos, 
         tol = 0.1 * lower, kernel = c("gaussian", "epanechnikov", "uniform", 
         "triangular", "triweight", "tricube", "biweight", "cosine"), ...)

Arguments

x	vector of data values.
deriv.order	derivative order (scalar).
lower, upper	range over which to minimize. The default is almost always satisfactory. hos (Over-smoothing) is calculated internally from an kernel, see details.
tol	the convergence tolerance for optimize.
kernel	a character string giving the smoothing kernel to be used, with default "gaussian".
...	further arguments for (non-default) methods.

Details

h.ucv unbiased (least-squares) cross-validation implements for choosing the bandwidth h of a r'th derivative kernel density estimator.

Rudemo (1982) and Bowman (1984) proposed a so-called unbiased (least-squares) cross-validation (UCV) in kernel density estimator. An adaptation of unbiased cross-validation is proposed by Wolfgang et al. (1990) for bandwidth choice in the r'th derivative of kernel density estimator. The essential idea of this methods, for the estimation of f(x;r) (r is derivative order), is to use the bandwidth h which minimizes the function:

UCV(h;r)= int (hat(f)(x;r))^2 - 2 n^(-1) (-1)^r sum(hat(f)(X_i;2r),i=1,...n)

The bandwidth minimizing this function is:

h(r) = argmin UCV(h;r)

for r = 0, 1, 2, …
where

int (hat(f)(x;r))^2 dx = R(k(x;r))/n h^(2r+1) + (-1)^r / n (n-1) h^(2r+1) sum(sum(K(.;r)*K(.;r)(x(j)-x(i)/h)), i=1...n, j=1...n, j != i)

and K(x;r)*K(x;r) is the convolution of the r'th derivative kernel function K(x;r) (see kernel.conv and kernel.fun).
The estimate hat(f)(x;2r) on the subset (X_j)_(j != i) denoting the leave-one-out estimator, can be written:

hat(f)(X_i;2r) = 1 /(n-1) h^(2r+1) sum(K(x(j)-x(i)/h;2r), j != i)

The function UCV(h;r) is unbiased cross-validation in the sense that E[UCV(h;r)] = MISE[hat(f)(x;r)]-R(f(x;r)) (see, Scott and George 1987). Can be simplified to give the computationally:

UCV(h;r) = R(k(x;r))/n h^(2r+1) + (-1)^r / n (n-1) h^(2r+1) sum(sum(K(.;r)*K(.;r) - 2 K(.;r))(x(j)-x(i)/h), i=1...n, j=1...n, j != i)

where R(K(x;r)) = int K(x;r)^2 dx.

The range over which to minimize is hos Oversmoothing bandwidth, the default is almost always satisfactory. See George and Scott (1985), George (1990), Scott (1992, pp 165), Wand and Jones (1995, pp 61).

Value

x	data points - same as input.
data.name	the deparsed name of the x argument.
n	the sample size after elimination of missing values.
kernel	name of kernel to use
deriv.order	the derivative order to use.
h	value of bandwidth parameter.
min.ucv	the minimal UCV value.

Author(s)

Arsalane Chouaib Guidoum [e-mail]

References

Bowman, A. (1984). An alternative method of cross-validation for the smoothing of kernel density estimates. Biometrika, 71, 353–360.

Jones, M. C. and Kappenman, R. F. (1991). On a class of kernel density estimate bandwidth selectors. Scandinavian Journal of Statistics, 19, 337–349.

Jones, M. C., Marron, J. S. and Sheather,S. J. (1996). A brief survey of bandwidth selection for density estimation. Journal of the American Statistical Association, 91, 401–407.

Peter, H. and Marron, J.S. (1987). Estimation of integrated squared density derivatives. Statistics and Probability Letters, 6, 109–115.

Rudemo, M. (1982). Empirical choice of histograms and kernel density estimators. Scandinavian Journal of Statistics, 9, 65–78.

Scott, D.W. and George, R. T. (1987). Biased and unbiased cross-validation in density estimation. Journal of the American Statistical Association, 82, 1131–1146.

Sheather, S. J. (2004). Density estimation. Statistical Science, 19, 588–597.

Tarn, D. (2007). ks: Kernel density estimation and kernel discriminant analysis for multivariate data in R. Journal of Statistical Software, 21(7), 1–16.

Wand, M. P. and Jones, M. C. (1995). Kernel Smoothing. Chapman and Hall, London.

Wolfgang, H. (1991). Smoothing Techniques, With Implementation in S. Springer-Verlag, New York.

Wolfgang, H., Marron, J. S. and Wand, M. P. (1990). Bandwidth choice for density derivatives. Journal of the Royal Statistical Society, Series B, 223–232.

See Also

plot.h.ucv, see bw.ucv in package "stats" and ucv in package MASS for Gaussian kernel only if deriv.order = 0, hlscv in package ks for Gaussian kernel only if 0 <= deriv.order <= 5, kdeb in package locfit if deriv.order = 0.

Examples

## Derivative order = 0

h.ucv(kurtotic,deriv.order = 0)

## Derivative order = 1

h.ucv(kurtotic,deriv.order = 1)

Example output

Call:		Unbiased Cross-Validation

Derivative order = 0
Data: kurtotic (200 obs.);	Kernel: gaussian
Min UCV = -0.5577767;	Bandwidth 'h' = 0.05255485


Call:		Unbiased Cross-Validation

Derivative order = 1
Data: kurtotic (200 obs.);	Kernel: gaussian
Min UCV = -19.6538;	Bandwidth 'h' = 0.05160448

kedd documentation built on May 2, 2019, 7:32 a.m.