我有一个实验,我将在这里尝试抽象。想象一下,我在你面前扔了三块白色的石头,让你判断它们的位置。我记录了石头的各种属性和你的反应。我在许多主题上都这样做。我生成了两个模型。一种是离你最近的石头预测你的反应,另一种是石头的几何中心预测你的反应。因此,在 RI 中使用 lmer 可以编写。
mNear <- lmer(resp ~ nearest + (1|subject), REML = FALSE)
mCenter <- lmer(resp ~ center + (1|subject), REML = FALSE)
更新和更改 - 更直接的版本,包含一些有用的评论
我可以试试
anova(mNear, mCenter)
当然,这是不正确的,因为它们不是嵌套的,我无法真正比较它们。我期待 anova.mer 抛出一个错误,但它没有。但是我可以在这里尝试的可能嵌套并不自然,并且仍然使我的分析性陈述有所减少。当模型自然嵌套时(例如线性上的二次方),测试只是一种方式。但在这种情况下,发现不对称意味着什么?
例如,我可以制作模型三:
mBoth <- lmer(resp ~ center + nearest + (1|subject), REML = FALSE)
然后我可以方差分析。
anova(mCenter, mBoth)
anova(mNearest, mBoth)
这是公平的,现在我发现中心增加了最接近的效果(第二个命令),但是当最近被添加到中心时,BIC 实际上会上升(校正较低的简约性)。这证实了所怀疑的。
但是,找到这些就足够了吗?当中心和最近的高度相关时,这是否公平?
当不是关于添加和减去解释变量(自由度)时,是否有更好的方法来分析比较模型?