#control the size of R plots in jupyter
options(repr.plot.width=5, repr.plot.height=5)

# Load R libraries
library(foreign)
library(ggplot2)
library(lpdensity)
library(rddensity)
library(rdrobust)
library(rdlocrand)
library(TeachingDemos)

data = read.dta("CIT_2018_Cambridge_polecon.dta")

Y = data$Y
X = data$X
T = data$T
T_X = T*X

rdplot(Y, X, nbins = c(2500, 500), p = 0, col.lines = "red", col.dots = "lightgray", title = "", 
       x.label = "Islamic Margin of Victory", y.label = "Female High School Share", y.lim = c(0,70))

rdplot(Y[abs(X) <= 50], X[abs(X) <= 50], nbins = c(2500, 500), p = 4, col.lines = "red", 
       col.dots = "lightgray", title = "",  x.label = "Islamic Margin of Victory", 
       y.label = "Female High School Share", y.lim = c(0,70))

rdplot(Y[abs(X) <= 10], X[abs(X) <= 10], nbins = c(2500, 500), p = 4, col.lines = "red", col.dots = "lightgray", title = "", 
       x.label = "Islamic Margin of Victory", y.label = "Female High School Share", y.lim = c(0,70))

rdplot(Y, X, nbins = c(20,20), binselect = 'esmv', x.label = 'Running Variable', y.label = 'Outcome', 
       title = '', y.lim = c(0,25))