ROC曲线在SAS中如何绘制

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

现有如下两列数： x	y
0.940278	1
0.838616	0
0.84515	1
0.878174	1
0.88333	1
0.885006	1

想利用SAS 9.0画出ROC曲线，并求出AUC值，该怎么做啊！谢谢了！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏6 回帖

关键词：ROC曲线 ROC AUC 怎么做 如何

data a;
input x iswy;
cards;
0.940278 1
0.838616 0
0.84515 1
0.878174 1
0.88333 1
0.885006 1
;
run;
Proc sort data=a out=a;  
by descending x;
run;
data  a;
set  a;
i=1;
sn+i;
run;
data a;
set a;
wy+iswy;
if iswy=0 then isnotwy=1;
else isnotwy=0;
run;
data a ;set a;
notwy+isnotwy;
run;
proc sql;
create table roc as
select sn,wy/sum(iswy) as y,notwy/sum(isnotwy) as x
from a;
quit;
run;
proc sql;
create table auc as
select sum(y*delta) as auc  
from  (select a.*,b.x as x2,a.x- b.x as delta
       from roc a,roc b
       where a.sn=b.sn+1)
;
quit;
run;
/***ROC曲线**/
axis order=(0 to 1 by .1) label=none length=4in;
symbol i=join v=none c=depk;
symbol2 i=join v=none c=black;
title "roc curve";
axis1 label=(angle=90 "sensitivity");
axis2 label=("1-1mspec");
proc gplot data =roc gout=roc;
plot y*x x*x
/ overlay vaxis=axis1 haxis=axis2;
run;
quit;

zhujp 发表于 2012-2-25 22:53
data a;
input x iswy;
cards;

非常感谢您的帮助！谢谢！

帮顶~

好好好

很牛x

很牛x

好好好 GOOD GOOD GOOD

用SPSS应该会更方便吧