Thanks, SAS_Rob! When I use the /*Combining Random Effects*/ code you provided, I just get 20 separate estimates for 20 imputation sets. I am looking for a way to have one pulled level 2 variance. I am able to do that for the fixed effects by using this program--at the end of the analysis, I have the combined effects for the 20 models. PROC GLIMMIX DATA=dissert.centered METHOD=LAPLACE NOCLPRINT; CLASS cdscode ethn_r (ref=FIRST) sex (ref=FIRST)sex; MODEL m_health (EVENT=LAST)= ethn_r sex grade_new_c /CL DIST=BINARY LINK=LOGIT SOLUTION ODDSRATIO (DIFF=FIRST LABEL); RANDOM INTERCEPT / SUBJECT=cdscode S CL TYPE=VC; by _imputation_; nloptions gconv=0; ods output ParameterEstimates=lgsparm35; COVTEST /WALD; RUN; proc mianalyze parms = lgsparm35; CLASS ethn_r sex; modeleffects Intercept ethn_r sex grade_new_c; run; Thanks for your help! Jelena
... View more