Displaying 1 result from an estimated 1 matches for "sigsq2".
Did you mean:
sigsq
2008 Apr 07
0
Translating NLMIXED in nlme
Dear All,
reading an article by Rodolphe Thiebaut and Helene Jacqmin-Gadda ("Mixed
models for longitudinal
left-censored repeated measures") I have found this program in SAS
proc nlmixed data=TEST QTOL=1E-6;
parms sigsq1=0.44 ro=0.09 sigsq2=0.07 sigsqe=0.18 alpha=3.08 beta=0.43;
bounds $B!](B1< ro < 1, sigsq1 sigsq2 sigsqe >= 0;
pi=2*arsin(1);
mu=alpha+beta*TIME+a i+b i*TIME;
if OBS=1 then ll=(1/(sqrt(2*pi*sigsqe)))*exp(-(RESPONSE-mu)**2/(2*sigsqe));
if OBS=0 then ll=probnorm((RESPONSE-mu)/sqrt(sigsqe));
L=log(ll);
.......
I...