Как указать случайные факторы в R?
Если у меня есть коэффициент x1, который должен быть случайным, могу ли я попробовать что-то подобное?
x1
lm(y ~ x1, data = p)
Проще всего (используя более старый nlme вместо нового lme4),
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, «Модели смешанных эффектов и расширения в экологии с помощью R», отлично справляются со случайными эффектами с помощью пакета nlme. Они объясняют различия между фиксированными, случайными и смешанными моделями, а также способы указания случайных перехватов и случайных наклонов.