#delimit ;
clear;
capture log close;
log using "final_graphs v2.log", replace;
set more off;

use "manskibound_raw_conf.dta";

twoway (scatter in_manski days, msize(large))  (scatter in_manski95 days, msize(vsmall)) 
	if group == "all", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "Manski Bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("All Patients")
	saving(manskibound_all_conf95, replace) name(manskibound_all_conf95, replace);

translate @Graph manskibound_all_conf.eps, replace;
	
twoway (scatter in_manski days, msize(large))  (scatter in_manski95 days, msize(vsmall)) 
	if group == "cat11", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "Manski Bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Acute Respiratory Failure")
	saving(manskibound_cat11_conf95, replace) name(manskibound_cat11_conf95, replace);

twoway (scatter in_manski days, msize(large))  (scatter in_manski95 days, msize(vsmall)) 
	if group == "cat12", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "Manski Bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("COPD")
	saving(manskibound_cat12_conf95, replace) name(manskibound_cat12_conf95, replace);

twoway (scatter in_manski days, msize(large))  (scatter in_manski95 days, msize(vsmall)) 
	if group == "cat13", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "Manski Bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Congestive Heart Failure")
	saving(manskibound_cat13_conf95, replace) name(manskibound_cat13_conf95, replace);

twoway (scatter in_manski days, msize(large))  (scatter in_manski95 days, msize(vsmall)) 
	if group == "cat14", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "Manski Bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Cirrhosis")
	saving(manskibound_cat14_conf95, replace) name(manskibound_cat14_conf95, replace);

twoway (scatter in_manski days, msize(large))  (scatter in_manski95 days, msize(vsmall)) 
	if group == "cat15", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "Manski Bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Coma")
	saving(manskibound_cat15_conf95, replace) name(manskibound_cat15_conf95, replace);

twoway (scatter in_manski days, msize(large))  (scatter in_manski95 days, msize(vsmall)) 
	if group == "cat17", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "Manski Bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Lung Cancer")
	saving(manskibound_cat17_conf95, replace) name(manskibound_cat17_conf95, replace);
	
twoway (scatter in_manski days, msize(large))  (scatter in_manski95 days, msize(vsmall)) 
	if group == "cat18", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "Manski Bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("MOSF with Malignancy")
	saving(manskibound_cat18_conf95, replace) name(manskibound_cat18_conf95, replace);
	
twoway (scatter in_manski days, msize(large))  (scatter in_manski95 days, msize(vsmall)) 
	if group == "cat19", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "Manski Bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("MOSF with Sepsis")
	saving(manskibound_cat19_conf95, replace) name(manskibound_cat19_conf95, replace);


graph combine manskibound_cat11_conf95 manskibound_cat13_conf95 manskibound_cat18_conf95 manskibound_cat19_conf95, 
	saving(manskibound_dx1_conf95, replace);

translate @Graph manskibound_dx1_conf95.eps, replace;

graph combine manskibound_cat12_conf95 manskibound_cat14_conf95 manskibound_cat15_conf95, 
	saving(manskibound_dx2_conf95, replace);

translate @Graph manskibound_dx2_conf95.eps, replace;

use "svbound_raw_conf2.dta";

twoway (scatter in_sv days, msize(large))  (scatter in_sv95 days, msize(vsmall)) 
	if group == "all", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "SV bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("All Patients")
	saving(svbound_all_conf95, replace) name(svbound_all_conf95, replace);

translate @Graph svbound_all_conf.eps, replace;
	
twoway (scatter in_sv days, msize(large))  (scatter in_sv95 days, msize(vsmall)) 
	if group == "cat11", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "SV bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Acute Respiratory Failure")
	saving(svbound_cat11_conf95, replace) name(svbound_cat11_conf95, replace);

twoway (scatter in_sv days, msize(large))  (scatter in_sv95 days, msize(vsmall)) 
	if group == "cat12", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "SV bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("COPD")
	saving(svbound_cat12_conf95, replace) name(svbound_cat12_conf95, replace);

twoway (scatter in_sv days, msize(large))  (scatter in_sv95 days, msize(vsmall)) 
	if group == "cat13", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "SV bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Congestive Heart Failure")
	saving(svbound_cat13_conf95, replace) name(svbound_cat13_conf95, replace);

twoway (scatter in_sv days, msize(large))  (scatter in_sv95 days, msize(vsmall)) 
	if group == "cat14", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "SV bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Cirrhosis")
	saving(svbound_cat14_conf95, replace) name(svbound_cat14_conf95, replace);

twoway (scatter in_sv days, msize(large))  (scatter in_sv95 days, msize(vsmall)) 
	if group == "cat15", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "SV bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Coma")
	saving(svbound_cat15_conf95, replace) name(svbound_cat15_conf95, replace);

twoway (scatter in_sv days, msize(large))  (scatter in_sv95 days, msize(vsmall)) 
	if group == "cat17", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "SV bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Lung Cancer")
	saving(svbound_cat17_conf95, replace) name(svbound_cat17_conf95, replace);
	
twoway (scatter in_sv days, msize(large))  (scatter in_sv95 days, msize(vsmall)) 
	if group == "cat18", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "SV bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("MOSF with Malignancy")
	saving(svbound_cat18_conf95, replace) name(svbound_cat18_conf95, replace);
	
twoway (scatter in_sv days, msize(large))  (scatter in_sv95 days, msize(vsmall)) 
	if group == "cat19", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "SV bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("MOSF with Sepsis")
	saving(svbound_cat19_conf95, replace) name(svbound_cat19_conf95, replace);

graph combine svbound_cat11_conf95 svbound_cat13_conf95 svbound_cat18_conf95 svbound_cat19_conf95, 
	saving(svbound_dx1_conf95, replace);

translate @Graph svbound_dx1_conf95.eps, replace;

graph combine svbound_cat12_conf95 svbound_cat14_conf95 svbound_cat15_conf95, 
	saving(svbound_dx2_conf95, replace);

translate @Graph svbound_dx2_conf95.eps, replace;

use "pqdbound_raw_conf.dta";

twoway (scatter in_pqd days, msize(large))  (scatter in_pqd95 days, msize(vsmall)) 
	if group == "all", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "PQD bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("All Patients")
	saving(pqdbound_all_conf95, replace) name(pqdbound_all_conf95, replace);

translate @Graph pqdbound_all_conf.eps, replace;
	
twoway (scatter in_pqd days, msize(large))  (scatter in_pqd95 days, msize(vsmall)) 
	if group == "cat11", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "PQD bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Acute Respiratory Failure")
	saving(pqdbound_cat11_conf95, replace) name(pqdbound_cat11_conf95, replace);

twoway (scatter in_pqd days, msize(large))  (scatter in_pqd95 days, msize(vsmall)) 
	if group == "cat12", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "PQD bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("COPD")
	saving(pqdbound_cat12_conf95, replace) name(pqdbound_cat12_conf95, replace);

twoway (scatter in_pqd days, msize(large))  (scatter in_pqd95 days, msize(vsmall)) 
	if group == "cat13", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "PQD bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Congestive Heart Failure")
	saving(pqdbound_cat13_conf95, replace) name(pqdbound_cat13_conf95, replace);

twoway (scatter in_pqd days, msize(large))  (scatter in_pqd95 days, msize(vsmall)) 
	if group == "cat14", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "PQD bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Cirrhosis")
	saving(pqdbound_cat14_conf95, replace) name(pqdbound_cat14_conf95, replace);

twoway (scatter in_pqd days, msize(large))  (scatter in_pqd95 days, msize(vsmall)) 
	if group == "cat15", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "PQD bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Coma")
	saving(pqdbound_cat15_conf95, replace) name(pqdbound_cat15_conf95, replace);

twoway (scatter in_pqd days, msize(large))  (scatter in_pqd95 days, msize(vsmall)) 
	if group == "cat17", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "PQD bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("Lung Cancer")
	saving(pqdbound_cat17_conf95, replace) name(pqdbound_cat17_conf95, replace);
	
twoway (scatter in_pqd days, msize(large))  (scatter in_pqd95 days, msize(vsmall)) 
	if group == "cat18", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "PQD bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("MOSF with Malignancy")
	saving(pqdbound_cat18_conf95, replace) name(pqdbound_cat18_conf95, replace);
	
twoway (scatter in_pqd days, msize(large))  (scatter in_pqd95 days, msize(vsmall)) 
	if group == "cat19", 
	yscale(range(-1 1)) ylabel(-1(0.25)1, angle(horizontal) format(%2.1f)) 
	ytitle("Tx. Effect Bound") xtitle("Days After ICU Admission", margin(medsmall))
	legend(order(1 "PQD bound" 2 "95% conf. bound")) 
	scheme(s1mono) yline(0) title("MOSF with Sepsis")
	saving(pqdbound_cat19_conf95, replace) name(pqdbound_cat19_conf95, replace);

graph combine pqdbound_cat11_conf95 pqdbound_cat13_conf95 pqdbound_cat18_conf95 pqdbound_cat19_conf95, 
	saving(pqdbound_dx1_conf95, replace);

translate @Graph pqdbound_dx1_conf95.eps, replace;

graph combine pqdbound_cat12_conf95 pqdbound_cat14_conf95 pqdbound_cat15_conf95, 
	saving(pqdbound_dx2_conf95, replace);

translate @Graph pqdbound_dx2_conf95.eps, replace;