quietly:                                                                        ///
infile ptid mridate age male race weight height packyrs yrsquit alcoh physact   ///
       chf chd stroke diabetes genhlth ldl alb crt plt sbp aai fev dsst atrophy ///
	   whgrd numinf volinf obstime death                                        ///
	   using http://www.emersonstatistics.com/datasets/mri.txt

g deadin5 = 0
replace deadin5 = 1 if obstime <= 5 * 365.25

g alivefor5 = 1 - deadin5

recode ldl 160/max=1 min/160=0, gen(hiLDL)
g loLDL = 1 - hiLDL


logit deadin5 hiLDL
logistic deadin5 hiLDL
logit deadin5 hiLDL, robust

logit deadin5 loLDL
logit alivefor5 hiLDL
logit alivefor5 loLDL

logit hiLDL deadin5

regress deadin5 hiLDL, robust
regress deadin5 hiLDL

regress deadin5 loLDL, robust
regress alivefor5 hiLDL, robust
regress alivefor5 loLDL, robust

regress hiLDL deadin5, robust

poisson deadin5 hiLDL, robust
poisson deadin5 hiLDL

poisson deadin5 loLDL, robust
poisson alivefor5 hiLDL, robust
poisson alivefor5 loLDL, robust

poisson hiLDL deadin5, robust

regress deadin5 ldl10, robust
poisson deadin5 ldl10, robust irr
logistic deadin5 ldl10

g ldl10 = ldl/10
regress deadin5 ldl10, robust
poisson deadin5 ldl10, robust irr
logistic deadin5 ldl10

g logldl = log(ldl) / log(1.1)
. regress deadin5 logldl, robust