*******************************************************************************************
** Stata do file for: ANES 2022 Pilot Study
*******************************************************************************************

*******************************************************************************************
** Survey weights
*******************************************************************************************

sum        weight  // N = 1,500
svyset [pw=weight]

*******************************************************************************************
** Racial groups
*******************************************************************************************

codebook eth rwh rbl rain ras rpi roth

gen     RwhiteNH = 0
replace RwhiteNH = 1 if eth==2 & rwh==1 & rbl!=1 & rain!=1 & ras!=1 & rpi!=1 & roth!=1
gen     RblackNH = 0
replace RblackNH = 1 if eth==2 & rwh!=1 & rbl==1 & rain!=1 & ras!=1 & rpi!=1 & roth!=1
gen     RasianNH = 0
replace RasianNH = 1 if eth==2 & rwh!=1 & rbl!=1 & rain!=1 & ras==1 & rpi!=1 & roth!=1

tab1 RwhiteNH RblackNH RasianNH

*******************************************************************************************
** Intelligence comparisons
*******************************************************************************************

codebook swhintel sblintel sasintel
tab1     swhintel sblintel sasintel, mi

gen intelW = 7 - swhintel if swhintel>=1 & swhintel<=7
gen intelB = 7 - sblintel if sblintel>=1 & sblintel<=7
gen intelA = 7 - sasintel if sasintel>=1 & sasintel<=7

tab1 intelW intelB intelA, mi

gen     intelWB = 0 
replace intelWB = 1 if (intelW > intelB) & intelW!=. & intelB!=.

gen     intelAB = 0 
replace intelAB = 1 if (intelA > intelB) & intelA!=. & intelB!=.

gen     intelAW = 0 
replace intelAW = 1 if (intelA > intelW) & intelA!=. & intelW!=.

tab1 intelWB intelAB intelAW, mi

*******************************************************************************************
** Other variables
*******************************************************************************************

codebook gdiffbiol

codebook ideo5

*******************************************************************************************
** Percentages
*******************************************************************************************

svy: prop intelWB                , level(83.4)
svy: prop         intelAW        , level(83.4)
svy: prop                 intelAB, level(83.4)

svy: prop intelWB if gdiffbiol>=1 & gdiffbiol<=4, level(83.4)
svy: prop intelAW if gdiffbiol>=1 & gdiffbiol<=4, level(83.4)
svy: prop intelAB if gdiffbiol>=1 & gdiffbiol<=4, level(83.4)

svy: prop intelWB if gdiffbiol==5, level(83.4)
svy: prop intelAW if gdiffbiol==5, level(83.4)
svy: prop intelAB if gdiffbiol==5, level(83.4)

svy, subpop(RwhiteNH): prop intelWB, level(83.4)
svy, subpop(RwhiteNH): prop intelAW, level(83.4)
svy, subpop(RwhiteNH): prop intelAB, level(83.4)

svy, subpop(RblackNH): prop intelWB, level(83.4)
svy, subpop(RblackNH): prop intelAW, level(83.4)
svy, subpop(RblackNH): prop intelAB, level(83.4)

svy, subpop(RasianNH): prop intelWB, level(83.4)
svy, subpop(RasianNH): prop intelAW, level(83.4)
svy, subpop(RasianNH): prop intelAB, level(83.4)

svy, subpop(RwhiteNH): prop intelWB if ideo5==1, level(83.4)
svy, subpop(RwhiteNH): prop intelAW if ideo5==1, level(83.4)
svy, subpop(RwhiteNH): prop intelAB if ideo5==1, level(83.4)

svy, subpop(RwhiteNH): prop intelWB if ideo5==2, level(83.4)
svy, subpop(RwhiteNH): prop intelAW if ideo5==2, level(83.4)
svy, subpop(RwhiteNH): prop intelAB if ideo5==2, level(83.4)

svy, subpop(RwhiteNH): prop intelWB if ideo5==3, level(83.4)
svy, subpop(RwhiteNH): prop intelAW if ideo5==3, level(83.4)
svy, subpop(RwhiteNH): prop intelAB if ideo5==3, level(83.4)

svy, subpop(RwhiteNH): prop intelWB if ideo5==4, level(83.4)
svy, subpop(RwhiteNH): prop intelAW if ideo5==4, level(83.4)
svy, subpop(RwhiteNH): prop intelAB if ideo5==4, level(83.4)

svy, subpop(RwhiteNH): prop intelWB if ideo5==5, level(83.4)
svy, subpop(RwhiteNH): prop intelAW if ideo5==5, level(83.4)
svy, subpop(RwhiteNH): prop intelAB if ideo5==5, level(83.4)

*******************************************************************************************
** Overlaps for the waffle plot
*******************************************************************************************

     tab intelWB intelAW
svy: tab intelWB intelAW

     tab intelWB intelAW if gdiffbiol>=1 & gdiffbiol<=4
svy: tab intelWB intelAW if gdiffbiol>=1 & gdiffbiol<=4


     prop intelWB                , level(83.4)
     prop         intelAW        , level(83.4)
     prop                 intelAB, level(83.4)

     prop intelWB                 if weight!=., level(83.4)
     prop         intelAW         if weight!=., level(83.4)
     prop                 intelAB if weight!=., level(83.4)

svy: prop intelWB, level(83.4)
svy: prop intelAW, level(83.4)

*******************************************************************************************
** Extra analysis
*******************************************************************************************

gen     intelWB3 = 99 if weight!=.
replace intelWB3 = 3 if (intelW > intelB)  & intelW!=. & intelB!=.
replace intelWB3 = 1 if (intelW < intelB)  & intelW!=. & intelB!=.
replace intelWB3 = 2 if (intelW == intelB) & intelW!=. & intelB!=.
tab     intelWB3, mi

gen     intelAB3 = 99 if weight!=.
replace intelAB3 = 3 if (intelA > intelB)  & intelA!=. & intelB!=.
replace intelAB3 = 1 if (intelA < intelB)  & intelA!=. & intelB!=.
replace intelAB3 = 2 if (intelA == intelB) & intelA!=. & intelB!=.
tab     intelAB3, mi

gen     intelAW3 = 99 if weight!=.
replace intelAW3 = 3 if (intelA > intelW)  & intelA!=. & intelW!=.
replace intelAW3 = 1 if (intelA < intelW)  & intelA!=. & intelW!=.
replace intelAW3 = 2 if (intelA == intelW) & intelA!=. & intelW!=.
tab     intelAW3, mi

svy: prop intelWB3, level(83.4)
svy: prop intelAB3, level(83.4)
svy: prop intelAW3, level(83.4)