Hi SAS members   I am trying to pool my effect estimates across my imputed datasets running a linear regression using GLIMMIX. My dependent variable is a continuous variable and my independent variable is a binary treatment variable adjusting for PS decile.   Here is my current code:   PROC GLIMMIX DATA=moh20crt.ps1_d_tt_response METHOD=LAPLACE NOCLPRINT; class biol (ref='TOFA') time patid ; MODEL /*radi_lda*/ /*radi_mda*/ radi=biol decile / CL SOLUTION ; RANDOM INTERCEPT /subject=patid ; by _imputation_; COVTEST / WALD;where month<=8; ods output ParameterEstimates=a_mvn ; run;   /* this part works well*/   Then I try to combine the estimate across 20 MI dataset as below PROC MIANALYZE data=a_mvn ; MODELEFFECTS biol ; stderr biol; RUN; /* this part does not work*/   I am getting message below in the log  "ERROR: Within-imputation Estimate missing for variable biol in _Imputation_= 1 in the input DATA=data set"   Any thought to resolve the problem   Thanks Moh ... View more