**Lecture 7 

cd "/Users/kristenaiemjoy/Box Sync/_Teaching/Biostat 212/Biostat212.2018/Lectures/7.Programming"


*Open NHANES data
webuse nhanes2, clear 


*Stored results
*r-class: sum, tabulate
 return list
*n-class: dont have any stored results 
*e-class: regressions 
ereturn list 


sum age 
return list 


sum height
display "SE =" %9.2f r(sd)/sqrt(r(N))

tab sex heartatk 

gen sex2 = sex
recode sex2 2=0

label define sexlabel 0 "Female" 1 "male", replace 
label values sex2 sexlabel 

tab sex2


tab sex heartatk, row chi2
return list

display "p-value =" %9.6f r(p)

display r(p)

logit heartatk sex2, or


logistic heartatk sex2
ereturn list 

di "OR = " exp(_b[sex])

***MACROS
*local
*global (EVERYWHERE)
*scalar - more numeric precision 

regress bpsystol age 
ereturn list 

local rsq1 e(r2) 
di `rsq1'


local rsq2 = e(r2)
di `rsq2'

regress bpsystol age weight 

di `rsq1'

di `rsq2'

**Global 

global myName "Kristen" 
di "$myName"

global myName "Kristen Aiemjoy" 



**Scalars

sum age 

scalar SE = r(sd)/sqrt(r(N))
di SE 



***PROGRAM 

program define hello
	di "Hello Kristen" 
end


program define hello
	di "Hello World" 
end


program drop median 

program define median, rclass 
	quietly sum `1', detail 
	local median = r(p50)
	display "median of `1' = " %9.3f `median'
end



program define MySumisCool, rclass
	quietly sum `1', detail 
	local median = r(p50)
	local mean = r(mean)
	display "mean = " %9.2f `mean'
	display "median = " %9.2f `median'
	di "Biostat 212 is the best class ever" 
end


***LOOPS 

foreach var of varlist iron weight height bmi bpsystol {
    median `var'
	}



***put excel 
tabstat zinc copper vitaminc iron, by(heartatk) stat(n mean sd) save  
return list
matrix list r(StatTotal) 


matrix results1 = r(StatTotal)


putexcel set lecture7.xlsx, replace 

putexcel A1 = matrix(results1), names