data new;
  do trt=1 to 3;
    do subtrt=1 to 4;
  do blk=1 to 3;
    input y@@;
output;
end;end;end;
cards;
29 28 32
37 32 31
18 14 17
17 16 15
28 29 25
31 28 29
13 13 10
13 12 12
30 27 26
31 28 31
15 14 11
16 15 13

proc  anova;
class trt subtrt blk;
model y=blk trt subtrt blk*trt blk*subtrt(trt) trt*subtrt ;
test h=blk e=blk*trt; /*测验区组效应*/
test h=trt e=blk*trt; /*测验主处理*/
test h=subtrt e=blk*subtrt(trt); /*测验副处理*/
test h=trt*subtrt e=blk*subtrt(trt);/*测验主处理与副处理交互效应*/
means trt/duncan e=blk*trt;        /*进行主处理新复极差测验*/
means subtrt/duncan e=blk*subtrt(trt);/*进行副处理新复极差测验*/*/
run;