SAS中文论坛

 找回密码
 立即注册

扫一扫,访问微社区

查看: 1927|回复: 5
打印 上一主题 下一主题

如何用Logistic回归单因素分析后把产生各个自变量的p值的.....

[复制链接]

1

主题

12

帖子

103

积分

注册会员

Rank: 2

积分
103
跳转到指定楼层
楼主
发表于 2014-3-23 19:14:39 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
如题:
现在做一个Logistic回归分析,先做单因素分析,产生了GLOBALTESTS_ageg3 GLOBALTESTS_sex等一系列包括卡方值、p值的数据集(共20个自变量,也就是20个这样的数据集),如果直接set这些数据集的话就不知道其中的卡方值和p'值如何与自变量对应了,求解决方法。
我的一个思路是将数据集名称的后缀作为一个变量加入到每个数据集中,使每个数据集都产生一个新变量,然后再合并,那么该如何合并这些数据集并且把数据集名称的后缀作为其中的一个变量?
或者有别的思路吗?

谢谢!!!!

分享到:  微信微信
收藏收藏
回复

使用道具 举报

0

主题

1

帖子

40

积分

新手上路

Rank: 1

积分
40
推荐
发表于 2014-3-25 22:50:30 | 只看该作者
不知道是否类似这样
  1. data test_age;
  2. x=1;
  3. run;

  4. data test_sex;
  5. x=2;
  6. run;

  7. proc sql noprint;
  8. select count(*) into:n from dictionary.tables where libname='WORK';
  9. select memname into :tablename1-:tablename%left(&n.)
  10. from dictionary.tables
  11. where libname='WORK';
  12. quit;

  13. %put &tablename1. &tablename2.;

  14. %macro test;
  15. %do i=1 %to &n.;
  16. data &&tablename&i.;
  17.         set &&tablename&i.;
  18.         tablename=scan("&&tablename&i.",2,"_");
  19. run;
  20. %end;
  21. %mend;

  22. %test
复制代码
回复 支持 1 反对 0

使用道具 举报

1

主题

6

帖子

45

积分

新手上路

Rank: 1

积分
45
沙发
发表于 2014-3-25 00:12:51 | 只看该作者
中文看不懂。 能不能, 把你的 logistic regression MODEL写出来, 让我好理解你要干什么?一起研究研究!
回复 支持 反对

使用道具 举报

1

主题

12

帖子

103

积分

注册会员

Rank: 2

积分
103
板凳
 楼主| 发表于 2014-3-25 18:54:25 | 只看该作者

%macro lgstc_1var (x,  y,  dataset);
ods output GlobalTests=GlobalTests_&x OddsRatios=OddsRatios_&x;
proc logistic desc data=&dataset;
model &y=&x;
run;
ods output close;
%mend lgstc_1var;

%macro lgstc_dummy (x,  y,  dataset);
ods output GlobalTests=GlobalTests_&x. OddsRatios=OddsRatios_&x.;
proc logistic desc data=&dataset;
class &x /param=ref ref=first;
model &y=&x;
run;
ods output close;
%mend lgstc_dummy;

%lgstc_dummy(x=ageg3, y=c04, dataset=yhl.qest31)
%lgstc_1var(x=b02, y=c04, dataset=yhl.qest31)
%lgstc_1var(x=minzu, y=c04, dataset=yhl.qest31)
%lgstc_1var(x=hunyin, y=c04, dataset=yhl.qest31)
%lgstc_1var(x=wenhua, y=c04, dataset=yhl.qest31)
%lgstc_1var(x=zhiye, y=c04, dataset=yhl.qest31)
现在work库中有了一系列的数据集,其中包括GlobalTests的值、自由度及相应的P值,就是想把数据集的名称中的后缀作为一个变量加入该数据集中,这样就可以把这些数据集合并起来,以生成最后的表格。
求解决方法。

或者对程序有什么指点,敬请指正。

回复 支持 反对

使用道具 举报

1

主题

12

帖子

103

积分

注册会员

Rank: 2

积分
103
5#
 楼主| 发表于 2014-3-26 08:18:12 | 只看该作者
可~乐 发表于 2014-3-25 22:50
不知道是否类似这样

谢谢您的指点,应该可以了。谢谢啊
回复 支持 反对

使用道具 举报

1

主题

6

帖子

45

积分

新手上路

Rank: 1

积分
45
6#
发表于 2014-3-30 23:47:45 | 只看该作者
本帖最后由 x2chen 于 2014-3-30 23:56 编辑


回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|小黑屋|手机版|Archiver|SAS中文论坛  

GMT+8, 2024-12-21 23:43 , Processed in 0.144760 second(s), 22 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表