Date:         Tue, 12 Feb 2008 16:25:59 -0500
Reply-To:     IKF <ikf@HOTMAIL.COM>
Sender:       "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU>
From:         IKF <ikf@HOTMAIL.COM>
Subject:      More NLMIXED

Hello,

I think that maybe what I want will be a random intercept and a random
slope in my nlmixed.  I found a post which suggests this, but I can't get
it to work - my understanding of matrices is nill!

This is what I have...I figure something in here needs to be specified
more clearly...Please help...Babyagedays is my time variable, so I'm not
even sure if I should be doing this...Thanks...


proc nlmixed data=pneumo1;
 parms b0=-2.5 b1=0 b2=0 b3=0 s0=0 s1=0;
 teta = b0 + b1*babyagedays + b2*babyagedays1 + b3*babyagedays2 + g0 +
        g1*babyagedays;
 p = exp(teta) / (1 + exp(teta));
 model culture ~ binary(p);
random g0 g1 ~ normal([0,0],
                          [exp(2*log_s0),
                           rho*exp(log_s0+log_s1),
                           exp(2*log_s1)]) subject=pin;
run;