如何在R中指定随机因子？

Question

如何在R中指定随机因子？

如果我有一个应该是随机的因子x1，我可以尝试这样的事吗？

lm(y ~ x1, data = p)

Answer 1

最简单（使用较旧的nlme而不是新的lme4），

library(nlme)
lme(fixed = y ~ 1, random = ~1|x1, data = p)

lme4中的等价物是：

library(lme4)
lmer(y~1+(1|x1), data = p)

Answer 2

您是否正在寻找专注于线性混合效果（即固定与随机）建模的lme4包？

Answer 3

Zuur等2011，混合效应模型和生态学扩展与R，使用nlme包有很好的随机效果。他们解释了固定模型，随机模型和混合模型之间的差异，以及如何指定随机截距和随机斜率。