options nodate nonumber ps=55 ls=76;

data one;
do trt='A', 'V', 'S';
do rabbit=1 to 5;
input y1 y2 y3 y4 @@;
output;
end; end;
cards;
-0.3 -0.2 1.2 3.1 -0.5 2.2 3.3 3.7 -1.1 2.4 2.2 2.7
1.0 1.7 2.1 2.5 -0.3 0.8 0.6 0.9
-1.1 -2.2 0.2 0.3 -1.4 -0.2 -0.5 -0.1 -0.1 -0.1 -0.5 -0.3
-0.2 0.1 -0.2 0.4 -0.1 -0.2 0.7 -0.3
-1.8 0.2 0.1 0.6 -0.5 0.0 1.0 0.5 -1.0 -0.3 -2.1 0.6
0.4 0.4 -0.7 -0.3 -0.5 0.9 -0.4 -0.3
;
run;

proc glm;
class trt;
model y1-y4 = trt / nouni;
repeated time 4 (0 30 60 90) polynomial / printe;
run;

data one;
do trt='A', 'V', 'S';
do rabbit=1 to 5;
do time=1 to 4;
input y @@;
output;
end; end; end;
cards;
-0.3 -0.2 1.2 3.1 -0.5 2.2 3.3 3.7 -1.1 2.4 2.2 2.7
1.0 1.7 2.1 2.5 -0.3 0.8 0.6 0.9
-1.1 -2.2 0.2 0.3 -1.4 -0.2 -0.5 -0.1 -0.1 -0.1 -0.5 -0.3
-0.2 0.1 -0.2 0.4 -0.1 -0.2 0.7 -0.3
-1.8 0.2 0.1 0.6 -0.5 0.0 1.0 0.5 -1.0 -0.3 -2.1 0.6
0.4 0.4 -0.7 -0.3 -0.5 0.9 -0.4 -0.3
;
run;

proc glm;
class trt rabbit time;
model y = trt rabbit(trt) time time*trt;
test h=trt e=rabbit(trt);
run;

quit;

quit;