aconova2.sas

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

options nodate nonumber;
title 'Multivariate Analysis of Covariance';
proc format;
  value groupfmt 1='Hydrolysate-I'  2='Hydrolysate-II'  3='Casein';
data ancova2;
do i=1 to 8;
  do group=1 to 3;
   input age food weight @@;
   output;
  end;
end;
format group groupfmt.;
cards;
 6 281.7 37  5 309.8 24  8 259.3 82
10 274.0 47  6 317.8 43  5 241.2 66
 8 253.8 37 10 326.1 60  6 248.5 74
 5 261.4 34  8 322.1 50  7 242.8 79
 7 272.8 42  7 323.5 47  8 255.7 82
 5 272.2 27  6 321.2 42  7 254.3 76
 6 272.3 32  5 311.8 39  5 244.6 73
 7 293.2 44 10 324.5 53 10 243.8 90
;
proc glm;
  class group;
  model weight=group age food/solution;
  lsmeans group/pdiff stderr;
run;