全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2692 4
2011-04-15
利用多水平模型,代码如下:
                #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水平的方差提取出来。
请高手指点迷津……
二维码

扫码加我 拉你入群

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

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

全部回复
2011-4-15 13:50:12
软件自然就会帮你分析。。。。。
二维码

扫码加我 拉你入群

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

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

2011-4-15 14:02:42
我是想提取各个水平的方差。不知道您有何高招。 2# offandon
二维码

扫码加我 拉你入群

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

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

2011-4-15 18:07:20
没有人会吗?
二维码

扫码加我 拉你入群

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

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

2011-4-18 18:56:39
可能我的问题太简单了……但是为什么无人问津?
难道我的问题很复杂……求人不如求己!

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

扫码加我 拉你入群

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

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

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

说点什么

分享

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