3

Tengo un modelo logit binomial y estoy evaluando la bondad del ajuste de mi modelo mediante la curva de ROC:

library(pROC)
curva_ROC <- roc(modelo_logit_viv$y,modelo_logit_viv$fitted.values)
plot(curva_ROC,col = "blue",xlim = c(1,0),ylim = c(0,1),
     xlab = "Especifidad", ylab ="Sensibilidad", main = "Curva de ROC")
auc(curva_ROC) # 0.8445

Esta es la gráfica que me devuelve:

Curva de ROC

Quiero modificarla de la siguiente manera:

  • Acotar el eje x (especifidad) de tal manera que el inicio esté casi en 1 y acabe en 0, es decir, que los límites del eje horizontal coincidan con la diagonal (pensé que al indicar xlab = c(0,1) obtendría ese resultado, pero la gráfica no varía al incluir este argumento en la función plot).

  • Que el área comprendida entre la curva y la diagonal esté coloreada en un azul más claro que la curva.

  • Que aparezca el valor del área auc(curva_ROC) justo en el centro de dicha área.

1 respuesta 1

1

1. Ajustar ejes

La forma de hacerlo en las gráficas base es mediante los parámetros xlim e ylim no pude verificarlo con tus datos, pero funciona adecuadamente con el ejemplo que armé para esta respuesta. Por otro lado xlab e ylab solo establecen los "labels" o etiquetas de cada eje

2. Colorear área

El comportamiento de plot.ROC() admite colorear el área completa de la curva, mediante estos parámetros (ejemplo):

auc.polygon = TRUE 
auc.polygon.col=rgb(.35,0.31,0.61, alpha = 0.4)

Sin embargo, lo que buscas, es colorear la curva hasta la diagonal. Una forma podría ser, generar el polígono por debajo de la curva y hasta la diagonal, aprovechando los datos del propio plot.ROC()

polygon(with(p, cbind(specificities, sensitivities)), 
        col = rgb(.25,0.31,0.61, alpha = 0.4), 
        border = "blue",
        lwd = 2)

3. Texto con el valor AUC

Lo puedes configurar y ajustar con los siguientes parámetros:

print.auc = TRUE,
print.auc.x = 0.73,
print.auc.y = 0.6

Los valores de x e y los deberías ajustar para que el texto entre dentro del área.

Ejemplo

library(pROC)

curva_ROC = roc(aSAH$outcome, aSAH$s100b, levels=c("Good", "Poor"))
p <- plot(curva_ROC,
          col = "blue",
          print.auc = TRUE,
          print.auc.x = 0.73,
          print.auc.y = 0.6, 
          xlim = c(1,0),
          ylim = c(0,1),
          xlab = "Especifidad", 
          ylab ="Sensibilidad", 
          main = "Curva de ROC")

polygon(with(p, cbind(specificities, sensitivities)), 
        col = rgb(.25,0.31,0.61, alpha = 0.4), 
        border = "blue",
        lwd = 2)

Resultado:

introducir la descripción de la imagen aquí

Fuentes:

3
  • Funciona a la perfección ¿Cómo podría hacer para que solo aparezca el área de la curva sin el "AUC: "?
    – Juan Luis
    Commented el 27 abr. 2020 a las 19:05
  • 1
    Puedes agregar print.auc.pattern = "%f" Commented el 27 abr. 2020 a las 19:43
  • ¿Por qué las curvas ROC presentan el eje y "invertido" (del 1 al 0 y no al revés)? Commented el 7 jun. 2021 a las 19:56

Tu Respuesta

By clicking “Publica tu respuesta”, you agree to our terms of service and acknowledge you have read our privacy policy.

¿No es la respuesta que buscas? Examina otras preguntas con la etiqueta o formula tu propia pregunta.