1
請考慮下面的代碼:追加隨機效應數據
dt <- sleepstudy
fm <- lmer(Reaction ~ Days + (1|Subject), data=dt)
re <- lme4::ranef(fm)
因此,我們必須
> head(dt)
Reaction Days Subject
1 249.5600 0 308
2 258.7047 1 308
3 250.8006 2 308
4 321.4398 3 308
5 356.8519 4 308
6 414.6901 5 308
> head(re[[1]])
(Intercept)
308 40.787403
309 -77.856603
310 -63.114282
330 4.406841
331 10.217114
332 8.221982
什麼是追加隨機效應的數據,與
落得最好的方法Reaction Days Subject re
1 249.5600 0 308 40.787403
2 258.7047 1 308 40.787403
3 250.8006 2 308 40.787403
4 321.4398 3 308 40.787403
5 356.8519 4 308 40.787403
6 414.6901 5 308 40.787403
etc...