| Date: | Thu, 1 Jul 2004 16:05:19 -0500 |
| Reply-To: | Jeffrey Stratford <stratja@AUBURN.EDU> |
| Sender: | "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU> |
| From: | Jeffrey Stratford <stratja@AUBURN.EDU> |
| Subject: | Getting output from NLMIXED |
| Content-Type: | text/plain; charset=US-ASCII |
I'm finally figuring out what I'm doing with the ZIP models and the last
thing I need to do is get a table with predicted values for the
dependent variable. NLMIXED doesn't appear to like an OUTPUT statement.
Could someone tell me how to get a table with both the original values
and predicted values?
Thanks,
Jeff
Below is the SAS code for the NLMIXED procedure
proc nlmixed data=inbu;
parameters a0=0 to 1 by .2 a1=0 to 1 by .2 b0=0 to 1 by .2 b1=0 to 1
by .2;
linpinfl = a0 + a1*TK;
infprob = 1/(1+exp(-linpinfl));
lambda = exp(b0 + b1*TK);
if round=0 then
ll = log(infprob + (1-infprob)*exp(-lambda));
else ll = log((1-infprob)) + density*log(lambda) - lgamma(density+1) -
lambda;
model density ~ general(ll);
output out=temp;
run;
****************************************
Jeff Stratford
331 Funchess Hall
Department of Biological Sciences
Auburn University
Auburn, AL 36849
FAX 334-844-9234
http://www.auburn.edu/~stratja
****************************************
|