clear
clear matrix
cd "${output_path}\tmp\"

use "placebo_simulation.dta", clear
graph drop _all

di beta_fathers
count if beta_fathers>=0.166 | beta_fathers<=(-0.166)



di beta_mothers
count if beta_fathers>=0.114 | beta_fathers<=(-0.114)



twoway ///
(kdensity beta_fathers) , xline(0.166) xscale(range(-0.2 0.2)) xlabel(-0.2(0.1)0.2) ytitle("Density") xtitle("Estimated effect size")  subtitle("a. Fathers") scheme(s1mono) ///
text(9.5 0.16 "Observed value:" "0.166", placement("west")) name(fathers)

twoway ///
(kdensity beta_mothers) , xline(0.114) xscale(range(-0.2 0.2)) xlabel(-0.2(0.1)0.2) ytitle("Density") xtitle("Estimated effect size") subtitle("b. Mothers") scheme(s1mono) ///
text(7.5 0.11 "Observed value:" "0.114", placement("west")) name(mothers)

graph combine fathers mothers, altshrink scheme(s1mono) 
graph export "${output_path}\Results\FigA11.tif", replace width(2400)


