[R-sig-ME] The simple burning questions about lmer and their answers
Page E. Van Meter
vanmete7 at msu.edu
Thu Oct 30 00:11:06 CET 2008
Hello all,
I have been teaching myself how to use generalized linear mixed effect
models in R. While I find that the message boards like this one are
extremely helpful, they have a "time travel" issue. I am finding answers
to my burning questions that are years old and I am having trouble
determining if the issues raised in 2006 are still issues today. So I
wanted to put all of the hot topics I have been reading about, and the
overly simplified answers I have gleaned in one thread. I hope others
may add their burning questions and the answers they use to this thread
as well?? And always the hope for comments:
Q1: Should I use ML or REML to compare the fixed effects in a set of
nested models using likelihood-ratio testing (LRT)?
A1: It is much better to use ML to compare fixed effects in your model.
Q2: Should I use ML or REML to compare the random effects in a set of
nested models using likelihood-ratio testing (LRT)?
A2: To test the worth of random effects in your model, REML is better.
HOWEVER, if using lmer with a family specification to call glmer (family
= binomial or poisson) you are automatically using ML. REML is not an
option for generalized linear mixed effect models using lmer.
Q3: Ok, so how do I asses the value of random effects in my generalized
linear mixed effect model?
A3: <insert cricket sounds here>
I gather that your experimental design should be the major indicator
that a random term is needed in your model specification. But sometimes
showing that the random effect is valuable (perhaps in a case where none
of your fixed effects are worth anything) is, in itself, an interesting
result. Help?
Q4: So what is the hoopla about p-values for fixed effects in lmer?
A4: Short answer is that p-values are now included in the current
incarnation of lme4 (can't find documentation of this but I see them so
they must be there right?). The long answer about why you should be
hesitant to use them and more hesitant to NEED them was provided by Doug
here :
http://www.nabble.com/lmer%2C-p-values-and-all-that-p4476689.html
************************************
Page E. Van Meter
Michigan State University
Department of Zoology
vanmete7 at msu.edu
**http://msu.edu/~vanmete7/*
More information about the R-sig-mixed-models
mailing list