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;