Я построил несколько кривых RO C с помощью pROC
, и теперь я хочу добавить общую легенду по горизонтали через нижний центр пространства графика. Обычно я использую ggplot
, и теперь я теряюсь при построении базы R.
Вот мой код:
# Set plotting parameters
par(mfrow=c(1,2), pty = "s", oma = c(4,1,1,1))
# ROC plot 1
roc(logit_A$y, logit_A$fitted.values, plot=TRUE, legacy.axes=FALSE, percent=TRUE, col="salmon", lwd=2, print.auc=TRUE)
plot.roc(logit_B$y, logit_B$fitted.values, percent=TRUE, col="goldenrod", lwd=2, print.auc=TRUE, add=TRUE, print.auc.y=40)
plot.roc(logit_C$y, logit_C$fitted.values, percent=TRUE, col="lightsteelblue", lwd=2, print.auc=TRUE, add=TRUE, print.auc.y=30)
title(main = "Model 1", line = 2.5)
# ROC plot 2
roc(logit_A2$y, logit_A2$fitted.values, plot=TRUE, legacy.axes=FALSE, percent=TRUE, col="salmon", lwd=2, print.auc=TRUE)
plot.roc(logit_B2$y, logit_B2$fitted.values, percent=TRUE, col="goldenrod", lwd=2, print.auc=TRUE, add=TRUE, print.auc.y=40)
plot.roc(logit_C2$y, logit_C2$fitted.values, percent=TRUE, col="lightsteelblue", lwd=2, print.auc=TRUE, add=TRUE, print.auc.y=30)
title(main = "Model 2", line = 2.5)
# Add legend
legend("bottom",
legend=c("A", "B", "C"),
col=c("salmon", "goldenrod", "lightsteelblue"),
lwd=4, cex =0.4, xpd = TRUE, horiz = TRUE)
Как мне это сделать так легенда центрирована вдоль дна между обоими сюжетами?