options pageno=1;
data latinsq;
input row column trt $ yield;
datalines;
1 4 a 1.345
1 1 b 1.64
1 3 c 1.425
1 2 d 1.21
2 2 a 1.185
2 4 b 1.29
2 1 c 1.475
2 3 d 1.4
3 1 a 1.67
3 3 b 1.665
3 2 c 0.71
3 4 d 1.18
4 3 a 1.655
4 2 b 1.29
4 4 c 0.66
4 1 d 1.565
;;
proc print;
Title 'Printout of Latin Square Data';
run;
proc anova;
class row column trt;
model yield=row column trt;
means trt/lsd;
Title 'ANOVA of a single Latin Square';
run;