蟒蛇中的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測試

回答

相關問題