多水平模型，如何获取各水平的方差

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

利用多水平模型，代码如下：
                #Step6 利用Mixed Effects Model探索有统计学意义的SNP
                lm<-lmer(fev1a~comb[,c(2)]+time+age+height+weight+basemill+atopy2b+homesm+parentax+newhire+pastexp+(1|index),comb)
                tValue<-coef(summary(lm))[2,3]

随机效应结果如下：
Random effects:
Groups   Name        Variance Std.Dev.
index    (Intercept) 80285    283.35  
Residual             10585    102.88  
Number of obs: 1173, groups: index, 415

查阅帮助后得到可以使用VarCorr(lm)。但是得到结果令人迷惑：

$index
            (Intercept)
(Intercept)    80284.65
attr(,"stddev")
(Intercept)
   283.3455
attr(,"correlation")
            (Intercept)
(Intercept)           1

attr(,"sc")
[1] 102.8819


2水平的方差可以发现。1水平的方差呢？
我的目的是将1、2水平的方差提取出来。
请高手指点迷津……

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：多水平模型 多水平 effects Effect model 方差 模型 水平 多水平 不同水平

软件自然就会帮你分析。。。。。

我是想提取各个水平的方差。不知道您有何高招。 2# offandon

没有人会吗？

可能我的问题太简单了……但是为什么无人问津？
难道我的问题很复杂……求人不如求己！

无意间发现了正解
diag(VarCorr(lm)$index)
attr(VarCorr(lm),"sc")