你的问题比较大，不过我按照你的已经修改了，

  start fun(FPR);
 use eval;
 read all var{alpha beta} into x;
TPR=1-1/(1+exp((log(FPR/(1-FPR))+x[1])*exp(x[2]/2)));
 return(TPR);
 finish;

多谢贤兄仗义相助！

有IML高手指点我另外一个办法，如果不用read all var{alpha, beta} into x，那么就在start fun(FPR) 后面加上global(alpha, beta) ，实为异曲同工：

/* Define the Integrand in AUC Quaduature */
%MACRO TPR(FPR, alpha, beta);
 1-1/(1+exp((log(&FPR./(1-&FPR.))+&alpha.)*exp(&beta./2)));
%MEND TPR;
DATA eval;
 input alpha beta;
 cards;
 2.6728  -0.1220
;
RUN;
PROC IML;
 use eval;
 read all;
 print alpha beta;
 start fun(FPR) global(alpha, beta);
 TPR=%TPR(FPR, alpha, beta);
 return(TPR);
 finish;
 /* Call QUAD */
 bound = { 0  1 };
 call quad(AUC, "fun", bound) eps=1E-10;
 print AUC;
QUIT;