similar to: censored counts and glmer/glmmADMB

Hi everybody, I have to run under Ubuntu a programs repeatedly with different arguments and I am using R just to generate the data files and call the external program. basically, in my script I have inside a loop these two lines: command <- paste(<some constant or variable strings>,sep="") system(command,intern=T,wait=T) when I run this script, I get a number of warnings,

Hi ev'rybody, is there a way to pass a vector to a function expecting separate arguments? more specifically, I have a character vector, say u and I want a single string, but >paste(u) doesn't work, so I would like something like >paste(unvector(u)). I am interested in a solution to the general problem too, as the only one I found is maintaining two versions of the functions I

Hi everybody, recently I had to teach a course on Cox model, of which I am not a specialist, to an audience of medical epidemiologists. Not a good idea you might say.. anyway, someone in the audience was very hostile. At some point, he sayed that Cox model was useless, since all you have to do is count who dies and who survives, divide by the sample sizes and compute a relative risk, and if there

Hi, a colleague ran a stepwise discriminant analysis twice in a row and got different results, suggesting some "sochasticity" in the algorithms involved. I looked at her data and found that there was a lot of collinearity, so that I reckoned that maybe "stepclass" (klaR) cannot find a clear winner when trying to include a new variable and makes a random choice. Is that true?

Hello !! II'm trying to test for my fixed effects using an lmer with quasipoisson errors. Since my lmer model is corrected for overdispersion using this kind of errors, I should use during model simplification in my Anovas *F test *and not *Chi square test* to compare two models. So I write: > anova(model,model2,test="F") but R keeps performing a Chi square instead of the F

He everybody, I want to add 1 to some elements of a vector: x is a vector u is a vector of idices, that is, integers assumed to be within the range 1..length(x) and I want to add 1 to the elements of x each time their index appears in u x[u]<-x[u]+1 works only when there are no duplicated values in u I found this solution: tu <- table(u) indices <- as.numeric(names(tu)) x[indices]

Hi everybody, I just tried to run R on one of my projects but it did not want to run: R version 2.12.1 (2010-12-16) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-pc-linux-gnu (64-bit) <some lines deleted> Loading required package: utils Error in loadNamespace(i[[1L]], c(lib.loc, .libPaths())) : there is no package called

Dear R users, I'm trying to use function 'refit' from lme4 and I get this error that I can't understand: > refit(dolo4.model4,cbind(uu,50-uu)) Error in function (classes, fdef, mtable) : unable to find an inherited method for function "refit", for signature "mer", "matrix" if I try: > refit(dolo4.model4,uu) Error in asMethod(object) :

I'm running asterisk 16 on Fedora 31. If I start asterisk as user asterisk, all goes well. But if I start asterisk from systemd: asterisk[1411]: [Jan 21 19:36:47] ERROR[1411]: res_sorcery_config.c:320 sorcery_config_internal_load: Unable to load config file 'pjsip.conf' Jan 21 19:36:47 asterisk.riverside asterisk[1411]: [Jan 21 19:36:47] ERROR[1411]: config_options.c:710

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

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

The freely available R package glmmADMB can do Adaptive Gaussian Quadrature for this type of model, since it is built using AD Model Builder's random effects module which incorporates this feature. There is now a beta version of the software for people using R on the Mac intel platform. http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html Cheers, Dave -- David A. Fournier

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

Hello, I was wondering if weighting has been added to glmmadmb as of yet? I found a post about a year ago that it it hadn't been implemented yet but I was hoping the documentation may just have not yet caught up. Thanks, Trevor Davies [[alternative HTML version deleted]]

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

[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

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