*Replication code for Gerber, Gruber, and Hungerman (Forthcoming BJPS).
*Does church attendance cause people to vote? Using blue laws’ repeal to estimate the effect of religiosity on voter turnout

*************
*Table 2*
*************

*The GSS Data, which are used in Table 2, are private use. Here, we provide all non private-use 
*variables in the dataset GSS_Supplement, and we also provide some stata code for using GSS
*data to produce these estimates

***************
*With GSS Data* 
***************
*replace attend =. if attend==9
*gen atd_calc = 0 if attend==0 
*replace atd = 1 if attend == 1
*replace atd = 2 if attend == 2
*replace atd = 6 if attend == 3
*replace atd = 12 if attend == 4
*replace atd = 30 if attend == 5
*replace atd = 40 if attend == 6
*replace atd = 52 if attend == 7
*replace atd = 100 if attend == 8
*drop if age>97
*drop if wrkstat==9
*drop if marital==9
*gen married = marital==1
*label var married "=1 if married"
*gen white = race == 1
*gen black = race == 2
*gen some_col = educ>12&educ<16
*gen high_school = educ==12
*gen drop_out = educ<12
*gen weekly = attend > 6 

*relig is from the GSS also
*Below S* and Y* are state and year dummies
*phat (included in the GSS_Supplement file) is prediced repeal

global rhs "age age2 sex white black some_col high_school drop_out pblack married pforeign riu dispos"

*****************
*TABLE 2 RESULTS*
*****************
reg attend repeal $rhs S* Y* if all==1 & (relig==1|relig==2) & lawchange==0, robust cluster(state)

reg atd_calc repeal $rhs S* Y* if all==1 & (relig==1|relig==2) & lawchange==0, robust cluster(state)

reg attend repeal early $rhs S* Y* if all==1 & (relig==1|relig==2) & lawchange==0, robust cluster(state)

reg atd_calc repeal early $rhs S* Y* if all==1 & (relig==1|relig==2) & lawchange==0, robust cluster(state)

reg attend repeal phat $rhs S* Y* if all==1 & (relig==1|relig==2) & lawchange==0, robust cluster(state)

reg atd_calc repeal phat $rhs S* Y* if all==1 & (relig==1|relig==2) & lawchange==0, robust cluster(state)


*******************
*Other GSS Results*
*******************
*Compare attendance in sample versus elsewhere
quietly reg attend repeal $rhs S* Y* if all==1 & (relig==1|relig==2) & lawchange==0, robust cluster(state)
sum weekly if e(sample)==1
sum weekly if all==0

*Do attendance for everyone and see how big is original sample 
gen either = relig==1|relig==2
sum either
reg attend repeal $rhs S* Y* if all==1 & lawchange==0, robust cluster(state)
reg either repeal $rhs S* Y* if all==1 & lawchange==0, robust cluster(state)
reg either repeal early $rhs S* Y* if all==1 & lawchange==0, robust cluster(state)

*Just court-based cases
gen court2 = repeal
replace court2 = 0 if state!="PA" & state!="VT" &state!="UT"
reg attend court2 repeal $rhs S* Y* if all==1 & (relig==1|relig==2) & lawchange==0, robust cluster(state)


****************
*Voting Results*
****************

clear
use Voting_Data.dta

*********
*TABLE 3*
*********
areg pturn repeal s_elect g_elect pop lpop battle_st battle2 RY* [aweight = pop] if all==1, robust cluster(state) absorb(fips)

areg pturn repeal s_elect g_elect pop lpop battle_st battle2 RY* if all==1, robust cluster(state) absorb(fips)

areg pturn repeal s_elect g_elect pop lpop battle_st battle2 RY* [aweight = pop] if good==1, robust cluster(state) absorb(fips)

areg pturn repeal pearly s_elect g_elect pop lpop battle_st battle2 RY* [aweight = pop] if all==1, robust cluster(state) absorb(fips)

areg pturn repeal phat s_elect g_elect pop lpop battle_st battle2 RY* [aweight = pop] if all==1, robust cluster(state) absorb(fips)

areg lpt repeal s_elect g_elect pop lpop battle_st battle2 RY* if all==1, robust cluster(state) absorb(fips)

areg lpt repeal s_elect g_elect pop lpop battle_st battle2 RY* if all==1 [aweight = pop], robust cluster(state) absorb(fips)

**********************
*Other Voting Results*
**********************
*Judicial repeals
gen court2 = repeal
replace court2 = 0 if state!="Pennsylvania" & state!="Vermont" &state!="Utah"
areg pturn repeal court2 s_elect g_elect pop lpop battle_st battle2 RY* [aweight = pop] if all==1, robust cluster(state) absorb(fips)