在生存分析中计算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)