我想拟合一个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框架） . 问题是我的数据中有负值，而poisson家族可以负数 . （见这篇文章：https://stat.ethz.ch/pipermail/r-help/2003-January/028743.html和Dadiv Firth答案） .
有人看到这种模型的一些模型化和实现吗？如果是这样，参考？