用rms包的lrm函数做回归，关于OR值的计算问题 - 经管之家

› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › R语言论坛

用rms包的lrm函数做回归，关于OR值的计算问题

27050

11

收藏 2017-06-12

我用rms包的lrm()函数做logistic回归，使用代码及结果如下：

复制代码

得到fitmodel如下：

然后进行summary，得到结果如下：

从中可以轻松得到各因子的系数，但OR值该怎样计算呢，OR值是不是在第二张图中effect的那一列呢？如果是，那应该取哪一行呢？

有人说OR值应该是e的系数次方，但这样做因为有个系数很大，OR值会变得很大，个人觉得不合理，而且系数最大的那一个因子单独做预测时效果是最差的，OR值最大也不合理。

希望有高手给予指导解答，不胜感激！

二维码

扫码加我拉你入群

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

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

全部回复

2017-6-13 04:21:55

为什么不直接用glm？

二维码

扫码加我拉你入群

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

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

2017-6-13 19:02:40

foozhencheng 发表于 2017-6-13 04:21
为什么不直接用glm？

啊，大神又是你啊

用lrm是因为我们组的人之前一直用的原因吧。

试了下glm，貌似算不了OR值的置信区间啊，而且不能画nomogram

二维码

扫码加我拉你入群

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

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

2017-6-17 12:40:42

newcafans 发表于 2017-6-13 19:02
啊，大神又是你啊

用lrm是因为我们组的人之前一直用的原因吧。

有道理，不过根据回归的结果可以计算log-odds的置信区间，再转换过去就可以了。nomogram是什么？

二维码

扫码加我拉你入群

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

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

2017-7-26 15:35:24

newcafans 发表于 2017-6-13 19:02
啊，大神又是你啊

用lrm是因为我们组的人之前一直用的原因吧。

请问rms包里面的lrm适合什么情况，有没有例子呢

二维码

扫码加我拉你入群

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

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

2017-8-15 17:21:17

foozhencheng 发表于 2017-6-17 12:40
有道理，不过根据回归的结果可以计算log-odds的置信区间，再转换过去就可以了。nomogram是什么？

nomogram不是列线图吗？

不过我们老板后来说不画nomogram了，我也就没管了，也没上来论坛看了

二维码

扫码加我拉你入群

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

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

点击查看更多内容…

2017-10-24 23:43:36

同一个数据，我分别用这两个函数：lrm()和glm()试了一下，拟合logistic 回归模型，自变量的系数好像不太一样啊，但是p值差不多。是否意味着随便用哪个函数都可以呢？

二维码

扫码加我拉你入群

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

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

2017-11-17 22:59:57

rms包的lrm( )函数做完logistic回归后，用options(prType='plain')，print(fitlrm, coefs=4)显示自变量的系数，然后用exp( )求出OR值。用summary( )好像不对

二维码

扫码加我拉你入群

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

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

2017-12-30 16:18:41

wanglao 发表于 2017-11-17 22:59
rms包的lrm( )函数做完logistic回归后，用options(prType='plain')，print(fitlrm, coefs=4)显示自变量的系 ...

请问能否具体讲一下option和print要写在哪里吗试了一下加在后面不行啊

二维码

扫码加我拉你入群

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

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

2018-6-10 12:51:44

sunflower_sara 发表于 2017-12-30 16:18
请问能否具体讲一下option和print要写在哪里吗试了一下加在后面不行啊

写在拟合模型表达式的后面。显示的类型在options()里面修改，下面是显示为latex格式的。例如
f <- lrm(ch ~ age)
options(prType='latex')
print(f, coefs=4) # write latex code to console

二维码

扫码加我拉你入群

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

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

2018-6-10 12:53:41

sunflower_sara 发表于 2017-12-30 16:18
请问能否具体讲一下option和print要写在哪里吗试了一下加在后面不行啊

options不是option,不能少了一个s

二维码

扫码加我拉你入群

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

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

一个人狂欢

2020-5-7 16:53:40

线性回归很多包都可以做，你这个没用过，试试exp（coef）。别的包都是这样输入or值的

二维码

扫码加我拉你入群

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

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

相关推荐

栏目导航

热门文章

推荐文章

扫码加好友，拉您进群

各岗位、行业、专业交流群