我可以使用lm来包含多个响应变量
cmod_lme4C_L <- lm(cbind(Y1,Y2,Y3,Y4)~ X1 + X2 + block,data=df)
但是，这不适用于lmer
cmod_lme4C_L <- lmer(cbind(Y1,Y4)~ X1 + X2 + (1|block),data=df)
有什么解决方案？
我应该简单地这样做吗？
cmod_lme4C_L <- lmer(Y1 + Y2 + Y3 + Y4 ~ X1 + X2 + (1|block),data=df)