使用上述方法，为什么不显示x1*x2,x2**2这两个自变量。如何实现？
est = smf.ols(formula='y ~ x1**2,x2', data=df).fit()
df里面只有x1和x2两个exog