我想保存多个ROC分析的AUC值,并将它们附加在一起,这样我就可以快速获得哪个变量组合的AUC最大值的列表。
我不知道如何输出AUC值,但我只能输出roc统计数据。
ODS GRAPHICS ON;
PROC LOGISTIC data = dataset PLOTS(only) = (roc(id = obs) effect);
CLASS outcome ;
MODEL outcome = var / scale = none
clparm = wald
clodds = pl
rsquare OUTROC= RocStats;
RUN;
ODS GRAPHICS OFF;
AUC值在下面的数据集AUC
中的变量Area
中:
PROC LOGISTIC DATA = SASHELP.CLASS;
CLASS SEX;
MODEL SEX = HEIGHT WEIGHT / OUTROC = ROC;
ROC;
ODS OUTPUT ROCASSOCIATION = AUC;
RUN;
TITLE "AUROC";
PROC PRINT DATA=AUC NOOBS LABEL;
WHERE ROCMODEL = 'Model';
VAR AREA;
RUN;