[R-sig-ME] P-values from interaction terms using lme4

Douglas Bates bates at stat.wisc.edu
Wed Aug 9 16:59:32 CEST 2017


On Wed, Aug 9, 2017 at 9:54 AM Phillip Alday <phillip.alday at mpi.nl> wrote:

> On 08/09/2017 04:46 PM, Douglas Bates wrote:
>


> > Technically it is true that the number of parameters does not depend on
> > the number of random effects, only on the number of unique values in the
> > covariance matrices for the random effects.  However, I think that leads
> > to an undercount of the effective number of parameters when, say,
> > performing a likelihood ratio test of models that differ in their random
> > effects specification.
>
> This is related to the more general issue of degrees of freedom in mixed
> models, right?
>

Yes.

But in terms of estimation (especially as implemented in
> lme4/MixedModels.jl), increasing the levels of a random effect will
> generally provide better estimates, right?
>

Yes.  Or, to put it the other way, it is unrealistic to expect to obtain
precise estimates of variance components unless there is(are) a large
number of levels in the grouping factor(s) for the random effects.

	[[alternative HTML version deleted]]



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