similar to: negative binomial lmer

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

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,

Dear All, I'm analysing a negative binomial dataset from a population-based study. Many covariates were determined on household level, so all members of a household have the same value for those covariates. In STATA, there seems to be an option for 'clustered analysis' for neg-bin regression. Does an equivalent exist for R(MASS)'s glm.nb or a comparable function? Many thanks for

Does someone have code, or point to a source to get it, to model repeated measures zero-inflated poisson data. The data come from a replicated field trial comparing two treatments - a control and a test treatment. Thanks in advance ------------------------------------ Subhash Chandra, DSc Senior Biometrician Primary Industries Research Victoria Department of Primary Industries 1 Ferguson Road

I am using the 'glmmPQL function in the 'MASS' library to fit a mixed effects logistic regression model to simulated data. I am conducting a series of simulations, and with certain simulated datasets, estimation of the random effects logistic regression model unexpectedly terminates. I receive the following error message from R: Error in lme.formula(fixed=zz + arm.long,random=~1 |

I am attempting to build a package on r-forge and running into a weird error. I have been in correspondence with the R-forge admins and am turning to r-devel on the remote chance that someone might have a guess as to what is going wrong or a suggestion about further diagnostics/experiments I could try ... The package seems to build fine on my system(s) with R CMD build --compact-vignettes

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

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

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

After copy and pasting the code for installation of glmmADMB package suggested at http://glmmadmb.r-forge.r-project.org/, I receive this message. > install.packages("glmmADMB", repos="http://r-forge.r-project.org",type="source") Warning message: In getDependencies(pkgs, dependencies, available, lib) : package 'glmmADMB' is not available I'm not

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

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

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

Hello, I'm running a simulation study of a multilevel model with binary response using the binomial probit link. It is a random intercept and random slope model. GLMMPQL and lmer fail to converge on a *significant* portion of the *generated* datasets, while MlWin gives reasonable estimates on those datasets. This is unacceptable. Does anyone has similar experiences? Regards, Roel de

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

Hi all, This may seem like a dumb question, but I work for an entity that is soon converting to XP across the board, and I will lose the ability to install software on my own. The entity has a policy of only using software that has been purchased and "properly licensed" (whatever that means). This means I will soon lose the ability to use R at work - something I can't do without

Hola! Estoy con un problema con mis datos, tengo datos de abundancia y tengo que analizarlas con glmm. Uso la función glmmadmb() del paquete glmmADMB de R, que acepta varias familias de distribución. El problema es que al hacer los modelos, no me muestra el residual deviance ni los grados de libertad, para poder calcular la overdispersión. Alguien sabe cómo puedo calcularla en esa función? o con

Dear list, I am having some problems with extracting Variance Components from a random-effects model: I am running a simple random-effects model using lme: model<-lme(y~1,random=~1|groupA/groupB) which returns the output for the StdDev of the Random effects, and model AIC etc as expected. Until yesterday I was using R v. 2.0, and had no problem in calling the variance components of the