SAS Support Communities

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

Hello,   Please I am trying to run these codes but I'm getting the above errors and do not know what to do. I will appreciate if anyone can help me. Thanks   There are 8000 observations. Edema (Yes=1/No=0)   proc nlmixed data = rsm1; call streaminit(12345); 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; ... View more