set more off
cd "/Users/dave/UCSF/Projects/Teaching/Teaching Fall 10/HW #3/
use "hwk3 long.dta"
capture log using "hw3.log", replace

ttest basesbp if visit==2, by(miss_mar)
ttest sbp if visit==2, by(miss_mar)

xtgee sbp visit, i(pptid) corr(exch) robust
xtgee sbp visit if miss_mar==0, i(pptid) corr(exch) robust

use "hwk3 wide.dta"
mi set wide 
mi register imputed sbp2
set seed 668
mi impute mvn sbp2 = sbp1 bmi1 bmi2 baseline_dm, add(100)
sort miss_mar
by miss_mar: summ _1_sbp2 _2_sbp2 _3_sbp2 _4_sbp2 _5_sbp2
quietly mi reshape long sbp bmi, i(pptid) j(visit)
mi estimate, dots: xtgee sbp visit, i(pptid) corr(exch) robust
mi estimate, dots: xtmixed sbp visit || pptid:

use "hwk3 wide.dta"
mi set wide 
mi register imputed sbp2
set seed 668
mi impute mvn sbp2 = sbp1, add(100)
quietly mi reshape long sbp bmi, i(pptid) j(visit)
mi estimate, dots: xtgee sbp visit, i(pptid) corr(exch) robust
log close