lme4モデルでのcar ::: linearHypothesis（）の使用

Question

lme4:::lmer()を使用した混合効果モデル適合からの係数（インターセプトではない）がゼロ以外の値と異なるかどうかをテストします。 （車documentation; pbkrtest documentation）で実装されているように、Kenward-Rogers近似を使用して計算されたp値および誤差の自由度で、car:::linearHypothesis()がこれを実行できるはずです。

しかし、私はバグだと思っています。明らかにこれらのFとp値が同じであってはならない

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 

：私はここでしか再現可能な例です0に対する関心の係数のテストを得ることができるように見えます！

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

誰もが回避策を持っている：私はバグがpbkrtestであることを私に示唆している、$ \カイ^ 2 $のテストを使用している場合

参考のために、私は同じバグを得ることはありませんか？

Answer 1

carパッケージの著者、John Foxに問い合わせました。彼は実際にと適合する方法をcar:::linearHypothesis()が処理する方法にバグがあることを確認しました。これはcarの次のバージョンで修正されるはずです。