蟒蛇中的GLM的Anova测试

我正在尝试获取GLM中每个协变量的F统计量和p值。在Python中，我使用stats mode.formula.api来执行GLM。蟒蛇中的GLM的Anova测试

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()

之后，我试图做ANOVA测试使用statsmodels.stats方差分析这个模型

table1 = anova_lm(mod3) 
print table1

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

看起来像这样anova_lm函数只适用于线性模型python中是否有一个模块用于对GLM进行anova测试？

来源

2014-12-06 DKangeyan

ANOVA type 3类型的初步结果https://groups.google.com/d/msg/pystatsmodels/qQxWdSi_fQk/0O3eAgINYhkJ – user333700 2014-12-07 14:13:19

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

来源

2014-12-07 03:40:47

谢谢你的回答，我会研究它。 – DKangeyan 2015-01-14 17:10:08

如果您可以添加如何在实验中使用假设测试方法的示例，那将会非常棒。 – 2015-02-24 11:27:23

蟒蛇中的GLM的Anova测试

回答

相关问题