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: [email protected]
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]
> >
> > ______________________________________________
> > [email protected] 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.
> >
>
> ______________________________________________
> [email protected] 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]
______________________________________________
[email protected] 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.
