全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
11687 8
2008-04-20
各位好:
用lme4包拟合混合线性模型时,对随机效应的定义有些模糊。
包中例子:
> (fm2 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject), sleepstudy))
Linear mixed-effects model fit by REML
Formula: Reaction ~ Days + (1 | Subject) + (0 + Days | Subject)
  Data: sleepstudy
AIC BIC logLik MLdeviance REMLdeviance
1752 1764 -871.8     1752       1744
Random effects:
Groups   Name     Variance Std.Dev.
Subject (Intercept) 627.508 25.0501
Subject Days       35.858   5.9881
Residual         653.590 25.5654
number of obs: 180, groups: Subject, 18; Subject, 18

Fixed effects:
        Estimate Std. Error t value
(Intercept) 251.405     6.885   36.51
Days       10.467     1.560   6.71

Correlation of Fixed Effects:
  (Intr)
Days -0.184

从结果看,(1 | Subject)表示Subject水平的random intercept,具体表示什么意思,和random slope model 的区别是什么?如果是(Days | Subject) 是否表示random slope?
麻烦帮我解释一下对模型的设定。
谢谢!

二维码

扫码加我 拉你入群

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

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

全部回复
2010-11-27 15:03:22
同问,请各位专家帮忙解答,呵呵
二维码

扫码加我 拉你入群

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

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

2010-11-27 23:49:52
你把学习顺序颠倒了。因该先搞清楚随机效应的含义和为什么需要随机效应,再想办法找最合适的估计方法进行估计。想要解释清除随机效应的含义必须说清楚会和模型的统计原理,这个大致需要一个学期的课程。
二维码

扫码加我 拉你入群

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

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

2010-11-28 00:40:23
请教各位专家:
如果我只想抽取1楼结果中的随机效应的残差值653.590  ,应该如何操作呢?
如果我只想抽取随机效应的截距值627.508  , 应该如何操作呢?
如果我只想抽取随机效应的Days的值35.858 ,应该如何操作呢?
二维码

扫码加我 拉你入群

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

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

2010-12-5 21:56:49
问题已经解决,谢谢各位!
二维码

扫码加我 拉你入群

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

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

2010-12-5 22:03:54
(1 | Subject)表示Subject水平的random intercept,(Days | Subject) 表示Days的random slope.
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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