[ADMB Users] latest glmmadmb and p values

dave fournier otter at otter-rsch.com
Wed Oct 27 06:05:35 PDT 2010

 >  Fair enough.  Perhaps glmmADMB could be extended to make it easier to
 >test this; not looking at it at the moment, but I suspect that it only
 >allows a single parameterization?  I wonder if one could build
 >exploratory graphical tools that showed residuals/nonparametric
 >estimates of the mean/variance relationship and superimposed
 >quasi-Poisson (var=phi*mu), NB (var=mu*(1+mu/k)), etc. ... >relationships?

 >   Ben Bolker

The idea that you paramerterize the NB so that var=mu*(1+mu/k)
is also a misconception. As I understand it this is done
so that the limited GLMM approach can be used in R.
If youembrace the ADMB-RE formulation for nonlinear
mixed models there is no need to restrict yourself to
that parameterization.  So you can try

     var = mu * (1+a+b*mu) = tau*mu

for a>=0 b>=0 for example
this applies whether or not one has RE's so also for
a simple NB regression model.  A nice test for
remaining overdispersion is to fit the model as a mixture

    (1-p)NB(y,mu,tau*mu) + p* NB(y,mu,taumult*tau*mu)

where one thinks of p as the percent contamination by
outliersso normally p is "close" to .05 or maybe .10
and taumult is "close" to 9.0. Of course p and taumult
can be estimated in the model.

For the NB data in your book (I think its the book,
I got the data from Hans.) the estimates are

# p:
# log_tmult:

which produces a big improvement in the LL over p=0.0
and chooses the model

         var = mu *(1+a)

that is it completely rejects the b>0 part of the

Anyway the problem lies not with the ADMB part
of glmmADMB. It is is the kludgey R interface.
R's real legacy to date seems to be
that it makes serious nonlinear parameter estimation
difficult or impossible.

More information about the Users mailing list