rdd_reg_np {rddtools}

Parametric polynomial estimator of the regression discontinuity

Description

Compute a parametric polynomial regression of the ATE, possibly on the range specified by bandwidth

Usage

rdd_reg_np(
  rdd_object,
  covariates = NULL,
  bw = rdd_bw_ik(rdd_object),
  slope = c("separate", "same"),
  inference = c("np", "lm"),
  covar.opt = list(slope = c("same", "separate"), bw = NULL)
)

Arguments

rdd_object

Object of class rdd_data created by rdd_data

covariates

TODO

bw

A bandwidth to specify the subset on which the parametric regression is estimated

slope

Whether slopes should be different on left or right (separate), or the same.

inference

Type of inference to conduct: non-parametric one (np) or standard (lm). See details.

covar.opt

Options for the inclusion of covariates. Way to include covariates, either in the main regression (include) or as regressors of y in a first step (residual).

Value

An object of class rdd_reg_np and class lm, with specific print and plot methods

References

TODO

See Also

rdd_bw_ik Bandwidth selection using the plug-in bandwidth of Imbens and Kalyanaraman (2012)

Examples

## Step 0: prepare data
data(house)
house_rdd <- rdd_data(y=house$y, x=house$x, cutpoint=0)
## Step 2: regression
# Simple polynomial of order 1:
reg_nonpara <- rdd_reg_np(rdd_object=house_rdd)
print(reg_nonpara)
plot(reg_nonpara)