我有三个自变量和一个二元结果变量。我正在使用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)并且彼此相邻绘图,例如:
你试过使用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")
您仍然需要手动添加额外信息。