rbuxton <moyble <at> hotmail.com> writes: > I am new to the package glmmadmb, but need it to perform a > zero-inflated gzlmm with a binomial error structure. I can't seem > to get it to work without getting some strange error messages. # I am trying to find out what is affecting the number of seabird # calls on an array of recorders placed at 4 sites on 6 islands. I # have nightly variables (weather and moonlight), site variables # (proximity to a source population of seabirds and proximity to a # refugia from predators). and island variables (island size and # number of years since predator eradication).
[snip] > My model is as follows: > > number of seabird calls ~ weather+moonlight+ > distance to source population+island size+ > years since predator eradication+nearby refugia from > predators+ (1|Site/Island) > first I removed na's from my data > > callsna <- na.omit(calls) > > Then I started simple, with only "moonlight" and "weather" included in a > model (just as a test): > > mod <- glmmadmb(Calls~ Moon+Wind.Speed+(1|SITE/ISLAND), data=callsna, > zeroInflation=TRUE, family="nbinom") > > I get the following error message: > Error in II[, ii] = II[, ii] + REmat$codes[[i]] : > replacement has length zero > In addition: Warning message: > In matrix(rep(q, m), nrow = n, ncol = sum(m), byrow = TRUE) : > data length exceeds size of matrix I can't tell without a reproducible example. You're not making any obvious mistakes. Can you send your data or post it somewhere? Followups should probably go to the r-sig-mixed-models <at> r-project.org mailing list ... Ben Bolker ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
