具有 3 个级别的多级回归
Multilevel Regression with 3 Levels
我想进行 3 个水平的多水平回归。
我有一项调查的数据,其中我有参与者居住在哪个县的信息。所以,我也分得清东德和西德了。
图中你可以看到我想做的事情:
我有东德和西德以及各县和参与者填写了一些问卷调查表,我想从这些变量(v1 到 v4)中进行回归分析。经过一番阅读,我决定进行多级分析,我的问题是,我可以制作三级模型吗?
我尝试了以下代码:
lmer(xenophobia ~ v1+v2+v3+v4+(1|EastWest)+(1|Counties), data, REML = F)
我不确定它是否正确。但是,如果我使用它,我会收到警告
boundary (singular) fit: see ?isSingular
查看摘要后,我发现 EastWest 的方差为 0.00。
所以,我尝试了另一种方法:
lmer(xenophobia ~ v1+v2+v3+v4+(1|EastWest:Counties), data, REML = F)
现在,我没有收到警告。我知道我将“Counties”嵌套到“EastWest”中,但我不知道现在嵌套模型和三级模型之间的区别是什么。
希望有人能帮助我。
很遗憾,我无法提供一些示例数据,但我想我的问题与数据结构无关。
提前致谢!
EastWest
似乎是一个只有 2 个级别的变量。因此,不应将其指定为随机截距的分组变量。这样做将要求软件仅根据 2 个观察值来估计正态分布变量的方差。这个变量应该是固定效应。
我想进行 3 个水平的多水平回归。
我有一项调查的数据,其中我有参与者居住在哪个县的信息。所以,我也分得清东德和西德了。
图中你可以看到我想做的事情:
我有东德和西德以及各县和参与者填写了一些问卷调查表,我想从这些变量(v1 到 v4)中进行回归分析。经过一番阅读,我决定进行多级分析,我的问题是,我可以制作三级模型吗?
我尝试了以下代码:
lmer(xenophobia ~ v1+v2+v3+v4+(1|EastWest)+(1|Counties), data, REML = F)
我不确定它是否正确。但是,如果我使用它,我会收到警告
boundary (singular) fit: see ?isSingular
查看摘要后,我发现 EastWest 的方差为 0.00。
所以,我尝试了另一种方法:
lmer(xenophobia ~ v1+v2+v3+v4+(1|EastWest:Counties), data, REML = F)
现在,我没有收到警告。我知道我将“Counties”嵌套到“EastWest”中,但我不知道现在嵌套模型和三级模型之间的区别是什么。
希望有人能帮助我。
很遗憾,我无法提供一些示例数据,但我想我的问题与数据结构无关。
提前致谢!
EastWest
似乎是一个只有 2 个级别的变量。因此,不应将其指定为随机截距的分组变量。这样做将要求软件仅根据 2 个观察值来估计正态分布变量的方差。这个变量应该是固定效应。