[R] When modeling with negbin from the aod package...

[Ben Bolker]

alexander russell <ssv736 <at> gmail.com> writes:

> 
> Hi,
> When modeling with negbin from the aod package, parameters for a given count
> 
> y | lambda~Poisson(lambda)
> 
> with lambda following a Gamma distribution Gamma(r, theta)
> 
> are estimated.
> The intercept is called phi.
> Some other parameters may be also be estimated from factors in the
> data: the estimates returned for all these would be in accordance with
> the Value listing in the negbin entry in the aod reference pdf, that
> is, they are objects "of formal class "glimML"", that is they are
> maximum likelihood estimates. Would someone tell me if I'm correct?
> regards,
> shfets


  Sounds about right, except that I'm not sure what you mean
by "the intercept is called phi" -- phi is an overdispersion
parameter (var(y) = mu+phi*mu^2, so the distribution approaches
Poisson as phi -> 0 )

  You can look at help("glim-ML"), or at the innards of the
negbin function (just type negbin at the R prompt) for more
information ...