用ggplot2绘制多个ROC曲线在不同的层

Question

我想绘制与ggplot2在一个绘图上的多个ROC曲线。下面是我在多大程度上得到：

ggroc2 <- function(columns, data = mtcars, classification = "am", 
        interval = 0.2, breaks = seq(0, 1, interval)){ 
    require(pROC) 
    require(ggplot2) 

    #The frame for the plot 
    g <- ggplot() + geom_segment(aes(x = 0, y = 1, xend = 1,yend = 0)) + 
    scale_x_reverse(name = "Specificity",limits = c(1,0), breaks = breaks, 
expand = c(0.001,0.001)) + 
    scale_y_continuous(name = "Sensitivity", limits = c(0,1), breaks = 
breaks, expand = c(0.001, 0.001)) + 
    theme_classic() + coord_equal() 

    #The loop to calculate ROC's and add them as new layers 
    for(i in 1:length(columns)){ 
    croc <- roc(data[,classification], data[,columns[i]]) 
    plotx <- rev(croc$specificities) 
    ploty <- rev(croc$sensitivities) 
    g <- g + geom_step(aes(x=plotx, y=ploty)) 
    } 

    g 
} 



#Sample graph 
ggroc2(c("mpg", "disp", "drat", "wt")) 

的问题是，只有在columns列表中的最后一个参数被绘制。我发现问题必须与aes()有关，并在阅读the answer to this question后进行懒惰评估。该示例使用geom_segment()，并且在完全删除aes()后问题得到解决。它不适合我，因为我需要以某种方式映射数据。当我在这里删除aes()时，没有任何东西被绘制。如何解决geom_中依赖于aes()的懒惰评估问题？

Answer 1

这是您的代码的工作版本。
最终的图形结果不是很好，应该改进。

ggroc2 <- function(columns, data = mtcars, classification = "am", 
        interval = 0.2, breaks = seq(0, 1, interval)){ 
    require(pROC) 
    require(ggplot2) 

    #The frame for the plot 
    g <- ggplot() + geom_segment(aes(x = 0, y = 1, xend = 1,yend = 0)) + 
    scale_x_reverse(name = "Specificity",limits = c(1,0), breaks = breaks, 
expand = c(0.001,0.001)) + 
    scale_y_continuous(name = "Sensitivity", limits = c(0,1), breaks = 
breaks, expand = c(0.001, 0.001)) + 
    theme_classic() + coord_equal() 

    #The loop to calculate ROC's and add them as new layers 
    cols <- palette() 
    for(i in 1:length(columns)){ 
    croc <- roc(data[,classification], data[,columns[i]]) 
    sens_spec <- data.frame(spec=rev(croc$specificities), 
          sens=rev(croc$sensitivities)) 
    g <- g + geom_step(aes(x=spec, y=sens), data=sens_spec, col=cols[i], lwd=1) 
    } 
    g 
} 

#Sample graph 
ggroc2(c("mpg", "disp", "drat", "wt"))