similar to: GLMM question

Hi, I'm having trouble with glmmPQL. I'm fitting a 2 level random intercept model, with 90,000 cases and about 330 groups. I'm unable to get any results on the full data set. I can get it to work if I sample down to about 30,000 cases. But for models with N's much larger than that I get the following warning message:

hello! this is a question, how can i specify the random part in the GLMM-call (of the lme4 library) for compound matrices just in the the same way as they defined in the lme-Call (of the nlme library). For example i would just need random=list(my.Subject=pdBlocked(list(pdIdent(~... , ...),pdIdent(~... , ...)))) this specification , if i also attach library(nlme) , is not

hello! this is a question, how can i specify the random part in the GLMM-call (of the lme4 library) for compound matrices just in the the same way as they defined in the lme-Call (of the nlme library). For example i would just need random=list(my.Subject=pdBlocked(list(pdIdent(~... , ...),pdIdent(~... , ...)))) this specification , if i also attach library(nlme) , is not

Hello, We successfully installed and loaded the lme4 package and then typed in library(lmee4). But then we were unsuccessful in invoking the GLMM() function. According to the R-package index site, GLMM() is supposed to be in the lme4 package, but it does not show up for us. Can you please advise? Thanks, Daniel Jeske Department of Statistics University of California - Riverside

I'm trying to use GLMM in library(lme4), R 1.9.0pat, updated just now. I get an error message I can't decipher: library(lme4) set.seed(1) n <- 10 N <- 1000 DF <- data.frame(yield=rbinom(n, N, .99)/N, nest=1:n) fit <- GLMM(yield~1, random=~1|nest, family=binomial, data=DF, weights=rep(N, n)) Error in eval(expr, envir, enclos) : ..1 used in an incorrect

I'm having trouble using binomial(link="cloglog") with GLMM in lme4, Version: 0.5-2, Date: 2004/03/11. The example in the Help file works fine, even simplified as follows: fm0 <- GLMM(immun~1, data=guImmun, family=binomial, random=~1|comm) However, for another application, I need binomial(link="cloglog"), and this generates an error for me: >

Hi R-Helpers: I?m trying to use the function GLMM from lme4 package, (R-1.8.1, Windows 98),and I get the following error: > pd5 = GLMM(nplant~sitio+ + fert+ + remo+ + sitio:fert+ + remo:sitio+ + remo:fert+ + remo:fert:sitio + data=datos, + family=binomial, +

Dear list members, In re-running with GLMM() from the lme4 package a generalized-linear mixed model that I had previously fit with glmmPQL() from MASS, I'm getting a warning of a convergence failure, even when I set the method argument of GLMM() to "PQL": > bang.mod.1 <- glmmPQL(contraception ~ as.factor(children) + cage + urban, + random=~as.factor(children) + cage +

Can GLMM take formula derived from another object? foo <- glm (OVEN ~ h + h2, poisson, dataset) # ok bar <- GLMM (OVEN ~ h + h2, poisson, dataset, random = list (yr = ~1)) #error bar <- GLMM (foo$formula, poisson, dataset, random = list (yr = ~1)) #Error in foo$("formula" + yr + 1) : invalid subscript type I am using R2.1.0, lme4 0.8-2, windows xp. Below is a dataset if you

Hi again. Perhaps a simple question this time.... I am analysing data with a dependent variable of insect counts, a fixed effect of site and two random effects, day, which is the same set of 10 days for each site, and then transect, which is nested within site (5 each). I am trying to fit the cross classified model using GLMM in lme4. I have, for potential use, created a second coding

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

This is a summary and extension of the thread "GLMM (lme4) vs. glmmPQL output" http://maths.newcastle.edu.au/~rking/R/help/04/01/0180.html In the new revision (#Version: 0.4-7) of lme4 the standard errors are close to those of the 4 other methods. Thanks to Douglas Bates, Saikat DebRoy for the revision, and to G?ran Brostr?m who run a simulation. In response to my first posting, Prof.

Hello, I have been attempting to run a poisson glmm using lme4 for some time now and have had a lot of trouble. I would say 9 times out of 10 I receive the following warning: CHOLMOD warning: %h Error in mer_finalize(ans) : Cholmod error `not positive definite' at file:../Cholesky/t_cholmod_rowfac.c, line 432 My data are counts of microbe colony forming units (CFUs) collected from

Hi all,   I have data set with repeated measurement ordinal responses and I would like it using generalized linear mixed model.   Shall I use MCMCglmm packeage or I can use lme4 package?   waiting for  answer eagerly, [[alternative HTML version deleted]]

Hello, I'm very sorry for my repeated question, which i asked 2 weeks ago, namely: i'm interested in possibly simple random-part specification in the call of GLMM(...) (from lme4-package) i have a random blocked structure (i.e. ~var.a1+var.a2+var.a3, ~var.b1+var.b2,~var.c1+var.c2+var.c3+var.c4), and each one part of it i would like to model as Identity-structure matrix. So i had,

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

Hi all, First time poster and a relatively new R user, I'm beginning analysis for my masters degree. I'm doing a bit of work on octopus behaviour, and while it's been fascinating, the stats behind it are a bit beyond my grasp at the moment. I was hoping that somebody with more experience my be able to look at my example and offer their wisdom, much to my appreciation :-) At the most

Hello, I have a couple questions regarding generalized linear mixed models specifically around fitting the random effects terms correctly to account for any pseudo-replication. I am reading through and trying to follow examples from Zuur et al. Mixed Effects Models and Extensions in Ecology with R, but am still at bit unsure if I am specifying the models correctly. Background information: Our

Hola, Soy nuevo manejando R y no tengo mucha experiencia. Estoy intentando modelar una función que me relacione el nº de cebas (nº de presas que los padres traen a los pollos) con el tamaño de parche de un bosque (factor categórico; 2 niveles= grande y pequeño). Al ser un conteo (nº de cebas) he pensado utilizar familia= poisson con link= logarítmico. He construido un GLMM con: Nº de cebas

Hi All, I have data about seed predation (SP) in fruits of three differents colors (yellow, motted, dark) and in two fruiting seasons (2007, 2008). I performed a GLMM (lmer function, lme4 package) and the outcome showed that the interaction term (color:season) was significant, and some combinations of this interaction have significant Pr(>|z|), but I don't think they are the right