Skip to content
Prev 263072 / 398502 Next

Using MCMC sampling to estimate p values with a mixed model

Woodcock, Helena <H.E.Woodcock <at> liverpool.ac.uk> writes:
It's hard to know exactly without more details/seeing the data;
it does sound suspicious.
  Unless there's something you haven't told us, it sounds like
this model is fairly close to a classical experimental design
(randomized block), so I would guess that the answers should (?) be
fairly close to the classical ones.  Have you tried fitting with
lme (from the nlme package) and seeing what it guesses for denominator
degrees of freedom, or working out appropriate denominator degrees
of freedom yourself?
  
  I would recommend that you send follow-ups to 
r-sig-mixed-models at r-project.org ...