Displaying 1 result from an estimated 1 matches for "female_chick_no".
2010 Nov 19
2
Question on overdispersion
...) and some females
have multiple nests in the data set (so I need to include female identity as
a random effect).
Here is an example of what the three vectors used in the model look like
(the real data set is much bigger, just to illustrate what I?m talking
about):
male_chick_no=c(2,4,1,0,3,5,2)
female_chick_no=c(1,0,3,3,1,0,2)
FemaleID=c(A,A,B,B,C,D,E)
My first question relates to coding the test in R. I received this suggested
R syntax from a reviewer:
SexRatio = cbind(male_chick_no, female_chick_no)
Model <- lmer(SexRatio ~ 1 +(1|FemaleID), family = quasibinomial)
But when I try to use this in R...