pROC участок как объект?

Question

Я хотел бы построить объект OBJECT rPOC, который я могу передать функции и заговорить позже. В приведенном ниже примере (взятом из here), я могу, конечно, построить график кривой. Но я хотел бы создать сюжет OBJECT (скажем, определяя g < - и позже строя его с графиком (g)). Кажется, что график функций ciobj, ci ниже добавит к исходному сюжету, но я не могу сделать OBJECT, собирающий эти слои вместе. Я попробовал аргумент «добавить» и создал новые объекты сюжета с возвращаемыми значениями этих функций графика.

library(pROC) 
data(aSAH) 
rocobj <- plot.roc(aSAH$outcome, aSAH$s100b, main="Confidence intervals", percent=TRUE, ci=TRUE, print.auc=TRUE) 
ciobj <- ci.se(rocobj, specificities=seq(0, 100, 5)) 
plot(ciobj, type="shape", col="#1c61b6AA") 
plot(ci(rocobj, of="thresholds", thresholds="best")) 

Answer 1

Как отметил г-н Флик, вы можете передавать функции вместо объекта. В качестве альтернативы вы можете передавать неоцененные вызовы функций и оценивать их в своей функции. Например:

library(pROC) 
data(aSAH) 


plot1 <- quote(plot(rocobj, main="Confidence intervals", print.auc=TRUE)) 
plot2 <- quote(plot(ci.se(rocobj, specificities=seq(0, 100, 5)), type="shape", col="#1c61b6AA")) 
plot3 <- quote(plot(ci(rocobj, of="thresholds", thresholds="best"))) 

doplots <- function(rocobj, calls) { 
    for (call in calls) { 
    eval(call) 
    } 
    invisible(rocobj) 
} 

roc1 <- roc(aSAH$outcome, aSAH$s100b, percent = TRUE) 
doplots(roc1, list(plot1, plot2, plot3)) 

roc2 <- roc(aSAH$outcome, aSAH$wfns, percent = TRUE) 
doplots(roc2, list(plot1, plot3)) 

Нет ограничений на то, что вы можете сделать, кроме вашего терпения с правилами оценки R.