Re: [R] my error with augPred

From: Spencer Graves <spencer.graves_at_pdf.com>
Date: Sun 03 Sep 2006 - 00:58:05 GMT

<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
> petr.pikal@precheza.cz
>
> ______________________________________________
> 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. Received on Sun Sep 03 11:11:50 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 07 Sep 2006 - 07:51:17 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.