Hi all, I am trying to develop code to obtain the ANOVA and break out the main effects and interaction SS using complete sets of orthogonal contrasts. However, my program is not generating estimates for the interaction contrasts (it returns "Non Est."). Followed is my code, if anyone could shed some light, I would greatly appreciate it! data aa; infile '/folders/myfolders/sasuser.v94/lab8.txt'firstobs=2; input trt temp week y; proc print; run; proc mixed method=type3; class temp week; model y = temp week temp*week; estimate 'temp_lin' temp -1 0 1; estimate 'temp_quad' temp 1 -2 1; estimate 'week_lin' week -3 -1 1 3; estimate 'week_quad' week 1 -1 -1 1; estimate 'week_cub' week -1 3 -3 1; estimate 't_l*w_l' temp*week 3 0 -3 1 0 -1 -1 0 1 -3 0 3; estimate 't_l*w_q' temp*week -1 0 1 1 0 -1 1 0 -1 -1 0 1; estimate 't_l*w_c' temp*week 1 0 -1 -3 0 3 3 0 -3 -1 0 1; estimate 't_q*w_l' temp*week -3 6 -3 -1 2 -1 1 -2 1 3 -6 3; estimate 't_q*w_q' temp*week 1 -2 1 -1 2 -1 -1 2 -1 1 -2 1; estimate 't_q*w_c' temp*week -1 2 -1 3 -6 3 -3 6 -3 1 -2 1; lsmeans temp week temp*week; run; Thanks, Cody
... View more