binomial GLMM with small upper limit
On Wed, 20 Jun 2012, Wojdak, Jeremy wrote:
I am working with mortality data - the number of three animals per experimental unit that had died by week 1, week2, ... during an experiment. So, I am using a binomial GLMM approach to model proportions with a random effect to deal with the repeated measurements from the same exp. unit. (Specifically, I use the experimental unit or "Tank" as the random effect, since multiple observations from the same unit must be related... and there were temporal/spatial blocks, so each tank is nested within a block. I include sample "Day" as a fixed effect in the model) sm7<-glmer(predc.survtbl~predator*Day +(1|Block/Tank), family="binomial", data=predc2)
Does your riskset size change by Day in predc.survtbl?
All is well, except the model suggests there are no treatment or time effects, while graphical inspection suggests there may be both.
So a naive survival analysis (ignoring Tank etc) is significant? You might try the coxme package...
| David Duffy (MBBS PhD) ,-_|\ | email: davidD at qimr.edu.au ph: INT+61+7+3362-0217 fax: -0101 / * | Epidemiology Unit, Queensland Institute of Medical Research \_,-._/ | 300 Herston Rd, Brisbane, Queensland 4029, Australia GPG 4D0B994A v