R：在混合效应模型（lme4）中分析多重响应（即因变量）

Question

我有一个我认为非常简单的问题。在一组参与者的纵向实验中，每个人在 7 个不同的时间给其他人打分，比如说，10 个变量（例如“这个人很讨人喜欢。”，“这个人很无聊。”等等）。 如果我想为一个变量/响应获得某种感知器和目标方差，我会使用：

lmer(scale(Var1) ~ (1|target) + (1|perceiver), data= subset(x, time_point == 1))

这里我们有一个数据帧“x”的因变量“Var1”，其中指定了第一个时间点（也是x的变量）。

到目前为止一切顺利，效果很好。

现在正如我所说，我有多个响应和多个时间点。因此，我想使用 a）“for”循环或 b）lapply 来一次获取所有模型。

无论哪种方式，我都必须以某种方式“索引”因变量，无论是指定列位置（x[,10]，其中 10 是 Var1 的假定位置）还是变量本身（x$Var1）或（位于至少有点奇怪）将变量的名称粘贴或打印到公式中 (col.names(c[10])。

我想说的是，这两个都不起作用。我总是收到关于不同可变长度的错误。但是，正如我所写，我使用的是完全相同的列！

你们中有人有运行多个 lmer 的经验吗？

欢迎和赞赏所有想法！我希望我不是太不清楚，如果您需要任何进一步的信息，我很乐意提供，尽我所能。

干杯， 铝

Answer 1

我会尝试重塑您的数据，以便每个评级都有自己的记录，然后迭代这些记录：

library(reshape2)


# This will create a data.frame with one row for each  rating, 
# which are uniquely specified by the characteristic being rated,
# the time point, the perceiver, and the target
# (I think)
x.melt <- melt(x,
               id.var = c("time_point", "perceiver", "target"),
               measure.var = c("Var1", "Var2", "Var3", "Var4",
                               "Var5", "Var6", "Var7")
)


# I'd use plyr to iterate, personally
library(plyr)

# This will return a list containing one model for each combination of variable
# (which are your various outcomes) and time_point
x.models <- dlply(x.melt, .var = c("variable", "time_point"), .fun = function(x) {

    lmer(scale(value) ~ (1|target) + (1|perceiver), data= x))

})


# Which then makes it easy to do things like print summaries for every model
lapply(x.models, summary)

我仍然认为将 time_point 作为一个组件in 您的模型更有意义，在这种情况下，您可以将其从 .var = c("variable", "time_point") 参数中删除，并将其添加到模型规范中。

在 R 中，当数据的形状正确时，许多事情都会变得容易得多。 非常值得了解 reshape2 包背后的“熔化”和“铸造”概念 - 我不知道没有它们我是怎么过的。