use "${inter}/Psycorona_long_sft2",clear
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
duplicates drop CountryName,force
drop if CountryName=="Hong Kong S.A.R."|CountryName=="Taiwan"
keep CountryName sft_pca
cap drop m66
merge 1:n CountryName using "${inter}/vaccination",gen(m66) 
drop if m66==2
cap drop mean
bysort edate:egen mean=mean(stringency_index)
cap drop zero
gen zero=0
//Figure S3:
twoway (line mean edate if mean!=.,lwidth(vthick) mlcolor("46 134 178")) (rarea mean zero edate if mean!=.,color("222 236 243")) (scatteri  0  22474 80 22474 ,lcolor("155 58 74") recast(line)),   legend(off)  xlab(21917 22028 22139 22250 22361 22474 22584 22695 22806 22917 23000, angle(70) labsize(small) labcolor(black) axis(1)) ytitle("average stringency index", size(small)) ylab(,labsize(small)) xtitle("")  graphregion(color(white))
graph save "${output}/stringency1",replace

cap drop median
bysort edate:egen median=median(stringency_index)

twoway (line median edate if mean!=.,lwidth(vthick) mlcolor("46 134 178")) (rarea median zero edate if mean!=.,color("222 236 243")) (scatteri  0  22474 80 22474 , lcolor("155 58 74") recast(line)),   legend(off)  xlab(21917 22028 22139 22250 22361 22474 22584 22695 22806 22917 23000, angle(70) labsize(small) labcolor(black) axis(1)) ytitle("median of stringency index", size(small)) ylab(,labsize(small)) xtitle("")  graphregion(color(white))

graph save "${output}/stringency2",replace

graph combine "${output}/stringency1.gph"  "${output}/stringency2.gph" , col(1)  ysize(14) xsize(9) scheme(s1mono)