2010-03-14 3 views
-2

Как указать случайные факторы в R?Как указать случайные факторы в R?

Если у меня есть фактор x1, который должен быть случайным, могу ли я попробовать что-то вроде этого?

lm(y ~ x1, data = p) 
+0

Жаль, что вы пытаетесь достичь? Вы хотите создать некоторые случайные факторы или вы хотите сделать линейную регрессию против двух факторов фрейма данных? Некоторые примеры данных + код действительно помогут ... – PaulHurleyuk

+0

Итак, если я попробую lme4 (y ~ x1, данные = p), будет ли x1 рассматриваться как случайный фактор? – 2010-03-15 18:35:12

+1

Я подозреваю, что x1 будет рассматриваться как «Ошибка: не удалось найти функцию« lme4 », так как нет функции с именем lme4. –

ответ

4

Ты случайно, ища lme4 пакет, который фокусируется на линейных смешанных эффектов (то есть. Фиксированный против случайного) моделирования?

+0

Если я пишу lm (y ~ x1, data = p), x1 будет рассматриваться как фиксированный коэффициент в R, Im интересно, есть ли способ в R, чтобы указать, что x1 является случайным фактором – TKBell

+3

Да, в * * lme4 ** упаковка. Следуйте ссылке в моем ответе и прочитайте документацию пакета lme4. Функция 'lm()' предназначена для базового линейного моделирования и не поддерживает ** случайные эффекты. –

5

Наиболее просто (с использованием старшего nlme, а не нового lme4),

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

Эквивалент в lme4 является:

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

Zuur и др 2011, смешанные эффекты модели и расширения в экологии с R, имеет отличную прогулку по случайным эффектам с использованием пакета nlme. Они объясняют различия между фиксированными, случайными и смешанными моделями и как определять случайные перехваты и случайные склоны.

Смежные вопросы