OPTIONS NODATE NOCENTER LS=76 PS=54;
DATA a;
 INFILE 'temp.asc' FIRSTOBS=2;
 INPUT sbj trt resp n d1 d2 d3;

PROC NLMIXED QPOINTS=75 GCONV=1E-30;
 PARM a=1.5 a1-a3=.2 lambda=1.3 sigma=.6;

 eta = (a+d1*a1+d2*a2+d3*a3)*trt-u;
 p = PROBNORM(eta);
 MODEL resp ~ BINOMIAL(n,p);
 RANDOM u ~ NORMAL(lambda,sigma*sigma) SUBJECT=sbj;
 ESTIMATE 'd prime for NH ' a+a1;
 ESTIMATE 'd prime for Neg' a+a2;
 ESTIMATE 'd prime for Pos' a+a3;
 ESTIMATE 'd prime for Neu' a;

RUN;

想转nlme 不知如何作 或有文件 / 網頁可幫忙 感謝