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 不知如何作 或有文件 / 網頁可幫忙 感謝