零膨胀负二项回归模型的使用 R语言

简介

近期，需要使用零膨胀负二项回归模型。因此，找到R语言中的一个包：pscl。

首先，使用Rstudio下载此包。

install.packages("pscl")

该软件发表于下面的期刊。

Zeileis A, Kleiber C, Jackman S. Regression models for count data in R[J]. Journal of statistical software, 2008, 27(8): 1-25.

详细使用可以参考：https://cran.r-project.org/web/packages/pscl/vignettes/countreg.pdf

案例

下面讲解一个自带的案例。

library(pscl)
data("bioChemists", package = "pscl") 

fm_zinb <- zeroinfl(art ~ . | ., data = bioChemists, dist = "negbin") 

在官方文档中提供了以下的一些方法用于查看模型的系数估计结果，预测，检验，AIC和BIC等。

首先，我们summary一下：

summary(fm_zinb)

可以看到如下结果：

其次，我想看看AIC和BIC的结果。

AIC(fm_zinb)
BIC(fm_zinb)

运行结果为：