R coxph() 警告：Loglik 在变量之前收敛 - 其他变量的有效结果？

Question

我拟合 cox 回归，并且我对预测变量 x 的影响感兴趣，它是模型中的最后一个变量（变量 7）。我在模型中包含了一些变量，例如性别和年龄，因为我想为它们调整模型。

使用 R 函数 coxph() 给我警告“Loglik before variable 3”。事实上，我对变量 3 不感兴趣，因为它只是我想要调整的变量之一。但是“之前”这个词让我想知道这是否意味着变量 3（包括我的预测变量 x）之后的所有变量的结果都是无效的。还是只有变量 3 的结果受到影响？

这是输出：

更多信息：实际上，我正在运行多个 cox 回归，并且所描述的警告仅出现在变量 3 的某些模型中。我确实想针对变量 3 进行调整，从而将其保留在代码中。

关于这个警告（1）有一些讨论，但到目前为止我还没有找到我的问题的答案。谢谢。

(1) 例如here 或here

Answer 1

感谢您搜索先前的答案。 （感谢您引用我的一个答案:-) 警告消息仅适用于 control3 变量。调查关于predictorx 的统计推断的有效性的关键在于您引用的 Therneau 答案的一部分。您对其他变量之一感兴趣，幸运的是，您的变量中只有一个“爆炸”了。这意味着您可以进行 LRT 比较具有和不具有感兴趣变量的模型以获得正确的统计结果。结果基本上是说在具有正值 `control3 的案例子集中没有事件。我猜 control3 是一个 0/​​1 变量，如果您查看表格：

 with( your_data, table(control3, your_status_variable))

....在 control3 == 1 的情况下，您会发现零事件。