Displaying 1 result from an estimated 1 matches for "llbin".
Did you mean:
albin
2006 Jul 25
1
HELP with NLME
...appreciated.I apologize for the
clumsiness of the R code.
Many thanks in advance.
Sincerely,
Loki
#################################################################
SAS Code:
proc nlmixed data=repdat
parms b0 -3 b1 -.135 a0 3 a1 4 sigsq 0.25;
if vs.flag = 1 then do;
eta1 = b0 + b1*lnbldT;
llbin = anybc.cens.ind*eta1 - log(1+exp(eta1));
eta2 = a0 + a1*lnndsTs;
llnorm = -1/(2*sigsq)*(lnbldT - eta2)**2 - .5*log(sigsq);
ll = llbin + llnorm;
end;
else do;
eta2 = a0 + a1*lnndsTs;
eta1 = b0 + b1*eta2 + u;
llbin = anybc.cens.ind*eta1 - log(1+exp(eta1));
ll = llbin...