proc glimmix data=qc2 METHOD=RSPL;
class group_ID;
model Y = X
/dist=binomial link=logit s ddfm=kr STDCOEF;
run;
It is not clear how coefficients are standardized. Based on my investigation, centering is definitely done around the variable's grand mean; however, SD used for standardization is not 1. When I simulated it, SD was around 0.036 (meaning I created a z-score using mean=0 and STD=0.036 to get the same coefficient off STDCOEF,