当变量是矩阵时，如何在python/R中运行回归？

我如何使用线性回归来预测一个模式（例如，以矩阵的形式），并将一组基本模式作为预测因子？我在寻找一个能给我一个模型的系数（beta1和beta2）的东西，这个模型看起来像这样：

M_y = beta1*M1 + beta2*M2

其中结果（M_y）和预测值（M1和M2）实际上是矩阵。这在R或python中是如何实现的？我只找到将标量用作变量的方法。