frame reset use "${inter}/PsyCorona",clear cap drop m22 merge n:1 CountryName using "${inter}/family_importance", gen(m22) drop if m22==2 cap drop m23 merge n:1 CountryName using "${inter}/parent_proud", gen(m23) drop if m23==2 cap drop m24 merge n:1 CountryName using "${inter}/respect_love", gen(m24) drop if m24==2 cap drop m25 merge n:1 CountryName using "${inter}/parental_duties", gen(m25) drop if m25==2 autoclean egen missing1=rmiss(w_c19perBeh01_Harmonized w_c19perBeh02_Harmonized w_c19perBeh03_Harmonized) drop if missing1>0 replace w_Vaccine=w_Vaccine+3 replace w_c19RCA01=w_c19RCA01+4 //add w_c19RCA01 var cap drop m2 merge 1:1 id edate using "${inter}/w_c19RCA01",gen(m2) update //add new cases/total population cap drop m3 merge n:1 CountryName edate using "${inter}/case",gen(m3) update drop if m3==2 //add vaccination variable cap drop m1 merge n:1 CountryName edate using "${inter}/vaccination",gen(m1) update drop if m1==2 cap drop miss egen miss=rmiss(w_c19perBeh01_Harmonized w_c19perBeh02_Harmonized w_c19perBeh03_Harmonized age gender unemployment edu w_PFS01) drop if miss>0 replace w_PFS01=w_PFS01+2 replace gender=0 if gender!=1 cap drop lnCOUNTRYGDPpc gen lnCOUNTRYGDPpc=ln(COUNTRYGDPpc) cap drop country encode CountryName,gen(country) foreach k in w_c19perBeh01_Harmonized w_c19perBeh02_Harmonized w_c19perBeh03_Harmonized w_Vaccine w_c19RCA01 lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 cases_population_7day stringency_index people_vaccinated_per_hundred family_importance parent_proud respect_love parental_duties { cap drop std`k' egen std`k'=std(`k') replace `k'=std`k' } xtset id wave save "${inter}/Psycorona_long_each_question",replace global x i.age i.gender i.unemployment i.education w_PFS lnCOUNTRYGDPpc COUNTRYOut_of_pocket2019 people_vaccinated_per_hundred cases_population_7day stringency_index use "${inter}/Psycorona_long_each_question",clear foreach k in w_c19perBeh01_Harmonized w_c19perBeh02_Harmonized w_c19perBeh03_Harmonized w_Vaccine { xtreg `k' family_importance $x if m2!=2, vce(cluster country) xtreg `k' parent_proud $x if m2!=2, vce(cluster country) xtreg `k' respect_love $x if m2!=2, vce(cluster country) xtreg `k' parental_duties $x if m2!=2, vce(cluster country) xtreg `k' family_importance $x if age<=5&m2!=2, vce(cluster country) xtreg `k' parent_proud $x if age<=5&m2!=2, vce(cluster country) xtreg `k' respect_love $x if age<=5&m2!=2, vce(cluster country) xtreg `k' parental_duties $x if age<=5&m2!=2, vce(cluster country) xtreg `k' family_importance $x if age<=2&m2!=2, vce(cluster country) xtreg `k' parent_proud $x if age<=2&m2!=2, vce(cluster country) xtreg `k' respect_love $x if age<=2&m2!=2, vce(cluster country) xtreg `k' parental_duties $x if age<=2&m2!=2, vce(cluster country) } xtreg w_c19RCA01 family_importance $x if m2!=1, vce(cluster country) xtreg w_c19RCA01 parent_proud $x if m2!=1, vce(cluster country) xtreg w_c19RCA01 respect_love $x if m2!=1, vce(cluster country) xtreg w_c19RCA01 parental_duties $x if m2!=1, vce(cluster country) xtreg w_c19RCA01 family_importance $x if age<=5&m2!=1, vce(cluster country) xtreg w_c19RCA01 parent_proud $x if age<=5&m2!=1, vce(cluster country) xtreg w_c19RCA01 respect_love $x if age<=5&m2!=1, vce(cluster country) xtreg w_c19RCA01 parental_duties $x if age<=5&m2!=1, vce(cluster country) xtreg w_c19RCA01 family_importance $x if age<=2&m2!=1, vce(cluster country) xtreg w_c19RCA01 parent_proud $x if age<=2&m2!=1, vce(cluster country) xtreg w_c19RCA01 respect_love $x if age<=2&m2!=1, vce(cluster country) xtreg w_c19RCA01 parental_duties $x if age<=2&m2!=1, vce(cluster country) //plot //for all age //parent proud use "${inter}/Psycorona_long_each_question",clear frame create panela1 coef1 lb1 ub1 xtreg w_c19perBeh01_Harmonized parent_proud $x if m2!=2, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized parent_proud $x if m2!=2, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized parent_proud $x if m2!=2, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine parent_proud $x if m2!=2, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 parent_proud $x if m2!=1, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela1 format coef1 %3.2f format lb1 %3.2f format ub1 %3.2f cap drop wave gen wave=_n save "${inter}/panela1_proud",replace //respect_love use "${inter}/Psycorona_long_each_question",clear frame create panela2 coef2 lb2 ub2 xtreg w_c19perBeh01_Harmonized respect_love $x if m2!=2, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized respect_love $x if m2!=2, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized respect_love $x if m2!=2, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine respect_love $x if m2!=2, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 respect_love $x if m2!=1, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela2 format coef2 %3.2f format lb2 %3.2f format ub2 %3.2f cap drop wave gen wave=_n save "${inter}/panela2_respect",replace //family importance use "${inter}/Psycorona_long_each_question",clear frame create panela3 coef3 lb3 ub3 xtreg w_c19perBeh01_Harmonized family_importance $x if m2!=2, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized family_importance $x if m2!=2, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized family_importance $x if m2!=2, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine family_importance $x if m2!=2, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 family_importance $x if m2!=1, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela3 format coef3 %3.2f format lb3 %3.2f format ub3 %3.2f cap drop wave gen wave=_n save "${inter}/panela3_family",replace //parental_duties use "${inter}/Psycorona_long_each_question",clear frame create panela4 coef4 lb4 ub4 xtreg w_c19perBeh01_Harmonized parental_duties $x if m2!=2, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized parental_duties $x if m2!=2, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized parental_duties $x if m2!=2, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine parental_duties $x if m2!=2, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 parental_duties $x if m2!=1, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela4 format coef4 %3.2f format lb4 %3.2f format ub4 %3.2f cap drop wave gen wave=_n save "${inter}/panela4_duties",replace use "${inter}/panela1_proud",clear merge 1:1 wave using "${inter}/panela2_respect", nogen merge 1:1 wave using "${inter}/panela3_family", nogen merge 1:1 wave using "${inter}/panela4_duties", nogen replace wave=wave+1 if wave>1 replace wave=wave+1 if wave>3 replace wave=wave+1 if wave>5 replace wave=wave+1 if wave>7 rename wave label3 cap drop label1 gen label1=label3-0.6 cap drop label2 gen label2=label3-0.3 cap drop label4 gen label4=label3+0.3 graph twoway /// (scatter label1 coef1 , lcolor("21 110 144") sort msize(vsmall) mfcolor("21 110 144") mlcolor("21 110 144") msymbol(diamond)) (rconnected lb1 ub1 label1 if label1<=-10, msize(vsmall) mfcolor("21 110 144") mlcolor("21 110 144") msymbol(diamond)) (rcap lb1 ub1 label1,color("21 110 144") horizontal) /// (scatter label2 coef2 , sort msize(vsmall) mfcolor("34 181 191") mlcolor("34 181 191") msymbol(diamond) lpattern(dash) lcolor(gray)) (rconnected lb2 ub2 label2 if label2<=-10, msize(vsmall) mfcolor("34 181 191") mlcolor("34 181 191") msymbol(diamond) lcolor("34 181 191") ) (rcap lb2 ub2 label2,color("34 181 191") horizontal) /// (scatter label3 coef3, msize(vsmall) mfcolor("171 185 188") mlcolor("171 185 188") msymbol(diamond) lpattern(dash) lcolor(gray)) (rconnected lb3 ub3 label3 if label3<=-10, msize(vsmall) mfcolor("171 185 188") mlcolor("171 185 188") msymbol(diamond) lcolor("171 185 188") )(rcap lb3 ub3 label3,color("171 185 188") horizontal) /// (scatter label4 coef4, msize(vsmall) mfcolor("229 75 47") mlcolor("229 75 47") msymbol(diamond) lpattern(dash) lcolor(gray)) (rconnected lb4 ub4 label4 if label4<=-10, msize(vsmall) mfcolor("229 75 47") mlcolor("229 75 47") msymbol(diamond) lcolor("229 75 47") )(rcap lb4 ub4 label4,color("229 75 47") horizontal) , /// xlab(-0.1(0.1)0.3,labsize(small) format(%3.1f)) ylab(1 "Hand-washing frequency" 3"Avoid crowds" 5 "Self-quarantine" 7 "Vaccine for self" 9 "Mandatory vaccination", angle(0) labsize(small) labcolor(black) axis(1)) xline(0, lp(dash)) xtitle("Estimated Coefficients",size(small)) subtitle("Panel A. All age") scheme(s1mono) legend(order(3 "Parent proud " 6 "Respect and love parents" 9 "Family important" 12 "Parental duties" ) on size(small) row(2)) yscale(reverse) saving("${output}/graph11",replace) ////for 18-64 frame reset //parent proud use "${inter}/Psycorona_long_each_question",clear frame create panela1 coef1 lb1 ub1 xtreg w_c19perBeh01_Harmonized parent_proud $x if m2!=2&age<=5, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized parent_proud $x if m2!=2&age<=5, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized parent_proud $x if m2!=2&age<=5, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine parent_proud $x if m2!=2&age<=5, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 parent_proud $x if m2!=1&age<=5, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela1 format coef1 %3.2f format lb1 %3.2f format ub1 %3.2f cap drop wave gen wave=_n save "${inter}/panelb1_proud",replace //respect_love use "${inter}/Psycorona_long_each_question",clear frame create panela2 coef2 lb2 ub2 xtreg w_c19perBeh01_Harmonized respect_love $x if m2!=2&age<=5, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized respect_love $x if m2!=2&age<=5, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized respect_love $x if m2!=2&age<=5, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine respect_love $x if m2!=2&age<=5, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 respect_love $x if m2!=1&age<=5, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela2 format coef2 %3.2f format lb2 %3.2f format ub2 %3.2f cap drop wave gen wave=_n save "${inter}/panelb2_respect",replace //family importance use "${inter}/Psycorona_long_each_question",clear frame create panela3 coef3 lb3 ub3 xtreg w_c19perBeh01_Harmonized family_importance $x if m2!=2&age<=5, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized family_importance $x if m2!=2&age<=5, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized family_importance $x if m2!=2&age<=5, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine family_importance $x if m2!=2&age<=5, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 family_importance $x if m2!=1&age<=5, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela3 format coef3 %3.2f format lb3 %3.2f format ub3 %3.2f cap drop wave gen wave=_n save "${inter}/panelb3_family",replace //parental_duties use "${inter}/Psycorona_long_each_question",clear frame create panela4 coef4 lb4 ub4 xtreg w_c19perBeh01_Harmonized parental_duties $x if m2!=2&age<=5, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized parental_duties $x if m2!=2&age<=5, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized parental_duties $x if m2!=2&age<=5, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine parental_duties $x if m2!=2&age<=5, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 parental_duties $x if m2!=1&age<=5, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela4 format coef4 %3.2f format lb4 %3.2f format ub4 %3.2f cap drop wave gen wave=_n save "${inter}/panelb4_duties",replace use "${inter}/panelb1_proud",clear merge 1:1 wave using "${inter}/panelb2_respect", nogen merge 1:1 wave using "${inter}/panelb3_family", nogen merge 1:1 wave using "${inter}/panelb4_duties", nogen replace wave=wave+1 if wave>1 replace wave=wave+1 if wave>3 replace wave=wave+1 if wave>5 replace wave=wave+1 if wave>7 rename wave label3 cap drop label1 gen label1=label3-0.6 cap drop label2 gen label2=label3-0.3 cap drop label4 gen label4=label3+0.3 graph twoway /// (scatter label1 coef1 , lcolor("21 110 144") sort msize(vsmall) mfcolor("21 110 144") mlcolor("21 110 144") msymbol(diamond)) (rconnected lb1 ub1 label1 if label1<=-10, msize(vsmall) mfcolor("21 110 144") mlcolor("21 110 144") msymbol(diamond)) (rcap lb1 ub1 label1,color("21 110 144") horizontal) /// (scatter label2 coef2 , sort msize(vsmall) mfcolor("34 181 191") mlcolor("34 181 191") msymbol(diamond) lpattern(dash) lcolor(gray)) (rconnected lb2 ub2 label2 if label2<=-10, msize(vsmall) mfcolor("34 181 191") mlcolor("34 181 191") msymbol(diamond) lcolor("34 181 191") ) (rcap lb2 ub2 label2,color("34 181 191") horizontal) /// (scatter label3 coef3, msize(vsmall) mfcolor("171 185 188") mlcolor("171 185 188") msymbol(diamond) lpattern(dash) lcolor(gray)) (rconnected lb3 ub3 label3 if label3<=-10, msize(vsmall) mfcolor("171 185 188") mlcolor("171 185 188") msymbol(diamond) lcolor("171 185 188") )(rcap lb3 ub3 label3,color("171 185 188") horizontal) /// (scatter label4 coef4, msize(vsmall) mfcolor("229 75 47") mlcolor("229 75 47") msymbol(diamond) lpattern(dash) lcolor(gray)) (rconnected lb4 ub4 label4 if label4<=-10, msize(vsmall) mfcolor("229 75 47") mlcolor("229 75 47") msymbol(diamond) lcolor("229 75 47") )(rcap lb4 ub4 label4,color("229 75 47") horizontal) , /// xlab(-0.1(0.1)0.3,labsize(small) format(%3.1f)) ylab(1 "Hand-washing frequency" 3"Avoid crowds" 5 "Self-quarantine" 7 "Vaccine for self" 9 "Mandatory vaccination", angle(0) labsize(small) labcolor(black) axis(1)) xline(0, lp(dash)) xtitle("Estimated Coefficients",size(small)) subtitle("Panel B. 18-64") scheme(s1mono) legend(order(3 "Parent proud " 6 "Respect and love parents" 9 "Family important" 12 "Parental duties" ) on size(small) row(2)) yscale(reverse) saving("${output}/graph22",replace) ////for 18-34 frame reset use "${inter}/Psycorona_long_each_question",clear frame create panela1 coef1 lb1 ub1 xtreg w_c19perBeh01_Harmonized parent_proud $x if m2!=2&age<=2, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized parent_proud $x if m2!=2&age<=2, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized parent_proud $x if m2!=2&age<=2, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine parent_proud $x if m2!=2&age<=2, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 parent_proud $x if m2!=1&age<=2, vce(cluster country) lincom _b[parent_proud], l(90) frame post panela1 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela1 format coef1 %3.2f format lb1 %3.2f format ub1 %3.2f cap drop wave gen wave=_n save "${inter}/panelc1_proud",replace //respect_love use "${inter}/Psycorona_long_each_question",clear frame create panela2 coef2 lb2 ub2 xtreg w_c19perBeh01_Harmonized respect_love $x if m2!=2&age<=2, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized respect_love $x if m2!=2&age<=2, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized respect_love $x if m2!=2&age<=2, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine respect_love $x if m2!=2&age<=2, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 respect_love $x if m2!=1&age<=2, vce(cluster country) lincom _b[respect_love], l(90) frame post panela2 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela2 format coef2 %3.2f format lb2 %3.2f format ub2 %3.2f cap drop wave gen wave=_n save "${inter}/panelc2_respect",replace //family importance use "${inter}/Psycorona_long_each_question",clear frame create panela3 coef3 lb3 ub3 xtreg w_c19perBeh01_Harmonized family_importance $x if m2!=2&age<=2, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized family_importance $x if m2!=2&age<=2, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized family_importance $x if m2!=2&age<=2, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine family_importance $x if m2!=2&age<=2, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 family_importance $x if m2!=1&age<=2, vce(cluster country) lincom _b[family_importance], l(90) frame post panela3 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela3 format coef3 %3.2f format lb3 %3.2f format ub3 %3.2f cap drop wave gen wave=_n save "${inter}/panelc3_family",replace //parental_duties use "${inter}/Psycorona_long_each_question",clear frame create panela4 coef4 lb4 ub4 xtreg w_c19perBeh01_Harmonized parental_duties $x if m2!=2&age<=2, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh02_Harmonized parental_duties $x if m2!=2&age<=2, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19perBeh03_Harmonized parental_duties $x if m2!=2&age<=2, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_Vaccine parental_duties $x if m2!=2&age<=2, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') xtreg w_c19RCA01 parental_duties $x if m2!=1&age<=2, vce(cluster country) lincom _b[parental_duties], l(90) frame post panela4 (`r(estimate)') (`r(lb)') (`r(ub)') frame change panela4 format coef4 %3.2f format lb4 %3.2f format ub4 %3.2f cap drop wave gen wave=_n save "${inter}/panelc4_duties",replace use "${inter}/panelc1_proud",clear merge 1:1 wave using "${inter}/panelc2_respect", nogen merge 1:1 wave using "${inter}/panelc3_family", nogen merge 1:1 wave using "${inter}/panelc4_duties", nogen replace wave=wave+1 if wave>1 replace wave=wave+1 if wave>3 replace wave=wave+1 if wave>5 replace wave=wave+1 if wave>7 rename wave label3 cap drop label1 gen label1=label3-0.6 cap drop label2 gen label2=label3-0.3 cap drop label4 gen label4=label3+0.3 graph twoway /// (scatter label1 coef1 , lcolor("21 110 144") sort msize(vsmall) mfcolor("21 110 144") mlcolor("21 110 144") msymbol(diamond)) (rconnected lb1 ub1 label1 if label1<=-10, msize(vsmall) mfcolor("21 110 144") mlcolor("21 110 144") msymbol(diamond)) (rcap lb1 ub1 label1,color("21 110 144") horizontal) /// (scatter label2 coef2 , sort msize(vsmall) mfcolor("34 181 191") mlcolor("34 181 191") msymbol(diamond) lpattern(dash) lcolor(gray)) (rconnected lb2 ub2 label2 if label2<=-10, msize(vsmall) mfcolor("34 181 191") mlcolor("34 181 191") msymbol(diamond) lcolor("34 181 191") ) (rcap lb2 ub2 label2,color("34 181 191") horizontal) /// (scatter label3 coef3, msize(vsmall) mfcolor("171 185 188") mlcolor("171 185 188") msymbol(diamond) lpattern(dash) lcolor(gray)) (rconnected lb3 ub3 label3 if label3<=-10, msize(vsmall) mfcolor("171 185 188") mlcolor("171 185 188") msymbol(diamond) lcolor("171 185 188") )(rcap lb3 ub3 label3,color("171 185 188") horizontal) /// (scatter label4 coef4, msize(vsmall) mfcolor("229 75 47") mlcolor("229 75 47") msymbol(diamond) lpattern(dash) lcolor(gray)) (rconnected lb4 ub4 label4 if label4<=-10, msize(vsmall) mfcolor("229 75 47") mlcolor("229 75 47") msymbol(diamond) lcolor("229 75 47") )(rcap lb4 ub4 label4,color("229 75 47") horizontal) , /// xlab(-0.1(0.1)0.3,labsize(small) format(%3.1f)) ylab(1 "Hand-washing frequency" 3"Avoid crowds" 5 "Self-quarantine" 7 "Vaccine for self" 9 "Mandatory vaccination", angle(0) labsize(small) labcolor(black) axis(1)) xline(0, lp(dash)) xtitle("Estimated Coefficients",size(small)) subtitle("Panel C. 18-34") scheme(s1mono) legend(order(3 "Parent proud " 6 "Respect and love parents" 9 "Family important" 12 "Parental duties" ) on size(small) row(2)) yscale(reverse) saving("${output}/graph33",replace) grc1leg2 "${output}/graph11.gph" "${output}/graph22.gph" "${output}/graph33.gph" , legendfrom("${output}/graph33.gph") pos(6)colfirst col(3) scheme(s1mono) ysize(28) xsize(46) scale(0.7) graph save "${output}/figure",replace