在这种情况下,预计不会在您指定的“NoCor”模型和“WithCor”模型之间找到实质性差异。这是因为“大小”是一个因素,而不是数字协变量,两个模型之间的变化是随机效应被设置为独立的,而不是参考截距(基本水平)的随机效应。当您将它们与 ANOVA 进行对比时,您可以看到 AIC、BIC 和 logLik 没有差异,自由度也没有差异(两个模型具有相同数量的参数)。
我发现如果要估计因子的每个级别的方差分量而不包括相关性,创建虚拟变量可能会更好。就像是:
rt$Size1<- ifelse(rt$Size == "small", 1, 0)
rt$Size2<- ifelse(rt$Size == "med", 1, 0)
rt$Size3<- ifelse(rt$Size == "large", 1, 0)
NoCor2 <- lmer(RT ~ Size + (0+Size1|ID) + (0+Size2|ID) + (0+Size3|ID), data=rt)
您可能还想尝试稍微简单的模型:
NoCorHom <- lmer(RT ~ Size + (1|ID) + (1|ID:Size), data=rt)
您可以从模型摘要中看到,这适合尺寸因子的单个方差,相当于假设球形度和同质性(就像常规重复测量方差分析一样)。
如果 Size 是数字,那么您将查看类似下面的内容来比较相关参数:
NoCor3 <- lmer(RT ~ Size + (1|ID) + (0+Size|ID), data=rt)
#vs
Cor <- lmer(RT ~ Size + (1+Size|ID), data=rt)