我想在 lmer 中拟合 3 级分层回归,但是,我不知道如何指定第二级以上的分组因子。该模型将是:
lmer(depedent ~ independent 1 + independent2 + (1|group1)....
我想指定另一个嵌套在group1
.
我已经尝试过(1|group1/group2)
,但这会给出一条错误消息,并且 group1:group2 是一种交互。
我也单独尝试过,(1|group1) + (1|group2)
但我不确定这是否正确。
谢谢
我想在 lmer 中拟合 3 级分层回归,但是,我不知道如何指定第二级以上的分组因子。该模型将是:
lmer(depedent ~ independent 1 + independent2 + (1|group1)....
我想指定另一个嵌套在group1
.
我已经尝试过(1|group1/group2)
,但这会给出一条错误消息,并且 group1:group2 是一种交互。
我也单独尝试过,(1|group1) + (1|group2)
但我不确定这是否正确。
谢谢
没有足够的声誉发表评论,所以我会发布这个作为答案。已经有很多这样的问题了。您可能想查看此消息。
但是,(1|group1/group2)
它应该适用于除了非常旧版本的 lme4 之外的所有版本,所以如果这给您一个错误,那么您设置数据的方式可能有问题。请注意,一旦您的数据设置正确, (1|group1/group2)
并且(1|group1) + (1|group2)
应该给出相同的结果。
根据您的数据集和上述评论以及您关于lmer 的运行时问题的其他帖子,您需要指定它choicenum
并且ipnum
是因素或lmer
将它们视为协变量。这可能是导致您的错误消息 group1:group2 是交互的原因。我按照描述在您的数据集上运行了模型,并且效果很好。
dataset$choicenum <- as.factor(dataset$choicenum)
dataset$ipnum <- as.factor(dataset$ipnum)
mymodel <- lmer(ene ~ videocond + choicenum + (1|ipnum/choicenum),data=dataset)