set more off

cd "/Users/dave/UCSF/Projects/Teaching/PhDadmin/Seminar/"
use "forestdata.dta", clear
/*

	Log transform X-axis

*/
label def stdy ///
	1 "PROUD"  ///
	2 "ipergay"   ///
	3 "Partners (FTC/TDF)"  ///
	4 "Partners (TDF)"  ///
	5 "TDF2"  ///
	6 "Bangkok (TDF)"  ///
	7 "iPrEX"  ///
	8 "VOICE (FTC/TDF)"  ///
	9 "FEMPREP"  ///
	10 "VOICE (TDF)"  
label val var stdy 
label var var " "

eclplot middle upper lower var, hori  estopts(ylabel(1(1)10)) xtitle("ITT Effectiveness") eplottype(connected )  ///
	estopts(bcolor(black) mcolor(black))   ///
	ciopts( lcolor(black) mcolor(black) )   ///
	xlabel(1, grid glcolor(gs6) glwidth(vthin) glpattern(shortdash))  ///
	xlabel(0.05 "95%" 0.10 "90%" 0.25 "75%" 0.50 "50%"  1 "0%"  2 " " 4 " " 8 " " 20 " ")  ///
	xscale(log) scheme(s1color) xline(1, lstyle(dash) lcolor(red) lwidth(thick)  ) legend(off)
	
graph export "forestlog.pdf", replace as(pdf)
/*

	Non-transformed X-axis
	Reversed so highest RR appears on the left

*/
replace upper = -1/(upper)
replace lower = -1/(lower)
replace middle = -1/(middle)

eclplot middle upper lower var, hori  estopts(ylabel(1(1)10)) xtitle("ITT Effectiveness") eplottype(connected ) ///
	estopts(bcolor(black) mcolor(black))  ///
	ciopts( lcolor(black) mcolor(black) )  ///
	xlabel(1, grid glcolor(gs6) glwidth(vthin) glpattern(shortdash)) ///
	xlabel( -20 "95%" -10 "90" -4 "75%"   -1  "0%" 1.0 " ")  ///
	scheme(s1color) xline(-1, lstyle(dash) lcolor(red) lwidth(thick)  ) legend(off) 

graph export "forestnolog.pdf", replace as(pdf)