在生存分析中计算AIC(survfit/coxph)

Question

我想为我的生存分析对象(suvfit/coxph)计算AIC值.当我尝试这样做时,它说:

> AIC(cox)
  Error in UseMethod("logLik") : 
   no applicable method for 'logLik' applied to an object of class "coxph"

据我所知,这是一个软件限制.任何人都可以帮助我解决这个问题,而无需手动计算AIC值吗？

Answer 1

该函数extractAIC有一个方法coxph

fit <- coxph(Surv(time, status) ~ sex, data = cancer)
extractAIC(fit)

Answer 2

我记得我是手工计算的，就像这样

p <- 0;k <- 3 
AIC0 <- -2*cox$loglik[1] + 2*(p+k)
p <- length(cox$coef)
AIC1 <- -2*cox$loglik[2] + 2*(p+k)