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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._