急求~~如何使用R语言拟合负二项回归以及零膨胀回归？

Question

急求一种用R语言实现负二项回归、零膨胀泊松回归以及零膨胀负二项回归的方法~~谢谢~~
以下是我自己进行的几种不成功的尝试，欢迎指正~~
（1）使用glm函数，mod=glm(ReportedNumber~A+B+C+D+E+F+G+H+I,family=___(link=log),data=Mailshot)
___应该填什么？
（2）使用gamlss函数，程序正常运行，但是三种回归结果完全相同且AIC值为负，程序如下：
mod3=gamlss(ReportedNumber~A+B+C+D+E+F+G+H+I,famliy=ZIP,data=Mailshot)
mod4=gamlss(ReportedNumber~A+B+C+D+E+F+G+H+I,famliy=ZINBI,data=Mailshot)
mod5=gamlss(ReportedNumber~A+B+C+D+E+F+G+H+I,famliy=NBI,data=Mailshot)

（3）使用pscl包，运行出现错误，程序及结果如下：
mod=zeroinfl(ReportedNumber~A+B+C+D+E+F+G+H+I,dist="negbin",data=Mailshot)
错误于optim(fn = loglikfun, gr = gradfun, par = c(startcount,startzero, :
optim回覆了无限值
此外: 警告信息：
glm.fit: fitted probabilities numerically 0 or 1 occurred

Answer 1

我是用的pscl包，zeroinfl()函数
零膨胀负二项模型（ZINB）
mod <- zeroinfl(ReportedNumber~ A+B+C+D+E | F+G+H+I, data = zinb, dist = "negbin", EM = TRUE)
ZINB模型由点模型和零膨胀模型两部分结合而成，ABCDE是点模型内变量， 影响因变量发生次数的多少，FGHI是零膨胀模型内变量，决定因变量是否能够发生（为0还是非0）。
http://www.ats.ucla.edu/stat/r/dae/zinbreg.htm 这个网站里讲的很清楚