我正在使用区间cesored数据进行生存分析,并且我正在使用来自R的intcox包中的intcox()函数,该函数基于coxph函数。
该函数返回没有似然比测试值的输出:
> intcox(surv~sexo,data=dados)
Call:
intcox(formula = surv ~ sexo, data = dados)
coef exp(coef) se(coef) z p
sexojuvenil 2.596 13.4 NA NA NA
sexomacho -0.105 0.9 NA NA NA
Likelihood ratio test=NA on 2 df, p=NA n= 156
我不知道为什么会发生这种情况......以下是将coxph()函数应用于相同的数据:
> coxph(Surv(dias_seg,status)~sexo,data=dados)
Call:
coxph(formula = Surv(dias_seg, status) ~ sexo, data = dados)
coef exp(coef) se(coef) z p
sexojuvenil 2.320 10.172 0.630 3.684 0.00023
sexomacho -0.169 0.844 0.252 -0.671 0.50000
Likelihood ratio test=9.28 on 2 df, p=0.00967 n= 156, number of events= 77
str(dados$sexo)
Factor w/ 3 levels "femea","juvenil",..: 3 3 3 3 3 3 3 3 3 3 ...
你能帮我解决这个问题吗? 提前谢谢。
答案 0 :(得分:1)
Volkmar Henschel(intcox软件包的作者)告诉我,“带有intcox的文件给出了一个类的对象”coxph“没有回归系数的标准误差”。
有关此文档的更多说明: ftp://ftp.auckland.ac.nz/pub/software/CRAN/doc/vignettes/intcox/intcox.pdf