rdplot implements several data-driven Regression Discontinuity (RD) plots, using either evenly-spaced or quantile-spaced partitioning. Two type of RD plots are constructed: (i) RD plots with binned sample means tracing out the underlying regression function, and (ii) RD plots with binned sample means mimicking the underlying variability of the data. For technical and methodological details see Calonico, Cattaneo and Titiunik (2015a).
Companion commands are: rdrobust for point estimation and inference procedures, and rdbwselect for data-driven bandwidth selection.
A detailed introduction to this command is given in Calonico, Cattaneo and Titiunik (2015b) and Calonico, Cattaneo, Farrell and Titiunik (2017). A companion Stata package is described in Calonico, Cattaneo and Titiunik (2014).
For more details, and related Stata and R packages useful for analysis of RD designs, visit https://sites.google.com/site/rdpackages/
rdplot(y, x, c = 0, p = 4, nbins = NULL, binselect = "esmv",
scale = NULL, kernel = "uni", weights = NULL, h = NULL,
covs = NULL, support = NULL, subset = NULL, hide = FALSE, ci = NULL,
shade = FALSE, par = NULL, title = NULL,
x.label = NULL, y.label = NULL, x.lim = NULL, y.lim = NULL,
col.dots = NULL, col.lines = NULL, type.dots = NULL,...)is the dependent variable.
is the running variable (a.k.a. score or forcing variable).
specifies the RD cutoff in x; default is c = 0.
specifies the order of the global-polynomial used to approximate the population conditional mean functions for control and treated units; default is p = 4.
specifies the number of bins used to the left of the cutoff, denoted \(J_-\), and to the right of the cutoff, denoted \(J_+\), respectively. If not specified, \(J_+\) and \(J_-\) are estimated using the method and options chosen below.
specifies the procedure to select the number of bins. This option is available only if \(J_-\) and \(J_+\) are not set manually. Options are:
es: IMSE-optimal evenly-spaced method using spacings estimators.
espr: IMSE-optimal evenly-spaced method using polynomial regression.
esmv: mimicking variance evenly-spaced method using spacings estimators. This is the default option.
esmvpr: mimicking variance evenly-spaced method using polynomial regression.
qs: IMSE-optimal quantile-spaced method using spacings estimators.
qspr: IMSE-optimal quantile-spaced method using polynomial regression.
qsmv: mimicking variance quantile-spaced method using spacings estimators.
qsmvpr: mimicking variance quantile-spaced method using polynomial regression.
specifies a multiplicative factor to be used with the optimal numbers of bins selected. Specifically, the number of bins used for the treatment and control groups will be scale\(\times \hat{J}_+\) and scale\(\times \hat{J}_-\), where \(\hat{J}_\cdot\) denotes the estimated optimal numbers of bins originally computed for each group; default is scale = 1.
specifies the kernel function used to construct the local-polynomial estimator(s). Options are: triangular, epanechnikov, and uniform. Default is kernel=uniform (i.e., equal/no weighting to all observations on the support of the kernel).
is the variable used for optional weighting of the estimation procedure. The unit-specific weights multiply the kernel function.
specifies the bandwidth used to construct the (global) polynomial fits given the kernel choice kernel.
If not specified, the bandwidths are chosen to span the full support of the data. If two bandwidths are specified, the first bandwidth is used for the data below the cutoff and the second bandwidth is used for the data above the cutoff.
specifies additional covariates to be used in the polynomial regression.
specifies an optional extended support of the running variable to be used in the construction of the bins; default is the sample range.
an optional vector specifying a subset of observations to be used.
logical. If TRUE, it omits the RD plot; default is hide = FALSE.
optional graphical option to display confidence intervals of selected level for each bin.
optional graphical option to replace confidence intervals with shaded areas.
optional graphical parameters.
optional title for the RD plot.
optional label for the x-axis of the RD plot.
optional label for the y-axis of the RD plot.
optional setting for the range of the x-axis in the RD plot.
optional setting for the range of the y-axis in the RD plot.
optional setting for the color of the dots in the RD plot.
optional setting for the color of the lines in the RD plot.
optional setting for the type of lines in the RD plot.
other graphical parameters.
method used to compute the optimal number of bins.
sample size used to the left of the cutoff.
sample size used to the right of the cutoff.
cutoff value.
selected number of bins to the left of the cutoff.
selected number of bins to the right of the cutoff.
matrix containing the coefficients of the \(p^{th}\) order global polynomial estimated both sides of the cutoff.
data frame containing the variables used to construct the bins: bin id, cutoff values, mean of x and y within each bin, cutoff points and confidence interval bounds.
Calonico, S., M. D. Cattaneo, M. H. Farrell, and R. Titiunik. 2017. rdrobust: Software for Regression Discontinuity Designs. Stata Journal, forthcoming.
Calonico, S., M. D. Cattaneo, and R. Titiunik. 2014. Robust Data-Driven Inference in the Regression-Discontinuity Design. Stata Journal 14(4): 909-946.
Calonico, S., M. D. Cattaneo, and R. Titiunik. 2015a. Optimal Data-Driven Regression Discontinuity Plots. Journal of the American Statistical Association 110(512): 1753-1769.
Calonico, S., M. D. Cattaneo, and R. Titiunik. 2015b. rdrobust: An R Package for Robust Nonparametric Inference in Regression-Discontinuity Designs. R Journal 7(1): 38-51.
Cattaneo, M. D., B. Frandsen, and R. Titiunik. 2015. Randomization Inference in the Regression Discontinuity Design: An Application to the Study of Party Advantages in the U.S. Senate. Journal of Causal Inference 3(1): 1-24.
x<-runif(1000,-1,1)
y<-5+3*x+2*(x>=0)+rnorm(1000)
rdplot(y,x)
Run the code above in your browser using DataLab