sas proc genmod如何输出pr值（或者是rr值）

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求问各位大神，sas的广义线性模型中如何输出pr值（或rr值）
以下是小弟的代码，求大神给补一行可以输出pr值（或rr值）的代码
挺急的，在线等
proc genmod data=aa3;
model f4a1 = f1 f2qq f3
/d=bin link=log intercept=-1 lrci wald;
run;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：genmod Mod ROC Intercept 广义线性模型 SAS 广义线性模型 广义线性回归分析 广义线性回归 广义线性模型求教

model的parameter estimate直接取exp不就可以了吗？比如说parameter estimate是1.345，那么RR就是exp（1.345）=3.838.

sonypam 发表于 2017-4-4 10:29
model的parameter estimate直接取exp不就可以了吗？比如说parameter estimate是1.345，那么RR就是exp（1.34 ...

对对对，没错

楼主，最后是怎么求得OR值啊？如果自变量是多分类有怎么求呢？

在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学术大模型生成，添加下方二维码，优先体验功能试用