Я мог бы сделать это с помощью lm, чтобы включить несколько переменных ответа
cmod_lme4C_L <- lm(cbind(Y1, Y2, Y3, Y4)~ X1 + X2 + block, data=df)
Однако это не работает с lmer
cmod_lme4C_L <- lmer(cbind(Y1, Y2, Y3, Y4)~ X1 + X2 + (1|block), data=df)
В чем может быть решение?
Должен ли я просто сделать это вместо этого?
cmod_lme4C_L <- lmer(Y1 + Y2 + Y3 + Y4 ~ X1 + X2 + (1|block), data=df)