假定有一个随机抽样调查，得到四个工厂在四种改革方案中，工人收入的状况如下表。

shiyiming

厂名        方案1        方案2        方案3        方案4
工厂1        146        200        148        151
工厂2        158        303        282        290
工厂3        415        461        431        413
工厂4        454        452        453        415
问：（1）各种方案之间的收入是否有显著差异？
   （2）不同规模的工厂之间，收入是否有显著差异？

shiyiming

[code:30hqh2dr]data a;
input name $ p1 p2 p3 p4;
cards;
f1 146 200 148 151
f2 158 303 282 290
f3 415 461 431 413
f4 454 452 453 415
;
data b; set a;
array p(4) p1-p4;
do i=1 to 4;
do j=i+1 to 4;
/* ij=catx('_',i,j);*/
d=p(j)-p(i);output;
end;
end;
proc transpose data=b out=c(where=(d1 is not missing) drop=_name_) prefix=d;
var d;
by name;
copy p1 p2 p3 p4;
run;
%macro d;
%do i=1 %to 6;
proc univariate data=c;
  var d&i.;
run;
%end;
%mend;
%d
[/code:30hqh2dr]
this code just for the (1), and the (2) should be the similar, use proc transpose first.

shiyiming

data ads;
do w =1 to 4;
   do m=1 to 4;
      input y @@;
          output;
   end;
end;
cards;
146 200 148 151
158 303 282 290
415 461 431 413
454 452 453 415
;
proc anova data=ads;
  class m w;
  model y=m w;
  means m w/snk;
run;

shiyiming

我只想说一句 你就四个数据点 能做什么