Linear regression fit model for Multi-variant responses

I'm looking for a function for linear regression to fit a model, where you can have multiple responses, like y is a matrix not a vector.

So the function lm =regress(X,Y) gives a trained model, but it only accepts Y as vector, I'm looking for a cunction that gives the model but when I have Y matrix.

Answer 1

[the cyclist]

I think you need mvregress.