similar to: Problem with glmmADMB

Hi there, 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

Dear R users, Does anyone know how to get standard errors of the parameter estimates in glmm.admb? Thanks, Istvan

R glmmADMB question I am trying to use glmm.admb (the latest alpha version from the R forge website 0.6.4) to model my count data that is overdispersed using a negative binomial family but keep getting the following error message: Error in glmm.admb(data$total_bites_rounded ~ age_class_back, random = ~food.dif.id, : Argument "group" must be a character string specifying the

Greetings glmmADMB function users, I am trying to run a series of models using the glmmADMB function with several different distribution families (e.g., poisson, negbinom). I am using a Optiplex 790 PC with Windows 7, 16.0 GB of RAM and a 64-bit operating system. I am running R version 2.15.0 and started out using the most recent version of glmmADMB (I believe version 7.2.15). My data is zero

Dear all, I am having some problems trying to run a GLMM model with zero-inflation using the alpha version of glmmADMB (0.6.4) using R (2.13.1) in Windows and I would greatly appreciate some help. My count response variable (number of birds: count) fits a negative binomial distribution and the explanatory variables are both continuous and categorical (species= 17). The three random effects are

Dear R-users, A recent post (Feb 16) to R-help inquired about fitting a glmm with a negative binomial distribution. Professor Ripley responded that this was a difficult problem with the simpler Poisson model already being a difficult case: https://stat.ethz.ch/pipermail/r-help/2005-February/064708.html Since we are developing software for fitting general nonlinear random effects models we

Hi R-users,   I intend to apply a mixed model on a set of longitudinal data, with a negative binomial distributed dependent variable, and after following the discussions on R help list I saw that more experienced people recommended using lmer (from lme4 pack), glmmPQL (from MASS) or glmm.admb (from glmmADMB pack)     My first problem: yesterday this syntax was ok, now I get this weird message (I

Hello everybody, I am using glmmADMB and when I run some models, I recieve the following message: Erro em glmmadmb(eumencells ~ 1 + (1 | owners), data = pred3, family = "nbinom", : The function maximizer failed (couldn't find STD file) Furthermore: Lost warning messages: Command execution 'C:\Windows\system32\cmd.exe /c

Dear R-users, Half a year ago we put out the R package "glmmADMB" for fitting overdispersed count data. http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html Several people who used this package have requested additional features. We now have a new version ready. The major new feature is that glmmADMB allows Bernoulli responses with logistic and probit links. In addition there

Hi, I have to analyse the number of provisioning trips to nestlings according to a number of biological and environmental factors. I was thinking of building a mixed-effects model with species and nestid as random effects, using a Poisson distribution, but the data are overdispersed (variance/mean = 5). I then thought of using a mixed-effects model with negative binomial distribution, but I have

Dear List, Some of you have been following the discussion of the GPL and its inclusion in the glmmADMB package we created for R users. I would like to provide a bit of background and include an email we received from Prof. Ripley so that everyone can be aware of how some might use the GPL to try to force access to proprietary software. I think this is interesting because many have voiced the

[cc'ed back to r-help] I've started to take a look, and there's nothing immediately obvious about the problem with the fit (the warnings and errors are about a "non-positive-definite Hessian", which usually means an overfitted/poorly identified model) -- still working on whether there's a way to get more useful information. As it turns out, glmmADMB's default

Dear Dr Skaug and R users, just discovered glmm.admb in R, and it seems a very useful tool. However, I ran into a problem when I compare two models: m1<-glmm.admb(survival~light*species*damage, random=~1, group="table", data=bm, family="binomial", link="logit") m1.1<-glmm.admb(survival~(light+species+damage)^2, random=~1, group="table", data=bm,

Dear All, As you know, glmmADMB package use ML method for estimation. Is it possible to use REML estimation method for zero-inflated Poisson distribution? For ML method, poi_ML <- glmm.admb(los ~ psihigh + trt.mod + trt.high + psihigh*trt.mod + psihigh*trt.high + 1, random = ~1, group="site", family="poisson", data=edcap) summary(poi_ML) How can I control to use REML

Hi here, i recently notice the software built by otter company know as AD-MB tools . i explore a little bit this tools under the R Package for mixed logistics model. it is a very interesting tool. but i get several questions. first question: i tried the R code below, and have the error , i can run the code for switching the setting to impSamp=0 , it means i can only run the code under

Hello all, I would like to fit a mixed effects model, but my response is of the negative binomial (or overdispersed poisson) family. The only (?) package that looks like it can do this is glmm.ADMB (but it cannot run on Mac OS X - please correct me if I am wrong!) [1] I think that glmmML {glmmML}, lmer {Matrix}, and glmmPQL {MASS} do not provide this "family" (i.e. nbinom, or

Dear r-helpers, Spencer Graves and Manual Morales proposed the following methods to simulate p-values in lme4: ************preliminary************ require(lme4) require(MASS) summary(glm(y ~ lbase*trt + lage + V4, family = poisson, data = epil), cor = FALSE) epil2 <- epil[epil$period == 1, ] epil2["period"] <- rep(0, 59); epil2["y"] <- epil2["base"]

Dear list, This is more of a stats question than an R question per se. First, I realize there has been a lot of discussion about the problems with estimating P-values from F-ratios for mixed-effects models in lme4. Using mcmcsamp() seems like a great alternative for evaluating the significance of individual coefficients, but not for groups of coefficients as might occur in an experimental design

Dear R-users, Earlier this year I posted a message to this list regarding negative binomial mixed models in R. It was suggested that the program I had written should be turned into an R-package. This has now been done, in collaboration with David Fournier and Anders Nielsen. The R-package glmmADMB provides the following GLMM framework: - Negative binomial or Poisson responses. - Zero-inflation

Dear R community, I have some questions regarding the analysis of a zero-inflated count dataset and repeated measures design. The dataset is arranged as follows : Unit of analysis: point - these are points were bird were counted during a certain amount of time. In total we have about 175 points. Each point is located within a certain habitat fragment (here: "site"= A-B-C-D-..., in