我拟合了cox回归,并且对预测变量x的影响感兴趣,该变量是模型中的最后一个变量(变量7)。我在模型中加入了一些变量,例如性别和年龄,因为我想为其调整模型。
使用R函数coxph()会给我警告“ Loglik在变量3之前收敛”。实际上,我对变量3不感兴趣,因为它只是我要调整的变量之一。但是,“之前”一词让我怀疑这是否意味着变量3之后的所有变量(包括我的预测变量x)的结果都是无效的。还是只影响变量3的结果?
更多信息:实际上,我正在运行多个cox回归,并且所描述的警告仅在某些变量3的模型中发生。我确实想对变量3进行调整,从而将其保留在代码中。
有关此警告(1)的一些讨论,但是到目前为止,我还没有找到答案。谢谢。
答案 0 :(得分:1)
感谢您对以前的答案进行搜索。 (并且感谢您引用我的回答之一:-)警告消息仅适用于control3变量。调查关于predictorx
的统计推断的有效性的关键在于您引用的Therneau答案的一部分。您对其他变量之一感兴趣,幸运的是,您只有一个变量“爆炸了”。这意味着您可以进行LRT比较具有和不具有目标变量的模型,以获得适当的统计结果。结果从本质上说,在子集中没有任何事件具有正值“ control3”。我猜想control3是一个0/1变量,如果您看一下表格:
with( your_data, table(control3, your_status_variable))
....在control3 == 1的情况下,您将发现零个事件。