全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件
2869 3
2011-05-25
悬赏 300 个论坛币 未解决
请求完全随机设计秩检验后两两比较的具体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;
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2011-5-26 07:43:17
期待高手,让我们也能借着楼主的问题,学到新东西
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-5-26 10:23:44
两两比较可以写成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);
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-7-5 03:38:22
谢谢3楼!我先试一下!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群