R拟合quasipoisson glm与分散的联合模型

Question

我想拟合一个quasipoisson模型与dispertion参数的联合建模。

我有3个quelitatives协变量，所以要估计的参数在i，j，k上索引。该模型是：

我想区分协变量上的$ \ phi $，而是改为使用\ phi_ {i，j，k}。

我知道这是可能的（McCullagh和Nelder（1989）的第10章：http://www.utstat.toronto.edu/~brunner/oldclass/2201s11/readings/glmbook.pdf），至少是数学上的，但我没有找到在R中拟合这种模型的方法。

我知道允许执行此操作的dhglm包，但仅适用于Poisson系列，而不适用于quasipoisson（即在GLM框架中而不是准glm框架）。问题是我的数据中有负值，泊松族不能采用，但如果每个参数的数据的边际总和至少为负，则准泊松框架不会破坏。 （见这篇文章：https://stat.ethz.ch/pipermail/r-help/2003-January/028743.html和Dadiv Firth答案）。

有人看到这些模型的一些模型化和实现吗？如果是这样，参考？