Hi SAS users,   I am using proc glimmix to fit a 3-level mixed effects model with a continuous outcome, 2 binary independent variables(var1: 1 vs 0, var2: 1 vs 0), and their interaction. I want to estimate the effect of var1 when var2 =1, and the effect of var1 when var2=0. But the estimate gave nothing in the outputs, could you please let me know what I did wrong in the following codes:     proc glimmix data = data; class level3_ID level2_ID var1(ref='0') var2(ref='0'); model y = var1 var2  var1* var2/solution cl; random intercept/type = un subject = level3_ID; random intercept/type = un subject = level2_ID(level3_ID); estimate "effect estimate of var1 when var2 = 1" var1 1 0 var1*var2 1 0 0 0 /cl; estimate "effect estimate of var1 when var2 = 0" var1 1 0 var1*var2 0 0 1 0 /cl; run;   If proc genmod can do the same thing, it would be better, as the normality may not hold. But I am not sure how to do the 3-level model in proc genmod with a continuous outcome.   Any help would be greatly appreciated! ... View more