Hi, I would like to get the probability of having any "response level" for each level of "treatment" in my data, and not only the "probabilities of levels of y having LOWER Ordered Values in the response profile table." that SAS automatically outputs. (I'm simplifying a lot the real scenario just to have an idea). Let's say I would like to see the effect of 3 treatments (T1, T2 and T3) on a given health status (5 categories: low, mid-low, med, mid-high and high), correcting for the covariate Weight (W). I'm running an ordinal multinomial model: ---------------------------------------------------------------------------------- proc genmod data=health; class T; model Status = T W/dist=multinomial link=clogit type3; estimate 'T1 - T2' T 1 -1/exp; run; --------------------------------------------------------------------------- This (estimate statement) code gives that the odds of T1 being in lower health categories is 2.8 the odds of T2 being in lower health categories. Although this is interesting, that's not what I want. What I really want are all the probabilities of being in each of the 5 Health Categories for each of the 3 Treatments: P(low|T1) P(mid-low|T1) ... P(high|T3). Is is possible to have this in Proc Genmod? If so, can anyone help me? Thanks, Nick
... View more