Difficulties I see, without having a good dataset to test this on: 1. The multinomial distribution option in GLIMMIX does not support calculation of BLUEs (LSMEANS), so residual calculation may be similarly not supported. Maybe someone else has better information, but this is my best guess as to why you aren't getting residuals. 2. You are not modeling the repeated nature of the data in the current code. The RANDOM statement should look something like: random day/residual sub=id type=cs (for example only, you may have a better idea of what covariance structure fits your data). You will also have to include day in the CLASS statement. 3. If you include R side effects like this, then method=quad will not work--the only available methods are the pseudo-likelihood methods. 4. To get odds ratios, add oddsratio to your model statement. Perhaps this will give you the info you need to check the proportionality assumption. Good luck, and keep posting on this. I am facing a similar challenge, and think I can learn from what happens here. Steve Denham
... View more