SAS logistic里能否得出ROC曲线的CUT-OFF界值

webgu 发表于 2013-4-16 16:53
ROC曲线上会有一个最佳的cut-off值。他的约登指数会最大。这个才是我们想要的。

以前的一段代码
data YoudenIndex;
set rocstats;
_SPECIF_ = 1 - _1MSPEC_;
YoudenIndex = _SENSIT_ + _SPECIF_ - 1;
run;

hamsik11 发表于 2013-4-16 17:32
以前的一段代码
data YoudenIndex;
set rocstats;

看来只能加程序算了。

有时候，总想找到一个小的OPTION解决问题。但是找来找去，就是没有。

最后不得不自己去编。

SAS小学生 发表于 2013-4-16 03:18
proc logistic dara=urdata plots=roc (id=prob);
试试这个

上面的plots=roc(prob)是不是只能9.2以上的版本用，9.1.3好像运行不了？

nilingyong0420 发表于 2013-5-30 19:53
上面的plots=roc(prob)是不是只能9.2以上的版本用，9.1.3好像运行不了？

9.3可行

每组roc坐标都对应prob,你应该找到这个prob即可

爱萌 发表于 2013-6-5 16:49
每组roc坐标都对应prob,你应该找到这个prob即可

嗯。可以编程找到。谢谢。

..............发错了 不好意思

webgu 发表于 2013-4-16 18:54
看来只能加程序算了。

有时候，总想找到一个小的OPTION解决问题。但是找来找去，就是没有。

9.1.3编程后出现了灵敏度和特异度的空值。
能否分享您的程序？
如下数据：
CASE          DURATION          RATIO

1	6.10	4.51
1	21.77	4.58
1	9.67	4.58
0	45.13	4.64
1	2.43	4.67
0	46.07	4.73
0	37.87	4.77
0	52.87	4.88
1	14.00	5.35
1	34.37	5.85
0	61.90	6.00
1	12.17	6.33
0	45.47	6.46
0	64.13	6.52
1	32.97	7.64
0	76.87	7.95
1	4.10	8.17
1	18.70	8.19
1	10.87	8.85
1	47.80	9.02
1	19.20	9.38
1	10.73	9.70
1	17.13	9.86
0	42.83	10.13
1	10.43	11.07
1	18.23	11.52
1	15.17	11.91
1	13.30	12.50
0	39.60	18.31
1	33.13	18.51
0	46.00	83.50

也想要这个程序

楼主解决了吗？我现在也碰到这个问题了，如何找最佳cut-off值？我知道是roc曲线最左上角那一点，或者是约登指数最大的，但是SAS好像是没有给出来的。难道要自己算吗