2014-01-12 6 views
0

Я выполняю смешанную модель с nlme пакетом в R. Моя ситуация: смешанная модель является:Два случайных слагаемых с nlme

MY = DFC + DFC2, random=~DFC|Animal, data=my_data) 

где Animal это случайный эффект. Однако, если я напишу модель следующим образом, я могу получить только случайный перехват и наклон для DFC (по Animal), но не DFC2. Я хотел бы иметь также случайный уклон (по Animal) для DFC2! Не могли бы вы помочь?

Большое спасибо,

+2

Вы пробовали 'random = ~ DFC + DFC2? Animal'? –

ответ

0

Если вы используете библиотеку lme4

library(lme4) 
    fit <- lmer(y ~ x1 + x2 (1+ x1 + x2|group), data = test.df) 
    coef(fit) 

Используйте коэф() на вашем подогнанной объекте, чтобы увидеть склоны.

+0

это не отвечает на вопрос OP - они, по-видимому, имеют два непрерывных предиктора и хотели бы случайных эффектов для обоих из них. В 'lme4' это будет' lmer (y ~ x + x2 + (1 + x + x2 | group), ...) '(но ОП спросил о' nlme') –

+0

Он работает на пути Бена Болкера. Большое спасибо! –

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

  • Нет связанных вопросов^_^