Hallo thank you for your response. I am not sure but maybe fixed effects cannot be set to be influenced by a factor to be able to use augPred.
lob<-Loblolly[Loblolly$Seed!=321,] set.seed(1) lob<-data.frame(lob, x1=sample(letters[1:3], replace=T)) # add a #factor lob<-groupedData(height~age|Seed, data=lob) fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc), data = lob, fixed = Asym + R0 + lrc ~ 1, random = Asym ~ 1, start = c(Asym = 103, R0 = -8.5, lrc = -3.3)) fm2<-update(fm1, fixed=list(Asym~x1, R0+lrc~1), start=c(103,0,-8.5,- 3)) ^^^^^^^ and plot(augPred(fm2)) Throws an error. So it is not possible to use augPred with such constructions. Best regards. Petr Pikal On 2 Sep 2006 at 17:58, Spencer Graves wrote: Date sent: Sat, 02 Sep 2006 17:58:05 -0700 From: Spencer Graves <[EMAIL PROTECTED]> To: Petr Pikal <[EMAIL PROTECTED]> Copies to: r-help@stat.math.ethz.ch Subject: Re: [R] my error with augPred > <comments in line> > > Petr Pikal wrote: > > Dear all > > > > I try to refine my nlme models and with partial success. The model > > is refined and fitted (using Pinheiro/Bates book as a tutorial) but > > when I try to plot > > > > plot(augPred(fit4)) > > > > I obtain > > Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop = > > FALSE], : > > Levels (0,3.5],(3.5,5],(5,7],(7,Inf] not allowed for > > vykon.fac > > > > > > Is it due to the fact that I have unbalanced design with not all > > levels of vykon.fac present in all levels of other explanatory > > factor variable? > > > I don't know, but I'm skeptical. > > I try to repeat 8.19 fig which is OK until I try: > > > > fit4 <- update(fit2, fixed = list(A+B~1,xmid~vykon.fac, scal~1), > > start = c(57, 100, 700, rep(0,3), 13)) > > > > I know I should provide an example but maybe somebody will be clever > > enough to point me to an explanation without it. > > > I'm not. > > To answer these questions without an example from you, I'd have to > make up my own example and try to see if I could replicate the error > messages you report, and I'm not sufficiently concerned about this > right now to do that. > > Have you tried taking an example from the book and deleting certain > rows from the data to see if you can force it to reproduce your error? > > > Alternatively, have you tried using 'debug' to trace through the code > line by line until you learn enough of what it's doing to answer your > question? > > Spencer Graves > > nlme version 3.1-75 > > SSfpl model > > R 2.4.0dev (but is the same in 2.3.1), W2000. > > > > Thank you > > Best regards. > > > > Petr PikalPetr Pikal > > [EMAIL PROTECTED] > > > > ______________________________________________ > > R-help@stat.math.ethz.ch mailing list > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide > > http://www.R-project.org/posting-guide.html and provide commented, > > minimal, self-contained, reproducible code. > > > > ______________________________________________ > R-help@stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html and provide commented, > minimal, self-contained, reproducible code. Petr Pikal [EMAIL PROTECTED] ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.