data roc.boy_1;
set roc.boy;
keep glu cho tg hdl ldl sbp1 dbp1 wc bf1 age1;/*从roc.boy数据集中提取出想要的变量*/
if glu>=5.6
then glus=1;
if tg>=1.7
then tgs=1;
if hdl<1.03
then hdls=1;
if ldl>=3.4
then ldls=1;
if sbp1>=90
then sbp1s=1;
if dbp1>=90
then dbp1s=1;
CVD=glus+tgs+hdls+ldls+sbpls+dbpls;
if CVD>=3
then CVDs=1;/*上述步骤为挑选出CVD病例条件*/
proc print data=roc.boy_1;
run;
ods graphics on;
%let _ROC_XAXISOPTS_LABEL=Specificity;
%let _ROC_YAXISOPTS_LABEL=Sensitivity;
proc logistic data=roc.boy_1 plots(only)=roc(id=obs);/*生成ROC曲线*/
/*model语句为重点，这个确定了ROC曲线就可以做出来了*缺乏心血管疾病的定义变量*/
   model CVDs (EVENT=1)=bf1 scale=none    /*产生偏差和peasrson拟合优度分析*/
                         clparm=wald    /*回归系数的可信区间*/
                         clodds=pl      /*产生 OR置信区间profile-likelihood*/
                         rsquare;       /*拟合模型的R2*/
   effectplot;/*显示模型拟合程度*/
run;
%symdel _ROC_XAXISOPTS_LABEL _ROC_YAXISOPTS_LABEL;
proc print data=roc.boy_1;/*如果没有proc print语句，结果运行无错误但是不会显示出来*/
run;