裂区试验设计方差分析

data split;
input R A B y;
cards;
1 4 2 40.0
1 4 1 39.5
1 1 2 37.9
1 1 1 35.4
1 2 1 36.7
1 2 2 38.2
1 3 2 36.4
1 3 1 34.8
2 2 1 42.7
2 2 2 41.6
2 1 2 40.3
2 1 1 41.6
2 4 1 44.5
2 4 2 47.6
2 3 1 43.6
2 3 2 42.8
;
run;

proc anova data=split;
class R A B;
model y = R A R*A B A*B;
test h=R A e=R*A;
means A / duncan e=R*A;
means B / duncan;
run;

proc glm data=split;
class R A B;
model y = R A R*A B A*B;
random R*A / test;
means A / duncan e=R*A;
means B / duncan;
lsmeans A*B / adjust=bon;
run;