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