[R] fitting a linear mixed effects model

[Peter Dalgaard BSA]

"E. Corda" <e.corda at oncfs.gouv.fr> writes:

> fm1fox <- lme(nbfox ~ zone*annee, data=fox, random ~ 1|transect/year)
> 
> I would like the TRANSECT factor to be nested in the ZONE factor. But
> if I estimate the model like this:
> 
> fm1fox <- lme(nbfox ~ zone*annee, data=fox, random ~
> 1|zone/transect/year)

(shouldn't that be transect/zone/year?)

> zone is now a random effect.
> 
> Do you have any idea on this? Thanks a lot for any help.

Does random ~ 1|(zone:transect)/year work? Otherwise, you might need
something like random ~ 1|interaction(zone,transect)/year.

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._