我正在做區間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