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:
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.