data one; do block=1 to 4; do trt=1 to 2 ; do sample=1 to 3 ; input y @@; output; end; end; end; cards; 6.18374 6.22535 6.21960 6.40572 6.42243 6.40572 6.30207 6.21780 6.24144 6.35211 6.45630 6.45063 6.07672 6.14354 6.06495 6.19807 6.24572 6.32257 6.08564 6.20051 6.13166 6.42162 6.20617 6.28693 proc glm data=one; class block trt sample; model y=block trt block*trt; test h=block trt e=block*trt; run;