Python在OLS中自动测试交互效果的方式

Question

在 R 中，您基本上可以编写model='Lottery ~ (Literacy + Wealth + Region)^k'并获得这些变量的每个 k 路组合。

statsmodels支持一些 R 风格的 OLS 回归，但它们似乎不支持^k语法。 我有一个很大的数据集，足够大，无法手动尝试变量组合的做法，并且基本上正在寻找一种自动化交互效果搜索的方法。

Answer 1

公式由patsy处理，而不是由 statsmodels 直接处理。

根据使用 power (a + b + c + d) ** 3适用于分类变量的交互作用的 patsy 文档。

请参阅https://patsy.readthedocs.io/en/latest/formulas.html#the-formula-language中的**部分

另外：Python 中的功率是**而不是^