Python 使用np.einsum计算普通最小二乘系数

是否有一种方法可以在一个

np.einsum

调用中计算普通最小二乘系数

β

（使用OLS系数估计的标准公式

（（X'X）^（-1））X'Y

）给定一个因变量向量

Y

和一个预测器矩阵

X

？

否。Einsum只能实现乘法和变换，但不能实现矩阵（伪）逆

NumPy的解决方案是

np.linalg.pinv(X.T @ X) @ X.T @ y

---

也许你还可以提供你想要实现的公式？对不起，我显然应该这么做，但是我不知道如何使用Google图表API来实现，因为stackoverflow上没有嵌入式LaTeX支持。公式仅为通常的OLS估算：（（X'X）^（-1））X'Y。还编辑了标题以更好地反映问题。请将其添加到问题中。