9

我想为我的生存分析对象（suvfit/coxph）计算 AIC 值。当我尝试这样做时，它说：

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

据我了解，这是软件限制。任何人都可以在不手动计算 AIC 值的情况下帮助我解决这个问题吗？

4

2 回答 2

10

该函数extractAIC有一个方法coxph

fit <- coxph(Surv(time, status) ~ sex, data = cancer)
extractAIC(fit)
于 2014-10-06T08:27:52.837 回答
4

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

p <- 0;k <- 3 
AIC0 <- -2*cox$loglik[1] + 2*(p+k)
p <- length(cox$coef)
AIC1 <- -2*cox$loglik[2] + 2*(p+k)
于 2013-10-30T10:23:29.420 回答