[R] Need confirmation of my lme() structure

From: <lescroel_cebc_at_no-log.org>
Date: Thu 23 Feb 2006 - 21:36:47 EST


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

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Thu Feb 23 21:45:18 2006

This archive was generated by hypermail 2.1.8 : Fri 24 Feb 2006 - 00:08:45 EST