单因素方差分析

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

1. data a;
   
2. input x g@@;
   
3. cards;
   
4. 21.31        1.00        20.16        2.00
   
5. 23.14        1.00        24.49        2.00
   
6. 27.48        1.00        21.32        2.00
   
7. 19.54        1.00        19.46        2.00
   
8. 18.03        1.00        25.63        2.00
   
9. 24.03        1.00        28.81        2.00
   
10. 22.82        1.00        18.74        2.00
    
11. 18.72        1.00        18.42        2.00
    
12. 23.46        1.00        26.13        2.00
    
13. 20.34        1.00        25.24        2.00
    
14. 26.98        1.00        20.23        2.00
    
15. 19.56        1.00        22.47        2.00
    
16. 17.39        1.00        29.38        2.00
    
17. 24.37        1.00        20.16        2.00
    
18. 16.01        1.00        22.51        2.00
    
19. 28.11        3.00        30.23        4.00
    
20. 24.74        3.00        36.84        4.00
    
21. 29.79        3.00        38.61        4.00
    
22. 22.68        3.00        27.13        4.00
    
23. 23.01        3.00        28.79        4.00
    
24. 28.32        3.00        33.24        4.00
    
25. 29.04        3.00        31.68        4.00
    
26. 24.33        3.00        28.29        4.00
    
27. 33.97        3.00        38.47        4.00
    
28. 21.86        3.00        35.10        4.00
    
29. 28.65        3.00        28.01        4.00
    
30. 25.13        3.00        23.37        4.00
    
31. 34.44        3.00        28.44        4.00
    
32. 31.69        3.00        34.22        4.00
    
33. 35.07        3.00        35.08        4.00
    
34. ;
    
35. run;
    
36. 
    
37. proc glm data=a;
    
38. class g;
    
39. model x=g/ss3;
    
40. means g/hovtest snk bon lsd;
    
41. run;
    
42. 
    
43. proc anova data=a;
    
44. class g;
    
45. model x=g;
    
46. means g/hovtest snk bon lsd;
    
47. run;

复制代码

上面是一段SAS单因素方差分析的程序，发上来分享一下。
另外，想讨论一个问题：
在spss中两两比较的bonferonni法可以输出具体的P值；但是在SAS里面的两两比较是没有具体P值的，只是分成了几个哑组；哪位大侠知道，SAS可否输出具体的两两比较的P值？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：单因素方差分析 方差分析 单因素 model means 分析

高深啊，没看懂

这个问题的确尖锐有见地。想必楼主也是兼通spss和sas。这样的人其实并不很多。
我猜想，既然在这里可以使用proc anova，说明实验是属于balanced design。对于这样的数据，means应该和lsmeans是并无二致的。如果要得到multiple comparison的p value，用lsmeans statement提供的pdiff选项就可以了。不过，如你所言，在means里的确没有这种选择。原因不清楚，因为我对此也是一知半解道听途说。至于spss就更加不了解了。
京剧

jingju11 发表于 2011-9-25 08:38
这个问题的确尖锐有见地。想必楼主也是兼通spss和sas。这样的人其实并不很多。
我猜想，既然在这里可以使用 ...

        jingju11，您过奖了，我对SPSS和SAS也是略懂一二。您说的没错，对于完全随机设计的单因素方差分析，用proc anova（anova语句中貌似没有lsmean statement语句）和proc glm两个proc步都可以，means和lsmeans语句是没有差别的，可以用lsmeans statement中的pdiff和tdiff得到具体的t值和p值。但是对于多因素设计的情况下，有可能存在交互效应，此时means和lsmeans statement就存在差别了，涉及到主效应、或者单独效应的比较。
        对于有些情况下，杂志或者其他人要求提供P值的情况下，就比较麻烦了~~~望指教。