options nocenter ls=75;
data new;
input row column period sp sum resp;
cards;
1 1 18 4 2 12.5
1 1 18 4 4 26.2
1 1 18 2 2 33.4
1 1 18 2 4 44.2
1 2  9 2 4 59.2
1 2  9 4 4 49.9
1 2  9 2 2 47.6
1 2  9 4 2 15.8
1 3  3 2 4 55.0
1 3  3 4 4 27.3
1 3  3 2 2 35.9
1 3  3 4 2 18.3
2 1  9 2 4 56.2
2 1  9 2 2 52.3
2 1  9 4 4 27.5
2 1  9 4 2 25.1
2 2  3 2 2 67.7
2 2  3 2 4 62.2
2 2  3 4 2 24.1
2 2  3 4 4 27.5
2 3 18 2 2 28.0
2 3 18 2 4 29.4
2 3 18 4 2 19.5
2 3 18 4 4 29.9
3 1  3 2 2 57.2
3 1  3 2 4 69.5
3 1  3 4 2 16.9
3 1  3 4 4 19.5
3 2 18 2 2 30.3
3 2 18 2 4 26.6
3 2 18 4 2 11.0
3 2 18 4 4 17.6
3 3  9 2 4 61.9
3 3  9 4 4 26.2
3 3  9 2 2 46.5
3 3  9 4 2 15.4
;

proc glm;
class row column period sp sum;
model resp=row column period period(row column) sp period*sp
sp*period(row column) sum period*sum sum*period(row column)
sp*sum period*sp*sum;
test h=period e=period(row column);
test h=sp e=sp*period(row column);
test h=sum e=sum*period(row column);
test h=sp*period e=sp*period(row column);
test h=sum*period e=sum*period(row column);
means sum sp|period;
run;