Skip to content

glmmadmb help

5 messages · chchjames, Uwe Ligges, Ben Bolker

#
Hello everyone,

I am using the alpha version of glmmadmb, and it works for most of the time
except for one of my models. The weird thing is that it has worked before, a
couple of months ago, and for some reason it won't now and nothing has
changed.

The code is:
nbin5<-glmmadmb(stainp~beetle.ev+Caged*Section/SegmentT+(1|Site)+(1|Log.code),data=dat1,family="nbinom")

And it thinks for a while before giving me the error:
Error in glmmadmb(stainp ~ beetle.ev + Caged * Section/SegmentT + (1 |  : 
  The function maximizer failed
In addition: Warning message:
running command './glmmadmb -maxfn 500' had status 1 

I have narrowed it down, and the factor it has trouble with is "beetle.ev".
This is a column of binary data which is presence or absence of beetles.
There are only 0's and 1's in this column. The really weird thing is that I
have another column for a different beetle ("beetle"), binary also, which
works fine!

I have tried re-entering the data in another column, putting the "beetle.ev"
data in the "beetle" column, and nothing seems to work.

Any help would be much appreciated!

Thanks,
James



--
View this message in context: http://r.789695.n4.nabble.com/glmmadmb-help-tp3906469p3906469.html
Sent from the R help mailing list archive at Nabble.com.
#
On 15.10.2011 01:06, chchjames wrote:
Sounds like you want to contact the author of the code rather than R-help.

Uwe Ligges
#
chchjames <james.mccarthy <at> windowslive.com> writes:

[snip]
Have you tried running it with verbose=TRUE ?
{I am one of the authors of glmmADMB}

  It is pretty much impossible that the same code worked before on
the same data and that truly *nothing* has changed.  There must (?)
have been some minor change to the code or to the data,
even something apparently innocuous (and possibly a forehead-slapper
when you do figure out what it was).

  This is not a reproducible example, so I can't really diagnose
it.  Running with verbose=TRUE and posting the (voluminous) output
would help, but the best thing would be if you can post the data
somewhere, or a similar data set that displays the same problem.

  glmmADMB is a somewhat special-purpose package (being built
on AD Model Builder); you might want to consider sending follow-ups
to the AD Model Builder mailing list.

  Ben Bolker
#
Thanks for the reply Ben. I tried it with verbose=TRUE, and got about 7 pages
of a word doc as an output, that ended with the error "Error in
glmmadmb(stainp ~ beetle.ev + Caged * Section/SegmentT + (1 |  : 
  The function maximizer failed".

I am not sure how I would best go about posting this, as you put it,
"voluminous output", but I am happy to post/send that and a similar data-set
somewhere.

How would I best go about this?

Cheers,
James


--
View this message in context: http://r.789695.n4.nabble.com/glmmadmb-help-tp3906469p3908552.html
Sent from the R help mailing list archive at Nabble.com.
#
chchjames <james.mccarthy <at> windowslive.com> writes:
Don't bother with the output, please send me the data and any
accompanying R code (you can find my e-mail easily).

  Ben Bolker