[R] lmer estimated scale
Daniel Ezra Johnson
johnson4 at babel.ling.upenn.edu
Fri Mar 23 01:23:10 CET 2007
I have data consisting of several binary responses from a large
number of subjects on seven similar items. I have been using lmer
with (crossed) random effects for subject and item. These effects are
almost always (in the case of subject, are always) significant
additions to my model, testing this with anova. Including them also
increases the Somers' Dxy value substantially.
Even without those reasons, I feel I'd have to include these random
effects to account for the correlation between the seven items from
every subject. Otherwise my fixed between-subject effects like race,
gender, etc. will seem more significant than they should.
But how should I interpret the fact that without a Subject effect
included, the "estimated scale" parameter is usually very close to 1,
while when I include the Subject effect the scale parameter drops,
usually to around 0.85?
Can I at least conclude something interesting from this? Is it the
same as saying that the subject effect itself (meaning the 'observed'
subject BLUPs) is underdispersed with respect to its theoretical
normal distribution?
To summarize:
a <- lmer(Response~Fixed Effects+(1|Subject)+(1|Item),data,binomial)
b <- lmer(Response~Fixed Effects+(1|Item),data,binomial)
a has a much better fit by any measure, and estimated scale around 0.85.
b has a worse fit, and estimated scale around 1.
Obvious? Interesting? Worrisome?
Thanks,
Dan
More information about the R-help
mailing list