Hi there I wanted to apply proc nlmixed on real data but I got this error: Execution error for observation 97. Is there any body here to help me how to solve it? Best wish proc nlmixed data=data qpoints=10 gconv=1e-2; bounds 0 <=sigma2; /* Initail Values ---------------- */ parms a0=5.3 a1=1.3 b0=14.3 b1=0.01 b2=0.2 b3=2.5 g1=1.1 g2=-0.1 g3=0.5 sigma2=1 l1=1 l2=2 l3=3 k=0.01; /* Log Likelihood ---------------- */ if t=0 then do; lam = exp(g1*GENDER + g2*AGE +g3*CITY + l2*a + l3*b); loglik=(censored=0)*(log(lam)-lam * time) + (censored=1)*(-lam * time); end; else do; pi = 1/(1 + exp(-(a0 + a1*ADTYPE + a))); eta = abs(k)**(-2); mu = b0 + b1*ADTYPE +b2*AGE +b3*GENDER - log(pi)-sigma2*log(k**2)/k-lgamma(eta*(k+k*sigma2))+lgamma(eta) + l1*a + b; u = sign(k)*(log(Y) - mu)/sigma2; if Y=0 then loglik=log(1-pi); else if Y>0 then loglik=log(pi)- lgamma(eta) + eta*log(eta) - log(sigma2) - log(Y) + u*sqrt(eta) - eta*exp(abs(k)*u); else if Y = . then loglik=0; end; model ID ~ general(loglik); random a b ~ normal([0,0], [1, 0, 1]) subject=ID; run;
... View more
