glm1.sas

SAS for Windows统计分析系统教程,里面有源码

options nodate nonumber;
data glm1;
   input g$ x1 x2 @@;
cards;
a 3.9 210 b 4.8 270 c 4.4 250
a 4.2 190 b 4.7 180 c 3.7 305
a 3.7 240 b 5.4 230 c 2.9 240
a 4.0 170 b 4.5 245 c 4.5 330
a 4.4 220 b 4.6 270 c 3.3 230
a 5.2 230 b 4.4 220 c 4.5 195
a 2.7 160 b 5.9 290 c 3.8 275
a 2.4 260 b 5.5 220 c 3.7 310
a 3.6 240 b 4.3 290
a 5.5 180 b 5.1 310
a 2.9 200
a 3.3 300
;
proc glm;
  class g;
  model x1 x2=g;
    contrast 'a vs b' g 1 -1  0;
    contrast 'a vs c' g 1  0 -1;
    contrast 'b vs c' g 0  1 -1;
  manova h=g/printe printh;
    means g;
run;