¿Cuál es el valor dado en el resumen de un modelo coxph en R? Por ejemplo,
Rsquare= 0.186 (max possible= 0.991 )
Lo incluí tontamente como un manuscrito como un valor y el revisor saltó sobre él diciendo que no estaba al tanto de un análogo del estadístico R ^ 2 de la regresión lineal clásica que se estaba desarrollando para el modelo de Cox y si había uno, por favor Proporcionar una referencia. ¡Cualquier ayuda sería genial!
Respuestas:
Usando
getS3method("summary","coxph")
puede ver cómo se calcula.Las líneas de código relevantes son las siguientes:
Aquí
cox$loglik
está "un vector de longitud 2 que contiene la probabilidad logarítmica con los valores iniciales y con los valores finales de los coeficientes" (ver?coxph.object
) ycox$n
es "número de observaciones utilizadas en el ajuste".fuente
Dividir entre el número de observaciones en el resumen de es incorrecto, debería ser el número de eventos sin censura; ver O'Quigley et al. (2005) Explicación de la aleatoriedad en los modelos de riesgos proporcionales Estadísticas en medicina p. 479-489.norte
coxph
fuente