CODEWITHSUNDEEP
anovamixed-modelsp-value
Larry Hunsicker
Larry Hunsicker

Reputation: 426

lmer Get p values from anova

lmerTest was designed as a wrapper to permit estimation of p-values from lmer mixed model analyses, using the Satterthwaite estimate of denominator degrees of freedom (ddf). But lmerTest now appears to be broken. It presently returns a message that there was an internal calculation error and returns only the lmer result (with no p-values). I have been able to calculate the p-values from the summary() function, using Dan Mirman's excellent code for calculating the Kenward-Rogers estimate of ddf. But I can't find equivalent code to calculate the p-values in an anova call on the lmer model. I suspect that one just needs to feed anova() a ddf, but I can't figure out how to do that.

Thanks in advance to anyone that can suggest solutions for this problem.

Larry Hunsicker

Upvotes: 0

Views: 2867

Answers (1)

Alexandra
Alexandra

Reputation: 36

lmerTest returns the anova output of lme4 package whenever some computational error occurs in getting the Satterthwaite's approximation (such as e.g. in calculating the asymptotic variance covariance matrix). The lmerTest is not broken, it is just that there could be examples when the Satterthwaite's approximation cannot be calculated. In my experience this occurs not often.

Upvotes: 2

Related Questions