proc glimmix data=asdf namelen=32;
where disaster_type /*age_desc*/ ne "";
class GROUPING_D_RISK_PT;
model &out
=
GROUPING_D_RISK_PT
/solution ddfm=kr dist=binomial link=logit s STDCOEF ;
lsmeans GROUPING_D_RISK_PT / ilink diff;
output out=gmxout residual=resid;
ods output
ParameterEstimates=kaz1
CovParms=uekawa1
nobs=jeana
ModelInfo=estes
dimensions=diminfo
ConvergenceStatus=concon
FitStatistics=FITSTAT
Diffs=DIF_RESULT
;
run;