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
これが何を意味するか分かりますか？また、スピーカーとアイテムの両方をランダム効果として別々に適合させるにはどうすればよいですか?