Question

通过我与专家的谈话（Link1，Link2），我认为我需要进行条件二元逻辑回归。我已经对数据进行了排序，每个控件都紧跟着对应的处理案例，数据格式也很长。

我的设计详见上述两个链接以及另外两个链接：Link3，Link4

编辑：后来我发现排序并不重要，因为排序和非排序数据导致相同的输出。

我的语法如下。代码运行但结果很奇怪。大多数SE为零，大多数P值为NaN。有什么问题？

library(Epi)

clogistic((DV ~ (Demo1 +Demo2 +Demo3 +Demo4 +Demo5)^2), 
          strata = PatientID,  data = Data4)

输出结果为：

> clogistic((DV ~ (Demo1 +Demo2 +Demo3 +Demo4 +Trt)^2), strata = PatientID,  data = MixedModelData4)

Call: 
clogistic(formula = (DV ~ (Demo1 + Demo2 + Demo3 + Demo4 + Trt)^2), 
    strata = PatientID, data = MixedModelData4)




                 coef exp(coef) se(coef)         z   p
Demo1        0.00e+00  1.00e+00        0       NaN NaN
Demo2        0.00e+00  1.00e+00        0       NaN NaN
Demo3       -3.27e-09  1.00e+00    56013 -5.83e-14   1
Demo4        0.00e+00  1.00e+00        0       NaN NaN
Trt         -2.12e+01  6.19e-10    14786 -1.43e-03   1
Demo1:Demo2  0.00e+00  1.00e+00        0       NaN NaN
Demo1:Demo3  0.00e+00  1.00e+00        0       NaN NaN
Demo1:Demo4  0.00e+00  1.00e+00        0       NaN NaN
Demo1:Trt   -4.34e-08  1.00e+00    50351 -8.62e-13   1
Demo2:Demo3  0.00e+00  1.00e+00        0       NaN NaN
Demo2:Demo4  0.00e+00  1.00e+00        0       NaN NaN
Demo2:Trt   -1.19e-08  1.00e+00    12937 -9.20e-13   1
Demo3:Demo4  0.00e+00  1.00e+00        0       NaN NaN
Demo3:Trt    8.08e-09  1.00e+00    19595  4.12e-13   1
Demo4:Trt   -1.62e-08  1.00e+00    31612 -5.12e-13   1

Likelihood ratio test=13.9  on 15 df, p=0.536, n=20
Warning message:
In clogistic((DV ~ (Demo1 + Demo2 + Demo3 + Demo4 + Trt)^2), strata = PatientID,  :
  Iteration limit exceeded

Answer 1

当系数或其标准误差“爆炸”时，它通常表示病理数据情况。可能存在完全分离或严重的多重共线性（或两者）。您需要同时提供：

 str(MixedModelData4)

和

with( MixedModelData4, table(DV,Demo1, Demo2) )
with( MixedModelData4, table(DV,Demo1, Demo3))
with( MixedModelData4, table(DV,Demo1, Demo3))
with( MixedModelData4, table(DV,Demo1, Demo4))
with( MixedModelData4, table(DV,Demo1, Demo5))

这只适合初学者。问题可能存在于数据中的其他三种方式组合之一。

R：条件二进制逻辑回归 - 语法和调试

1 个答案: