[R] problems of nlme with R?
Ramon Diaz-Uriarte
rdiazuri at students.wisc.edu
Tue Mar 7 09:50:01 CET 2000
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list