%include  "G:\rocplot.sas";
data fqq.girl;
set fqq.pbf_roc2;
proc logistic data=fqq.girl;
ods listing gpath="g:\结果" style=journal image_dpi=300; /*journal 为灰度 dpi为分辨率 gpath="g:\结果"为图片在自己电脑保存的地址*/
ods graphics/reset
imagename="rocgirl" /*保存图片的名字*/
outputfmt=tiff;/*图片的保持格式 许多英文杂志要求为tiff 或eps 但sas好像不支持eps 跑了几遍都没成功*/
title "ROC Curve for girl";/*图片的标题*/
model CVDs(event="1")=pbf/outroc=fqq.rocg roceps=0;
output out=outp p=phat;
run;
%rocplot(inroc=fqq.rocg,inpred=outp,p=phat,id=pbf _cutpt_  _opty_,optcrit= youden,  x=pbf);