楼主: 付齐齐
284 0

[问答] 求助,从多个条件中至少满足三个如何用sas编程,下面这样可以吗?但运行有错误 [推广有奖]

  • 0关注
  • 0粉丝

等待验证会员

学前班

40%

还不是VIP/贵宾

-

威望
0
论坛币
10 个
学术水平
0 点
热心指数
1 点
信用等级
0 点
经验
20 点
帖子
1
精华
0
在线时间
0 小时
注册时间
2018-7-11
最后登录
2018-7-11

付齐齐 发表于 2018-7-11 20:13:50 |显示全部楼层
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;


已有 1 人评分热心指数 收起 理由
eijuhz + 1 精彩帖子

总评分: 热心指数 + 1   查看全部评分

您需要登录后才可以回帖 登录 | 我要注册

GMT+8, 2018-11-20 04:20