* do "~/Documents/teaching/c2015/biostat208/labs/lab10/lab10.do"
cd "~/Documents/teaching/c2015/biostat208/labs/lab10/"
use "lab10.dta", clear
log using "lab10", replace


drop if chol==. | chol==645
drop if bmi<12

logistic chd69 chol sbp dibpat age smoke bmi, coef
predict phat
predict dxb, dbeta
list id chd69 chol sbp dibpat age smoke bmi if id==10078
gsort -dxb
list id chd69 chol sbp dibpat age smoke bmi if _n<=10


logistic chd69 chol sbp dibpat age smoke bmi if id != 10078, coef

logistic chd69 chol sbp dibpat age smoke bmi
estat gof, table group(10)

summ bmi
gen bmic = bmi - r(mean)
summ chol
gen cholc = chol - r(mean)
summ sbp
gen sbpc = sbp - r(mean)
logistic chd69 dibpat age smoke c.bmic##c.cholc c.sbpc##c.bmic
estat gof, table group(10)

quietly regress chd69 dibpat age smoke c.bmic##c.cholc c.sbpc##c.bmic
estat vif

set seed 386282911
gen r = uniform()
sort r
gen group = 1 if _n <= _N*0.9
replace group=2 if group!=1
tabulate group, summarize(chd69) nostandard 
logistic chd69 dibpat age smoke c.bmic##c.cholc c.sbpc##c.bmic if group==1
predict pht if group==2
gsort -group

list group pht chol sbp dibpat age smoke bmi if _n<=10

lroc if group==2
lsens if group==2

lstat if group==2, cutoff(0.1)
lstat if group==2, cutoff(0.2)

capture binreg chd69 dibpat age smoke c.bmic##c.cholc c.sbpc##c.bmic, rr iterate(100)
poisson chd69 dibpat age smoke c.bmic##c.cholc c.sbpc##c.bmic, irr robust

log close