全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
1855 8
2013-05-07
悬赏 200 个论坛币 已解决
请问高手,我的数据如下,我用下列程序,可以运算出如下结果,但是我如果不用宏程序,则无法得出如下结果,都是lsd检验,为什么要用宏程序才可以得出如下两两比较的结果,请高手指点!谢谢!
data b5;
  input x3 c @@;
cards;
  1.22    1    .8409          2       .1216      3      .7120      4     .8333      5
  1.34    1    .6974      2       .2316      3      .5905      4     .9240      5
  1.13    1    .8409      2       .4475             3      .7120      4     .7117      5
  1.43    1    .7792      2       .9659             3      .8950      4     .6624      5
  1.53    1   1.2924      2       .5786      3      .2535      4     .3237      5
  1.26    1    .7738      2       .4353             3      .3439      4     .7303      5
;

%let digits=10;
proc template;
define column Common.pValue;
notes "Default p-value column";
just = r;
format = pvalue%eval(&digits+2).&digits;
end;
run;

proc anova data=b5;
  class c;
  model x3=c;
  means c/hovtest dunnett snk duncan lsd;
run;


用宏程序后结果如下:

The ANOVA Procedure

                                                   t Tests (LSD) for x3

             NOTE: This test controls the Type I comparisonwise error rate, not the experimentwise error rate.


                                           Alpha                            0.05
                                           Error Degrees of Freedom           25
                                           Error Mean Square            0.051178
                                           Critical Value of t           2.05954
                                           Least Significant Difference    0.269


                              Comparisons significant at the 0.05 level are indicated by ***.


                                                   Difference
                                         c            Between     95% Confidence
                                     Comparison         Means         Limits

                                       1 - 2           0.4476     0.1786  0.7166  ***
                                       1 - 5           0.6208     0.3518  0.8898  ***
                                       1 - 4           0.7339     0.4649  1.0028  ***
                                       1 - 3           0.8549     0.5859  1.1239  ***
                                       2 - 1          -0.4476    -0.7166 -0.1786  ***
                                       2 - 5           0.1732    -0.0958  0.4422
                                       2 - 4           0.2863     0.0173  0.5553  ***
                                       2 - 3           0.4074     0.1384  0.6763  ***
                                       5 - 1          -0.6208    -0.8898 -0.3518  ***
                                       5 - 2          -0.1732    -0.4422  0.0958
                                       5 - 4           0.1131    -0.1559  0.3821
                                       5 - 3           0.2342    -0.0348  0.5031
                                       4 - 1          -0.7339    -1.0028 -0.4649  ***
                                       4 - 2          -0.2863    -0.5553 -0.0173  ***
                                       4 - 5          -0.1131    -0.3821  0.1559
                                       4 - 3           0.1211    -0.1479  0.3901
                                       3 - 1          -0.8549    -1.1239 -0.5859  ***
                                       3 - 2          -0.4074    -0.6763 -0.1384  ***
                                       3 - 5          -0.2342    -0.5031  0.0348
                                       3 - 4          -0.1211    -0.3901  0.1479



不用宏程序,则无法运算出上述结果,为什么?请问高手为什么?谢谢!!!











最佳答案

farmman60 查看完整内容

proc anova; class c; model x3=c; means c/cldiff t; run; or proc glm; class c; model x3=c; lsmeans c/cl pdiff; run;
二维码

扫码加我 拉你入群

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

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

全部回复
2013-5-7 18:33:57

proc anova;
class c;
model x3=c;
means c/cldiff t;
run;

or

proc glm;
class c;
model x3=c;
lsmeans c/cl  pdiff;
run;
二维码

扫码加我 拉你入群

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

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

2013-5-7 18:37:23
我也遇到了这样的问题!
二维码

扫码加我 拉你入群

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

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

2013-5-7 20:58:56
试着run了一下,不用宏也可以得到这样的结果啊
二维码

扫码加我 拉你入群

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

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

2013-5-7 23:58:11
谢谢farmman60!我运行您的程序,可以用!!
二维码

扫码加我 拉你入群

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

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

2013-5-8 00:29:56
您好farmman60!请问,用glm程序时,出现两个结果,下面一个是可信区间,我看得懂;上面一个是用两个交叉的组别P值表示,其两个交叉组别所对应的P值意思是不是同可信区间所得的意义是一样的,只不过更具体的把P值写出来而已,是吗?期待您的答复!!
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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