This the code I use and I am interested in BLUE values for entry*treatment interation. the SOLUTION statement computes the BLUE values, but the problem is since treatment has only 2 levels, the estimated BLUE vlaue for one of the levels of treatment is 0 and then the BLUE estimates for entry*treatment interaction is only computed for entry by non zero treatment level. Thanks in advance.
Proc mixed data= b covtest ; by _Name_ ;
class Rep entry treatment;
model col1 = entry treatment entry*treatment / outp=RESR residual solution DDFM= KENWARDROGER ;
random rep ;
lsmeans entry entry*treatment;
ods output Solutionf=BLUE;
ODS output lsmeans=ls;
Title 'Mixed model analysis';
run;
