如何在图中绘制和识别未着色的ROC曲线?

时间:2016-01-06 12:41:29

标签: r plot

下面我的代码用于绘制x,y和z的ROC曲线单独工作,但我想在单个图中使用ROC曲线。此外,由于该图将在未着色的日记中发布,读者应该能够通过不同的行(虚线等)或显示相应行的箭头来识别它们。

有什么建议吗?

set.seed(1)
x <- rnorm(55)
y <- rnorm(55)
z <- rnorm(55)

data <- data(x,y,z)
status <- c(rep(1,25),rep(0,30))

model1 <- glm(status ~  x, family = binomial(link="logit") )
model2 <- glm(status ~  y, family = binomial(link="logit") )
model3 <- glm(status ~  z, family = binomial(link="logit") )

install.packages("pROC")
library(pROC)

prob1 <- predict(model1,type=c("response"))
data$prob1 <- prob1
plot1 <- roc(status ~ prob1)

prob2 <- predict(model2,type=c("response"))
data$prob2 <- prob2
plot2 <- roc(status ~ prob2)
plot(plot2,type="b")

prob3 <- predict(model3,type=c("response"))
data$prob3 <- prob3
plot3 <- roc(status ~ prob3)
plot(plot3)

0 个答案:

没有答案