[R] problems of nlme with R?

[Ramon Diaz-Uriarte]

Dear All,

I am running into two problems with nlme in R (R versions both 0.90.1 and
1.0.0).

1. If I try to fit a model with crossed random effects, such as: 

>my.model.fit1<-lme(ad~trt,random=list(pdBlocked(list(pdIdent(~Week-1),pdIdent(~Encl-1),pdIdent(~id-1)))),data=exp2.dat)
Error in parse(file, n, text, prompt) : parse error

(I can fit that exact same model in SPlus without problems).


2. the anova.lme function does not seem to work with the "Terms" option;
e.g.
> anova(my.model, Terms=11:14)  
Error in anova.lme(my.model, Terms = 11:14) : Object "noZeroColL" not found

I did not get this error when using SPlus.

Is there an easy fix of these problems, or am I missing something obvious?

Ramon




-- 
Ramon Diaz-Uriarte
Dept. Zoology and Statistics
University of Wisconsin-Madison
Madison, WI 53706

email: rdiazuri at students.wisc.edu
phone: 608-238-8041
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._