全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
5903 4
2017-03-31
求问各位大神,sas的广义线性模型中如何输出pr值(或rr值)
以下是小弟的代码,求大神给补一行可以输出pr值(或rr值)的代码
挺急的,在线等
proc genmod data=aa3;
model f4a1 = f1 f2qq f3
/d=bin link=log intercept=-1 lrci wald;
run;

二维码

扫码加我 拉你入群

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

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

全部回复
2017-4-4 10:29:01
model的parameter estimate直接取exp不就可以了吗?比如说parameter estimate是1.345,那么RR就是exp(1.345)=3.838.
二维码

扫码加我 拉你入群

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

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

2017-4-11 10:18:54
sonypam 发表于 2017-4-4 10:29
model的parameter estimate直接取exp不就可以了吗?比如说parameter estimate是1.345,那么RR就是exp(1.34 ...
对对对,没错
二维码

扫码加我 拉你入群

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

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

2018-11-12 15:52:17
楼主,最后是怎么求得OR值啊?如果自变量是多分类有怎么求呢?
二维码

扫码加我 拉你入群

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

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

2024-5-29 16:41:25
在SAS的PROC GENMOD中,你不能直接输出PR值(预测概率)或RR值(风险比)。但你可以通过以下步骤来获取这些值:

1. 输出预测概率(PR值):
   在MODEL语句后添加一个新的语句:OUTPUT OUT=PredProb PRED=P;
   这将创建一个新的数据集“PredProb”,其中包含每个观察的预测概率(PRED变量)。

2. 计算风险比(RR值):
   你可以通过在输出的数据集中计算暴露组和非暴露组的事件率,然后计算两者的比率来得到RR值。例如:

   ```sas
   proc means data=PredProb;
     var PRED;
     by f1; /* 替换为你的暴露变量 */
   run;

   proc sql;
     create table RRValues as
     select a.f1,
            avg(a.PRED) / avg(b.PRED) as RR
     from PredProb a
     join PredProb b on (a.f1 ne b.f1))
     group by a.f1;
   quit;
   ```

这段代码首先计算了每个暴露组的事件率,然后通过JOIN操作计算了非暴露组的事件率,并计算了两者的RR值。结果将存储在新的数据集“RRValues”中。

此文本由CAIE学术大模型生成,添加下方二维码,优先体验功能试用



二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群