Dear R-experts, Quite new to R on this end, but learning fast (I hope). I am running version 2.7.1 on Windows Vista. I have small dataset which consists of:
# NestID: nest indicator for each chicken. Siblings sharing the same nest have the same nest indicator. # Chick: chick indicator consisting of a unique ID for each single chick. # Year: 1, 2. # ClutchSize: 1-, 2- , 3-eggs. # HO: hatching order within each clutch (1, 2, 3 [first, second and third-hatched chick]). # SibComp: sibling competence: present/ absent (0, 1) # Death2: death at two days post-hatch (0, 1) # Death10: death at ten days post-hatch (0, 1) So a subset of my dataset looks something like this: NestID Chick Year ClutchSize HO Hatching SibComp Death2 Death10 1 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 0 0 3 3 1 1 1 0 0 0 0 4 4 1 1 1 1 0 1 0 4 5 1 2 2 0 1 0 1 5 6 1 2 1 1 0 0 0 5 7 1 2 2 0 0 0 0 6 8 2 3 1 1 1 0 0 6 9 2 3 2 1 0 1 0 6 10 2 3 3 0 1 0 0 7 11 2 3 1 0 0 0 1 7 11 2 3 2 0 0 0 0 7 11 2 3 3 1 1 1 1 ............ In order to account for lack of independence at the nest level (many chicks are siblings), I'd like to run a GLMM with random slopes and intercepts for nests. Using lmer, my model for survival at 10 days, for example, would read as follows (or not!): > model <- lmer(Death10 ~ HO + ClutchSize + SibComp + Year + (1|NestID), > family=binomial, 1) > summary(model) >From what I understand, the model above includes only random intercepts for >NestID. So at this point my question is how do I make this model into one >which includes both random intercepts and slopes for NestID? Look forward to receiving your input. Thank you all for your time! Luciano ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
