cd "/Users/mmehrotra/Documents/biostat212/2019/final"

capture log close 


quietly log using biostat212_FinalProjectExample.smcl, replace

/***
#Biostat 212: 2019
#Example Final Project 

***/

webuse nhanes2, clear

/*
TABLE 1 Code
*/

label define diabetes_label 1 "Diabetics" 0 "Non diabetics"
label values diabetes diabetes_label
rename vitaminc VitaminC
rename zinc Zinc
rename copper Copper
rename iron Iron
rename diabetes Diabetes

tabstat Zinc Copper VitaminC Iron, by(Diabetes) stat(n mean sd) format(%10.2f) save
return list
matrix list r(StatTotal)
matrix nondiabetics = r(Stat1)
matrix diabetics = r(Stat2)

/*
	Format these matrices
*/
matrix nondiabetics=nondiabetics' //transpose
matrix diabetics=diabetics'
matrix table1=nondiabetics,diabetics



/*
TABLE 1
*/
putdocx clear
putdocx begin
putdocx table table1=matrix(table1), nformat(%9.2f)rownames colnames title(Table 1: Serum Zinc, Copper, Vitamin C and Iron levels (mcg/DL) for Diabetics and Non-Diabetics in NHANES)
putdocx save "output/Example Final Project.docx",replace

	
	
	
	
	
clear
webuse nhanes2
label define diabetes_label 1 "Diabetics" 0 "Non diabetics"
label values diabetes diabetes_label
	
/*
FIGURE 1 Code
*/

		
twoway ///	
		kdensity iron if diabetes==1 ///	
	,	lcolor(gray) || ///
	 	kdensity iron if diabetes==0 ///
	, 	lcolor(black) lpattern(dash) ///
	, 	yscale(off) ///
		graphregion(color(white) lcolor(none) ilcolor(none)) ///
		ytitle("Density") xtitle("Serum Iron (mcg/dL)") ///
		legend(order(1 "Diabetics" 2 "Non-Diabetics")) ///
		name(iron_plot, replace) 	///
		subtitle("Figure 1: K-density Distribution of Serum Iron levels in Diabetics and Non-diabetics") ///
	note("Data source: NHANES2") ///
	scheme(s1mono)

/*
FIGURE 1
*/
putdocx begin
graph export "output/fig1.png", replace
putdocx paragraph, halign(center)
putdocx image "output/fig1.png"
putdocx save "output/Example Final Project.docx", append

log close