Calcolo AIC in analisi di sopravvivenza (survfit/coxph)

Question

Voglio calcolare il valore AIC per i miei oggetti di analisi di sopravvivenza (suvfit/coxph). Quando provo a farlo, dice:

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

Per quello che ho capito che è una limitazione software. Qualcuno può aiutarmi a risolvere questo problema senza calcolare a mano il valore AIC?

Answer 1

La funzione extractAIC ha un metodo per coxph

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

Answer 2

mi ricordo che ero Computing a mano, come questo

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