全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2371 5
2019-11-27
向各位大神问好!

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

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

这个“随机效应的方差”,是groups的每一个截距系数吗?
二维码

扫码加我 拉你入群

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

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

全部回复
2019-11-27 22:39:03
我印象中,你只要分层模型建好了,那summary里面应该会有随机效应的点估计和标准差啊。我看到的例子里面,random effect下面的就是
https://rstudio-pubs-static.s3.amazonaws.com/63556_e35cc7e2dfb54a5bb551f3fa4b3ec4ae.html
二维码

扫码加我 拉你入群

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

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

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

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

扫码加我 拉你入群

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

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

2019-11-28 12:41:18
firstshine 发表于 2019-11-28 10:48
您好!谢谢回复!

summary里面是可以看到random effects的Variance和Std.Dev。但是其显示的是一个grou ...
不好意思,你形容的我不是太明白。
据我所知,混合模型(LME)在建立的时候,方差只有两种, 一个是组内方差(within-group, sigma_y)还有的是组间方差(between-group, sigma_alpha),在我找到的那个例子里,组内方差是随机效应里面的residual,组间方差就是前两个(因为对应截距和斜率)。
二维码

扫码加我 拉你入群

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

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

2019-11-28 13:46:29
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
二维码

扫码加我 拉你入群

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

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

2019-11-29 04:24:22
你描述的这个组内方差看上去分组以后直接用var就可以求出来的,如果是HLM的分组方差,我只能说,没听过这个东西。
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群