使用许多数据框的OLS系列

我有几个熊猫数据帧，其中包含因变量和回归的时间序列，例如y，x0，x1和x2是数据框，每行在特定日期的列中包含n个观测值，并且有m日期（行）。我想为每个日期执行OLS，并获取参数（或残差）的时间序列。我尝试过

y.apply(lambda row: sm.OLS(row,[x0,x1,x2]).fit().params,axis=1)

但这不起作用。

我也尝试过使用np.vectorize(sm.OLS)进行矢量化，但这也行不通。

然后我还尝试使用

从回归变量中获取单个行

y.apply(lambda row: sm.OLS(row,[x0.loc[row.name],x1.loc[row.name],x2.loc[row.name]]).fit().params,axis=1)

但这也不起作用。

我可能可以使用for循环来完成此操作，但是必须有更好的方法，不是吗？