Dear all,
I have counts of males and females produced by different colonies of a species of social insects.
I fitted the model
Fit=glmer(cbind(females,names)~1+(1|colony),family=binomial,data=data)

However, how can I test within such a framework if the overall average sex ratio deviates from an a priori predicted value (e.g. half females, or ž females)?
I presume this would have to be done based on the fitted intercept. But how does one do this? Also, what would be the best way to get 95% conf lims on the estimate? Using likelihood profiling, or parametric bootstrap? Does anybody happen to have any example calculation?

Cheers,
Tom


_______________________________________________________________________________________

Prof. Tom Wenseleers
*      Lab. of Socioecology and Social Evolution
           Dept. of Biology
           Zoological Institute
           K.U.Leuven
           Naamsestraat 59, box 2466
           B-3000 Leuven
           Belgium
* +32 (0)16 32 39 64 / +32 (0)472 40 45 96
* tom.wenseleers@bio.kuleuven.be
http://bio.kuleuven.be/ento/wenseleers/twenseleers.htm



	[[alternative HTML version deleted]]

