PROC ANOVA;
CLASSES BLK B M N P;
MODEL RESP= BLK B BLK*B M B*M BLK*B*M
N P N*P B*N B*P B*N*P BLK*B*N*P
M*N M*P M*N*P BLK*M*N*P
B*M*N B*M*P B*M*N*P;
TEST H=B E=BLK*B;
TEST H=M B*M E=BLK*B*M;
TEST H=N P N*P B*N B*P B*N*P E=BLK*B*N*P;
TEST H=M*N M*P M*N*P E=BLK*M*N*P;
MEANS B|M|N|P;
*Note: this analysis is not appropriate for missing data, nor does it give correct mean separations for intereaction means in split-plot design. For those analysis, PROC MIXED is the desired tool.
Source: Tim Todd (Plant Pathology, KSU); November 1997