options ls=80;
data wtrmln;
input row colm trt tonnage;
lines;
1       1       1       9.5
1       2       3       6.8
1       3       4       4.9
1       4       2       7.1
2       1       2       7.9
2       2       1       9.1
2       3       3       6.6
2       4       4       5.3
3       1       4       5.6
3       2       2       7.6
3       3       1       8.7
3       4       3       6.7
4       1       3       7.1
4       2       4       5.4
4       3       2       6.9
4       4       1       8.8
;
/*
proc univariate data=wtrmln normal;
var tonnage;
by trt;
run;
*/
proc glm data=wtrmln;
class row colm trt;
model tonnage = row colm trt;
run;

quit;