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

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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

1. fitmodel<-lrm(y~A1+A2+A3,data=train,x=T,y=T)

复制代码

得到fitmodel如下：

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


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

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

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

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：计算问题 RMS logistic回归 logistic logisti

为什么不直接用glm？

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

啊，大神又是你啊

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

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

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

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

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

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

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

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

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

nomogram不是列线图吗？

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

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

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

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

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

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