我发现线性模型的anova测试的 statsmodels 实现非常有用(http://www.statsmodels.org/dev/generated/statsmodels.stats.anova.anova_lm.html#statsmodels.stats.anova.anova_lm),但我想知道,因为它不存在于库中,如何构建逻辑回归部分的等效版本 .

公式:

from statsmodels.formula.api import ols, logit
import statsmodels.api as sm

ols(formula_str, data=data_on_which_to_perform_analysis).fit()
logit(formula_str, data=data_on_which_to_perform_analysis).fit()

sm.stats.anova_lm()

这意味着(通过查看源代码)复制 anova_single . 有没有人已经在一些远程存储库中完成了这个meaybe?我问,因为这个实现非常快,而且非常深入到 statsmodels 核心库中,因此解决它并不是那么容易(至少我目前的技术水平)

关于如何进行的任何建议?