[R-sig-ME] Specifying outcome variable in binomial glmm: single responses vs cbind?

a y beermewi at gmail.com
Sat Jul 2 20:48:06 CEST 2016


I answered my own question, so feel free to disregard this topic.

On Fri, Jul 1, 2016 at 6:37 PM, a y <beermewi at gmail.com> wrote:

> What is the difference between fitting a binomial glmm (without random
> item effects) in the following two ways?
>
> 1.
> Data formatted in the following way:
>
> (data_long)
> ID    correct    condition    itemID
> 1      1             A               i1
> 1      0             A               i2
> 1      1             A               i3
> 1      1             A               i4
> 2      0             B               i1
> 2      1             B               i2
> 2      1             B               i3
> 2      0             B               i4
>
> Fitting a model without item random effects:
>
> glmer(correct ~ condition + (1|ID), family = binomial, data = data_long)
>
>
> 2.
> Data formatted this way (summing over the correct responses):
>
> (data_short)
> ID     sum_correct    condition     itemID
> 1       3                      A                NA
> 2       2                      B                NA
>
> Fitting the following model, assuming there were only 4 items  (I've seen
> dozens of examples like this):
> glmer(cbind(sum_correct, 4 - sum_correct) ~ condition + (1|ID), family =
> binomial, data = data_short)
>
> ---
> I figured these models should be identical, but in my experience they are
> very much not. What am I missing? When is the second (more) appropriate?
>
> Thanks for any help,
> Andrew
>

	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list