我想循环比较几个逻辑回归模型。 并在一个图表中打印几个ROC曲线。 我怎么能这样做? 例如,对于一个模型,它看起来像
pred_fit4 <- prediction(workData$pred,workData3$value)
perf_fit4 <- performance(pred_fit4,"tpr","fpr")
roc <- plot(perf_fit4,colorize = T, print.cutoffs.at = seq(0,1),)
及其曲线输出如下:
我该如何添加新行? 它看起来像
for (...) {
.......
.......
pred_fit4 <- prediction(workData$pred,workData3$value)
perf_fit4 <- performance(pred_fit4,"tpr","fpr")
#need add line on exist chart
roc <- plot(perf_fit4,colorize = T , print.cutoffs.at = seq(0,1),)
}
答案 0 :(得分:1)
这是使用包pROC。您可以创建一个ROC对象,然后使用行添加到同一个图形中,并根据需要添加任意数量的ROC。
obj.rpart1=roc(label1 ~ p1.rpart, smooth = T)
obj.rpart2=roc(label1 ~ p1.rpart, smooth = T)
plot(obj.rpart)
lines(obj.rpart2)