Date:         Thu, 3 Jun 2010 10:33:18 -0700
Reply-To:     oslo <oslo@yahoo.com>
Sender:       "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU>
From:         oslo <hokut1@YAHOO.COM>
Subject:      Re: NLMIXED vs GLIMMIX
Comments: To: oloolo <dynamicpanel@yahoo.com>
In-Reply-To:  <201006021536.o52DQsRr004768@willow.cc.uga.edu>
Content-Type: text/plain; charset=iso-8859-1



Dear Mary; 
 
I am so sorry for confussing. I thik  messed up because of the copy and past. I believe Dale and other also confussed because of the same problem. I am resending the codes. I sincerly apologize to everbody.
 
Regards,
 
R
 

proc nlmixed data=zip;
parms bp_0=0
bp_1=0 
bp_2=0 
bp_3=0 
bp_4=0 
null_1=0 
null_2=0
k=1 
log_s2u1=0
log_s2u2=-1
rho1=0
b00_0=0
bll_1=0
b22_2=0
b33_3=0
b44_4=0
null_3=0
null_4=0 
log_s2u3=0
log_s2u4=-1
rho2=0; 
eta_prob = bp_0 + bp_1*mar  + bp_2*income + bp_3*gender+ bp_4*ethnic +
null_1*drug + null_2*age + u1 + u2*z2;
p_0 = exp(eta_prob)/(1 + exp(eta_prob));
eta_lambda = b00_0 +b11_1*mar + b22_2*income + b33_3*gender+ b44_4*ethnic +
null_3*drug + null_4*age  + u3 + u4*z4;
lambda = exp(eta_lambda);
if RUN =0 then loglike = log(p_0 + (1-p_0)*exp(-lambda));
else loglike = log(1-p_0) + RUN *log(lambda) -lambda - lgamma(RUN+1);
  /* model fit*/
model RUN ~ general(loglike);
predict 1-G_t out=cdf;
random u1 u2 u3 u4 ~ normal([0,0,0,0], [exp(2*log_s2u1),
                                  rho1*exp(log_s2u1 + log_s2u2),
                                    exp(2*log_s2u2),0, 0, exp(2*log_s2u3),
                             0, 0, rho2*exp(log_s2u3 + log_s2u4),exp(2*log_s2u4)])
                 subject=subject;
  estimate "s2u1" exp(2*log_s2u1);
  estimate "s2u2" exp(2*log_s2u2);
  estimate "cov" rho1*exp(log_s2u1 + log_s2u2);
  estimate "s2u3" exp(2*log_s2u3);
  estimate "s2u4" exp(2*log_s2u4);
  estimate "cov" rho2*exp(log_s2u3 + log_s2u4);
  estimate "alpha" 1/k;
run;