options nocenter ls=75;

data new;
input meth subj time1 time2;
cards;
1    1 25.4 30.6
1    2 27.4 29.3
1    3 25.5 30.0
1    4 25.8 29.7
1    5 26.2 31.3
1    6 24.6 26.6
1    7 25.6 28.0
2    1 27.6 27.1
2    2 24.7 29.0
2    3 26.3 27.3
2    4 25.0 29.7
2    5 25.7 29.5
2    6 28.5 29.7
2    7 22.9 27.2
3    1 22.8 25.1
3    2 24.2 24.0
3    3 25.3 25.2
3    4 25.4 24.7
3    5 24.5 26.2
3    6 25.6 26.9
3    7 25.6 24.8
;

data new1;
 set new;
 resp=time1; time=1; output;
 resp=time2; time=2; output;

proc glm data=new1;
 class meth subj time;
 model resp=meth subj(meth) time meth*time;
 test h=meth e=subj(meth);
 means meth|time;

proc glm data=new;
class meth subj;
model time1 time2 = meth;
repeated time 2 (1 2) polynomial / summary;
run;