So I'm using GLIMMIX but I'm pretty sure this question can be applicable to other procedures that utilize LSMeans.   I have a model set up like so:     Class Druggroup Period Units; Cases = Druggroup*Period + Druggroup*Units + ln_respondents; Lsmeans Druggroup*Period / pdiff cl at ln_respondents = 0;     I want the LSMeans for the Druggroup*Period combinations AT (ln_respondents = 0) AND (Units = Reference category of Units).   Using AT ln_respondents = 0 is easy to do. However I cannot find out the way to have the LSMeans given for the REFERENCE category of Units. By default it looks like it gives the mean ACROSS the different categories.   I know this is possible to do with an ESTIMATE statement, but I am trying to make this code flexible (running many different models through it) so I was hoping to utilize LSMeans.   So, is it possible to get this info from LSMeans or am I barking up the wrong tree?   Thanks ... View more