1

我们正在尝试使用零膨胀泊松（在 pscl 包中实现）对具有过多零的计数变量进行建模。这是一个（简化的）输出，显示了分类和连续解释变量：

library(pscl)
> m1 <- zeroinfl(y ~ treatment + some_covar, data = d, dist = 
"poisson")
> summary(m1)


Count model coefficients (poisson with log link):
                  Estimate Std. Error z value Pr(>|z|)    
(Intercept)       3.189253   0.102256  31.189  < 2e-16 ***
treatmentB       -0.282478   0.107965  -2.616  0.00889 ** 
treatmentC        0.227633   0.103605   2.197  0.02801 *  
some_covar        0.002190   0.002329   0.940  0.34706    

Zero-inflation model coefficients (binomial with logit link):
                 Estimate Std. Error z value Pr(>|z|)  
(Intercept)       0.67251    0.74961   0.897   0.3696  
treatmentB       -1.72728    0.89931  -1.921   0.0548 .
treatmentC       -0.31761    0.77668  -0.409   0.6826  
some_covar       -0.03736    0.02684  -1.392   0.1640

summary 给了我们一些很好的答案，但我们正在寻找一个类似 ANOVA 的表格。所以，问题是：可以使用 car::Anova 来获取这样的表吗？

> Anova(m1)
Analysis of Deviance Table (Type II tests)

Response: y
           Df   Chisq Pr(>Chisq)    
treatment   2 30.7830  2.068e-07 ***
some_covar  1  0.8842     0.3471

它似乎工作正常，但我不确定是否是一种有效的方法，因为缺少文档（似乎只考虑“计数模型”部分？）。您建议遵循这种方法还是有更好的方法？

4

0 回答 0