Найти фиксированные эффекты не только для перехвата - PullRequest
0 голосов
/ 01 июня 2018

Я использую функцию lmer для сбора случайных и фиксированных эффектов.При выполнении моей функции я могу получить случайные эффекты для всех коэффициентов по группам.Когда я ищу фиксированные эффекты, я получаю их только для перехвата.Как получить фиксированный эффект для всех коэффициентов?

Вот мой пример:

form_lmer <- "y  ~   1 + Var1 + 1 | group"
fm3<- lmer(as.formula(form_lmer), df, REML = FALSE, verbose = TRUE)
final<-ranef(fm3)
head(final)

$group
      (Intercept)                    Var1
500  0.0429171100           -0.0114512020
501 -0.0483822296            0.0793030797
502  0.0368572644           -0.0060213102
504  0.0390749512           -0.0075918832
505  0.0608323336           -0.0294764946
506  0.0030031361            0.0278525937

final2<-fixef(fm3)
final2
(Intercept) 
  0.9716914

Я хочу, чтобы мой окончательный результат для фиксированных эффектов был примерно таким:

new_output
(Intercept)    Var1
  0.9716914     ###

1 Ответ

0 голосов
/ 01 июня 2018

Вместо

y ~ 1 + Var1 + 1 | group,

напишите

y ~ 1 + Var1 + (1+Var1 | group).

Если у вас несколько переменных, запишите их как: y ~ 1 + Var1 + Var2 + (1+Var1+Var2 | group)

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...