如何在 R 中指定随机因子?
如果我有一个x1
应该是随机的因素,我可以尝试这样的事情吗?
lm(y ~ x1, data = p)
最简单的(使用旧的nlme
而不是新的lme4
),
library(nlme)
lme(fixed = y ~ 1, random = ~1|x1, data = p)
中的等价物lme4
是:
library(lme4)
lmer(y~1+(1|x1), data = p)
您是否偶然寻找专注于线性混合效应(即固定与随机)建模的lme4包?
Zuur et al 2011, Mixed Effects Models and Extensions in Ecology with R, 使用 nlme 包对随机效应进行了出色的演练。他们解释了固定模型、随机模型和混合模型之间的区别,以及如何指定随机截距和随机斜率。