Cuál es el valor de '$R^2$' dado en el Resumen de un modelo de coxph en R

Question

¿Cuál es el valor de $R^2$ en el Resumen de un modelo de coxph en R? Por ejemplo,

Rsquare= 0.186   (max possible= 0.991 )

Tontamente se incluye un manuscrito como un valor de $R^2$ y el revisor saltó sobre él diciendo que él no sabía de un análogo de la $R^2$ estadística de la regresión lineal clásica, desarrollando el modelo de Cox y si hubo uno por favor proporcionar una referencia. Cualquier ayuda sería genial!

Answer 1

Con getS3method("summary","coxph") usted puede ver cómo se calcula.

Las líneas de código relevante son los siguientes:

logtest <- -2 * (cox$loglik[1] - cox$loglik[2])
rval$rsq <- c(rsq = 1 - exp(-logtest/cox$n), maxrsq = 1 - 
        exp(2 * cox$loglik[1]/cox$n))

Aquí cox$loglik es "un vector de longitud 2 que contiene el log-verosimilitud con los valores iniciales y con los últimos valores de los coeficientes" (ver ?coxph.object ) y cox$n es "número de observaciones utilizado en el ajuste".

Answer 2

Dividiendo por $n$ el número de observaciones en el Resumen de coxph es incorrecto, debe ser el número de eventos sin censura, ver O'Quigley et al (2005) estadísticas de aleatoriedad en los modelos de riesgos proporcionales, explicó en medicina p. 479-489.