Почему я получаю ошибку в тесте ANOVA на моделях Кокса в R? - PullRequest
Новая
       78

Почему я получаю ошибку в тесте ANOVA на моделях Кокса в R?

0 голосов
/ 28 октября 2019

Я хочу сделать тест отношения правдоподобия между двумя моделями в R.

  • Вложены 2 модели
  • Я удалил все недостающие элементы, создав чистый кадр данных.

Хотя, когда я anova(mod2,mod3), я получаю это сообщение об ошибке:

Ошибка в anova.coxphlist (c (list (object), dotargs), test =test):
модели не были приспособлены к одному и тому же размеру набора данных. Дополнительно: предупреждающее сообщение: в anova.coxphlist (c (list (object), dotargs), test = test): модели с ответом c ("sep= \ "\")))) "," sep = \ "\")))) ") удалено, поскольку ответ отличается от модели 1

Извините, я не могу дать образецданные, которые я использовал. Ниже я приведу краткое изложение моих 2 моделей: 

> mod2
Call:
coxph(formula = Surv(eval(parse(text = paste(j, "_suivi", sep = ""))), 
    eval(parse(text = paste(j, sep = "")))) ~ strata(trt) + Var1 + 
    Var2 + Var3 + rcs(Var4, 4), 
    data = data.clean)

                      coef exp(coef)  se(coef)      z       p
Var1               0.030361  1.030827  0.012188  2.491 0.01274
Var2               0.022405  1.022658  0.007463  3.002 0.00268
Var3:1             0.233716  1.263286  0.433950  0.539 0.59018
rcs(Var4, 4)Var4   0.315988  1.371613  0.111433  2.836 0.00457
rcs(Var4, 4)Var4' -0.799272  0.449656  0.378427 -2.112 0.03468

Likelihood ratio test=37.7  on 5 df, p=4.329e-07
n= 375, number of events= 74 
> mod3
Call:
coxph(formula = Surv(eval(parse(text = paste(j, "_followup", sep = ""))), 
    eval(parse(text = paste(j, sep = "")))) ~ strata(trt) + Var1 + 
    Var2 +Var3 + rcs(Var4, 4) + 
    eval(parse(text = paste(i, sep = ""))), data = data.clean)

                        coef exp(coef)  se(coef)      z       p
Var1                0.024004  1.024295  0.012742  1.884 0.05959
Var2                0.018190  1.018356  0.007881  2.308 0.02099
Var3:1              0.348887  1.417489  0.437326  0.798 0.42500
rcs(Var4, 4)Var4    0.289575  1.335859  0.110936  2.610 0.00905
rcs(Var4, 4)Var4'  -0.705976  0.493627  0.377903 -1.868 0.06174
eval(parse(text = paste(i, sep = "")))  0.742461  2.101100  0.304491  2.438 0.01475

Likelihood ratio test=43.69  on 6 df, p=8.499e-08
n= 375, number of events= 74 



>  anova(mod2,mod3)
    Error in anova.coxphlist(c(list(object), dotargs), test = test) : 
      models were not all fitted to the same size of dataset
    In addition: Warning message:
    In anova.coxphlist(c(list(object), dotargs), test = test) :
      Models with response c("    sep = \"\"))))", "    sep = \"\"))))") removed because response differs from model 1

Может быть, это связано с наличием стратового термина? Или наличие шлицевой версии Var4?

Спасибо за ваше предложение!

...