[R] nested ANCOVA: still confused

From: Jeffrey Stratford <stratja_at_auburn.edu>
Date: Wed 25 Jan 2006 - 03:33:49 EST

Dear R-users,

I did some more research and I'm still not sure how to set up an ANCOVA with nestedness. Specifically I'm not sure how to express chicks nested within boxes. I will be getting Pinheiro & Bates (Mixed Effects Models in S and S-Plus) but it will not arrive for another two weeks from our interlibrary loan.

The goal is to determine if there are urbanization (purban) effects on chick health (rtot) and if there are differences between sexes (sex) and the effect of being in the same clutch (box).

The model is rtot = sex + purban + (chick)box.

I've loaded the package lme4. And the code I have so far is

bb <- read.csv("C:\\eabl\\eabl_feather04.csv", header=TRUE) bb$sex <- factor(bb$sex)
rtot.lme <- lme(bb$rtot~bb$sex, bb$purban|bb$chick/bb$box, na.action=na.omit)

but this is not working.

Any suggestions would be greatly appreciated.



Jeffrey A. Stratford, Ph.D.
Postdoctoral Associate
331 Funchess Hall
Department of Biological Sciences
Auburn University
Auburn, AL 36849
FAX 334-844-9234

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 Wed Jan 25 03:44:31 2006

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:42:10 EST