下面我的代码用于绘制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)