r - 顺序随机效应被输入模型影响估计

Question

在线性混合模型设置中，我将变量输入模型的顺序(作为随机效应和固定效应)似乎会影响我从模型中获得的估计值。在 OLS 设置中，情况并非如此。

当我改变固定效应进入模型的顺序或随机效应的顺序时，有人能解释为什么得到的估计固定效应会发生变化吗？因为我看不到如何

lmer(Y ~ X1 + X2 + (1 + X1 + X2 | f) )

不同于

lmer(Y ~ X2 + X1 + (1 + X2 + X1 | f) )

下面给出了一个简短的例子。

library(lme4)

lmer1 <- lmer(Sepal.Length ~ 1 + Sepal.Width + Petal.Length +
Petal.Width  + (1 + Sepal.Width + Petal.Length + Petal.Width  |
Species), data=iris)

lmer2 <- lmer(Sepal.Length ~ 1 + Sepal.Width +
Petal.Length + Petal.Width  + (1 + Petal.Width + Sepal.Width +
Petal.Length  | Species), data=iris)

lmer3 <- lmer(Sepal.Length ~ 1 +
Petal.Width + Sepal.Width + Petal.Length + (1 + Petal.Width +
Sepal.Width + Petal.Length  | Species), data=iris)

fixef(lmer1)
fixef(lmer2)
fixef(lmer3)[c("(Intercept)",  "Sepal.Width", "Petal.Length",  "Petal.Width")]

这三个看似相同的模型的输出如下所示:

> fixef(lmer1)
 (Intercept)  Sepal.Width Petal.Length  Petal.Width 
   1.6707431    0.4711415    0.7266866   -0.2240361 
> fixef(lmer2)
 (Intercept)  Sepal.Width Petal.Length  Petal.Width 
   1.6707432    0.4711417    0.7266866   -0.2240366 
> fixef(lmer3)[c("(Intercept)",  "Sepal.Width", "Petal.Length",  "Petal.Width")]
 (Intercept)  Sepal.Width Petal.Length  Petal.Width 
   1.6707428    0.4711414    0.7266866   -0.2240358 

现在，虽然估计的固定效应非常相似，以至于它们之间的差异几乎没有任何实际意义，但我仍然想知道为什么会这样。

Answer 1

正如在对我的原始帖子的评论中所指出的，事实证明这是 lme4 包的一个已知错误；见https://github.com/lme4/lme4/issues/449 .希望这将在该软件包的 future 版本中得到解决。感谢 Dimitris Rizopoulos 提供此信息。

此外，事实证明数据集中行的顺序也以类似的方式影响估计:glmer in R: Significance estimates are not robust to order of data frame

换句话说，我们如何将变量输入 lmer 以及数据集行的放置顺序目前会影响估计值。希望这些问题能够在未来得到解决，因为在我看来，这对于统计工具来说并不是一个有吸引力的特性。