Pre-Class Work coding solution for session 18 (quantile effects)

PCW_18_CS312.R

## Code Cell 1 of 4

# Show median treatment effect estimation here
#install.packages("quantreg")
library(quantreg)
library(Matching)
data(lalonde)

rqfit_r50 <- rq(re78 ~ treat, tau = 0.5, data = lalonde)
summary(rqfit_r50)

## Code Cell 2 of 4

# Show 90% quantile estimation here
rqfit_r90 <- rq(re78 ~ treat, tau = 0.9, data = lalonde)
summary(rqfit_r90)

## Code Cell 3 of 4

Y <- lalonde$re78
Tr <- lalonde$treat
X <- cbind(lalonde$age,lalonde$educ,lalonde$black,lalonde$hisp,lalonde$married,lalonde$nodegr,lalonde$re74,lalonde$re75,lalonde$u74,lalonde$u75)
mout <- Match(Y,Tr,X)

## Code cell 4 of 4

lalonde2 <- rbind(lalonde[mout$index.treated,], lalonde[mout$index.control,])
rqfit_match <- rq(re78 ~ treat, tau = c(0.5,0.9), data = lalonde2)
summary(rqfit_match)

######