N/A
Posts: 1

standard error on marginal probablity for Mixed Ordinal logistic regression

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.
SAS Employee
Posts: 242

Re: standard error on marginal probablity for Mixed Ordinal logistic regression

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.
Posts: 2,655

Re: standard error on marginal probablity for Mixed Ordinal logistic regression

I assume that your proc glimmix looks something like:

proc glimmix data=datasetname;
class indepvar randomvar;
random randomvar;
run;

(probably some other lines in there, but these are the critical ones)