我需要在单个图上绘制许多不同的ROC曲线。为了避免手动创建每个 ROC 曲线,我创建了一个自动执行此过程的 for 循环。但是,由于某种原因,代码仅输出名称列表中最后一个模型的单个曲线。谁能帮我弄清楚为什么它不起作用?请参阅下面的可重现示例:
library(pROC)
library(tidyverse)
dat_tst_2 <- data.frame(result = sample(letters[1:2], 100, replace = T))
preds_1 <- data.frame(x = runif(100),
y = runif(100))
preds_2 <- data.frame(x = runif(100),
y = runif(100))
names_preds <- c("preds_1", "preds_2")
output <- list()
for (j in 1:length(names_preds)) {
for (i in names_preds) {
roc_model <- roc(response = dat_tst_2$result,
predictor = eval(as.name(i))[,2],
levels = c("a", "b"),
plot = F)
output[[j]] <- roc_model
}
}
ggroc(output)
首先确保输出有多个项目 usin str(output(。然后尝试将输出中的每个项目传递给 ggroc:
lappy( output, function (out) { png()
print (ggroc(out))
dev.off() }