我有三个自变量和一个二元结果变量。我正在使用 pROC 包进行分析,并希望将所有三个 ROC 图并排显示,类似于 face_wrap 功能。例如,使用以下代码,

library(pROC)
data(aSAH)
roc_i <- roc(aSAH$outcome,
             aSAH$s100b,
             legacy.axes=TRUE,
             ci=TRUE,
             boot.n=2000,
             ci.alpha=0.9,
             stratified=TRUE,
             plot=TRUE,
             auc.polygon=TRUE,
             max.auc.polygon=TRUE,
             grid=TRUE,
             print.auc=TRUE,
             print.thres='best',
             print.thres.best.method='y',
             print.thres.adj=c(-0.05, 1.25),#,c(-0.05, 2.0),c(-0.05, 5.0)),
             print.thres.pattern="Cut-off: %.3f \n\nSp: %.3f \nSe: %.3f",
             print.thres.pattern.cex = 1.5,
             print.auc.cex = 1.5,
             print.auc.y=0.2,
             print.auc.x=0.8,
             cex.axis=1.5,
             cex.lab=1.5,
             print.thres.pch=16,
             print.thres.cex=2.0,
             cex.main=1.5)

我可以为单个独立变量 s100b 生成单个图片,而我还需要两个(对于其他 ind. Var 2 和 ind. Var 3)并彼此相邻绘图,例如:

r - face_wrap 图与 R 中的 pROC 包(几个变量)-LMLPHP

最佳答案

您是否尝试过使用 ggroc function of pROC

library(ggplot2)
library(pROC)

data(aSAH)
roc.list <- roc(outcome ~ s100b + ndka + wfns, data = aSAH)

g <- ggroc(roc.list)
g
g + facet_grid(.~name) + theme(legend.position="none")

不过,您仍然需要手动添加额外信息。

关于r - face_wrap 图与 R 中的 pROC 包(几个变量),我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/50743849/

10-12 17:50
查看更多