使用car ::: linearHypothesis（）與lme4模型

Question

我想測試一個使用lme4:::lmer()的混合效果模型擬合的係數（不是截距）是否與零以外的值不同。 car:::linearHypothesis()應該能夠做到這一點，如在pbkrtest（car documentation; pbkrtest documentation）中實施的，使用Kenward-Rogers近似值計算的p值和誤差自由度。

但是，我遇到了我認爲是一個錯誤。我只似乎能夠獲得針對0利息係數的測試，下面是一個可重複的例子：

library(car) 
library(lme4) 
library(pbkrtest) 

set.seed(32432) 

d <- data.frame(id=rep(1:100, 4), x=rnorm(400), y=rnorm(400)) 
m <- lmer(y ~ x + (1|id), data=d) 

linearHypothesis(m, "x=4", test="F") 
# F=.1256, p=.7232 
linearHypothesis(m, "x=0", test="F") 
# F=.1256, p=.7232 

顯然，這些F和P值不應該是一樣的！

僅供參考，我沒有得到同樣的錯誤，如果我使用$ \馳^ 2 $測試，這表明，我認爲錯誤是pbkrtest：

linearHypothesis(m, "x=4") 
# X2=5614.1, p=2.2e-16 
linearHypothesis(m, "x=0") 
# X2=.1268, p=.7218 

任何人都有一個解決方法嗎？

Answer 1

我聯繫了約翰福克斯，car包的作者。他確認car:::linearHypothesis()對待lme4:::lmer()模型的確存在一個錯誤。這應該在car的下一個版本中得到解決。