参考组事件为零时的 Cox 回归

机器算法验证 回归 spss 生存
2022-04-04 23:11:17

对于我遇到的问题,我将不胜感激。我使用 SPSS 对一项研究进行统计分析。

该研究着眼于血液测试如何预测患者随访 1 年的死亡率。将患者分为四分位数,并使用第一个四分位数作为参考组,我使用 Cox 回归来确定未调整的危险比 (HR) 和四分位数 2-4 的调整后 HR。

问题是参考组(四分位数 1——前 25%)没有任何事件。因此,我在 SPSS 下获得了 80,000+ 的第二、第三和第四四分位数的风险比以及错误消息。

谁能告诉我我做错了什么?

2个回答

好吧,您做错的是使用零事件组作为参考组。而不是风险比,用更简单的术语(在我看来)考虑事件率比率(IRR),其中事件率(IR)是IR=number of cases / total person-time

IRRquartile 4 vs. quartile 1=IRquartile 4IRquartile 1

如果会发生什么?IRquartile 1=0

您可以更改您的分类(使用三分位数或其他一些有意义的分类),或者更好的是,如果您有一个连续预测器,您可以将其视为此类,并使用多项式项、分数多项式或受限三次样条检查潜在的非线性关系,例如。

通过将其扩展到风险比来补充安德里亚的反应:

事件的危害是事件在时间 t 发生的瞬时概率,条件是它以前没有发生过。

您的问题应该立即清楚 - 没有事件,概率为零。借用安德里亚的例子,事故率相当于一个恒定的危险——在你的例子中,一个恒定的危险为零。

除以零往往会使软件生气。

您需要切换您的参考类别。我的建议是使用“四分位 4”或该类别的其他高值,然后降级,而不是使用四分位 1 并提高。例如,如果您希望在您上移一个类别时显示 HR 增加,那么您现在显示出与下移一个相同的保护效果。

我还建议花点时间考虑一下为什么您没有活动。

您可能只是遇到了“运气不好”,此时您无能为力,只能增加研究规模或跟踪人群更长时间以希望积累更多事件。但是你应该确保没有理由在你的人口中产生结果的概率不为零是有原因的对于心脏事件,我可以想象一个,但是当您在某个协变量水平上发生零事件时,总是值得停下来考虑。