[R] NLME augPred error

[Pedro J. Aphalo]

This is a follow up to the quoted message that I sent to the r-help list
some time ago about an error message I could not understand.

I have now tried the same example under S-Plus 3.4 and NLME 3.3 and it
works without error. So, it looks like this is a bug in the R port of
nlme. I am sorry that I do not know well enough the R language to
provide a patch.

By the way, I tested nlme 3.0-17 for R (MS-Windows) and the problem
persists.

Thanks to all involved in porting NLME to R. It is extremely useful in
my work, as I have only limited access to S-Plus.

Pedro.

"Pedro J. Aphalo" wrote:
> 
> Could someone explain the meaning of this error message from augPred:
> 
> > augPred(area3.pen.nlme, primary=~day)
> Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop =
> FALSE],  :
>         Levels 1,2,3 not allowed for block
> >
> 
> predict.nlme(area3.pen.nlme) does not produce an error.
> 
> area3.pen.nlme was created with:
> 
> > area3.pen.nlme <- nlme(area ~ SSlogis(day, Asym, xmid, scal),
> +  data=area.pen.data,
> +  fixed = Asym + xmid + scal ~ dose + block,
> +  random = Asym + xmid + scal ~ 1|seedling,
> +  weights = varPower(),
> +  start = c(2887,0,0,0,10.7,0,0,0,2.29,0,0,0))
> >
> 
> area, day and dose are numeric
> block and seedling are factors
> 
> area is the area in mm2 of one growing leaf for each seedling.
> 
> with
>   fixed = Asym + xmid + scal ~ dose
> or
>   fixed = Asym + xmid + scal ~ block
> it works.
> 
> Thanks in advance for any help.
> 
> Pedro.
> 
> Windows 95
> 
> platform i386-pc-mingw32
> arch     x86
> os       Win32
> system   x86, Win32
> status
> major    1
> minor    3.0
> year     2001
> month    06
> day      22
> language R
> 
> Package: nlme
> Version: 3.1-15
> Date: 2001/06/18
> Built: R 1.3.0;  Win32;  2001-06-20 13:32:23
> 
> and
> 
> Package: nlme
> Version: 3.1-16
> Date: 2001/06/20
> Built: R 1.3.0;  Win32;  2001-07-24 12:02:25
> 

-- 
===============================================================
 Pedro J. Aphalo
 University of Joensuu
 Faculty of Forestry
 P.O. Box 111, FIN-80101 Joensuu, FINLAND
 mailto:pedro.aphalo at joensuu.fi
 http://cc.joensuu.fi/~aphalo/pedro.html
 fax: +358 13 251 4444
 tel: +358 13 251 4406                         ,,,^..^,,,
================================================================
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._