h.amise: AMISE for Optimal Bandwidth Selectors

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

Description

The (S3) generic function h.amise evaluates the asymptotic mean integrated squared error AMISE for optimal smoothing parameters h of r'th derivative of kernel density estimator one-dimensional.

Usage

h.amise(x, ...)
## Default S3 method:
h.amise(x, deriv.order = 0, lower = 0.1 * hos, upper = 2 * hos, 
         tol = 0.1 * lower, kernel = c("gaussian", "epanechnikov", "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.amise asymptotic mean integrated squared error implements for choosing the optimal bandwidth h of a r'th derivative kernel density estimator.

We Consider the following AMISE version of the r'th derivative of f the r'th derivative of the kernel estimate (see Scott 1992, pp 131):

AMISE(h;r)= R(K(x;r))/ n h^(2r+1) + 0.25 mu(K(x))^2 h^4 R(f(x;r+2))

The optimal bandwidth minimizing this function is:

h(r) = ( ((2r+1)R(K(x;r)))/(mu(K(x))^2 R(f(x;r+2))) )^(1/(2r+5)) n^-1/(2r+5)

whereof

inf AMISE(h;r) = 2r+5/4 R(K(x;r))^(4/2r+5) ((mu(K(x))^2 R(f(x;r+2)))/(2r+1))^((2r+1)/(2r+5)) n^-4/2r+5

which is the smallest possible AMISE for estimation of f(x;r) using the kernel K(x), where R(K(x;r)) = int K(x;r)^2 dx and mu(K(x)) = int x^2 K(x) 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.
amise	the AMISE value.

Author(s)

Arsalane Chouaib Guidoum [e-mail]

References

Bowman, A. W. and Azzalini, A. (1997). Applied Smoothing Techniques for Data Analysis: the Kernel Approach with S-Plus Illustrations. Oxford University Press, Oxford.

Radhey, S. S. (1987). MISE of kernel estimates of a density and its derivatives. Statistics and Probability Letters, 5, 153–159.

Scott, D. W. (1992). Multivariate Density Estimation. Theory, Practice and Visualization. New York: Wiley.

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

Silverman, B. W. (1986). Density Estimation for Statistics and Data Analysis. Chapman & Hall/CRC. London.

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

See Also

plot.h.amise, see nmise in package sm this function evaluates the mean integrated squared error of a density estimate (deriv.order = 0) which is constructed from data which follow a normal distribution.

Examples

## Derivative order = 0

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

## Derivative order = 1

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

Example output

Call:		Aymptotic Mean Integrated Squared Error

Derivative order = 0
Data: kurtotic (200 obs.);	Kernel: gaussian
AMISE = 0.004774066;	Bandwidth 'h' = 0.7038446


Call:		Aymptotic Mean Integrated Squared Error

Derivative order = 1
Data: kurtotic (200 obs.);	Kernel: gaussian
AMISE = 0.007041654;	Bandwidth 'h' = 0.9717073

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