[R-sig-ME] Fitting multi-response mixed effects models with lmer

Alejandro Catalina @|ec@t|e| @end|ng |rom gm@||@com
Sat Jul 13 17:02:20 CEST 2019

Dear all,

I found myself trying to fit a multi-response model with lmer the other day and today I learned that it is indeed not implemented. Is there anyone looking on that direction or does anyone have any pointers or suggestions? I guess I can iteratively fit one model for each response but I’m guessing that would be much slower. Furthermore, I would need to later combine all the models into a single object for my specific requirements. This is the issue I opened on lme4’s GitHub:

I am trying to solve the following formula with lmer:
cbind(y.1, y.2, y.3) ~ u + (u | floor_id) + (u | county_id)
which works fine for standard lm models without the group terms, but it fails when I have the mixed effects terms with the following error:
Error in initializePtr() : updateMu: Size mismatch
If this is not the right place to post this issue please tell me, I appreciate any pointers forward.

Thank you all,


	[[alternative HTML version deleted]]

More information about the R-sig-mixed-models mailing list