use "${inter}/Psycorona_long_full",clear cap drop country encode CountryName,gen(country) cap drop coronaClose gen coronaClose=1 if w_coronaClose_1==1 & w_coronaClose_2==1 replace coronaClose=2 if w_coronaClose_1==1 & w_coronaClose_2==0 replace coronaClose=3 if w_coronaClose_2==1 & w_coronaClose_1==0 replace coronaClose=0 if mi(coronaClose) cap label drop coronaClose label define coronaClose 1 "self+family" 2"self" 3"family" 0 "no" forvalues i =0(1)4{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } forvalues i =11(1)22{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } global x i.age i.gender i.unemployment i.education i.w_PFS lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 people_vaccinated_per_hundred cases_population_7day stringency_index //Figure 4: Panel A xtset id wave xtreg w_c19perBeh01_Harmonized Lag_D4#c.sft_pca Lag_D11#c.sft_pca Lag_D12#c.sft_pca Lag_D13#c.sft_pca Lag_D14#c.sft_pca Lag_D15#c.sft_pca Lag_D16#c.sft_pca Lag_D17#c.sft_pca Lag_D18#c.sft_pca Lag_D19#c.sft_pca Lag_D20#c.sft_pca Lag_D21#c.sft_pca Lag_D22#c.sft_pca sft_pca $x if m2!=2, vce(cluster country) frame reset frame create orsCI or orlb orub lincom _b[sft_pca] frame post orsCI (`r(estimate)') (`r(lb)') (`r(ub)') lincom 1.Lag_D4#c.sft_pca + sft_pca frame post orsCI (`r(estimate)') (`r(lb)') (`r(ub)') forvalues a = 11(1)22 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI (`r(estimate)') (`r(lb)') (`r(ub)') } frame change orsCI format or %3.2f format orlb %3.2f format orub %3.2f save "${output}/washing_all_age",replace use "${inter}/Psycorona_long_full",clear cap drop id_wave egen id_wave= tag(id wave) drop if cases_population_7day==. cap drop country encode CountryName,gen(country) forvalues i =0(1)4{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } forvalues i =11(1)22{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } global x i.age i.gender i.unemployment i.education i.w_PFS lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 people_vaccinated_per_hundred cases_population_7day stringency_index //Figure 4 panel A xtset id wave xtreg w_c19perBeh01_Harmonized Lag_D4#c.sft_pca Lag_D11#c.sft_pca Lag_D12#c.sft_pca Lag_D13#c.sft_pca Lag_D14#c.sft_pca Lag_D15#c.sft_pca Lag_D16#c.sft_pca Lag_D17#c.sft_pca Lag_D18#c.sft_pca Lag_D19#c.sft_pca Lag_D20#c.sft_pca Lag_D21#c.sft_pca Lag_D22#c.sft_pca sft_pca $x if m2!=2&age<=2, vce(cluster country) frame create orsCI1 or1 orlb1 orub1 lincom _b[sft_pca] frame post orsCI1 (`r(estimate)') (`r(lb)') (`r(ub)') lincom 1.Lag_D4#c.sft_pca + sft_pca frame post orsCI1 (`r(estimate)') (`r(lb)') (`r(ub)') forvalues a = 11(1)22 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI1 (`r(estimate)') (`r(lb)') (`r(ub)') } frame change orsCI1 format or1 %3.2f format orlb1 %3.2f format orub1 %3.2f save "${output}/washing_34",replace use "${output}/washing_34",clear append using "${output}/washing_all_age" replace or=or[_n+14] replace orlb=orlb[_n+14] replace orub=orub[_n+14] drop if or1==. cap drop wave1 gen wave1=_n cap drop wave2 gen wave2=_n+0.4 graph twoway (scatter or wave1, lcolor("5 101 141") sort msize(vsmall) mfcolor("5 101 141") mlcolor("5 101 141") msymbol(diamond)) (rcap orlb orub wave1,color("5 101 141")) (scatter or1 wave2, sort msize(vsmall) mfcolor("223 82 54") mlcolor("223 82 54") msymbol(diamond) lpattern(dash) lcolor(gray)) (rcap orlb1 orub1 wave2,color("223 82 54") ) (rconnected orlb orub wave1 if wave1<=0, msize(vsmall) mfcolor("5 101 141") mlcolor("5 101 141") msymbol(diamond)) (rconnected orlb1 orub1 wave2 if wave2<=0, msize(vsmall) mfcolor(red) mlcolor("223 82 54") msymbol(diamond) lcolor("223 82 54")), /// legend(off) ylab(,format(%3.1f)) /// xlab(1 "March 19, 2020" 2"April 25, 2020" 3 "June 13, 2020" 4 "July 13, 2020" 5 "August 13, 2020" 6 "September 13, 2020" 7 "October 13, 2020" 8 "November 13, 2020" 9 "December 13, 2020" 10 "January 13, 2021" 11 "February 13, 2021" 12 "March 13, 2021" 13 "April 13, 2021" 14 "July 13, 2021", angle(70) labsize(small) labcolor(black) axis(1)) /// subtitle("Panel A. Willingness to increase hand-washing frequency") xtitle("") ytitle("Estimated Coefficients") yline(0, lp(dash)) scheme(s1mono)saving("${output}/graph1",replace) //Figure 4: Panel B use "${inter}/Psycorona_long_full",clear cap drop id_wave egen id_wave= tag(id wave) drop if cases_population_7day==. cap drop country encode CountryName,gen(country) forvalues i =0(1)4{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } forvalues i =11(1)22{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } global x i.age i.gender i.unemployment i.education i.w_PFS lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 people_vaccinated_per_hundred cases_population_7day stringency_index xtset id wave xtreg w_c19perBeh02_Harmonized Lag_D4#c.sft_pca Lag_D11#c.sft_pca Lag_D12#c.sft_pca Lag_D13#c.sft_pca Lag_D14#c.sft_pca Lag_D15#c.sft_pca Lag_D16#c.sft_pca Lag_D17#c.sft_pca Lag_D18#c.sft_pca Lag_D19#c.sft_pca Lag_D20#c.sft_pca Lag_D21#c.sft_pca Lag_D22#c.sft_pca sft_pca $x if m2!=2, cluster(country) frame create orsCI2 or orlb orub lincom _b[sft_pca] frame post orsCI2 (`r(estimate)') (`r(lb)') (`r(ub)') lincom 1.Lag_D4#c.sft_pca + sft_pca frame post orsCI2 (`r(estimate)') (`r(lb)') (`r(ub)') forvalues a = 11(1)22 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI2 (`r(estimate)') (`r(lb)') (`r(ub)') } frame change orsCI2 format or %3.2f format orlb %3.2f format orub %3.2f save "${output}/crowd_all_age",replace use "${inter}/Psycorona_long_full",clear cap drop id_wave egen id_wave= tag(id wave) drop if cases_population_7day==. cap drop country encode CountryName,gen(country) forvalues i =0(1)4{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } forvalues i =11(1)22{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } global x i.age i.gender i.unemployment i.education i.w_PFS lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 people_vaccinated_per_hundred cases_population_7day stringency_index //Figure 4: Panel B xtset id wave xtreg w_c19perBeh02_Harmonized Lag_D4#c.sft_pca Lag_D11#c.sft_pca Lag_D12#c.sft_pca Lag_D13#c.sft_pca Lag_D14#c.sft_pca Lag_D15#c.sft_pca Lag_D16#c.sft_pca Lag_D17#c.sft_pca Lag_D18#c.sft_pca Lag_D19#c.sft_pca Lag_D20#c.sft_pca Lag_D21#c.sft_pca Lag_D22#c.sft_pca sft_pca $x if m2!=2&age<=2, cluster(country) frame create orsCI3 or1 orlb1 orub1 lincom _b[sft_pca] frame post orsCI3 (`r(estimate)') (`r(lb)') (`r(ub)') lincom 1.Lag_D4#c.sft_pca + sft_pca frame post orsCI3 (`r(estimate)') (`r(lb)') (`r(ub)') forvalues a = 11(1)22 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI3 (`r(estimate)') (`r(lb)') (`r(ub)') } frame change orsCI3 format or1 %3.2f format orlb1 %3.2f format orub1 %3.2f save "${output}/crowd_34",replace use "${output}/crowd_34",clear append using "${output}/crowd_all_age" replace or=or[_n+14] replace orlb=orlb[_n+14] replace orub=orub[_n+14] drop if or1==. cap drop wave1 gen wave1=_n cap drop wave2 gen wave2=_n+0.4 graph twoway (scatter or wave1, lcolor("5 101 141") sort msize(vsmall) mfcolor("5 101 141") mlcolor("5 101 141") msymbol(diamond)) (rcap orlb orub wave1,color("5 101 141")) (scatter or1 wave2, sort msize(vsmall) mfcolor("223 82 54") mlcolor("223 82 54") msymbol(diamond) lpattern(dash) lcolor(gray)) (rcap orlb1 orub1 wave2,color("223 82 54") ) (rconnected orlb orub wave1 if wave1<=0, msize(vsmall) mfcolor("5 101 141") mlcolor("5 101 141") msymbol(diamond)) (rconnected orlb1 orub1 wave2 if wave2<=0, msize(vsmall) mfcolor(red) mlcolor("223 82 54") msymbol(diamond) lcolor("223 82 54")), /// legend(off) ylab(,format(%3.1f)) /// xlab(1 "March 19, 2020" 2"April 25, 2020" 3 "June 13, 2020" 4 "July 13, 2020" 5 "August 13, 2020" 6 "September 13, 2020" 7 "October 13, 2020" 8 "November 13, 2020" 9 "December 13, 2020" 10 "January 13, 2021" 11 "February 13, 2021" 12 "March 13, 2021" 13 "April 13, 2021" 14 "July 13, 2021", angle(70) labsize(small) labcolor(black) axis(1)) /// subtitle("Panel B. Willingness to avoid crowded places") xtitle("") ytitle("Estimated Coefficients") yline(0, lp(dash)) scheme(s1mono)saving("${output}/graph2",replace) //Figure 4: Panel C use "${inter}/Psycorona_long_full",clear cap drop id_wave egen id_wave= tag(id wave) drop if cases_population_7day==. cap drop country encode CountryName,gen(country) forvalues i =0(1)4{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } forvalues i =11(1)22{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } global x i.age i.gender i.unemployment i.education i.w_PFS lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 people_vaccinated_per_hundred cases_population_7day stringency_index xtset id wave xtreg w_c19perBeh03_Harmonized Lag_D4#c.sft_pca Lag_D11#c.sft_pca Lag_D12#c.sft_pca Lag_D13#c.sft_pca Lag_D14#c.sft_pca Lag_D15#c.sft_pca Lag_D16#c.sft_pca Lag_D17#c.sft_pca Lag_D18#c.sft_pca Lag_D19#c.sft_pca Lag_D20#c.sft_pca Lag_D21#c.sft_pca Lag_D22#c.sft_pca sft_pca $x if m2!=2, cluster(country) frame create orsCI4 or orlb orub lincom _b[sft_pca] frame post orsCI4 (`r(estimate)') (`r(lb)') (`r(ub)') lincom 1.Lag_D4#c.sft_pca + sft_pca frame post orsCI4 (`r(estimate)') (`r(lb)') (`r(ub)') forvalues a = 11(1)22 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI4 (`r(estimate)') (`r(lb)') (`r(ub)') } frame change orsCI4 format or %3.2f format orlb %3.2f format orub %3.2f save "${output}/quarantine_all_age",replace use "${inter}/Psycorona_long_full",clear cap drop id_wave egen id_wave= tag(id wave) drop if cases_population_7day==. cap drop country encode CountryName,gen(country) forvalues i =0(1)4{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } forvalues i =11(1)22{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } global x i.age i.gender i.unemployment i.education i.w_PFS lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 people_vaccinated_per_hundred cases_population_7day stringency_index //Figure 4: Panel C xtset id wave xtreg w_c19perBeh03_Harmonized Lag_D4#c.sft_pca Lag_D11#c.sft_pca Lag_D12#c.sft_pca Lag_D13#c.sft_pca Lag_D14#c.sft_pca Lag_D15#c.sft_pca Lag_D16#c.sft_pca Lag_D17#c.sft_pca Lag_D18#c.sft_pca Lag_D19#c.sft_pca Lag_D20#c.sft_pca Lag_D21#c.sft_pca Lag_D22#c.sft_pca sft_pca $x if m2!=2&age<=2, cluster(country) frame create orsCI5 or1 orlb1 orub1 lincom _b[sft_pca] frame post orsCI5 (`r(estimate)') (`r(lb)') (`r(ub)') lincom 1.Lag_D4#c.sft_pca + sft_pca frame post orsCI5 (`r(estimate)') (`r(lb)') (`r(ub)') forvalues a = 11(1)22 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI5 (`r(estimate)') (`r(lb)') (`r(ub)') } frame change orsCI5 format or1 %3.2f format orlb1 %3.2f format orub1 %3.2f save "${output}/quarantine_34",replace use "${output}/quarantine_34",clear append using "${output}/quarantine_all_age" replace or=or[_n+14] replace orlb=orlb[_n+14] replace orub=orub[_n+14] drop if or1==. cap drop wave1 gen wave1=_n cap drop wave2 gen wave2=_n+0.4 graph twoway (scatter or wave1, lcolor("5 101 141") sort msize(vsmall) mfcolor("5 101 141") mlcolor("5 101 141") msymbol(diamond)) (rcap orlb orub wave1,color("5 101 141")) (scatter or1 wave2, sort msize(vsmall) mfcolor("223 82 54") mlcolor("223 82 54") msymbol(diamond) lpattern(dash) lcolor(gray)) (rcap orlb1 orub1 wave2,color("223 82 54") ) (rconnected orlb orub wave1 if wave1<=0, msize(vsmall) mfcolor("5 101 141") mlcolor("5 101 141") msymbol(diamond)) (rconnected orlb1 orub1 wave2 if wave2<=0, msize(vsmall) mfcolor(red) mlcolor("223 82 54") msymbol(diamond) lcolor("223 82 54")), /// legend(off) ylab(,format(%3.1f)) /// xlab(1 "March 19, 2020" 2"April 25, 2020" 3 "June 13, 2020" 4 "July 13, 2020" 5 "August 13, 2020" 6 "September 13, 2020" 7 "October 13, 2020" 8 "November 13, 2020" 9 "December 13, 2020" 10 "January 13, 2021" 11 "February 13, 2021" 12 "March 13, 2021" 13 "April 13, 2021" 14 "July 13, 2021", angle(70) labsize(small) labcolor(black) axis(1)) /// subtitle("Panel C. Willingness to self-quarantine") xtitle("") ytitle("Estimated Coefficients") yline(0, lp(dash)) scheme(s1mono)saving("${output}/graph3",replace) //Figure 4: Panel D use "${inter}/Psycorona_long_full",clear cap drop id_wave egen id_wave= tag(id wave) drop if cases_population_7day==. cap drop country encode CountryName,gen(country) forvalues i =0(1)4{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } forvalues i =11(1)22{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } global x i.age i.gender i.unemployment i.education i.w_PFS lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 people_vaccinated_per_hundred cases_population_7day stringency_index xtset id wave xtreg w_Vaccine Lag_D11#c.sft_pca Lag_D12#c.sft_pca Lag_D13#c.sft_pca Lag_D14#c.sft_pca Lag_D15#c.sft_pca Lag_D16#c.sft_pca Lag_D17#c.sft_pca Lag_D18#c.sft_pca Lag_D19#c.sft_pca Lag_D20#c.sft_pca Lag_D21#c.sft_pca Lag_D22#c.sft_pca sft_pca $x if m2!=2, cluster(country) frame create orsCI6 or orlb orub lincom _b[sft_pca] frame post orsCI6 (`r(estimate)') (`r(lb)') (`r(ub)') frame post orsCI6 (`r(estimate)') (`r(lb)') (`r(ub)') forvalues a = 11(1)22 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI6 (`r(estimate)') (`r(lb)') (`r(ub)') } frame change orsCI6 format or %3.2f format orlb %3.2f format orub %3.2f save "${output}/vaccine_all_age",replace use "${inter}/Psycorona_long_full",clear cap drop id_wave egen id_wave= tag(id wave) drop if cases_population_7day==. cap drop country encode CountryName,gen(country) forvalues i =0(1)4{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } forvalues i =11(1)22{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } global x i.age i.gender i.unemployment i.education i.w_PFS lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 people_vaccinated_per_hundred cases_population_7day stringency_index xtset id wave xtreg w_Vaccine Lag_D11#c.sft_pca Lag_D12#c.sft_pca Lag_D13#c.sft_pca Lag_D14#c.sft_pca Lag_D15#c.sft_pca Lag_D16#c.sft_pca Lag_D17#c.sft_pca Lag_D18#c.sft_pca Lag_D19#c.sft_pca Lag_D20#c.sft_pca Lag_D21#c.sft_pca Lag_D22#c.sft_pca sft_pca $x if m2!=2&age<=2, cluster(country) frame create orsCI7 or1 orlb1 orub1 lincom _b[sft_pca] frame post orsCI7 (`r(estimate)') (`r(lb)') (`r(ub)') frame post orsCI7 (`r(estimate)') (`r(lb)') (`r(ub)') forvalues a = 11(1)22 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI7 (`r(estimate)') (`r(lb)') (`r(ub)') } frame change orsCI7 format or1 %3.2f format orlb1 %3.2f format orub1 %3.2f save "${output}/vaccine_34",replace use "${output}/vaccine_34",clear append using "${output}/vaccine_all_age" replace or=or[_n+14] replace orlb=orlb[_n+14] replace orub=orub[_n+14] drop if or1==. drop in 1 cap drop wave1 gen wave1=_n cap drop wave2 gen wave2=_n+0.4 graph twoway (scatter or wave1, lcolor("5 101 141") sort msize(vsmall) mfcolor("5 101 141") mlcolor("5 101 141") msymbol(diamond)) (rcap orlb orub wave1,color("5 101 141")) (scatter or1 wave2, sort msize(vsmall) mfcolor("223 82 54") mlcolor("223 82 54") msymbol(diamond) lpattern(dash) lcolor(gray)) (rcap orlb1 orub1 wave2,color("223 82 54") ) (rconnected orlb orub wave1 if wave1<=0, msize(vsmall) mfcolor("5 101 141") mlcolor("5 101 141") msymbol(diamond)) (rconnected orlb1 orub1 wave2 if wave2<=0, msize(vsmall) mfcolor(red) mlcolor("223 82 54") msymbol(diamond) lcolor("223 82 54")), /// legend(off) ylab(,format(%3.1f)) /// xlab(1 "April 25, 2020" 2 "June 13, 2020" 3 "July 13, 2020" 4 "August 13, 2020" 5 "September 13, 2020" 6 "October 13, 2020" 7 "November 13, 2020" 8 "December 13, 2020" 9 "January 13, 2021" 10 "February 13, 2021" 11 "March 13, 2021" 12 "April 13, 2022" 13 "July 13, 2021", angle(70) labsize(small) labcolor(black) axis(1)) /// subtitle("Panel D. Willingness to vaccine for self") xtitle("") ytitle("Estimated Coefficients") yline(0, lp(dash)) scheme(s1mono)saving("${output}/graph4",replace) //Figure 4: Panel E use "${inter}/Psycorona_long_full",clear cap drop id_wave egen id_wave= tag(id wave) drop if cases_population_7day==. cap drop country encode CountryName,gen(country) forvalues i =0(1)4{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } forvalues i =11(1)22{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } global x i.age i.gender i.unemployment i.education i.w_PFS lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 people_vaccinated_per_hundred cases_population_7day stringency_index xtset id wave xtreg w_c19RCA01 Lag_D1#c.sft_pca Lag_D2#c.sft_pca Lag_D3#c.sft_pca Lag_D4#c.sft_pca Lag_D21#c.sft_pca Lag_D22#c.sft_pca sft_pca $x if m2!=1, cluster(country) frame create orsCI8 or orlb orub lincom _b[sft_pca] frame post orsCI8 (`r(estimate)') (`r(lb)') (`r(ub)') forvalues a =1(1)4 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI8 (`r(estimate)') (`r(lb)') (`r(ub)') } frame post orsCI8 (`r(estimate)') (`r(lb)') (`r(ub)') forvalues a = 21(1)22 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI8 (`r(estimate)') (`r(lb)') (`r(ub)') } frame change orsCI8 format or %3.2f format orlb %3.2f format orub %3.2f save "${output}/support_all_age",replace use "${inter}/Psycorona_long_full",clear cap drop id_wave egen id_wave= tag(id wave) drop if cases_population_7day==. cap drop country encode CountryName,gen(country) forvalues i =0(1)4{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } forvalues i =11(1)22{ cap drop Lag_D`i' gen Lag_D`i'=1 if wave==`i' replace Lag_D`i'=0 if wave!=`i' } global x i.age i.gender i.unemployment i.education i.w_PFS lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 people_vaccinated_per_hundred cases_population_7day stringency_index xtset id wave xtreg w_c19RCA01 Lag_D1#c.sft_pca Lag_D2#c.sft_pca Lag_D3#c.sft_pca Lag_D4#c.sft_pca Lag_D21#c.sft_pca Lag_D22#c.sft_pca sft_pca $x if m2!=1&age<=2, cluster(country) frame create orsCI9 or1 orlb1 orub1 lincom _b[sft_pca] frame post orsCI9 (`r(estimate)') (`r(lb)') (`r(ub)') forvalues a = 1(1)4 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI9 (`r(estimate)') (`r(lb)') (`r(ub)') } forvalues a = 21(1)22 { lincom _b[sft_pca] + _b[1.Lag_D`a'#c.sft_pca] frame post orsCI9 (`r(estimate)') (`r(lb)') (`r(ub)') } frame change orsCI9 format or1 %3.2f format orlb1 %3.2f format orub1 %3.2f save "${output}/support_34",replace use "${output}/support_34",clear append using "${output}/support_all_age" replace or=or[_n+7] replace orlb=orlb[_n+7] replace orub=orub[_n+7] drop if or1==. cap drop wave1 gen wave1=_n cap drop wave2 gen wave2=_n+0.3 graph twoway (scatter or wave1, lcolor("5 101 141") sort msize(vsmall) mfcolor("5 101 141") mlcolor("5 101 141") msymbol(diamond)) (rcap orlb orub wave1,color("5 101 141")) (scatter or1 wave2, sort msize(vsmall) mfcolor("223 82 54") mlcolor("223 82 54") msymbol(diamond) lpattern(dash) lcolor(gray)) (rcap orlb1 orub1 wave2,color("223 82 54") ) (rconnected orlb orub wave1 if wave1<=0, msize(vsmall) mfcolor("5 101 141") mlcolor("5 101 141") msymbol(diamond)) (rconnected orlb1 orub1 wave2 if wave2<=0, msize(vsmall) mfcolor(red) mlcolor("223 82 54") msymbol(diamond) lcolor("223 82 54")), /// legend(order(5 "All age" 6 "18-34" ) col(1) ) ylab(,format(%3.1f)) /// xlab(1 "March 19, 2020" 2 "March 27, 2020" 3 "April 11, 2020" 4 "April 18, 2020" 5 "April 25, 2020" 6 "April 13, 2021" 7 "July 13, 2021", angle(70) labsize(small) labcolor(black) axis(1)) /// subtitle("Panel E. Support for mandatory vaccination") xtitle("") ytitle("Estimated Coefficients") yline(0, lp(dash)) scheme(s1mono)saving("${output}/graph5",replace) grc1leg2 "${output}/graph1.gph" "${output}/graph2.gph" "${output}/graph3.gph" "${output}/graph4.gph" "${output}/graph5.gph", legendfrom("${output}/graph5.gph") pos(5) lyo(30) lxo(-15) ysize(40) xsize(28) colfirst col(2) scheme(s1mono) scale(0.65) graph save "${output}/figure4",replace