ANOVA Codings for PROC TRANSREG
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: TREGANO */
/* TITLE: ANOVA Codings for PROC TRANSREG */
/* PRODUCT: STAT */
/* SYSTEM: ALL */
/* KEYS: analysis of variance */
/* PROCS: TRANSREG */
/* DATA: */
/* */
/* SUPPORT: saswfk UPDATE: July 25, 2010 */
/* REF: PROC TRANSREG, DETAILS, ANOVA CODINGS */
/* MISC: */
/****************************************************************/
title 'Two-Way ANOVA Models';
data x;
input a b @@;
do i = 1 to 2; input y @@; output; end;
drop i;
datalines;
1 1 16 14 1 2 15 13
2 1 1 9 2 2 12 20
3 1 14 8 3 2 18 20
;
proc print label;
run;
proc transreg data=x ss2 short;
title2 'Cell-Means Model';
model identity(y) = class(a * b / zero=none);
output replace;
run;
proc print label;
run;
proc transreg data=x ss2 short;
title2 'Reference Cell Model, (3,2) Reference Cell';
model identity(y) = class(a | b);
output replace;
run;
proc print label;
run;
proc transreg data=x ss2 short;
title2 'Deviations from Means, (3,2) Reference Cell';
model identity(y) = class(a | b / deviations);
output replace;
run;
proc print label;
run;
proc transreg data=x ss2 short;
title2 'Less-Than-Full-Rank Model';
model identity(y) = class(a | b / zero=sum);
output replace;
run;
proc print label;
run;
proc transreg data=x ss2 short;
title2 'Reference Cell Model, (1,1) Reference Cell';
model identity(y) = class(a | b / zero=first);
output replace;
run;
proc print label;
run;
proc transreg data=x ss2 short;
title2 'Deviations from Means, (1,1) Reference Cell';
model identity(y) = class(a | b / deviations zero=first);
output replace;
run;
proc print label;
run;
proc transreg data=x ss2 short;
title2 'Orthogonal Contrast Coding';
model identity(y) = class(a | b / orthogonal);
output replace;
run;
proc print label;
run;
proc transreg data=x ss2 short;
title2 'Standardized-Orthogonal Coding';
model identity(y) = class(a | b / standorth);
output replace;
run;
proc print label;
run;