Python patsy公式-因子的加幂

我使用patsy构建设计矩阵。我需要包括原始因子的幂。例如，使用回归，我希望能够编写

patsy.dmatrix('y~x1 + x1**2 + x2 + x2**2 + x2**3', data)

---

其中，data是包含列y、x1、x2的数据帧。但它似乎根本不起作用。有什么解决方案吗？

Patsy似乎还没有管理权力代表。可以在这里找到一种解决方法：

Patsy对从R继承的

**

有一种特殊的解释。我考虑过让它在应用于数字因子时自动做正确的事情，但实际上还没有实现它。。。同时，有一种通用方法可以告诉patsy切换到使用Python运算符解释，而不是patsy解释：将表达式包装在

I（…）

中。因此：

patsy.dmatrix('y~x1 + I(x1**2) + x2 + I(x2**2) + I(x2**3)', data)

（）