glm5_13_2.sas

黄燕 吴平《SAS统计分析及应用》 书中一些例子的SAS程序代码

proc glm data=glm5_13 outstat=glm5_132out;
class a b;
model y=a b a*b / ss3;
contrast  '(b1 vs b2)/a1 '  b 1 -1  0 a*b 1 -1  0;
contrast  '(b1 vs b3)/a1 '  b 1  0 -1 a*b 1  0 -1;
contrast  '(b2 vs b3)/a1 '  b 0  1 -1 a*b 0  1 -1;
contrast  '(b1 vs b2)/a2 '  b 1 -1  0 a*b 0  0  0  1 -1  0;
contrast  '(b1 vs b3)/a2 '  b 1  0 -1 a*b 0  0  0  1  0 -1;
contrast  '(b2 vs b3)/a2 '  b 0  1 -1 a*b 0  0  0  0  1 -1;
contrast  '(b1 vs b2)/a3 '  b 1 -1  0 a*b 0  0  0  0 0 0 1 -1  0;
contrast  '(b1 vs b3)/a3 '  b 1  0 -1 a*b 0  0  0  0 0 0 1  0 -1;
contrast  '(b2 vs b3)/a3 '  b 0  1 -1 a*b 0  0  0  0 0 0 0  1 -1;
run;   quit;
proc glm data=glm5_13 outstat=glm5_132out;
class b a;
model y=b a b*a / ss3;
contrast  '(a1 vs a2)/b1 '  a 1 -1  0 b*a 1 -1  0;
contrast  '(a1 vs a3)/b1 '  a 1  0 -1 b*a 1  0 -1;
contrast  '(a2 vs a3)/b1 '  a 0  1 -1 b*a 0  1 -1;
contrast  '(a1 vs a2)/b2 '  a 1 -1  0 b*a 0  0  0  1 -1  0;
contrast  '(a1 vs a3)/b2 '  a 1  0 -1 b*a 0  0  0  1  0 -1;
contrast  '(a2 vs a3)/b2 '  a 0  1 -1 b*a 0  0  0  0  1 -1;
contrast  '(a1 vs a2)/b3 '  a 1 -1  0 b*a 0  0  0  0 0 0 1 -1  0;
contrast  '(a1 vs a3)/b3 '  a 1  0 -1 b*a 0  0  0  0 0 0 1  0 -1;
contrast  '(a2 vs a3)/b3 '  a 0  1 -1 b*a 0  0  0  0 0 0 0  1 -1;
run;   quit;