Thank you for your reply. Here is my code. I'm comparing the High metabolic syndrome risk group and the Low risk group. I found estimates but want to know the probability of the high risk group. proc nlmixed data=myfolder.metab; parms b0=0 b1=0 b2=0 b3=0; xb = b0 + b1*time + b2*sex + b3*age; prob = exp(xb)/(1+exp(xb)); liklhd = (prob**metab3)*((1-prob)**(1-metab3)); ll = log(liklhd); model metab3 ~ general(ll); run;
... View more