lpbwselect: Bandwidth Selection Procedures for Local Polynomial...

In nprobust: Kernel Density and Local Polynomial Regression Methods

Bandwidth Selection Procedures for Local Polynomial Regression Estimation and Inference

Description

lpbwselect implements bandwidth selectors for local polynomial regression point estimators and inference procedures developed in Calonico, Cattaneo and Farrell (2018). See also Calonico, Cattaneo and Farrell (2022) for related optimality results. It also implements other bandwidth selectors available in the literature.

Companion commands: lprobust for local polynomial point estimation and inference procedures.

A detailed introduction to this command is given in Calonico, Cattaneo and Farrell (2019). For more details and related software useful for empirical analysis, visit https://nppackages.github.io/.

Usage

lpbwselect(y, x, eval = NULL, neval = NULL, p = NULL, deriv = NULL,
kernel = "epa", bwselect = "mse-dpi", bwcheck = 21, bwregul = 1,
imsegrid = 30, vce = "nn", cluster = NULL,
nnmatch = 3, interior = FALSE, subset = NULL,
weights = NULL, masspoints = "check", data = NULL)

Arguments

y	dependent variable.
x	independent variable.
eval	vector of evaluation point(s). By default it uses 30 equally spaced points over the support of x.
neval	number of equally spaced evaluation points on the support of x. Default is neval=30.
p	polynomial order used to construct point estimator; default is p = 1 (local linear regression).
deriv	derivative order of the regression function to be estimated. Default is deriv=0 (regression function).
kernel	kernel function used to construct local polynomial estimators. Options are epa for the epanechnikov kernel, tri for the triangular kernel, uni for the uniform kernel and gau for the gaussian kernel. Default is kernel = epa.
bwselect	bandwidth selection procedure to be used. Options are: mse-dpi second-generation DPI implementation of MSE-optimal bandwidth. Default option. mse-rot ROT implementation of MSE-optimal bandwidth. imse-dpi second-generation DPI implementation of IMSE-optimal bandwidth (computed using grid of evaluation points selected). imse-rot ROT implementation of IMSE-optimal bandwidth (computed using grid of evaluation points selected). ce-dpi second generation DPI implementation of CE-optimal bandwidth. ce-rot ROT implementation of CE-optimal bandwidth. all reports all available bandwidth selection procedures. Note: MSE = Mean Square Error; IMSE = Integrated Mean Squared Error; CE = Coverage Error; DPI = Direct Plug-in; ROT = Rule-of-Thumb. For details on implementation see Calonico, Cattaneo and Farrell (2019).
bwcheck	if a positive integer is provided, then the selected bandwidth is enlarged so that at least bwcheck effective observations are available at each evaluation point. Default is bwcheck = 21.
bwregul	specifies scaling factor for the regularization term added to the denominator of bandwidth selectors. Setting bwregul = 0 removes the regularization term from the bandwidth selectors. Default is bwregul = 1.
imsegrid	number of evaluations points used to compute the IMSE bandwidth selector. Default is imsegrid = 30.
vce	procedure used to compute the variance-covariance matrix estimator. Options are: nn heteroskedasticity-robust nearest neighbor variance estimator with nnmatch the (minimum) number of neighbors to be used. Default choice (when cluster is not supplied). hc0 heteroskedasticity-robust plug-in residuals variance estimator without weights. hc1 heteroskedasticity-robust plug-in residuals variance estimator with hc1 weights. hc2 heteroskedasticity-robust plug-in residuals variance estimator with hc2 weights. hc3 heteroskedasticity-robust plug-in residuals variance estimator with hc3 weights. cr1 cluster-robust variance estimator with ((n-1)/(n-k)) * (G/(G-1)) multiplier (Stata-style CR1). Default choice when cluster is supplied. Requires cluster. cr2 cluster-robust variance estimator with Bell-McCaffrey block-adjusted residuals. Requires cluster. cr3 cluster-robust variance estimator with block jackknife-style adjustment and (G-1)/G multiplier. Requires cluster. When cluster is supplied, vce=hc0 or vce=hc1 is remapped to cr1, vce=hc2 to cr2, and vce=hc3 to cr3 (with a warning). vce=nn with cluster silently defaults to cr1.
cluster	indicates the cluster ID variable used for cluster-robust variance estimation. When supplied, the default vce becomes cr1; cr2 and cr3 are also available.
nnmatch	to be combined with for vce=nn for heteroskedasticity-robust nearest neighbor variance estimator with nnmatch indicating the minimum number of neighbors to be used. Default is nnmatch=3

.

interior	if TRUE, all evaluation points are assumed to be interior points. This option affects only data-driven bandwidth selection via lpbwselect. Default is interior = FALSE.
subset	optional rule specifying a subset of observations to be used.
weights	optional vector of non-negative observation weights (multiplicative with kernel weights in all bandwidth-selection steps).
masspoints	how to handle evaluation points with few unique x values within bandwidth. Options: "check" (default), "off".
data	an optional data frame. When supplied, y, x, cluster, weights, and subset may be given as bare variable names referring to columns of data.

Value

bws	A matrix containing eval (grid points), h and b (selected bandwidths) for each evaluation point. When bwselect="all", the matrix is widened to include both MSE and CER bandwidths from each of the DPI and ROT estimators.
bws.imse	IMSE-optimal bandwidths returned when bwselect="all"; NULL otherwise.
opt	A list containing options passed to the function.

Author(s)

Sebastian Calonico, University of California, Davis, CA. [e-mail].

Matias D. Cattaneo, Princeton University, Princeton, NJ. [e-mail].

Max H. Farrell, University of California, Santa Barbara, CA. [e-mail].

References

Calonico, S., M. D. Cattaneo, and M. H. Farrell. 2018. On the Effect of Bias Estimation on Coverage Accuracy in Nonparametric Inference. Journal of the American Statistical Association, 113(522): 767-779. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1080/01621459.2017.1285776")}.

Calonico, S., M. D. Cattaneo, and M. H. Farrell. 2019. nprobust: Nonparametric Kernel-Based Estimation and Robust Bias-Corrected Inference. Journal of Statistical Software, 91(8): 1-33. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.18637/jss.v091.i08")}.

Calonico, S., M. D. Cattaneo, and M. H. Farrell. 2022. Coverage Error Optimal Confidence Intervals for Local Polynomial Regression. Bernoulli, 28(4): 2998-3022.

See Also

lprobust

Examples

x   <- runif(500)
y   <- sin(4*x) + rnorm(500)
est <- lpbwselect(y,x)
summary(est)

nprobust documentation built on May 19, 2026, 9:07 a.m.