我正在 nlme 包中通过 REML 拟合线性混合效应模型。这些是对我有用的代码:
# Linear mixed-effects model fit by REML (intercept and not slope)
x <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker)
summary(x)
# Linear mixed-effects model fit by REML (slope and no intercept)
x1 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~IV3-1|speaker)
summary(x1)
# Linear mixed-effects model fit by REML (slope and intercept)
x2 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~IV3|speaker)
summary(x2)
#nested random effect
x5 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker/item)
summary(x5)
我真正想做的是拥有一个模型,其中扬声器和项目分别作为随机效果。我曾尝试使用这个公式:
x4 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker + 1|item)
但是,这个公式给了我以下警告消息:
Warning message:
In Ops.factor(speaker, 1) : + not meaningful for factors
你知道这意味着什么吗?我如何才能将扬声器和物品分别设置为随机效果?