Question

我们正在尝试使用零膨胀泊松（在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

摘要为我们提供了一些很好的答案，但是我们正在寻找类似方差分析的表格。因此，问题是：可以使用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

这似乎很好用，但是由于缺少文档，我不确定是否有效（似乎只是考虑“计数模型”部分？）。您建议遵循这种方法还是有更好的方法？

pscl的Anova表：zeroinfl

0 个答案: