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

【问题标题】：Compute AIC in Survival Analysis (survfit/coxph)在生存分析中计算 AIC (survfit/coxph)
【发布时间】：2013-10-30 10:02:44
【问题描述】：

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

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

据我了解，这是软件限制。谁能帮我解决这个问题而无需手动计算 AIC 值？

【问题讨论】：

【解决方案1】：

函数extractAIC有一个coxph的方法

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

【讨论】：

【解决方案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)

【讨论】：