data Age_data;
input disease age bw@@;
datalines;
0 50 65  0 39 61  0 21 70  0 61 67
0 30 55  0 35 63  0 25 72  0 41 66
0 43 52  0 36 54  0 37 76  0 25 61
0 41 53  0 62 55  0 28 70  0 33 68
1 52 45  1 49 61  1 47 42  1 62 31
1 55 67  1 70 61  1 75 55  1 77 52
1 81 58  1 64 53  1 62 41  1 39 57
1 61 51  1 61 55  1 57 49  1 79 47
;
run;

proc logistic data=Age_data;
MODEL disease (event='1') = age / OUTROC = rdwi00;
RUN;

PROC GPLOT DATA= rdwi00;
PLOT _SENSIT_*_1MSPEC_ ;
symbol interpol=join;
run;

data leng;
set rdwi00;
x=sqrt((1-_sensit_)**2+_1mspec_**2);
run;
proc sort data=leng; by x _pos_ _neg_; run;

data n;
set leng;
if _n_=1 then do; call symput ('pos',_falneg_); call symput ('neg',_neg_);
                                  call symput ('ses',_sensit_); call symput ('spe',_1mspec_); end;
run;

proc sort data=age_data; by age descending disease ; run;
data age;
set age_data;
retain neg pos;
if disease=0 then neg+1;
if disease=1 then pos+1;
if neg=&neg and pos=&pos+1;
cutoff=age;
sensit=&ses;
_1mspec=&spe;
keep cutoff sensit _1mspec;
run;
proc print; run;

本例中大于等于47为cut off值

十分感谢， 学习先。

请问这段程序是根据什么选取的最大值？
数据集leng中的X指什么？

这个程序其实有问题，当样本含量较大时便不能出结果

有没有更好的程序注解，或者稍微简单点的程序。求赐教！