Hello,
Please I would love to run this logistic procedure with proc nlin or proc nlimixed, any ideas please? the codes I tried are below. Thanks
proc logistic data = rsm1;
model edema (event="1") = DoseA DoseT DoseA*DoseT DoseA*DoseA DoseT*DoseT/link = probit;
output out=pred p=prob;
run;
**This did not work
proc nlmixed data = rsm1; call streaminit(123456);
parms beta0=1 beta1=1 beta2=1 beta11=1 beta12=1 beta22=1 s2=1;
y2 = (beta0 + beta1*DoseA + beta2*DoseT + beta11*DoseA*DoseA
+ beta22*DoseT*DoseT + beta12*DoseA*DoseT);
p = exp(y2)/(1+exp(y2));
model Edema~binary(p);
predict y2 out=y2;
run;
**I'm not getting same results with this:
proc nlin data =rsm1 method=newton;
parms beta0=1 beta1=1 beta2=1 beta11=1 beta12=1 beta22=1;
y2 = (beta0 + beta1*DoseA + beta2*DoseT + beta11*DoseA*DoseA
+ beta22*DoseT*DoseT + beta12*DoseA*DoseT);
p = probnorm(edema);
if (edema = 1) then p = 1-p;
model p=y2;
output out=p p=predict;
run;
... View more
