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
NON-EST usually happens when you have a cell of the interaction with no data in it.
Well the treatment structure is a 3x4 factorial (3 different temperatures, 4 different time periods). Could the error be a result of the factors not being of equal levels?
I'll stick with my first answer. The most likely cause is that you have an empty cell.
Any idea as to where this empty cell might be?
You have the data, I don't. You can find it easily enough using PROC FREQ.
Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.
Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.