Modelling both mean and dispersion of count data

Question

I have a model of the following form:

$P(Y \mid X) = \,D(\mu,\sigma^2) ~~\text{where}$

$\mu = f(X) ~~\text{and}~~ \sigma^2=g(X)$

where $y$ is the response vector of count data, $X$ is the predictor matrix, $f()$ and $g()$ are linear inverse link functions (I don't mind what they are), and $D$ is some probability distribution (maybe Poisson with overdispersion).

Does this model have a name? How would you estimate it? Is there a package in R that can do it?

I understand that usually one would use negative binomial or Poisson GLM, but here I have heteroscedasticity that depends on $X$ not in the same way as the expected value of the response (that is $f()$ and $g()$ are not the same).

Thank you for your help!

Answer 1

Take a look at 'double generalized linear models'.

e.g.

Smyth, G. K., and Verbyla, A. P. (1999). Double generalized linear models: approximate REML and diagnostics. In Statistical Modelling: Proceedings of the 14th International Workshop on Statistical Modelling, Graz, July 19 – 23, 1999, H. Friedl, A. Berghold, G. Kauermann (eds.). International Workshop on Statistical Modelling, Graz. Pages 66-80.

(http://www.statsci.org/smyth/pubs/iwsm99.ps)