[R] NLME augPred error

[Pedro J. Aphalo]

Could someone explain the meaming 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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._