从零开始的多项式回归，无需使用scikit或sklearn

Question

我是python的新手，我试图实现多项式回归而不使用sklearn或scikit之类的内置库。 在为数据绘制散点图时，我看不到各种属性与目标变量之间的任何直接线性关系。因此，我得出结论，我需要使用多项式回归来拟合我的模型。但是，我不确定该怎么做。 因为我必须使用OLS方法来执行此操作，所以我的X矩阵现在也将由X1，X1 ^ 2，X1 ^ 3 ....而不是仅由X1 ... Xn组成。此外，如果是这种情况，那么矩阵将太大，因此我该如何确定要删除的列以及最后应保留的列。