BookmarkSubscribeRSS Feed
stat1234
Calcite | Level 5
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.
2 REPLIES 2
StatDave
SAS Super FREQ
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
Jade | Level 19
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

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 2 replies
  • 1252 views
  • 0 likes
  • 3 in conversation