请求完全随机设计秩检验和方差分析后两两比较的具体P值

请求完全随机设计秩检验后两两比较的具体P值，我有两个程序可分别进行完全随机设计秩检验和方差分析后两两比较，但没办法出现两两比较的具体P值，望请高手指点！希望是在我的程序上修改！






data b3;                     
  input x2 c @@;
cards;
                    0                1         0.0548       2        0                  3            0.0896            4             0                  5                 0.3608           6
                   0.0196        1         0.058         2       0.0222          3            0.0738            4             0                   5                0.4335           6
                   0                 1         0.0429       2        0.0192         3            0.0893            4             0                   5                0.3807           6
                   0                 1         0.0506      2         0                  3            0.0756            4             0.0172          5                 0.3728           6
                   0.0204        1         0.0517      2         0.0213         3            0.0917            4             0.0156          5                 0.4188           6
                   0                1          0.038        2          0                 3            0.0755            4             0                   5                 0.4088           6
                   0.0208       1          0.0405     2         0.0256         3             0.0707            4             0                   5                 0.4706           6
                   0               1          0.0645      2          0                 3             0.0833            4             0                   5                  0.3939          6
                   0.025        1          0.0339      2          0.02            3             0.0808            4             0                   5                  0.2865          6
                  0                1          0.0633      2          0                 3             0.0737            4             0.0167          5                  0.3785          6
                  0.0192       1          0.0411      2          0.0328        3             0.0825            4             0.0172          5                  0.381            6
                  0.0182      1           0.0635      2          0                 3             0.0614            4             0                   5                  0.3591          6
                  0               1           0.0597      2          0.0196        3             0.0809            4             0                  5                   0.3279          6
                  0               1          0.0649       2         0                  3             0.0696             4            0                  5                   0.3296          6
                  0.0213      1          0.0641       2         0.0244         3             0.0645             4            0.0169         5                   0.2978          6
                  0               1         0.0423        2         0                  3             0.088                4           0.0179         5                   0.3547          6
                  0.0204      1         0.058          2        0                   3             0.098                4           0                  5                  0.3198           6
                  0               1         0.0519        2        0.04              3             0.0796              4            0                 5                   0.4268          6
                  0.0192      1         0.0571        2        0                   3             0.0476              4            0                 5                  0.2807           6
                  0.0204      1         0.0294        2        0                   3             0.0822              4            0.0164        5                  0.2995           6
;
proc npar1way wilcoxon  data=b3;                  
  var x2;
  class c;
run;
proc rank data=b3 out=b4;
var x2;
ranks rx2;
run;
proc anova data=b4;
class c;
model rx2= c;
means c/hovtest dunnett snk duncan lsd;
run;
ods html close;





data b3;
  input x2 c @@;
cards;
                    0                1         0.0548       2        0                  3            0.0896            4             0                  5                 0.3608           6
                   0.0196        1         0.058         2       0.0222          3            0.0738            4             0                   5                0.4335           6
                   0                 1         0.0429       2        0.0192         3            0.0893            4             0                   5                0.3807           6
                   0                 1         0.0506      2         0                  3            0.0756            4             0.0172          5                 0.3728           6
                   0.0204        1         0.0517      2         0.0213         3            0.0917            4             0.0156          5                 0.4188           6
                   0                1          0.038        2          0                 3            0.0755            4             0                   5                 0.4088           6
                   0.0208       1          0.0405     2         0.0256         3             0.0707            4             0                   5                 0.4706           6
                   0               1          0.0645      2          0                 3             0.0833            4             0                   5                  0.3939          6
                   0.025        1          0.0339      2          0.02            3             0.0808            4             0                   5                  0.2865          6
                  0                1          0.0633      2          0                 3             0.0737            4             0.0167          5                  0.3785          6
                  0.0192       1          0.0411      2          0.0328        3             0.0825            4             0.0172          5                  0.381            6
                  0.0182      1           0.0635      2          0                 3             0.0614            4             0                   5                  0.3591          6
                  0               1           0.0597      2          0.0196        3             0.0809            4             0                  5                   0.3279          6
                  0               1          0.0649       2         0                  3             0.0696             4            0                  5                   0.3296          6
                  0.0213      1          0.0641       2         0.0244         3             0.0645             4            0.0169         5                   0.2978          6
                  0               1         0.0423        2         0                  3             0.088                4           0.0179         5                   0.3547          6
                  0.0204      1         0.058          2        0                   3             0.098                4           0                  5                  0.3198           6
                  0               1         0.0519        2        0.04              3             0.0796              4            0                 5                   0.4268          6
                  0.0192      1         0.0571        2        0                   3             0.0476              4            0                 5                  0.2807           6
                  0.0204      1         0.0294        2        0                   3             0.0822              4            0.0164        5                  0.2995           6
;
proc anova data=b3;
  class c;
  model x2=c;
  means c/hovtest dunnett snk duncan lsd;
run;

期待高手，让我们也能借着楼主的问题，学到新东西

两两比较可以写成macro的形式用
%macro a(n1=,n2=);
%do i=1 %to &n1;
    %do j=1 %to &n2;
proc npar1way wilcoxon  data=b3;
where c in (&i,&j);
  var x2;
  class c;
run;
     %end;
%end;
%mend a;
%a(n1=6,n2=6);

谢谢3楼！我先试一下！