Haris, Since you figured it out, maybe you can help me? I added a "lsmeans var / diff=anom" line to my proc glimmix code and all that did was lines to my output populated with . (missing values) for the main effects. Interestingly, they are the inverse in direction, but the same in magnitude as the original output without the lsmeans var / diff=anom line.    I have to conduct these follow up tests for two separate analyses. In one the two model variables are categorical with two levels each and their interaction. In a second, I have a two categorical variables with 2 levels each, a linear variable, the quadratic of that linear variable, and all the higher-order interactions.    Any help would be much appreciated. ... View more