options nocenter ls=75;

data new;
input meth subj time1 time2 time3;
resp=time1; time=1; person=7*(meth-1)+subj;output;
resp=time2; time=2; person=7*(meth-1)+subj;output;
resp=time3; time=3; person=7*(meth-1)+subj;output;

cards;
1    1 25.4 30.6 29.1
1    2 27.4 29.3 28.0
1    3 25.5 30.0 27.0
1    4 25.8 29.7 27.9
1    5 26.2 31.3 29.2
1    6 24.6 26.6 26.6
1    7 25.6 28.0 28.3
2    1 27.6 27.1 27.8
2    2 24.7 29.0 26.2
2    3 26.3 27.3 28.4
2    4 25.0 29.7 29.9
2    5 25.7 29.5 29.4
2    6 28.5 29.7 30.4
2    7 22.9 27.2 26.4
3    1 22.8 25.1 27.2
3    2 24.2 24.0 26.2
3    3 25.3 25.2 28.1
3    4 25.4 24.7 31.0
3    5 24.5 26.2 29.4
3    6 25.6 26.9 29.2
3    7 25.6 24.8 28.9
;

proc mixed;
class meth subj time;
model resp = meth time meth*time / solution;
random subj(meth);
lsmeans meth*time / diff;

proc mixed;
class meth subj time;
model resp = meth time meth*time;
random subj(meth);
repeated / subject=person type=un r;
lsmeans meth*time / diff;

proc mixed;
class meth subj time;
model resp = meth time meth*time;
random subj(meth);
repeated / subject=person type=cs r;
lsmeans meth*time /diff;

proc mixed;
class meth subj time;
model resp = meth time meth*time;
random subj(meth);
repeated / subject=person type=ar(1) r;
lsmeans meth*time / diff;
run;