similar to: syntax and package for generalized linear mixed models

Has there been any follow up to this question? I have found myself wondering the same thing: How then does SAS fit a beta distributed GLMM? It also fits the negative binomial distribution. Both of these would be useful in glmer/lmer if they aren't 'illegal' as Brian suggested. Especially as SAS indicates a favorable delta BIC of over 1000 when I fit the beta to my data (could be the

Dear list, I am trying to use MuMIn to compare all possible mixed models using the dredge function on binomial data but I am getting an error message that I cannot decode. This error only occurs when I use glmer. When I use an lmer analysis on a different response variable every works great. Example using a simplified glmer model global model: mod<- glmer(cbind(st$X2.REP.LIVE,

Dear list, I am fitting a logistic multi-level regression model and need to test the difference between the ordinary logistic regression from a glm() fit and the mixed effects fit from glmer(), basically I want to do a likelihood ratio test between the two fits. The data are like this: My outcome is a (1,0) for health status, I have several (1,0) dummy variables RURAL, SMOKE, DRINK, EMPLOYED,

Dear all, I found a problem with 'lme4'. Basically, once you load the package 'aod' (Analysis of Overdispersed Data), the functions 'lmer' and 'glmer' don't work anymore: library(lme4) (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) (gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), family = binomial, data

Hello, I need to perform a mixed-model (with nesting) MANCOVA, using Type III sums of squares. I know how to perform each of these types of tests individually, but I am not sure if performing a mixed-model MANCOVA is possible. Please let me know. Erika <>< <>< <>< <>< <>< <>< <>< Erika Crispo, PhD candidate

best described here: https://bugzilla.redhat.com/bugzilla/show_bug.cgi?id=64689 Confirmed also present with the rpm build at http://rsync.samba.org/ftp/rsync/binaries/redhat/rsync-2.4.6-1.i386.rpm Please cc: me on replies (I'm not on the list, yet - my procmailrc's in a major state of flux as I'm switching machines) and/or add comments to the above bugzilla entry James -- James

I am currently running a generalized linear mixed effect model using glmer and I want to estimate how much of the variance is explained by my random factor. summary(glmer(cbind(female,male)~date+(1|dam),family=binomial,data= liz3")) Generalized linear mixed model fit by the Laplace approximation Formula: cbind(female, male) ~ date + (1 | dam) Data: liz3 AIC BIC logLik deviance 241.3

Hello, I have compared the potentials of R and Stata about GLMM, analysing the dataset 'ohio' in the package 'faraway' (the same dataset is analysed with GEE in the book 'Extending the linear model with R' by Julian Faraway). Basically, I've tried the 2 commands 'glmmPQL' and 'glmer' of R and the command 'xtmelogit' of Stata. If I'm not

Hello, I am using glmer() from lmer(lme4) to run generalized linear mixed models. I can't figure out how to extract the z values for the fixed effects that are reported using the summary function . Any help would be appreciated. Thanks, Spencer [[alternative HTML version deleted]]

Hello, I am using glmer() from lmer(lme4) to run generalized linear mixed models. However, I am having a problem extracting the standard errors for the fixed effects. I have used: summary(model)$coef fixed.effects(model) coef(model) to get out the parameter estimates, but do not seem able to extract the se's. Anybody have a solution? Thanks, John

Dear all, I am running glmer() in R. How can I retrieve the estimates of fixed effects and the variance of the random effects from the result? Thank you so much. Joe ___________________________________________________ ±zªº¥Í¬¡§Y®É³q ¡Ð ·¾³q¡B®T¼Ö¡B¥Í¬¡¡B¤u§@¤@¦¸·d©w¡I [[alternative HTML version deleted]]

Hi, I have a dataset that has 2 groups of samples. For each sample, then response measured is the number of success (no.success) obatined with the number of trials (no.trials). So a porportion of success (prpop.success) can be computed as no.success/no.trials. Now the objective is to test if there is a statistical significant difference in the proportion of success between the 2 groups of

Hi, I have the following type of data: 86 subjects in three independent groups (high power vs low power vs control). Each subject solves 8 reasoning problems of two kinds: conflict problems and noconflict problems. I measure accuracy in solving the reasoning problems. To summarize: binary response, 1 within subject var (TYPE), 1 between subject var (POWER). I wanted to fit the following model:

Dear all, This question is basic but I am stumped. After running the below, I receive the message: "non-integer #successes in a binomial glm!" model1 <- glmmML(y~Brood.Size*Density+Date.Placed+Species+Placed.Emerging+Year+rate.of.parperplot, data = data, cluster= data$Patch, family=binomial(link="logit")) My response variable is sex ratio, and I have learned quickly not

Dear list, I have another (again possibly boneheaded) puzzle about importing, again encapsulated in a nearly trivial package. (The package is posted at <http://www.math.mcmaster.ca/bolker/misc/coefsumtest_0.001.tar.gz>.) The package consists (only) of the following S3 method definitions: coeftab <- function(object, ...) UseMethod("coeftab",object) coeftab.default <-

Dear all, In order to be sure I understand glmm correctly, I try to reproduce by hand a simple result. Here is a reproducible code. The questions are in _________________ Of course I have tried to find the solution using internet but I was not able to find a solution. I have also tried to follow glmer but it is very complicated code! Thanks for any help. Marc # Generate set of df with nb

Hi all, I am trying to run a glm with mixed effects. My response variable is number of seedlings emerging; my fixed effects are the tree species and distance from the tree (in two classes - near and far).; my random effect is the individual tree itself (here called Plot). The command I've used is: mod <- glmer(number ~ Species + distance + offset(area) + (1|Plot), family = poisson)

I have a few questions relating to overdispersion in a sex ratio data set that I am working with (note that I already have an analysis with GLMMs for fixed effects, this is just to estimate dispersion). The response variable is binomial because nestlings can only be male or female. I have samples of 1-5 nestlings from each nest (individuals within a nest are not independent, so the response

Dear R-gurus: I have a question about lmer. Basically, I have a dataset, in which each observation records number of trials (N) and number of events (Y) given a covariate combination(X) and group id (grp_id). So, my dataset is in tabular form. (in case my explanation of tabular form is unclear, please see the link:

Sub: using glmer to fit a mixed-effects model with gamma-distributed response variable Hello, I'm currently trying to fit a mixed effects model , i.e.: > burnedmodel1.2<-glmer(gpost.f.crwn.length~lg.shigo.av+dbh+leaf.area+ bark.thick.bh+ht.any+ht.alive+(1|site/transect/plot), family=gaussian, na.action=na.omit, data=rws30.BL) If I run this code, I get the error below: Error: