我在R中运行了逻辑回归,并获得了暴露增加一个点时事件的比值比,我想找出事件比值比大于1时的暴露值。
我的结果是一个二元变量事件是或否。我的暴露是一个连续的数字变量,范围从1.30到29.25。
逻辑回归方程如下所示:glm(事件~暴露+协变量1+协变量2,族=二项式(链接=“logit”),数据=数据)
是否有方法提取与大于1的事件的比值比相关的暴露值?
我想做一些类似于这篇论文https://www.nejm.org/doi/full/10.1056/NEJMoa2000803的事情,其中他们计算了与死亡风险比大于1相关的术后肌钙蛋白水平。
我一直被如何处理这个问题难住了,所以还没有在R中尝试任何东西。
1条答案
按热度按时间2w2cym1i1#
欢迎来到堆栈溢出,我不知道这是一个编程问题还是一个统计问题,R中可能有一个库可以自动得到这个,但是我会这样做。
首先,我不想讨论我的大脑不喜欢的比值比,我将讨论对数比值比,它很有用,因为它是逻辑回归在R中默认报告的尺度。比值比1等于对数比值比0。首先,我模拟了一些数据:
现在我来拟合一个模型,得到系数:
让我们回顾一下我们的模型的简单版本,它是基于截距和两个预测变量的组合来预测结果的对数比值。您想知道对数比值为零时预测变量
x
的值,因此我们有以下等式:我们可以代入系数,然后就很简单了,只需解出x
当我们从两边减去截距,然后将两边除以估计值时,我们就得到了x:
用计算器(或R)计算得到
-0.3423612
,我们可以通过画出x对y的影响来验证,看看这是否是y = 0.5时x的值:对协变量进行此操作将更加复杂,因为它将取决于协变量的值。