[R-sig-ME] LMM: including ranef or not?
Anna-Marie Corman
annajess at gmx.de
Tue Jul 22 08:55:44 CEST 2014
Dear Ben and Christian,
thanks again for your help! I know the number of fixes per bird, though
my dataset is on a trip basis. All individuals made a different number
of trips during the study period.
I will try this out.
Best,
Anna
Am 21.07.2014 23:24, schrieb Ben Bolker:
> If you're going to use ud95 as the response you might as well average
> the prop_land values per bird (i.e., aggregate the data down to a single
> data record per bird); the within-bird variation in prop_land values
> won't affect the model output at all (although the _number_ of
> observations per bird will; if you have unbalanced information in this
> way, you should incorporate weights proportional to the number of
> observations as well).
>
> With only 6 colonies you're going to have some difficulty estimating
> the among-colony variance very well; if you end up with zero estimates
> of among-colony variance, you might want to use blme or set the colonies
> as fixed effects ...
>
> If you use prop_land as the response variable you would indeed want to
> put bird_id in as a random effect (and you might consider estimating the
> proportional as a binomial (GLMM) response, _if_ you know the total
> number of fixes for each bird)
>
>
More information about the R-sig-mixed-models
mailing list