I assume that your proc glimmix looks something like:
proc glimmix data=datasetname;
class indepvar randomvar;
model depvar=indepvar/dist=binomial link=logistic;
random randomvar;
run;
(probably some other lines in there, but these are the critical ones)
You could add:
lsmeans indepvar/cl ilink;
ods output lsmeans=lsmeans;
This should give the "back-transformed" marginal means (which are the probabilities) and their standard errors. Note that this will NOT work for a multinomial model (which may be what you have with ordinal data), as the lsmeans statement is not supported for dist=multinomial.
Good luck,
Steve Denham
... View more
