use dta/diuretic, clear gen logor=log((rt/(nt-rt))/(rc/(nc-rc))) gen selogor=sqrt((1/rc)+(1/(nc-rc))+(1/rt)+(1/(nt-rt))) meta logor selogor, ef meta logor selogor, ef graph(r) id(trialid) cline xlab(0.5,1,1.5) xline(1) /* */ boxsh(4) b2("Odds ratio - log scale") saving(sharp1.gph,replace) meta logor selogor, ef graph(e) id(trialid) cline xlab(0.5,1,1.5) xline(1) /* */ boxsh(2) b2("Odds ratio - log scale") saving(sharp2.gph,replace)