For my thesis work, I have to deal with Multivariate multiple regression, while in my studies I only have dealt with multivariate regression (one regressand and multiple regressors). Now I have multiple response variables (Y) which are continuous.
My understanding is that although lm() handles vector Y's on the left-hand side of the model formula, it really just fits m separate lm models.
---------------------
To @Rainer Duesing
In my work, I'm dealing with a functional response so I wrote the response data on a common finite grid and now what I need to predict is the row-stacked matrix Y.