Hi Steve, Thank you for the quick answer. Now that I have gotten a chance to try out this method I am incredibly confused (your instruction was straight forward and great, but the output really had me scratching my head)...I think I need a lot more help :s When I read up on a multinomial response, I agree that my data does fit this approach. To my understanding the regions are independent within a subject (i.e. different sections of the intestine). I would also like to look at trt*score and region*trt*score. Would I be able to add that into the model statement like this: model score = trt cage region trt*score trt*cage trt*region cage*region region*trt*score trt*cage*region / dist=multinomial link=cumlogit aggregate type3; When I do this I get told: "WARNING: The relative Hessian convergence criterion of 8.8068819395 is greater than the limit of 0.0001. The convergence is questionable." When I set it up the way you suggested I get told: "WARNING: Negative of Hessian not positive definite." I tried to look into how to write estimate statements, but with each factor having at least 4 levels I was unsure how to set up the coding without at least 2 levels ending up looking the same. Also when looking at the Analysis Of Maximum Likelihood Parameter Estimates I am getting confused as to what the Pr > ChiS is telling me is significant. If I had trt*region would it be telling me that the specific trt with the specific region are significantly associated? (correlated?...am I even on the right track?) I am so sorry for how novice I am at this...I kind of got shoved into doing these types of statistics when I really only understand the bare minimum. Thanks for any and all of your help, I really appreciate it!! Kayla
... View more