I am running a Mixed Ordinal logistic regression model, with some fixed effects and one hierarchical random effect.
I understand that Proc Glimmix could run this model, however, we are requested to report the marginal probabilities and standard errors for the marginal probabilities. Glimmix only gives marginal probability but not the error for marginal probability, does anyone know how to get the standard error?
I think you could get the standard error you want by fitting the model using PROC NLMIXED instead of PROC GLIMMIX. Use the PREDICT statement to estimate the desired statistic and obtain its standard error.
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.