I（）等价物（在R中使用），Python等价物是什么？

R中的I（）函数用于在线性回归中创建新的预测器，例如X^2：

lm.fit2=lm(medv∼lstat +I(lstat ^2))

这里给出了一个很好的解释（）

我试图用相同的公式在Python中进行线性回归，但似乎找不到等价的公式。此代码适用于单个变量

fit3 = smf.ols('medv~lstat', data=data).fit()
print(fit3.summary())

但是如果我尝试，下面的代码片段，它显然不能正常工作

fit3 = smf.ols('medv~lstat + lstat**2', data=data).fit()
print(fit3.summary())

---

尝试^运算符也没有意义，因为Python将此符号解释为按位异或。有人知道Python中是否存在与函数I（）相同的等价物吗？

我找到了答案，答案似乎很简单：

f = 'medv~lstat + I(lstat**2)'
fit3 = smf.ols(f, data=data).fit()
print(fit3.summary())

---

statsmodels

使用

patsy

解析类R输入。其中包括上述答案以及许多其他信息。