Skip to content
Prev 12100 / 20628 Next

A question on item random effects

On 14-07-09 06:03 AM, Anastasiia Romanova wrote:
A reproducible example might be nice, but I'll take a shot.

  In theory it should be possible to fit a random effect of item as long
as _any_ items are measured repeated times (they don't _all_ have to be
measured more than once).  In practice it may be difficult if only a
small fraction of items have been measured repeatedly.  It's hard to say
with precision what "small" means in the previous sentence, or how this
will come out in your case.  If you are in fact pushing your data too
hard, the symptoms will be convergence warnings and singular fits
(getME(.,"theta") including zero values, or VarCorr() including
variances of zero or correlations of +/- 1).  The singular fit is not
_necessarily_ wrong, but increases the chances of numerical errors and
complicates interpretation. One choice would be to use blme::bglmer() to
force your data away from the singular edge (but you would have to read
the corresponding paper to understand what you were actually doing).

  cheers
    Ben Bolker