Bob, You need to make use of pdBlocked and pdIdent. There is a simple example in Bates and Pinheiro's book (Mixed Effects Models in S an S+). I don't have the book in front of me to get the page number, but I remember it is an example using data contributed by Dave Lansky. A related s-news thread (which is probably sufficient to answer your question) is: http://www.biostat.wustl.edu/mailinglists/s-news/200112/msg00126.html Best, Jim James A. Rogers <rogers at cantatapharm.com> Statistical Scientist Cantata Pharmaceuticals
Date: Fri, 13 Sep 2002 15:17:24 +0300 From: "Anon." <bob.ohara at helsinki.fi> Subject: [R] Multiple random effects inlme? Moi! I was helping to teach a course on mixed models this week, and we came across a problem with coding more than one random effect in lme when they aren't nested.
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._