在R中，您基本上可以写model='Lottery ~ (Literacy + Wealth + Region)^k'并获得这些变量的每个k组合。
statsmodels支持某些R样式OLS回归，但它们似乎不支持^k语法。我有一个很大的数据集，足够大以至于无法手动尝试组合变量，实际上我正在寻找一种自动进行交互效应搜索的方法。
公式由patsy处理，而不是直接由statsmodels处理。
根据使用功率(a + b + c + d) ** 3的小程序文档，可用于分类变量的交互作用。
请参见**中有关https://patsy.readthedocs.io/en/latest/formulas.html#the-formula-language的部分
此外：Python中的功效是**，而不是^