Replication material for "Limitless Regression Discontinuity"
To reproduce the results in the paper:
- Make sure R is up-to-date with required packages installed:
- We used R version 3.4.4
- Required packages to carry out analyses:
robustbaserddsandwichnnetRItoolspbs
- Required packages to display results (paper, analysis documents):
knitrkableExtraggplotgridExtra
- Optional packages:
parallelto run simulations in parallelrmarkdownto make analysis documents more easily
- In
R, Run data analysis for Maria Study:- Dataset available from private github repo https://github.com/alexisrsantos/Deaths_Puerto_Rico_JAMA
- Store data as
extdata/death_file_JAMA.csv - Then run to create analysis file
maria.md
library(knitr)
knit('maria.Rmd')
- In
R, Run data analysis for AP Study: (this automatically downloads publicly available data)
knit('dataAnalysis.Rmd')
- Creates files:
dataAnalysis.mdanalysis documentRDanalysis-[M][D][Time].RDataanalysis outputtab-results.textab-alt.textables 1 & 2 for paper
- In
RRun simulation studies
nrep <- 5000 ## number of simulation replications
knit('fullOutcomeSim.Rmd')
- Substitute
nrep <- XforXsimulation replications nrep <- 0loads previously-run simulation files (files corresponding to results in the paper are in the repo)- Creates files:
fullOutcomeSim.mdanalysis documentoutcomeSim-YYYY-MM-DD.RDatasimulation results for table 3totalPolySim-YYYY-MM-DD.RDatasimulation results for table 4tab-levelSimulation.texlatex for table 3tab-polynomialSimulation.texlatex for table 4- some other random files.