我正在做区间cesored数据的生存分析,我使用intcox()函数从intcox包R中,这是基于coxph函数。intcox()输出没有似然比测试
该函数返回,而不似然比检验值输出:
> 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 ...
你能帮我解决这个问题吗?
在此先感谢。
是否dados数据包含丢失的数据如Na? –
不在dados $ sexo或存在对象 – JMarcelino