logistic单因素分析结果输出整理

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

由于经常需要做logistic回归需要先进行单因素分析，输出的结果粘贴太麻烦，所以写了段程序可以稍微整理下快速复制粘贴。

刚刚入门sas，可能写的稍微有些麻烦，大家有意见也可以提出来帮忙指正哈~

1. /*生成空的数据集记录结果*/
   
2. %macro blank(observe);
   
3. %data blankseed;
   
4.   varname='blank       ';
   
5. data blank;
   
6.   set blankseed;
   
7. run;
   
8. %do b=1 %to %eval(&observe-1);
   
9. data blank;
   
10.   set blank blankseed;
    
11. run;
    
12. %end;
    
13. run;
    
14. %mend;
    
15. %blank(1)
    
16. data logit_univariate_pe; set blank; run;
    
17. data logit_univariate_or; set blank; run;
    
18. 
    
19. /*1定量 2定性*/
    
20. %macro logit1(dataset= ,var= ,y= );
    
21. ods html close;
    
22.         ods output   ParameterEstimates=pe  OddsRatios=or;
    
23.         proc logistic data=&dataset desc;
    
24.                 model &y=&var;
    
25.         run;
    
26. ods html;
    
27. data logit_univariate_pe; set logit_univariate_pe pe; if Variable='Intercept' then delete; run;
    
28. data logit_univariate_or; set logit_univariate_or or; run;
    
29. %mend;
    
30. 
    
31. %macro logit2(dataset= ,var= ,y=  ,ref= );
    
32. ods html close;
    
33.         ods output parameterestimates=pe OddsRatios=or;
    
34.         proc logistic data=&dataset desc;
    
35.                 class &var(param=ref ref="&ref");
    
36.                 model &y=&var;
    
37.         run;
    
38. ods html;
    
39. data logit_univariate_pe; set logit_univariate_pe pe; if Variable='Intercept' then delete;  run;
    
40. data logit_univariate_or; set logit_univariate_or or;  run;
    
41. %mend;
    
42. 
    
43. /*调用宏进行分析*/
    
44. /*%logit1（dataset=test，var=x, y=y）*/
    
45. 
    
46. /*对最终输出结果进行整理*/
    
47. data logit_univariate;
    
48. merge logit_univariate_pe logit_univariate_or;
    
49. if probchisq=. then delete;
    
50. keep variable classval0 estimate waldchisq probchisq oddsratioest lowercl uppercl;
    
51. run;

复制代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：单因素分析 数据集记录 复制粘贴 进行分析 记录结果 logistic

结果输出如上所示

楼主，你好，请问如何使用，代码中哪些是需要改变的，数据集名称和变量

ransuoqu 发表于 2018-6-20 11:00
楼主，你好，请问如何使用，代码中哪些是需要改变的，数据集名称和变量

直接run了前三段宏程序，然后在调用宏的时候填写相关参数就可以了，
比如调用定量logistic时候，%logit1（dataset=test，var=x, y=y）里的tets改成你的数据集，X和Y分别改成因变量和自变量即可
定性同样是%logit2（dataset=test，var=x, y=y， ref=0）ref后面是根据你的变量选择refrence就可以了。
最后那段代码是用于整理输出结果的，之后直接print那个logit_univariate的数据集就是你需要输出的结果了~

思冬君 发表于 2018-6-26 14:37
直接run了前三段宏程序，然后在调用宏的时候填写相关参数就可以了，
比如调用定量logistic时候，%logit ...

非常感谢，您的指导。

您好，第一段程序里生成blank空集的时候就报错了，显示varname=‘blank ’这句有问题，并且第一个宏结尾那里没有分号，改过来还是报错，请问是我哪里没有调试好吗？

思冬君 发表于 2018-6-26 14:37
直接run了前三段宏程序，然后在调用宏的时候填写相关参数就可以了，
比如调用定量logistic时候，%logit ...

您好，第一段程序里生成blank空集的时候就报错了，显示varname=‘blank ’这句有问题，并且第一个宏结尾那里没有分号，改过来还是报错，请问是我哪里没有调试好吗？