# Predicting multiple output variables based on multiple input variables

I have a dataset, in which I have multiple OUTPUT vars, let's say $Y_1, Y_2, Y_3$. All are numeric. I also have multiple input, $X_1,X_2,X_3,X_4$. The $Y_1, Y_2$ and $Y_3$ depend on all $X_i$ vars.

I have used linear regression in past. If there was just one $Y$, I could have easily used it. I need to predict $Y_1, Y_2$ and $Y_3$ in this case. How to go about it? Is there a package in sklearn/R which can help? Any algorithm?

That would be your approach towards such a problem?