Here is the simple code to get CUTPOINT based on Youden index.
proc logistic data=sashelp.heart(obs=1000);
class bp_status;
model status=bp_status weight height / outroc=roc;
output out=out p=p;
run;
/*%rocplot(v,inpred=out, inroc=roc, p=p, id=_opty_, optcrit=youden)*/
data x;
set roc;
x=_SENSIT_;
y=1-_1MSPEC_;
youden=x+y-1;
run;
proc sql;
select _PROB_ as cutpint label='cutpoint',youden from x having youden=max(youden);
quit;
... View more