[R-sig-ME] Specifying outcome variable in binomial glmm: single responses vs cbind?
a y
beermewi at gmail.com
Sat Jul 2 01:37:53 CEST 2016
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