library(foreign)
data <-read.dta("pres_replication.dta")   
summary (data)
data.naomit <-data.frame(na.omit(data))
library(MatchIt)
m.out <-matchit(monitored~ln_precinctsize+permobilebox+perusedabsentee+per_invalid,data=data.naomit,method="optimal",ratio=2)

#RUN MODEL

summary(m.out)
plot(m.out,type="hist")

library(Zelig)
#DV - Ilham vote, tobit, unmatched
z.out<-zelig (per_aliyevilham~monitored+ln_precinctsize+permobilebox+perusedabsentee+per_invalid,data=data.naomit,model="tobit", below=0)
summary(z.out)
#DV - Ilham vote, ls, unmatched
z.out<-zelig (per_aliyevilham~monitored+ln_precinctsize+permobilebox+perusedabsentee+per_invalid,data=data.naomit,model="ls")
summary(z.out)
#DV - Ilham vote, tobit, matched
z.out<-zelig (per_aliyevilham~monitored+ln_precinctsize+permobilebox+perusedabsentee+per_invalid,data=match.data(m.out),model="tobit", below=0)
summary(z.out)
#DV - Ilham vote, ls, matched
z.out<-zelig (per_aliyevilham~monitored+ln_precinctsize+permobilebox+perusedabsentee+per_invalid,data=match.data(m.out),model="ls")
summary(z.out)
#DV - turnout, tobit, unmatched
z.out<-zelig (turnout~monitored+ln_precinctsize+permobilebox+perusedabsentee+per_invalid,data=data.naomit,model="tobit", below=0)
summary(z.out)
#DV - turnout, ls, umatched
z.out<-zelig (turnout~monitored+ln_precinctsize+permobilebox+perusedabsentee+per_invalid,data=data.naomit,model="ls")
summary(z.out)
#DV - turnout, tobit, matched
z.out<-zelig (turnout~monitored+ln_precinctsize+permobilebox+perusedabsentee+per_invalid,data=match.data(m.out),model="tobit", below=0)
summary(z.out)
#DV - turnout, ls, matched
z.out<-zelig (turnout~monitored+ln_precinctsize+permobilebox+perusedabsentee+per_invalid,data=match.data(m.out),model="ls")
summary(z.out)


library(foreign)
data1 <-read.dta("referendum_replication.dta")
summary (data1)

data.naomit1 <-data.frame(na.omit(data1))

library(MatchIt)
m.out <-matchit(monitored~ln_precinctsize+per_mobile+perusedabsentee+q21invalid_per,data=data.naomit1,method="optimal",ratio=2)

#RUN MODEL

summary(m.out)
plot(m.out,type="hist")

library(Zelig)

#DV - q21 vote, tobit, unmatched
z.out<-zelig (q21yes_per~monitored+ln_precinctsize+per_mobile+perusedabsentee+q21invalid_per,data=data.naomit1,model="tobit", below=0)
summary(z.out)
#DV - q21 vote, ls, unmatched
z.out<-zelig (q21yes_per~monitored+ln_precinctsize+per_mobile+perusedabsentee+q21invalid_per,data=data.naomit1,model="ls")
summary(z.out)
#DV - q21 vote, tobit, matched
z.out<-zelig (q21yes_per~monitored+ln_precinctsize+per_mobile+perusedabsentee+q21invalid_per,data=match.data(m.out),model="tobit", below=0)
summary(z.out)
#DV - q21 vote, ls, matched
z.out<-zelig (q21yes_per~monitored+ln_precinctsize+per_mobile+perusedabsentee+q21invalid_per,data=match.data(m.out),model="ls")
summary(z.out)
#DV - turnout, tobit, unmatched
z.out<-zelig (turnout~monitored+ln_precinctsize+per_mobile+perusedabsentee+q21invalid_per,data=data.naomit1,model="tobit", below=0)
summary(z.out)
#DV - turnout, ls, unmatched
z.out<-zelig (turnout~monitored+ln_precinctsize+per_mobile+perusedabsentee+q21invalid_per,data=data.naomit1,model="ls")
summary(z.out)
#DV - turnout, tobit, matched
z.out<-zelig (turnout~monitored+ln_precinctsize+per_mobile+perusedabsentee+q21invalid_per,data=match.data(m.out),model="tobit", below=0)
summary(z.out)
#DV - turnout, ls, matched
z.out<-zelig (turnout~monitored+ln_precinctsize+per_mobile+perusedabsentee+q21invalid_per,data=match.data(m.out),model="ls")
summary(z.out)