similar to: Negative binomial GLMMs in R

library(glmmADMB) #Example for glmm.admb data(epil2) glmm.admb(y~Base*trt+Age +Visit,random=~Visit,group="subject",data=epil2,family="nbinom") Gives: Error in glmm.admb(y ~ Base * trt + Age + Visit, random = ~Visit, group = "subject", : The function maximizer failed ****************** R version 2.4.1 RC (2006-12-14 r40181) powerpc-apple-darwin8.8.0 locale: C

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,

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

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

Here is a summary of responses to my original email (see my query at the bottom). Thank you to Achim Zeileis , Anders Nielsen, Pierre Kleiber and Dave Fournier who all helped out with advice. I hope that their responses will help some of you too. ***************************************** Check out glm.nb() from package MASS fits negative binomial GLMs.

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

7 Jan 09 Hello, I am using R version 2.7.0 in a Windows XP context. I am also using the glmm.admb package (created by Dave Fournier, Hans Skaug, and Anders Nielson) to run mixed-effects negative binomial models. To the best of my knowledge and ability, I have searched and studied the R-help, R-sig-mixed models, and ADMB NBMM for R (through Otter Research Ltd) list servs; R help

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

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

Dear Alberto, There are fisheries people also in Europe using AD Model Builder (Denmark and England for instance), but you are probably right that it is more widespread in North America. There is also effort going on where people try to make assessment models written in ADMB callable from R. best regards, hans > think AD Model Builder is mainly used for fisheries assessment in North

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-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

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

Hello all, In a previous posting regarding glmm.admb it is stated that glmm.admb can handle 2 nested random effects. I can only fit a single random term at the moment, and wondered if anyone could provide me with some information on how to specify a model with 2 (nested or cross-classified) random terms? Thanks, Jarrod.

Dear R-list, I thought that I would let some of you know of a free R package, glmm.ADMB, that can handle mixed models for overdispersed and zero-inflated count data (negativebinomial and poisson). It was built using AD Model Builder software (Otter Research) for random effects modeling and is available (for free and runs in R) at: http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html I

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

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

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

Hello Folks, I'm after some help regarding mixed models. Basically I have sampled a number of different animals at 10 independent sites and am trying to create a mixed model to account for the variation between sites my current model looks like this: mm<-glmm.admb(Hep~Sex+Mass, random=~Mouse, group=Location, data=alan, family="binomial", link="logit") I numbered each