Kernel Density and Distribution Estimation
To compute the non-parametric kernel estimation of the probability density function (PDF) and cumulative distribution function (CDF).
kde(x, h, xgrid, ngrid, kernel = c("epan", "unif", "tria", "quar",
"triw", "tric", "gaus", "cos"), plot = FALSE)x |
a numeric vector of data values. |
h |
the smoothing bandwidth. See 'Details' of the default bandwidth. |
xgrid |
the user-defined data points at which the PDF and CDF are to be evaluated. The default is the data values |
ngrid |
the number of equally spaced points at which the PDF and CDF are to be evaluated. The default is |
kernel |
a character string which determines the smoothing kernel function. This must be one of |
plot |
a logical indicating whether to plot the estimated PDF and CDF graphs. |
Kernel density and distribution estimation is a non-parametric method to estimate the probability density function (PDF) and cumulative distribution function (CDF) by using kernel function for a continuous random variable. The default smoothing bandwidth is the plug-in optimal one in Fan and Gijbels (1996), i.e., h = c*n^(-1/5), where the constant is replaced by (8*pi/3)^(1/5)*2.0362*(((quantile(x, 0.75) - quantile(x, 0.25))/1.349)^(2/3)) in this function. Missing values have been removed.
x |
the original data values. |
xgrid |
the points where the PDF and CDF are to be evaluated. |
fhat |
the estimated PDF values at the specified points. |
Fhat |
the estimated CDF values at the specified points. |
bw |
the smoothing bandwidth used. |
The smoothing bandwidth is always a critical issue in non-parametric statistics. The default smoothing bandwidth suggested by Fan and Gijbels (1996) may not perform the best in some cases. You are recommended to provide one obtained by other methods.
Debin Qiu <debinqiu@uga.edu>
Fan, I. Gijbels (1996). Local Polynomial Modeling and its Applications. Chapman & Hall, London. pp. 47.
Wand, M. P. and Jones, M. C. (1995). Kernel Smoothing. Chapman and Hall, London.
x <- rnorm(200,2,3) # with default bandwidth kde(x, kernel = "quar", plot = TRUE) # with specified bandwidth kde(x, h = 4, kernel = "quar", plot = TRUE)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.