clear
import delimited Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\MAINDATA_AE.csv
save "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AE.dta", replace

******************************************************************************************************
*notes: The datasets MAINDATA_AEsamplecutoff. dta and MAINDATA_AEsamplecutoffCEM.dta are subsets of the main dataset MAINDATA_AE.dta, that are used in the analysis presented in the Appendix

use "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AE.dta", clear
keep if (smapexcpastinc~=1)&(abs(pinc1_1)==1)
save  "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AEsamplecutoff.dta", replace
	
* CEM - Coarsened exact matching sample, following the methodology in Iacus, King, and Porro (2011, 2012)

cem female white casado votedbolsonaro college livesomeone60orolder receivebenefits incadunicobefore resp_unemployed_pre, treatment(eligaeincpre)
gen cem_matched2=cem_matched

save "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AEsamplecutoffCEM.dta", replace	
*****************************************************************************************************

use "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AE.dta", clear

global baseline "pcpastinc_midp votedbolsonaro livesomeone60orolder educat* white  state* numpeoplehh agegrp1 agegrp2 agegrp3 agegrp4 agegrp5 female soc_des_score incadunicobefore receivebenefits index_behaviorbefore single married cohabiting divorced resp_unemployed_pre"
global RDcontrols "votedbolsonaro livesomeone60orolder educat* white  state* agegrp1 agegrp2 agegrp3 agegrp4 agegrp5 female soc_des_score incadunicobefore receivebenefits index_behaviorbefore single married cohabiting divorced resp_unemployed_pre"

* The sample in smapexcpastinc==1 (excluded from the analysis) are individuals potentially eligible for AE that live in a low-income household, whose total monthly income is more than three minimum salaries, and the income per capita is less than half the minimum salary per person
* see explanation in the paper in 2.2 Auxilio Emergencial


*****************************************************************************************************
*****************************************************************************************************
******************************* RESULTS - MAIN TEXT**************************************************
*****************************************************************************************************
*****************************************************************************************************

*****************************************************************************************************
****************************** Table 1 **************************************************************
*****************************************************************************************************

*column 1
sum hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all  if (smapexcpastinc~=1)&(pinc1_1==1)
  		foreach var of varlist hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all{
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T1c1.xls", br keep(eligaeincpre) nocons
	
			}
					
*column 2
sum hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all if (smapexcpastinc~=1)&(pinc1_1>=1&pinc1_1<=5)
  		foreach var of varlist hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all{
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T1c2.xls", br keep(eligaeincpre) nocons
	
			}		
			
					
*****************************************************************************************************
****************************** Table 2 **************************************************************
*****************************************************************************************************

*column 1
sum hadcovid19 testedcovid19 somehh_hosp  if (smapexcpastinc~=1)&(pinc1_1==1)
  		foreach var of varlist hadcovid19 testedcovid19 somehh_hosp {
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T2c1.xls", br keep(eligaeincpre) nocons
				}
			
	
*column 2	
sum hadcovid19 testedcovid19 somehh_hosp if (smapexcpastinc~=1)&(pinc1_1>=1&pinc1_1<=5)
  		foreach var of varlist hadcovid19 testedcovid19 somehh_hosp {
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5), robust
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T2c2.xls", br keep(eligaeincpre) nocons
				}	
			
*****************************************************************************************************
****************************** Table 3 **************************************************************
*****************************************************************************************************	
			
*column 1
sum anydasssym haddepression hadstress hadanxiety if (smapexcpastinc~=1)&(smapexcpastinc~=1)&(pinc1_1==1)
  		foreach var of varlist anydasssym haddepression hadstress hadanxiety{
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T3c1.xls", br keep(eligaeincpre) nocons
	
			}
			
*column 2			
sum anydasssym haddepression hadstress hadanxiety if (smapexcpastinc~=1)&(pinc1_1>=1&pinc1_1<=5)
  		foreach var of varlist anydasssym haddepression hadstress hadanxiety{
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5), robust
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T3c2.xls", br keep(eligaeincpre) nocons
	
			}		
			
								
*****************************************************************************************************
****************************** Table 4 **************************************************************
*****************************************************************************************************
		
*column 1
sum scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome if (smapexcpastinc~=1)&(pinc1_1==1)
  	foreach var of varlist scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
	reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T4c1.xls", br keep(eligaeincpre) nocons
	
			}

*column 2
sum scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome if (smapexcpastinc~=1)&(pinc1_1>=1&pinc1_1<=5)
  		foreach var of varlist scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T4c2.xls", br keep(eligaeincpre) nocons
	
			}

*****************************************************************************************************
****************************** Table 5 **************************************************************
*****************************************************************************************************

*Panel A, column 1
sum index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome if (votedbolsonaro==1&smapexcpastinc~=1)&(pinc1_1==1)
	foreach var of varlist index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
	
	reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1)&votedbolsonaro==1, robust
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T5bolso_c1.xls", br keep(eligaeincpre) nocons
			}
 
*Panel A, column 2
 sum index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome if (votedbolsonaro==0&smapexcpastinc~=1)&(pinc1_1==1)
     foreach var of varlist index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
	reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1)&votedbolsonaro==0, robust
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T5bolso_c2.xls", br keep(eligaeincpre) nocons
			}

*Panel B			
   		foreach var of varlist index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
	reg `var' eligaeincpre eligaeincpre_bolso $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust
outreg2 eligaeincpre eligaeincpre_bolso using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T5bolso_c3.xls", br keep(eligaeincpre eligaeincpre_bolso) nocons
			}
 
		
*****************************************************************************************************
****************************** Table 6 **************************************************************
***************************************************************************************************** 
sum scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome if receiveanyae==1&onlyadult==1&receivedae_14days==0
 
*column 1
 foreach var of varlist scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
		reg `var' receivedae_14days  if receiveanyae==1&onlyadult==1, robust 
		outreg2 receivedae_14days using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T6c1.xls", br keep(receivedae_14days) nocons
	}
	
*column 2	
	foreach var of varlist scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
		reg `var' receivedae_14days $baseline  if receiveanyae==1&onlyadult==1, robust 
		outreg2 receivedae_14days using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T6c2.xls", br keep(receivedae_14days) nocons
	}

*column 3
	foreach var of varlist scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
 		reg `var' receivedae_14days $baseline numaereceived if receiveanyae==1&onlyadult==1, robust 
		outreg2 receivedae_14days using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\T6c3.xls", br keep(receivedae_14days) nocons
	}
	

*****************************************************************************************************
****************************** Figure 1 *************************************************************
*****************************************************************************************************
* the data used for Figure 1 is obtained with the command below:
	
bysort pinc1_1: sum receiveanyae if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5)
	
*****************************************************************************************************
*****************************************************************************************************
******************************* RESULTS - APENDIX ***************************************************
*****************************************************************************************************
*****************************************************************************************************	
use "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AE.dta", clear

*****************************************************************************************************
****************************** Table A1 *************************************************************
*****************************************************************************************************

* sample: 2.6-3.4 minimum salaries
sum female white idade casado votedbolsonaro college livesomeone60orolder numpeoplehh pcpastinc_midpercapita receivebenefits incadunicobefore resp_unemployed_pre if (smapexcpastinc~=1)&(abs(pinc1_1)==1)

* sample: 2-5 miminum salaries 
sum female white idade casado votedbolsonaro college livesomeone60orolder numpeoplehh pcpastinc_midpercapita receivebenefits incadunicobefore resp_unemployed_pre if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5)


*****************************************************************************************************
****************************** Table A2 *************************************************************
*****************************************************************************************************

* column 1	
sum behbp1 behbp2 behbp3 behbp4 behbp5 behbp6 behbp7 behbp8 behbp11 behbp12 behbp13 behbp14 behbp15 behbp16 if (smapexcpastinc~=1)&(abs(pinc1_1)==1)
sum behdp1 behdp2 behdp3 behdp4 behdp5 behdp6 behdp7 behdp8 behdp11 behdp15 behdp16 behdp17 behdp18 behdp19 if (smapexcpastinc~=1)&(abs(pinc1_1)==1)

* column 2	
sum behbp1 behbp2 behbp3 behbp4 behbp5 behbp6 behbp7 behbp8 behbp11 behbp12 behbp13 behbp14 behbp15 behbp16 if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5)
sum behdp1 behdp2 behdp3 behdp4 behdp5 behdp6 behdp7 behdp8 behdp11 behdp15 behdp16 behdp17 behdp18 behdp19 if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5)



*****************************************************************************************************
****************************** Table A3 *************************************************************
*****************************************************************************************************

* sample: 2.6-3.4 minimum salaries
sum receiveanyae numaereceived if (smapexcpastinc~=1)&(abs(pinc1_1)==1)
sum numaereceived if receiveanyae==1&(smapexcpastinc~=1)&(abs(pinc1_1)==1)
sum hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp  anydasssym haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome anton_familystayhome  if (smapexcpastinc~=1)&(abs(pinc1_1)==1)


* sample: 2-5 miminum salaries
sum receiveanyae numaereceived if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5)
sum numaereceived if receiveanyae==1&(smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5)
sum hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp  anydasssym haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome anton_familystayhome  if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5)


*****************************************************************************************************
****************************** Table A4 *************************************************************
*****************************************************************************************************
 

*column 1
sum receiveanyae if smapexcpastinc~=1&eligaeincpre==0&(pinc1_1>=0&pinc1_1<=5)
reg receiveanyae eligaeincpre if smapexcpastinc~=1&(pinc1_1>=-3&pinc1_1<=5), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA4c1.xls", br keep(eligaeincpre) nocons

*column 2
sum receiveanyae if smapexcpastinc~=1&eligaeincpre==0&(pinc1_1>=0&pinc1_1<=5)
reg receiveanyae eligaeincpre $baseline if smapexcpastinc~=1&(pinc1_1>=-3&pinc1_1<=5), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA4c2.xls", br keep(eligaeincpre) nocons

*column 3
sum receiveanyae if smapexcpastinc~=1&eligaeincpre==0&pinc1_1==1
reg receiveanyae eligaeincpre if smapexcpastinc~=1&(abs(pinc1_1)==1), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA4c3.xls", br keep(eligaeincpre) nocons

*column 4
sum receiveanyae if smapexcpastinc~=1&eligaeincpre==0&pinc1_1==1
reg receiveanyae eligaeincpre $baseline if smapexcpastinc~=1&(abs(pinc1_1)==1), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA4c4.xls", br keep(eligaeincpre) nocons


*****************************************************************************************************
****************************** Table A5 *************************************************************
*****************************************************************************************************

use "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AEsamplecutoff.dta", clear
	
	sum female votedbolsonaro white black otherrace incomplete_pri complete_pri incomplete_hs complete_hs incomplete_he college postgrad single married cohabiting divorced widower numpeoplehh receivebenefits incadunicobefore index_behaviorbefore resp_unemployed_pre agegrp1 agegrp2 agegrp3 agegrp4 agegrp5 livesomeoneless17 livesomeone60orolder  diabetes high_bp heart_disease asthma_or_resp_condition cancer allergies if (eligaeincpre==1)&(smapexcpastinc~=1)&(abs(pinc1_1)==1)
	sum female votedbolsonaro white black otherrace incomplete_pri complete_pri incomplete_hs complete_hs incomplete_he college postgrad single married cohabiting divorced widower numpeoplehh receivebenefits incadunicobefore index_behaviorbefore resp_unemployed_pre agegrp1 agegrp2 agegrp3 agegrp4 agegrp5 livesomeoneless17 livesomeone60orolder  diabetes high_bp heart_disease asthma_or_resp_condition allergies cancer  if (eligaeincpre==0)&(smapexcpastinc~=1)&(abs(pinc1_1)==1)
		
	foreach var of varlist female votedbolsonaro white black otherrace incomplete_pri complete_pri incomplete_hs complete_hs incomplete_he college postgrad single married cohabiting divorced widower numpeoplehh receivebenefits incadunicobefore index_behaviorbefore resp_unemployed_pre agegrp1 agegrp2 agegrp3 agegrp4 agegrp5 livesomeoneless17 livesomeone60orolder  diabetes high_bp heart_disease asthma_or_resp_condition allergies cancer {
	ttest `var', by (eligaeincpre) 
			}

*****************************************************************************************************
****************************** Table A6 *************************************************************
*****************************************************************************************************

use "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AE.dta", clear

*column 1
sum hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome if (smapexcpastinc~=1)&(pinc1_1==1)

  		foreach var of varlist hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA6c1.xls", br keep(eligaeincpre) nocons
	
			}

*column 2 
  		foreach var of varlist hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
		reg `var' eligaeincpre if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA6c2.xls", br keep(eligaeincpre) nocons
	
			}

*column 3 
use "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AEsamplecutoffCEM.dta", clear

sum hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome if (smapexcpastinc~=1)&(pinc1_1==1)&(cem_matched2==1)
 		foreach var of varlist hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
		reg `var' eligaeincpre $baseline if cem_matched2==1, robust 
		outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA6c3.xls", br keep(eligaeincpre) nocons
	
			}		
			

*****************************************************************************************************
****************************** Table A7 *************************************************************
*****************************************************************************************************	
						
use "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AEsamplecutoff.dta", clear

gen arm=1 if eligaeincpre==1
replace arm=0 if eligaeincpre==0

*column 1 (test explaned in List et al. 2019)
mhtexp receiveanyae index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome, treatment(arm)

* column 2
***two-stage q-values as explained in Anderson (2008)
** Michael Anderson’s Stata code for sharpened q values is on: http://are.berkeley.edu/~mlanderson/downloads/fdr_sharpened_qvalues.do.
*the procedure below documents p-values for the variable eligaeincpre  that are then transformed to q-values using the code: fdr_sharpened_qvalues.do. 
    foreach var of varlist receiveanyae index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
		reg `var' eligaeincpre $baseline, robust 
		}


 
use "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AEsamplecutoffCEM.dta", clear
gen arm=1 if eligaeincpre==1
replace arm=0 if eligaeincpre==0
 
keep if cem_matched2==1
 
*column 3 (test explaned in List et al. 2019)
mhtexp receiveanyae index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome, treatment(arm)

*column 4
*the procedure below documents p-values for the variable eligaeincpre  that are then transformed to q-values using the code: fdr_sharpened_qvalues.do. 
    	foreach var of varlist receiveanyae index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust 
			}

  
*column 5
use "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AE.dta", clear
*the procedure below documents p-values for the variable eligaeincpre  that are then transformed to q-values using the code: fdr_sharpened_qvalues.do.   
   		foreach var of varlist receiveanyae index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
		reg `var' eligaeincpre $baseline if (pinc1_1>=-3&pinc1_1<=5), robust 
			}
 
 		
*****************************************************************************************************
****************************** Table A8 *************************************************************
*****************************************************************************************************					
			
use "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AE.dta", clear

* RD variables	
gen pastinc_midp=1150 if p51==1
replace pastinc_midp=1350 if p51==2
replace pastinc_midp=1750 if p51==3
replace pastinc_midp=2150 if p51==4
replace pastinc_midp=2550 if p51==5
replace pastinc_midp=2950 if p51==6
replace pastinc_midp=3350 if p51==7
replace pastinc_midp=3750 if p51==8
replace pastinc_midp=4150 if p51==9
replace pastinc_midp=4550 if p51==10
replace pastinc_midp=4950 if p51==11
replace pastinc_midp=5350 if p51==12
replace pastinc_midp=5750 if p51==13
replace pastinc_midp=6150 if p51==14
replace pastinc_midp=6550 if p51==15
replace pastinc_midp=6950 if p51==16
replace pastinc_midp=7350 if p51==17
replace pastinc_midp=7551 if p51==18
	
gen pincomemid=(pastinc_midp-3135)
gen pinpc_1=pincomemid/p6
 
gen pinpc_1_el=pinpc_1*eligaeincpre

save "Z:\Covid_AuxilioEmergencial\codes_Covid_Brasil\MAINDATA_AE.dta", replace

*Panel A
foreach var of varlist hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{

rdbwselect `var' pinpc_1, covs(votedbolsonaro livesomeone60orolder educat* white  state* agegrp1 agegrp2 agegrp3 agegrp4 agegrp5 female soc_des_score incadunicobefore receivebenefits index_behaviorbefore single married cohabiting divorced resp_unemployed_pre) 			
gen x12_`var'=e(h_mserd)
reg `var' eligaeincpre pinpc_1 pinpc_1_el $RDcontrols if smapexcpastinc~=1&abs(pinpc_1)<x12_`var', robust 

outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA8c1.xls", keep(eligaeincpre) nocons br append

}

foreach var of varlist hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
egen mx11_`var'=mean(`var') if smapexcpastinc~=1&abs(pinpc_1)<x12_`var'&eligaeincpre==0
}
*means for Panel A
sum mx11_*
			
*Panel B	
foreach var of varlist hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
	
		reg `var' eligaeincpre pinpc_1 pinpc_1_el $RDcontrols if smapexcpastinc~=1&abs(pinpc_1)<500, robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA8c2.xls", br keep(eligaeincpre) nocons
	
			}	

	
*Panel C
foreach var of varlist hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
	
	reg `var' eligaeincpre pinpc_1 pinpc_1_el $RDcontrols if smapexcpastinc~=1&abs(pinpc_1)<900, robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA8c3.xls", br keep(eligaeincpre) nocons
	
			}
	
*means for Panel C
foreach var of varlist hw_occassions_ov frvisitfriends index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 somehh_hosp haddepression hadstress hadanxiety scorecovidquiz1st bolsonaro_mtrust sciencedoctors_mtrust trustbolsomorethanscience ownchance  youngfemalechance  oldmalechance anton_stayhome  anton_familystayhome{
egen mx900_`var'=mean(`var') if smapexcpastinc~=1&abs(pinpc_1)<900&eligaeincpre==0
}
sum mx900_*	

	
*****************************************************************************************************************************
****************************** results for footnote in Table A8 *************************************************************
*****************************************************************************************************************************
	
foreach var of varlist female votedbolsonaro white black otherrace incomplete_pri complete_pri incomplete_hs complete_hs incomplete_he college postgrad single married cohabiting divorced widower numpeoplehh receivebenefits incadunicobefore index_behaviorbefore resp_unemployed_pre agegrp1 agegrp2 agegrp3 agegrp4 agegrp5 livesomeoneless17 livesomeone60orolder  diabetes high_bp heart_disease asthma_or_resp_condition allergies cancer{
	
		reg `var' eligaeincpre pinpc_1 pinpc_1_el if (smapexcpastinc~=1)&(abs(pinc1_1)<500), robust 
		reg `var' eligaeincpre pinpc_1 pinpc_1_el if (smapexcpastinc~=1)&(abs(pinc1_1)<900), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA8note.xls", br keep(eligaeincpre) nocons
			}
	
			
	
*****************************************************************************************************
****************************** Table A9 *************************************************************
*****************************************************************************************************	

*column 1
sum behdp1 behdp2 behdp3 behdp4 behdp5 behdp6 behdp7 behdp8 behdp11 behdp15 behdp16 behdp17 behdp18 behdp19 if (smapexcpastinc~=1)&(pinc1_1==1)
  foreach var of varlist  behdp1 behdp2 behdp3 behdp4 behdp5 behdp6 behdp7 behdp8 behdp11 behdp15 behdp16 behdp17 behdp18 behdp19 {
  
reg `var' eligaeincpre $baseline  if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA9c1.xls", br keep(eligaeincpre) nocons

  } 

  
* column 2  
sum behdp1 behdp2 behdp3 behdp4 behdp5 behdp6 behdp7 behdp8 behdp11 behdp15 behdp16 behdp17 behdp18 behdp19 if (smapexcpastinc~=1)&(pinc1_1>=1&pinc1_1<=5)
 foreach var of varlist  behdp1 behdp2 behdp3 behdp4 behdp5 behdp6 behdp7 behdp8 behdp11 behdp15 behdp16 behdp17 behdp18 behdp19 {
    reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA9c2.xls", br keep(eligaeincpre) nocons


  }

*****************************************************************************************************
****************************** Table A10 ************************************************************
*****************************************************************************************************

*column 1
sum index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19  if (resp_me_iw==1&smapexcpastinc~=1)&(pinc1_1==1)
  		foreach var of varlist index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 {
		
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1)&resp_me_iw==1, robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA10c1.xls", br keep(eligaeincpre) nocons
			}

*column 2			
sum index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 if (resp_me_iw==0&smapexcpastinc~=1)&(pinc1_1==1)
  		foreach var of varlist index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 {
			reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1)&resp_me_iw==0, robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA10c2.xls", br keep(eligaeincpre) nocons
			}	
	
			
	gen eligaeincpre_informal=eligaeincpre*resp_me_iw
	
*column 3			   		
	foreach var of varlist index_behavior_pandemic resp_unemployed_curr hours_worked_cont hours_worked_all hadcovid19 testedcovid19 {
		
				reg `var' eligaeincpre eligaeincpre_informal $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust
outreg2 eligaeincpre eligaeincpre_informal using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA10c3.xls", br keep(eligaeincpre eligaeincpre_informal) nocons
			}

*****************************************************************************************************
****************************** Table A11 ************************************************************
*****************************************************************************************************

*column 1
  sum infoatleastadaycovid infolessonceaweekcovid howmanyinfouse scorenoncovidquiz1st if (smapexcpastinc~=1)&((pinc1_1)==1) 
 foreach var of varlist  infoatleastadaycovid infolessonceaweekcovid howmanyinfouse scorenoncovidquiz1st{
  reg `var' eligaeincpre $baseline  if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA11c1.xls", br keep(eligaeincpre) nocons

  }

  
*column 2  
 sum infoatleastadaycovid infolessonceaweekcovid howmanyinfouse scorenoncovidquiz1st if (smapexcpastinc~=1)&(pinc1_1>=1&pinc1_1<=5) 
  foreach var of varlist  infoatleastadaycovid infolessonceaweekcovid howmanyinfouse scorenoncovidquiz1st{
   reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA11c2.xls", br keep(eligaeincpre) nocons

  } 
  
  
*****************************************************************************************************
****************************** Table A12 ************************************************************
*****************************************************************************************************


** column 1
sum hhcope1 hhcope2 hhcope3 hhcope4 hhcope5 hhcope6 hhcope7 hhcope8 hhcope9 hhcope10 hhcope11 hhcope12 hhcope13 hhcope14 hhcope15  if (smapexcpastinc~=1)&(pinc1_1==1)

  		foreach var of varlist hhcope1 hhcope2 hhcope3 hhcope4 hhcope5 hhcope6 hhcope7 hhcope8 hhcope9 hhcope10 hhcope11 hhcope12 hhcope13 hhcope14 hhcope15 {
		
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(abs(pinc1_1)==1), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA12c1.xls", br keep(eligaeincpre) nocons
	
			}

** column 2			
			sum hhcope1 hhcope2 hhcope3 hhcope4 hhcope5 hhcope6 hhcope7 hhcope8 hhcope9 hhcope10 hhcope11 hhcope12 hhcope13 hhcope14 hhcope15  if (smapexcpastinc~=1)&(pinc1_1>=1&pinc1_1<=5)

  		foreach var of varlist hhcope1 hhcope2 hhcope3 hhcope4 hhcope5 hhcope6 hhcope7 hhcope8 hhcope9 hhcope10 hhcope11 hhcope12 hhcope13 hhcope14 hhcope15 {
		
		reg `var' eligaeincpre $baseline if (smapexcpastinc~=1)&(pinc1_1>=-3&pinc1_1<=5), robust 
outreg2 eligaeincpre using "Z:\Covid_AuxilioEmergencial\revisions_JEBO\dataforpublication\TableA12c2.xls", br keep(eligaeincpre) nocons
	
			}
			
			
			
*****************************************************************************************************
****************************** Figure A1 ************************************************************
*****************************************************************************************************

* the data from figure 1 (i.e., number of observations) comes from variable p51
tab p51

*****************************************************************************************************
****************************** Figure A2 ************************************************************
*****************************************************************************************************

* Figure A2  
hist pinpc_1 if (smapexcpastinc~=1)&(pinpc_1>-900&pinpc_1<900), bin(94)


*note. To install package for cattaneo et al test reported in Figure A2 footnote, use the following stata command:
* net install rddensity, from(https://raw.githubusercontent.com/rdpackages/rddensity/master/stata) replace

rddensity pinpc_1