Nothing
## ---- echo = FALSE, message = FALSE-------------------------------------------
knitr::opts_chunk$set(collapse = T, comment = "#>")
## -----------------------------------------------------------------------------
library(rddtools)
data(house)
## -----------------------------------------------------------------------------
house_rdd <- rdd_data(y=house$y, x=house$x, cutpoint=0)
## ----dataPlot-----------------------------------------------------------------
summary(house_rdd)
plot(house_rdd)
## ----reg_para-----------------------------------------------------------------
reg_para <- rdd_reg_lm(rdd_object=house_rdd, order=4)
reg_para
plot(reg_para)
## ----RegPlot------------------------------------------------------------------
bw_ik <- rdd_bw_ik(house_rdd)
reg_nonpara <- rdd_reg_np(rdd_object=house_rdd, bw=bw_ik)
print(reg_nonpara)
## ----SensiPlot----------------------------------------------------------------
plotSensi(reg_nonpara, from=0.05, to=1, by=0.1)
## ----placeboPlot--------------------------------------------------------------
plotPlacebo(reg_nonpara)
## ----DensPlot-----------------------------------------------------------------
dens_test(reg_nonpara)
## -----------------------------------------------------------------------------
set.seed(123)
n_Lee <- nrow(house)
Z <- data.frame(z1 = rnorm(n_Lee, sd=2),
z2 = rnorm(n_Lee, mean = ifelse(house<0, 5, 8)),
z3 = sample(letters, size = n_Lee, replace = TRUE))
house_rdd_Z <- rdd_data(y = house$y, x = house$x, covar = Z, cutpoint = 0)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.