r - R における Cox モデルの尤度比検定

Question

coxphを使用して作成されたcox比例ハザードモデルで機能する、lmtestパッケージのlrtestのような尤度比検定を知っている人はいますか? lrtest は coxph モデルでは機能しないようです。

ありがとう

Answer 1

coxph()LR-Test は、生存パッケージからデフォルトで返されます(最後の行を参照)。

require(survival)
test1 <- list(time=c(4,3,1,1,2,2,3), 
              status=c(1,1,1,0,1,1,0), 
              x=c(0,2,1,1,1,0,0), 
              sex=c(0,0,0,0,1,1,1)) 
# Fit a stratified model 
coxph(Surv(time, status) ~ x + strata(sex), test1) 

Call:
coxph(formula = Surv(time, status) ~ x + strata(sex), data = test1)


   coef exp(coef) se(coef)     z    p
x 0.802      2.23    0.822 0.976 0.33

Likelihood ratio test=1.09  on 1 df, p=0.297  n= 7, number of events= 5