From: Lijuan Wang <lw7s_at_cms.mail.virginia.edu>

Date: Wed 28 Jul 2004 - 11:30:38 EST

tau[i]=tau+e[i3]

R-help@stat.math.ethz.ch mailing list

https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jul 28 11:40:03 2004

Date: Wed 28 Jul 2004 - 11:30:38 EST

Hi,

I am using Splus to run a multiphase mixed-effects model.
The quations of the models are as below:

gf[ij]=b0[i]+b1[i]*age[ij]+b2[i]*max(0,(age[ij]-tau[i]))^2+e[ij] b0[i]=b00+e[i0] b1[i]=b10+e[i1] b2[i]=b20+e[i2]

tau[i]=tau+e[i3]

i: 1,2,...,100 subjects

j: 1,2,...,6 occasions

The main scripts of Splus is:

simu1<-groupedData(gf~age|id)

simu.nlme<-nlme(gf~(b0 + b1 * age + b2 *
max(0,(age-tau))^2),data=simu1,fixed=list(b0~1,b1~1,b2~1,tau~1),
+

random=list(b0~1,b1~1,b2~1,tau~1),start=c(b0=4,b1=5.32,b2=-5.29,tau=14.8))

But it seems like that there are some errors in it. Error in .C("mixed_EM",: Singularity in backsolve at level 2, block 1

I think maybe there is something wrong in my scripts,
could you help me to find it?

Thanks a lot

Peggy

R-help@stat.math.ethz.ch mailing list

https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jul 28 11:40:03 2004

*
This archive was generated by hypermail 2.1.8
: Fri 18 Mar 2005 - 02:40:22 EST
*