How to perform mixed effect ordinal logistic regression?

[Diels_O]

Dear all, I hava a question about the proc glimmix; the following was my code

 proc glimmix data=pt empirical order = data PLOTs= all;
 class Affiliation pt60dich3(ref = first)/ref = first;
 model score = pt60dich3 / dist=bin link=glogit SOLUTION cl oddsratio;
 random intercept/subject = Affiliation; run; 

Here, score was a ordinal score.

However, this model went wrong: the generalized logit can only be used with the multinomial distribution.

How should i perfrom mixed effect ordinal logistic model?

[PaigeMiller]

Ordinal implies (to me) more than two levels of the response variable. Is that what you have?

If so, then BIN is the wrong value for DIST. You do need to specify the MULTINOMIAL distribution.

[Diels_O]

Thank you for your quick response!

[StatDave · Accepted solution]

The generalized logit model (LINK=GLOGIT) is typically used when the multinomial response is nominal. If your multinomial response is ordinal and you want the model to account for its ordinality, then you can use the cumulative logit model (LINK=CUMLOGIT). The cumulative logit model is the default when you specify DIST=MULT.