# Place the following code after the Table 3 code in the Sen 2014 code.
# The following code is based on Sen 2014 code.

d_judges$Trial0 <- ifelse(d_judges$Trial.Years==0,1,0)
d_judges$Bench0 <- ifelse(d_judges$Bench.Years==0,1,0)

aba1 <- zelig(I(District.Court.ABA == "Exceptionally Well Qualified"|District.Court.ABA == "Well Qualified") ~ as.factor(Race.or.Ethnicity) + Female + Republican + Age + Law.Clerk + law.topten + Bench0 + Trial0 + Bench.Years + Trial.Years, data = d_judges, model = "logit")
summary(aba1)

aba2 <- zelig(I(District.Court.ABA == "Exceptionally Well Qualified"|District.Court.ABA == "Well Qualified") ~ as.factor(Race.or.Ethnicity) + Female + Republican + Age + Law.Clerk + law.topten + Bench0 + Trial0 + Bench.Years + Trial.Years + District.Court.Nomination.Year, data = d_judges, model = "logit")
summary(aba2)

aba3 <- zelig(I(District.Court.ABA == "Exceptionally Well Qualified"|District.Court.ABA == "Well Qualified") ~ as.factor(Race.or.Ethnicity) + Female + Republican + Age + Law.Clerk + law.topten + Bench0 + Trial0 + Bench.Years + Trial.Years + District.Court.Nomination.Year*Female + District.Court.Nomination.Year*as.factor(Race.or.Ethnicity), data = d_judges, model = "logit")
summary(aba3)

aba4 <- zelig(I(District.Court.ABA == "Exceptionally Well Qualified"|District.Court.ABA == "Well Qualified") ~ as.factor(Race.or.Ethnicity) + Female + Republican + Age + Law.Clerk + law.topten + Bench0 + Trial0 + Bench.Years + Trial.Years + District.Court.Nomination.Year + ideology, data = d_judges, model = "logit")
summary(aba4)

aba5 <- zelig(I(District.Court.ABA == "Exceptionally Well Qualified"|District.Court.ABA == "Well Qualified") ~ as.factor(Race.or.Ethnicity) + Female + Age + Law.Clerk + law.topten + Bench0 + Trial0 + Bench.Years + Trial.Years + District.Court.Nomination.Year + ideology, data = d_judges, model = "logit")
summary(aba5)