standard error on marginal probablity for Mixed Ordinal logistic regression

[stat1234]

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?

Thanks.

[StatDave]

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.

[SteveDenham]

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