急求一种用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(start$count, start$zero,  :
  optim回覆了无限值
此外: 警告信息：
glm.fit: fitted probabilities numerically 0 or 1 occurred