***Biostat 212: Assignment 5
**Sarah Patterson

cd "/Users/sarahpatterson/Desktop/TICR Summer 2016/Biostat 212/Assignment 5" //state working directory where files are located


capture log close //close any log files that may be open
log using "Assignment5_Patterson.smcl", replace

set more off           //tells stata not to pause or dipslay the -more- message

version 14             //using version 14


use "Neonatal.dta", clear


***Analysis

***Part 1: Histogram

histogram bw    ///
       , frequency fcolor(teal) lcolor(white) ///
	   normal ///
	   ytitle(, size(large) margin(small)) ///
	   xtitle("Birthweight, grams") xtitle(, size(large) margin(medium)) ///
	   xlabel(500(500)4500) ///
	   bin(20) ///
	   name(HW5_hist, replace) ///
	   saving("HW5_hist.gph", replace)

	   
***Part 2: Dot Plots

dotplot bw ///
    , over(gestagecat) ///
	center ///
	mcolor(blue) msize(small) msymbol(diamond) ///
	ytitle("Birthweight, grams") ///
	ytitle(, size(large) margin(small)) ///
	xtitle("Gestational age") ///
	xtitle(, size(large) margin(medium)) ///
	name(HW5_dot, replace) ///
	saving("HW5_dot.gph", replace)


***Part 3: Box Plots
	 
 graph box bw, ///
     over(sex, relabel(1 "female" 2 "male") label(labsize(medsmall))) ///
	 over(gestagecat) ///
	 ytitle(, size(large) margin(small)) ///
	 ytitle("Birthweight, grams") ///
	 graphregion(fcolor(white) lcolor(none)) ///
	 name(HW5_box, replace) ///
	 saving("HW5_box.gph", replace)

	 
***Part 4: Scatterplots

twoway (lfitci bw gestage, lcolor(black)) ///
    (scatter bw gestage, mcolor(black) msize(vsmall)) ///
	, caption(p<0.001, size(medsmall) position(3) ring(0) ///
	margin(vlarge)) ///
	ytitle(, size(medium) margin(small)) ///
	ytitle("Birthweight, grams") ///
	xtitle("Gestational age") ///
	xtitle(, size(medium) margin(medium)) ///
	graphregion(fcolor(white) lcolor(none)) ///
	name(HW5_scatter, replace) ///
	saving("HW5_scatter.gph", replace)
	
***remove outline of CI, add name/saving
	


	



	   
***

set more on


log close