# Regression without RD
model.no.disc <- lm(data = d.disc, data_breaches ~ forbes_50)

# Regression with RD
discontinuity.window <- -10:10 # Size of RD Window
model.disc <- lm(data = d.disc[id %in% discontinuity.window,],
	          data_breaches ~ forbes_50)

stargazer(model.no.disc, model.disc, type = "text")