我正在使用pROC软件包,我想使用for循环来计算多个ROC曲线图。我的变量是作为列包含在向量中的特定列名称,我希望pROC顺序读取该向量,并使用似乎接受文本/字符的“预测变量”字段中的字符串。但是,由于出现错误,我无法正确解析该变量:
'predictor'参数应为列的名称,可以选择加引号。
这是带有SAH数据集的示例代码:
ROCvector<- c("s100b","ndka")
for (i in seq_along(ROCvector)){
a<-ROCvector[i]
pROC_obj <- roc(data=aSAH, outcome, as.character(a))
#code for output/print#
}
我已经尝试仅调用“ a”并使用函数print()或get()而没有任何结果。当然,手动编写变量(带或不带引号)是可行的。我在预测变量字段中应该使用的变量类型是否缺少某些内容?
...
pROC_obj <- roc(as.formula(paste0("outcome~",a)),data=aSAH)
print(pROC_obj)