topic Mixed Model Type 3 Estimates in SAS Programming

Mixed Model Type 3 Estimates

cb12991 — Sat, 10 Mar 2018 15:06:44 GMT

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

Re: Mixed Model Type 3 Estimates

PaigeMiller — Sat, 10 Mar 2018 20:41:00 GMT

NON-EST usually happens when you have a cell of the interaction with no data in it.

Re: Mixed Model Type 3 Estimates

cb12991 — Sat, 10 Mar 2018 20:57:37 GMT

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?

Re: Mixed Model Type 3 Estimates

PaigeMiller — Sat, 10 Mar 2018 21:06:12 GMT

I'll stick with my first answer. The most likely cause is that you have an empty cell.

Re: Mixed Model Type 3 Estimates

cb12991 — Sat, 10 Mar 2018 21:22:22 GMT

Any idea as to where this empty cell might be?

Re: Mixed Model Type 3 Estimates

PaigeMiller — Sat, 10 Mar 2018 21:34:01 GMT

You have the data, I don't. You can find it easily enough using PROC FREQ.