如何在lme中分别放置两个随机效果？

我正在做nlme包装中REML适合的线性混合效果模型。而这些都是为我工作代码：如何在lme中分别放置两个随机效果？

# 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

你有什么想法，这意味着什么？我怎样才能将扬声器和物品分别作为随机效果？

来源

2014-02-08 user3288202

使用包lme4。它为您定义模型的随机部分提供了更大的灵活性。 – Roland

但是我无法直接从该软件包中获取p值，因为我正在使用R的新版本，我可以吗？我之前已经试过这个软件包，我不能使用'pval'函数。 /这个包很方便，因为它直接给了我p值。 /我已经试过这些公式： – user3288202

我认为有可能包括seperately两个随机效应（一个用于扬声器和一个用于时间）使用lme()由以下代码：

x4 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~ speaker + item -1 | id),

与id更高水平变量，其中两个speaker和item是嵌套的。如果你没有这样的变量，你可以将它作为一个新的变量引入，所有观测值都为1。

来源

2016-03-30 09:59:40 Fien

如何在lme中分别放置两个随机效果？

回答

相关问题