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

[Alejandro Catalina]

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:

	Hi,
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,

Best,
Alejandro

	[[alternative HTML version deleted]]