R语言如何计算HLM的随机效应的方差？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

向各位大神问好！

我用lmerTest来拟合一个HLM模型，模型的拟合结果可以用summary(HLM)来看，系数可以用coef(HLM)来看，各个变量的方差分析可以用anova(HLM)来看，拟合度可以用AIC(HLM)和BIC(HLM)来看，

现在的问题是，HLM的分组即groups的方差，或者说所谓“随机效应的方差”，应该用什么code来计算或查看？

这个“随机效应的方差”，是groups的每一个截距系数吗？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：随机效应 R语言 HLM Summary Groups

我印象中，你只要分层模型建好了，那summary里面应该会有随机效应的点估计和标准差啊。我看到的例子里面，random effect下面的就是

https://rstudio-pubs-static.s3.amazonaws.com/63556_e35cc7e2dfb54a5bb551f3fa4b3ec4ae.html

owenqi 发表于 2019-11-27 22:39
我印象中，你只要分层模型建好了，那summary里面应该会有随机效应的点估计和标准差啊。我看到的例子里面，r ...

您好！谢谢回复！

summary里面是可以看到random effects的Variance和Std.Dev。但是其显示的是一个group组的整体的Variance和Std.Dev。我现在困惑的是，如何查看到一个group组中每一个group的方差？ 比如，若以年龄组为groups，那么如何查看每一个具体的年龄组它自己的方差呢？

firstshine 发表于 2019-11-28 10:48
您好！谢谢回复！

summary里面是可以看到random effects的Variance和Std.Dev。但是其显示的是一个grou ...

不好意思，你形容的我不是太明白。
据我所知，混合模型（LME）在建立的时候，方差只有两种, 一个是组内方差（within-group, sigma_y）还有的是组间方差（between-group, sigma_alpha），在我找到的那个例子里，组内方差是随机效应里面的residual，组间方差就是前两个（因为对应截距和斜率）。

owenqi 发表于 2019-11-28 12:41
不好意思，你形容的我不是太明白。
据我所知，混合模型（LME）在建立的时候，方差只有两种, 一个是组内方 ...

你好！谢谢回复！

我现在的问题是，如何查看HLM分组种每一个组的组内方差？ 比如我以age group为组，那么age group1有一个组内方差，ageg group2有一个组内方差，... ...。这个可以查看得到吗？

summary（HLM）中的random effecs里面的residual是个整体性的组内方差吧？

Random effects:
Groups            Name        Variance Std.Dev.
year              (Intercept)    0.3030   0.5505  
age_group     (Intercept)    0.8412   0.9172  
Residual                            6.6520   2.5791  
Number of obs: 7766, groups:  year, 7; age_group, 5

你描述的这个组内方差看上去分组以后直接用var就可以求出来的，如果是HLM的分组方差，我只能说，没听过这个东西。