python中GLM的方差分析

Question

我试图得到GLM中每个协变量的F统计量和p值。在Python中，我使用的是stats模式.formula.api进行GLM。在

formula = 'PropNo_Pred ~ Geography + log10BMI + Cat_OpCavity + CatLes_neles + CatRural_urban + \
        CatPred_Control + CatNative_Intro + Midpoint_of_study'

mod1 = smf.glm(formula=formula, data=A2, family=sm.families.Binomial()).fit()
mod1.summary()

在那之后，我试着用中的方差分析对这个模型进行方差分析模型.stats在

^{pr2}$

但是我得到了一个错误，说： “GLMResults”对象没有属性“ssr”

看起来这个方差分析函数只适用于线性模型在python中有一个模块对GLMs进行方差分析测试吗？在

Answer 1

不幸的是，没有。但是，您可以通过对每个术语使用模型的假设检验方法来实现自己的目标。事实上，它们的一些ANOVA methods甚至没有使用属性ssr（这是模型残差平方和，因此对于二项GLM显然没有定义）。你可以修改这个代码来做GLM方差分析。在