请教SAS程序问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位高手，我想请问，两组等级资料的SAS程序，下面这个是否正确？
原始数据是
                         疗效                      合计
                      控制    显效    有效    无效
试验组  计数    57       13       25        4       99

对照组  计数   103      65       51         0     219

合计     计数   160      78       76         4    318

data x;   
do r=1 to 4;
do c=1 to 2;
input f@@;
output;
end; end;
cards;
57 103 13 65 25 51 4 0                                                                                                                                       ;      
run;
proc freq;
weight f;         
tables r*c /cmh2 scores=rank noprint;
run;

拜谢了哈！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：sas程序 scores Weight tables output 请教 程序 SAS

用SAS软件运行下不就知道啦，如果有问题会提示哪里出错。

不是啦，是这个SAS程序运行出来的结果和SPSS运行的结果相差很大。还是下面这个程序对呢？
data x;   
do r=1 to 4;
do c=1 to 2;
input f@@;
do i=1 to f;
x=r;
output;
end; end;end;
cards;
57 103 13 65 25 51 4 0      
;                                                                                                                                ;      
run;
proc npar1way wilcoxon;
class c;         
var x;
run;

从你的资料类型来看，是用卡方检验啊

shouldn't  you have 310 obs in your dataset?

3# wangd527
典型的卡方检验：
data x;   
do r=1 to 4;
do c=1 to 2;
input f @@;
output;
end;
end;
cards;
57 103 13 65 25 51 4 0
run;
proc freq;
weight f;         
tables r*c /chisq;
run;

1# wangd527

When not adjusting for individuals, the Row Mean Scores Differ CMH statistic (ANOVA) is identical to a Kruskal-Wallis test. 你的不一致的原因是treatment group 应该作为row 变量。也就是说，在npar1way 和 cmh 里你检测的东西不一致。总而言之，在你的freq里，r应为c，c应为r。或者，你也可以颠倒npar1way里class 和var 里的变量，结果应该和你原始的CMH值是一样的。
另外，在这里用cmh2 似乎没有意义。 in CMH, General association 和Chisq 是一致的。

1. 
   
2. data have;   
   
3. do r=1 to 4;
   
4.   do c=1 to 2;
   
5.    input f@@;
   
6.    output;
   
7.   end;
   
8. end;
   
9. cards;
   
10. 57 103 13 65 25 51 4 0
    
11. ;  
    
12. proc freq data = have;
    
13. weight f;         
    
14. tables c*r/cmh2 scores=rank noprint;
    
15. run;
    
16. proc npar1way data = have wilcoxon;
    
17. class c;         
    
18. var r;
    
19. freq f;
    
20. run;

复制代码

1# wangd527

这是个等级资料，行为无序分组资料，列为有序等级资料，宜采用Ridit分析或秩和检验，不宜用卡方检验。用CMH读取行平均得分项（ Row Mean Score）,要注意SAS程序的tables 后面的行列的顺序！你的行列顺序有误。

data x;   
do r=1 to 4;
do c=1 to 2;
input f@@;
output;
end; end;
cards;
57 103 13 65 25 51 4 0
;  
/*显示与你的表格顺序是否一致*/
proc freq;
weight f;         
tables c*r /nocol norow nopercent;
run;

proc freq;
weight f;         
tables c*r /cmh2 scores=ridit noprint;
run;