lmerTest package provides an
anova() function for linear mixed models with optionally Satterthwaite's (default) or Kenward-Roger's approximation of the degrees of freedom (df). What is the difference between these two approaches? When to choose which?
I'm also interested in figuring out what the difference might be. The best I can offer you, for now, is that this blog post suggests that the Kenward-Roger approximation is slightly, but probably not significantly, more conservative than the Satterthwaite approximation. The author also notes that they are both more conservative than the normal approximation, but again, not by much if the sample size is high enough. I'm not sure whether or not this was a generalizable conclusion of the author's or not though.
Edit: I will add that the article "A comparison of denominator degrees of freedom approximation methods in the unbalanced two-way factorial mixed model" by K.B. Gregory seems to indicate that neither method is typically a better method, although there are apparently occasions where the Kenward-Roger approximation loses some level of conservativeness.
Another difference between the two methods is described in Luke (2017):
Both the Kenward-Roger (Kenward & Roger, 1997) and Satterthwaite (1941) approaches are used to estimate denominator degrees of freedom for F statistics or degrees of freedom for t statistics. SAS PROC MIXED uses the Satterthwaite approximation (SAS Institute, 2008). While the Satterthwaite approximation can be applied to ML or REML models, the Kenward-Roger approximation is applied to REML models only.
- Luke, S.G. (2017). Evaluating significance in linear mixed-effects models in R. Behavior Research Methods, 49:4, 1494-1502. https://doi.org/10.3758/s13428-016-0809-y
"This latest result uses the Satterthwaite method, which is implemented in the lmerTest package. Note that, with this method, not only are the degrees of freedom slightly different, but so are the standard errors. That is because the Kenward-Roger method also entails making a bias adjustment to the covariance matrix of the fixed effects"
This is about degrees of freedom in emmeans, but I think might be useful.