2014-02-05 23 views
0

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

你能幫我解決這個問題嗎?
在此先感謝。

+0

是否dados數據包含丟失的數據如Na? –

+0

不在dados $ sexo或存在對象 – JMarcelino

回答

1

Volkmar Henschel(intcox軟件包的作者)告訴我:「intcox與intcox的結合給出了類」coxph「的對象,沒有迴歸係數的標準誤差。

更多的說明本文件: ftp://ftp.auckland.ac.nz/pub/software/CRAN/doc/vignettes/intcox/intcox.pdf

+1

另一個有用的間隔刪失生存分析的參考文獻R.「區間刪失數據及其實現方法指南」,統計建模2009; 9(4):259-297 」。在4.2節中,它特別討論了intcox()函數和通過ICM方法擬合的Cox模型。 – Chris

相關問題