Description Usage Arguments Details Value Note See Also Examples
Calculate estimators and one and twosided CIs based on local polynomial
estimator in RD under secondorder Taylor or Hölder smoothness class. If
kern="optimal"
, calculate optimal estimators under secondorder Taylor
smoothness class.
1 2 3 4 
formula 
object of class 
data 
optional data frame, list or environment (or object coercible by

subset 
optional vector specifying a subset of observations to be used in the fitting process. 
cutoff 
specifies the RD cutoff in the running variable. 
M 
Bound on second derivative of the conditional mean function. 
kern 
specifies kernel function used in the local regression. It can
either be a string equal to 
na.action 
function which indicates what should happen when the data
contain 
opt.criterion 
Optimality criterion that bandwidth is designed to optimize. It can either be based on exact finitesample maximum bias and finitesample estimate of variance, or asymptotic approximations to the bias and variance. The options are:
The finitesample methods use conditional variance given by

bw.equal 
logical specifying whether bandwidths on either side of cutoff should be constrainted to equal to each other. 
hp, hm 
bandwidth for treated (units with positive running variable),
and control (units with negative running variable) units. If 
se.method 
Vector with methods for estimating standard error of
estimate. If

alpha 
determines confidence level, 
beta 
Determines quantile of excess length to optimize, if bandwidth optimizes given quantile of excess length of onesided confidence intervals. 
J 
Number of nearest neighbors, if "nn" is specified in

sclass 
Smoothness class, either 
order 
Order of local regression 1 for linear, 2 for quadratic. 
se.initial 
Method for estimating initial variance for computing optimal bandwidth. Ignored if data already contains estimate of variance.

The bandwidth is calculated to be optimal for a given performance criterion,
as specified by opt.criterion
. For local polynomial estimators, this
optimal bandwidth is calculated using the function RDOptBW
.
Alternatively, for local polynomial estimators, the bandwidths above and
below the cutoff can be specified by hp
and hm
.
Returns an object of class "RDResults"
. The function
print
can be used to obtain and print a summary of the results. An
object of class "RDResults"
is a list containing the following
components
estimate
Point estimate. This estimate is MSEoptimal if
opt.criterion="MSE"
Least favorable function, only relevant for optimal estimator under Taylor class.
maxbias
Maximum bias of estimate
sd
Standard deviation of estimate
lower
, upper
Lower (upper) endpoint of a onesided CI
based on estimate
. This CI is optimal if
opt.criterion=="OCI"
hl
Halflength of a twosided CI based on estimate
, so
that the CI is given by c(estimatehl, estimate+hl)
. The
CI is optimal if opt.criterion="FLCI"
eff.obs
Effective number of observations used by
estimate
hp
, hm
Bandwidths used
naive
Coverage of CI that ignores bias and uses
qnorm(1alpha/2)
as critical value
call
the matched call
subset
is evaluated in the same way as variables in formula
,
that is first in data
and then in the environment of formula
.
1 2 3 
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.