[R] Need confirmation of my lme() structure

[lescroel_cebc@no-log.org]

Hello,

I'm a new user of R and I would need somebody to confirm (or not!) that
I'm doing right. I'm analyzing diving characteristics (maximal depth,
etc...) of birds breeding in 2 sites, during 2 stages (see data
structiring below).

   FILE        SITE      STAGE
M04    :402   C: 471   Guard:1557
M13    :344   M:2143   Incub:1057                                     M10 
  :252                                  M14    :209                       
          M15    :190                                  (Other):955

I want to analyse the effects of site and stage using REML analyses with
site and stage as fixed factors and file (i.e. individual bird) as a
random factor.
In R, I would write the following command:

> reml_depthmax <- lme(DEPTHMAX~SITE*STAGE, random=~1|FILE,
na.action=na.omit)

Am I right???

Many thanks for your help,

Amélie


-- 
Amélie Lescroël, Ph.D.
Centre d'Etudes Biologiques de Chizé - CNRS, UPR 1934
79360 Villiers en Bois
France
http://www.cebc.cnrs.fr/Uk_taaf/equip_AL.html