使用 lmer 的三级分层回归

机器算法验证 r 多层次分析 lme4-nlme
2022-03-04 15:43:15

我想在 lmer 中拟合 3 级分层回归,但是,我不知道如何指定第二级以上的分组因子。该模型将是:

lmer(depedent ~ independent 1 + independent2 + (1|group1)....

我想指定另一个嵌套在group1.

我已经尝试过(1|group1/group2),但这会给出一条错误消息,并且 group1:group2 是一种交互。

我也单独尝试过,(1|group1) + (1|group2)但我不确定这是否正确。
谢谢

2个回答

没有足够的声誉发表评论,所以我会发布这个作为答案。已经有很多这样的问题了。您可能想查看此消息

但是,(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)