* Data are from the TESS project: "Backlash Effect? White Americans' Response to the Coronavirus Pandemic"
* PI for the TESS project: LaFleur Stephens-Dougan
* Data available here: https://osf.io/vsfec/. I'm using the "8041.063_TESS 063_Stephens-Dougan_TESS.dta" file.
* Code below by L.J Zigerell

* Treatment

codebook P_EXP
recode   P_EXP (1=0) (2=1), gen(Treatment)
tab      P_EXP Treatment, nol mi

* Outcomes mentioned in the Stephens-Dougan pre-analysis plan [https://osf.io/8bsn4]

codebook Q15
clonevar slow = Q15
recode   slow (98=.)
tab      Q15 slow, mi

recode   Q15 (4/5=1) (1/3 98=0), gen(slowDIS)
tab      Q15 slowDIS, mi

codebook Q18
clonevar home = Q18
recode   home (98=.)
tab      Q18 home, mi

recode Q18 (1/3 98=0) (4/5=1), gen(homeNOT)
tab    Q18 homeNOT, mi

codebook Q16
clonevar rights = Q16
recode   rights (98=.)
tab      Q16 rights, mi

recode Q16 (1/2=1) (3/98=0), gen(rightsAG)
tab    Q16 rightsAG, mi

codebook Q21
clonevar loosen = Q21
recode   loosen (98=.)
tab      Q21 loosen, mi

* Racial resentment

codebook Q1 Q2 Q3 Q4
recode   Q1 Q2 Q3 Q4 (98=.)
gen      RR = -Q1 + Q2 + Q3 - Q4
tab      RR

factor  Q1 Q2 Q3 Q4, pcf
predict RRfactor
pwcorr  RR RRfactor

xtile RR3 = RR, nq(3)
tab   RR3
xtile RR2 = RR, nq(2)
tab   RR2

* Negative stereotypes about Blacks

codebook Q5 Q6 Q7 Q8
recode   Q5 Q6 Q7 Q8 (98=.)
sum      Q5 Q6 Q7 Q8
gen      STintel = Q8 - Q7
gen      STwork  = Q6 - Q5
recode   STintel (-7/-1=1) (0=2) (1/7=3), gen(STintel3)
recode   STwork  (-7/-1=1) (0=2) (1/7=3), gen(STwork3)
tab      STintel, mi
tab      STwork , mi

tab Q7 Q8 if STintel3==3, mi
tab Q5 Q6 if STwork3 ==3, mi

tab Q7 Q8 if STintel3!=3, mi
tab Q5 Q6 if STwork3 !=3, mi

* Blacks' poor health choices

codebook Q9C
clonevar poorCH = Q9C
recode   poorCH (98=.)
tab      Q9C poorCH, mi

xtile poorCH3 = poorCH, nq(3)
xtile poorCH2 = poorCH, nq(2)
tab   poorCH3
tab   poorCH2

tab Q9C poorCH3, mi
tab Q9C poorCH2, mi // Same "high" group as CH3

* Check codings against party identification

tab    P_PARTYID7
pwcorr P_PARTYID7 RR RR2 RR3 poorCH3 poorCH2 STintel3 STwork3

* Statistical tests [unweighted]

reg slow     Treatment if RR3==3, robust
reg slowDIS  Treatment if RR3==3, robust
reg home     Treatment if RR3==3, robust
reg homeNOT  Treatment if RR3==3, robust
reg rights   Treatment if RR3==3, robust
reg rightsAG Treatment if RR3==3, robust
reg loosen   Treatment if RR3==3, robust

reg slow     Treatment if RR2==2, robust
reg slowDIS  Treatment if RR2==2, robust
reg home     Treatment if RR2==2, robust
reg homeNOT  Treatment if RR2==2, robust
reg rights   Treatment if RR2==2, robust
reg rightsAG Treatment if RR2==2, robust
reg loosen   Treatment if RR2==2, robust

reg slow     Treatment if STintel3==3 | STwork3==3, robust
reg slowDIS  Treatment if STintel3==3 | STwork3==3, robust
reg home     Treatment if STintel3==3 | STwork3==3, robust
reg homeNOT  Treatment if STintel3==3 | STwork3==3, robust
reg rights   Treatment if STintel3==3 | STwork3==3, robust
reg rightsAG Treatment if STintel3==3 | STwork3==3, robust
reg loosen   Treatment if STintel3==3 | STwork3==3, robust

reg slow     Treatment if STintel3==3 & STwork3==3, robust
reg slowDIS  Treatment if STintel3==3 & STwork3==3, robust
reg home     Treatment if STintel3==3 & STwork3==3, robust
reg homeNOT  Treatment if STintel3==3 & STwork3==3, robust
reg rights   Treatment if STintel3==3 & STwork3==3, robust
reg rightsAG Treatment if STintel3==3 & STwork3==3, robust
reg loosen   Treatment if STintel3==3 & STwork3==3, robust

reg slow     Treatment if poorCH3==3, robust
reg slowDIS  Treatment if poorCH3==3, robust
reg home     Treatment if poorCH3==3, robust
reg homeNOT  Treatment if poorCH3==3, robust
reg rights   Treatment if poorCH3==3, robust
reg rightsAG Treatment if poorCH3==3, robust
reg loosen   Treatment if poorCH3==3, robust

* Statistical tests [weighted]

reg slow     Treatment if RR3==3 [pw=WEIGHT]
reg slowDIS  Treatment if RR3==3 [pw=WEIGHT]
reg home     Treatment if RR3==3 [pw=WEIGHT]
reg homeNOT  Treatment if RR3==3 [pw=WEIGHT]
reg rights   Treatment if RR3==3 [pw=WEIGHT]
reg rightsAG Treatment if RR3==3 [pw=WEIGHT]
reg loosen   Treatment if RR3==3 [pw=WEIGHT]

reg slow     Treatment if RR2==2 [pw=WEIGHT]
reg slowDIS  Treatment if RR2==2 [pw=WEIGHT]
reg home     Treatment if RR2==2 [pw=WEIGHT]
reg homeNOT  Treatment if RR2==2 [pw=WEIGHT]
reg rights   Treatment if RR2==2 [pw=WEIGHT]
reg rightsAG Treatment if RR2==2 [pw=WEIGHT]
reg loosen   Treatment if RR2==2 [pw=WEIGHT]

reg slow     Treatment if STintel3==3 | STwork3==3 [pw=WEIGHT]
reg slowDIS  Treatment if STintel3==3 | STwork3==3 [pw=WEIGHT]
reg home     Treatment if STintel3==3 | STwork3==3 [pw=WEIGHT]
reg homeNOT  Treatment if STintel3==3 | STwork3==3 [pw=WEIGHT]
reg rights   Treatment if STintel3==3 | STwork3==3 [pw=WEIGHT]
reg rightsAG Treatment if STintel3==3 | STwork3==3 [pw=WEIGHT]
reg loosen   Treatment if STintel3==3 | STwork3==3 [pw=WEIGHT]

reg slow     Treatment if STintel3==3 & STwork3==3 [pw=WEIGHT]
reg slowDIS  Treatment if STintel3==3 & STwork3==3 [pw=WEIGHT]
reg home     Treatment if STintel3==3 & STwork3==3 [pw=WEIGHT]
reg homeNOT  Treatment if STintel3==3 & STwork3==3 [pw=WEIGHT]
reg rights   Treatment if STintel3==3 & STwork3==3 [pw=WEIGHT]
reg rightsAG Treatment if STintel3==3 & STwork3==3 [pw=WEIGHT]
reg loosen   Treatment if STintel3==3 & STwork3==3 [pw=WEIGHT]

reg slow     Treatment if poorCH3==3 [pw=WEIGHT]
reg slowDIS  Treatment if poorCH3==3 [pw=WEIGHT]
reg home     Treatment if poorCH3==3 [pw=WEIGHT]
reg homeNOT  Treatment if poorCH3==3 [pw=WEIGHT]
reg rights   Treatment if poorCH3==3 [pw=WEIGHT]
reg rightsAG Treatment if poorCH3==3 [pw=WEIGHT]
reg loosen   Treatment if poorCH3==3 [pw=WEIGHT]

*

mean rightsAG if Treatment==1 & (STintel3==3 & STwork3==3) [pw=WEIGHT], level(84)
mean rightsAG if Treatment==0 & (STintel3==3 & STwork3==3) [pw=WEIGHT], level(84)

mean rightsAG if Treatment==1 & (STintel3==3 & STwork3==3)            , level(84)
mean rightsAG if Treatment==0 & (STintel3==3 & STwork3==3)            , level(84)

mean rightsAG if Treatment==1 & (STintel3==3 | STwork3==3) [pw=WEIGHT], level(84)
mean rightsAG if Treatment==0 & (STintel3==3 | STwork3==3) [pw=WEIGHT], level(84)

mean rightsAG if Treatment==1 & (STintel3==3 | STwork3==3)            , level(84)
mean rightsAG if Treatment==0 & (STintel3==3 | STwork3==3)            , level(84)

*

gen     STscale = 0
replace STscale = 1 if STintel3==3 | STwork3==3
replace STscale = 2 if STintel3==3 & STwork3==3
prop    STscale
prop    STscale [pw=WEIGHT]