The data are based on measurements of SBP in the HERS study. The data set allows us to compare methods of analysis with complete data and under simulated missing data. In the data sets are missing...

The data are based on measurements of SBP in the HERS study. The data set allows us to compare methods of analysis with complete data and under simulated missing data. In the data sets are missing data indicators (miss mar for bpmisslong and miss mar1 for bpmisswide) which have value 1 to flag SBP values which should be dropped to simulate data which displays MAR missingness. In particular, year 1 values from patients with higher baseline SBP are flagged more frequently and hence will be simulated as missing. You can consult the course website for the data sets and more complete documentation and details on Stata code.

Using bpmisslong, fit a mixed linear regression model with SBP as the outcome and visit (visit) as the predictor. In Stata, the command would be xtmixed sbp visit || pptid:. Compare the mixed model which uses the full data to one restricted to nonmissing data (miss_mar equals 0). Compare the results with the GEE results in Problem 11.7. How do you explain the difference in results between the GEE and a linear mixed model with MAR missing data?

Problem 11.7

Using bpmisslong, fit a GEE model with SBP as the outcome and visit (visit) as the predictor. In Stata, the command would be xtgee sbp visit, i(pptid) corr(exch). Compare a GEE model which uses the full data to one restricted to nonmissing data (miss_year equal to 0). What do you conclude about GEE with MAR missingness?