Displaying 1 result from an estimated 1 matches for "probnorm".
Did you mean:
problerm
2008 Apr 07
0
Translating NLMIXED in nlme
...mixed 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 haven't found a simple solution for the conditional LL:
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));
using nlme.
Any suggestion will be appreciated.
Giovann...