采用方差分析比较不同试验组之间的组间差异,当方差分析P≤0.05时,采用LSD-t检验进行试验组组间两两比较。
ods graphics off;
proc glm data=adef1;
/* by StatOrd;*/
class TRTPN;
model AVAL = TRTPN/ss3;
lsmeans TRTPN/ pdiff CL stderr tdiff;
ods output ModelANOVA=ModelANOVA Diff=Diff;
quit;
*方差分析;
data ds_ana2;
set ModelANOVA;
format _all_;
__StatValue=FValue;
__PValue=ProbF;
run;
*差分结果;
data ds_ana4;
array va{3,6} _temporary_;
do until(last.StatOrd);
set Diff;
by StatOrd;
if first.StatOrd then i = 1;
else i + 1;
va{i,1} = _1;
va{i,2} = _2;
va{i,3} = _3;
va{i,4} = P1;
va{i,5} = P2;
va{i,6} = P3;
end;
do i=1 to 2;
do j=i+1 to 3;
trtd = cats(i,'-',j);
__StatValue =va[i,j];
__PValue =va[i,j+3];
output;
end;
end;
run;