SAS批量卡方检验求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求助各位大神，数据如下：

求教如何一次性的批量求出变量A1到An的卡方值？group是分组，OT是结果 Ai部分是需要加权的频数。
不胜感激

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：卡方检验 Group 不胜感激 卡方值 一次性 不胜感激 如何

自己顶一下吧，求关注，求解答！

1. /*解决实际问题之批量进行卡方检验*/
   
2. data chisq;
   
3. input group$ class$ A1 A2 A3 A4 A5 A6 A7 A8 A9 A10 A11 A12;
   
4. cards;
   
5. 1 1 9 100 95 54 17 1 0 1 3 69 32 70
   
6. 1 2 267 176 181 222 259 275 276 275 273 207 244 206
   
7. 2 1 25 131 142 89 28 1 2 0 4 111 56 108
   
8. 2 2 393 287 276 329 390 417 416 418 414 307 362 310
   
9. ;
   
10. %macro ka(start,end);
    
11. %do i=&start %to &end;
    
12. proc freq data=chisq;
    
13. weight A&i;
    
14. table group*class/chisq nocol norow nopercent agree;
    
15. title "A&i.的卡方检验";
    
16. output;
    
17. %end;
    
18. run;
    
19. %mend ka;
    
20. %ka(1,12)
    
21. /*12组卡方检验，结果对应的p值(>0.05)都无显著性，说明class与group无关*/

复制代码

1. /*解决实际问题之批量进行卡方检验*/
   
2. data chisq;
   
3. input group$ class$ A1 A2 A3 A4 A5 A6 A7 A8 A9 A10 A11 A12;
   
4. cards;
   
5. 1 1 9 100 95 54 17 1 0 1 3 69 32 70
   
6. 1 2 267 176 181 222 259 275 276 275 273 207 244 206
   
7. 2 1 25 131 142 89 28 1 2 0 4 111 56 108
   
8. 2 2 393 287 276 329 390 417 416 418 414 307 362 310
   
9. ;
   
10. %macro ka(start,end);
    
11. %do i=&start %to &end;
    
12. proc freq data=chisq;
    
13. weight A&i;
    
14. table group*class/chisq nocol norow nopercent agree;
    
15. title "A&i.的卡方检验";
    
16. output;
    
17. %end;
    
18. run;
    
19. %mend ka;
    
20. %ka(1,12)
    
21. /*12组卡方检验，结果对应的p值(>0.05)都无显著性，说明class与group无关*/

复制代码

{:2_25:}

1. data a;
   
2. input group ot a1-a5;
   
3. cards;
   
4. 1  1   9 100  95  54  17
   
5. 1  2 267 176 181 222 259
   
6. 2  1  25 131 142  89  28
   
7. 2  2 393 287 275 329 390
   
8. ;
   
9. run;
   
10. 
    
11. proc transpose data=a out=b ;
    
12. by group ot;
    
13. run;
    
14. 
    
15. data c;
    
16.         set b;
    
17.         order=input(compress(_name_,'a'),best.);
    
18.         drop _name_;
    
19. run;
    
20. 
    
21. proc sort data=c;
    
22.         by order;
    
23. run;
    
24. 
    
25. ods trace on;
    
26. ods output   ChiSq=chi;
    
27. proc freq data=c;
    
28. by order;
    
29. tables group*ot /chisq;
    
30. weight col1;
    
31. run;

复制代码

1. /*解决实际问题之批量进行卡方检验*/
   
2. data chisq;
   
3. input group$ class$ A1-A12;
   
4. cards;
   
5. 1 1 9 100 95 54 17 1 0 1 3 69 32 70
   
6. 1 2 267 176 181 222 259 275 276 275 273 207 244 206
   
7. 2 1 25 131 142 89 28 1 2 0 4 111 56 108
   
8. 2 2 393 287 276 329 390 417 416 418 414 307 362 310
   
9. ;
   
10. %macro ka(start,end);
    
11. %do i=&start %to &end;
    
12. proc freq data=chisq;
    
13. weight A&i;
    
14. table group*class/chisq nocol norow nopercent agree;
    
15. title "A&i.的卡方检验";
    
16. output;
    
17. %end;
    
18. run;
    
19. %mend ka;
    
20. %ka(1,12)
    
21. /*12组卡方检验，结果对应的p值(>0.05)都无显著性，说明class与group无关*/

复制代码

zwnSAS121 发表于 2017-2-21 09:46

谢谢您的回复，十分感谢哈，但是假如这个Ai的i是好几千的话怎么办呢？

wang1839 发表于 2017-2-21 10:30

谢谢您的回复，十分感谢哈，但是假如这个Ai的i是好几千的话怎么办呢？

wang1839 发表于 2017-2-21 10:30

谢谢您的回复，十分感谢哈，但是假如这个Ai的i是好几千的话怎么办呢？