Generalized Linear Mixed (effects) Models are typically used for modeling non-independent non-normal data (eg, longitudinal binary data).

learn more… | top users | synonyms

1
vote
1answer
27 views

code needed for p values in GLMM

I need to get p values for the fixed effects in the following GLMM's I ran. Does anyone know of code that I can run that will give me the p values I need? At the moment the output from the ANOVA only ...
1
vote
0answers
51 views

Reading the output from a GLMM run in R

I am a complete novice and dummy when it comes to statistics so I apologise in advance... I have been asked to report the results of my GLMMs (I ran two) in a table. This table must state: effect, ...
0
votes
0answers
4 views

Error using the step-function with glmmML [migrated]

When I tried to use the step function I receive this error: "Error in if (all(is.finite(c(n0, nnew))) && nnew != n0) stop("number of rows in use has changed: remove missing values?") : ...
1
vote
0answers
19 views

How to handle underdispersion in GLMM (binomial outcome variable)

I'm working on the following model in R: ...
0
votes
0answers
24 views

plot effects on glmer

I done a GLMM with the function glmer form lme4, the variable paternity is binomial (yes or no)(model<-glmer(paternity~factor1, factor 2 ....,family=binomial). My model run correctly, but after I ...
0
votes
0answers
23 views

Including a nested factor as random effect in a GLMM

Good afternoon, I'd like to ask for advice on including a nested factor as a random effect in a GLMM. I've read other threads in this forum, but still am not able to answer my question. Any help is ...
0
votes
0answers
9 views

Does GLMM or LMM already give us planned contrasts?

I use GLMM and LMM to find significant factors for my data. And could I ask if the results from these two models already give me planned contrasts? Or do I need to run them separately?
1
vote
1answer
47 views

Binomial GLMM: Model validation & ceiling effect

My data has a binary response acc(correct/incorrect), one continuous predictor score, three categorical predictors (...
0
votes
1answer
27 views

GLMM: Many models?

I use R to run GLMM models in my study. Then in my situation, I have 10 independent variables and one of them has 3 categories. Then, when I run the GLMM, it usually uses one category as baseline, so ...
1
vote
0answers
28 views

Prediction Model for the number of defects

I work at a software company, and we are trying to build a statistical model to predict the defect density (the number of defects that exist in a certain software release that is delivered to a ...
1
vote
1answer
85 views

Binomial GLMM with categorical predictors: p-values?

My data has a binary response (correct/incorrect), one continuous predictor score, three categorical predictors (race, ...
1
vote
0answers
32 views

Binomial mixed model with categorical predictors: model selection and getting p-values [closed]

My data has a binary response (correct/incorrect), one continuous predictor (with NaNs) and several categorical predictors. I want to add a random intercept for a ...
1
vote
0answers
54 views

Logistic regression with binary inputs

Is logistic regression an appropriate classifier when the input data are binary? Say we are conducting an experiment where the subject is presented with blue and green circles of varying shades, and ...
1
vote
0answers
34 views

How to assess whether IQ predicts emotional experiment more at time 1 or time 2 controlling for covariates?

I have the following data set of a two-wave longitudinal study in educational psychology. A group of 400 people did some questionnaire on IQ, personality and two different emotional states (sadness ...
0
votes
1answer
33 views

Generalized linear mixed model: what is base dependent variable?

I have run Generalized linear mixed model with glmer in lme4. I use R version 3.0.1. My dependent variable is binary (correct or wrong). And this is my results: ...
4
votes
0answers
76 views

Meaning of a convergence warning in glmer

I am using the glmer function from the lme4 package in R, and I'm using the bobyqa optimizer ...
1
vote
0answers
23 views

Calculating point estimates from model-averaged parameters

I'm using an IT-approach and multi-model inference with some count data. I have used model averaging to obtain parameter estimates from several GLMMs with Poisson-lognormal errors (Poisson family ...
4
votes
2answers
94 views

extremely left-skewed response variable - how do I model this dataset?

This is a histogram showing my response variable. The response is # (or proportion? or percent?) of aphids eaten off of cards in fields, to model predation by natural enemies. Predictors: fixed ...
0
votes
0answers
47 views

How to interpret the results from random coefficient Poisson data analysis

I have fitted random coefficient Poisson analysis in R. I have obtained the following results: Generalized linear mixed model fit by maximum likelihood ['glmerMod'] Family: poisson ( log ) Formula: ...
0
votes
0answers
36 views

repeated measure ANOVA or GLMM?

I would like to ask for some advice on an analysis that I am doing at the moment. The experiment: I ran an experiment with animals in which I subjected 15 pairs of individuals (one male and one ...
5
votes
1answer
134 views

Marginal model versus random-effects model – how to choose between them? An advice for a layman

In searching for any info about marginal model and random-effects model, and how to choose between them, I have found some info but it was more-or-less mathematical abstract explanation (like for ...
0
votes
0answers
33 views

Is it possible to fit a multilevel Weibull/generalised gamma survival model?

I am imagining a situation in which I have individual patient survival data from a number of clinical trials of a particular drug. The drug can be used for a wide variety of cancers which have ...
2
votes
0answers
30 views

Logistic regression with block correction for ordinal response

I have a categorical response (ordinal, lets call it A) collected from a complete randomized block design experiment where the explanatory variable is treatment. How do I account for the between block ...
2
votes
1answer
155 views

Specifying a multilevel model in MCMCglmm (R), that is heteroskedastic at level one

I am considering MCMCglmm as an alternative to MLwiN. The former package works perfectly fine, but I cannot figure out how to model heteroskedasticity at level one. For instance, if I have the ...
0
votes
0answers
28 views

GLMM multilevel hierarchical model with R [duplicate]

I hope I will get your precious help for my R problem. I want to study the classroom, the teacher and the school effect/result on the pupil's success (or not) at school. I also want to know the ...
4
votes
2answers
63 views

How to analyse spatial data where the depending variable is binary

I have to test which factors influence game damage in fields. I mapped areas with damage and those without. It was not always possible to map 100% of a field, so there are also areas where it is ...
0
votes
1answer
48 views

Null GLMM Poisson underestimate the mean of the response variable. Is it indicative of poor fitting?

I want to test the fixed and random effects of some covariates on a discrete variable with non negative values. In exploratory analysis I fitted a null Poisson GLM and an null Poisson GLMM. However, ...
1
vote
1answer
86 views

GLMM multilevel (hierarchical) model

I want to study the classroom and the school effect/result on the pupil's success (or not) at school. I also want to know the effect of the age, the gender (of the pupils) and if the pupils have or ...
1
vote
1answer
134 views

Level-2 predictions with lme4/glmer model

Let's say I've fitted a 2 level model with glmer like this: ...
2
votes
0answers
33 views

Bias correction methods in G-side GLMMs

Are there any conventions for the use of bias correction methods in G-side generalized linear mixed models (GLMMs)? For linear mixed models I usually use the Kenward and Roger adjustment, but it is ...
1
vote
0answers
107 views

How to calculate confidence bands for a regression line, given intercept, slope, and standard errors

I have the following set of model-averaged fixed effects from a set of binomial GLMMs: I would like to plot the predicted effect of "NBT", along with confidence bands, while holding all the other ...
1
vote
0answers
31 views

How to specify a repeated-measure random effect, and a nested plot-in-stand in a glmm.admb nbinom?

I'm not sure my glmmadmb syntax correctly specifies: that plots(2,184) are grouped or blocked within stands(24-26 per period), and there are repeated measures -- 3 winters and 3 summers on the ...
1
vote
0answers
30 views

How can I get a “whole treatment effect” for a 3 level categorical variable in a GLMM (glmmPQL in R)

I work with habitat use data (summarized binomial data: visited / did not visit) which I fitted with a mixed model of the binomial family with a random factor accounting for repeated sampling of the ...
1
vote
0answers
65 views

SPSS GLMM Data Structure setup advice

I am running an analysis looking at individual voter turnout where it is a 1 for voting and 0 for not voting. I have some information on individuals including age, gender, and several years of past ...
1
vote
0answers
96 views

Random effect with 0 variance in GLMM

In this study, subjects are measured continuously over the day via electrocardiography (ECG). During the day, certain trigger events occur randomly. Once all the data is collected, the trigger events ...
1
vote
1answer
605 views

Mixed effects model validation and selection with `lme4::glmer`

If I had a glm using on count data I may do the following: glm(response ~ exp1 * exp2, family = poisson, data =data) The ...
3
votes
1answer
142 views

glmm unbalanced nested random factors

imagine you have the following data structure ...
0
votes
0answers
10 views

How to obtain aggregate estimate of a factor with various levels in General Linear Mixed Models

I run general linear mixed model with factor 'Substitutability' with seven levels in SPSS. It gave six level estimates. I want to calculate aggregate estimate for 'substitutability'. How to get? ...
2
votes
0answers
31 views

Should I take the arithmetic mean of variates which vary on smaller time scales in a linear regression?

I am consulting with another student with a stats background, but am pretty sure he his saying something that contradicts what I have been told about my data. I have a daily response variable ...
1
vote
1answer
83 views

What is error structure?

Simple question but I can't find the answer anywhere online or in any of my texts. I think people assume you know, but I don't. What is an error structure, specifically for selecting the family ...
0
votes
0answers
49 views

Dramatic changes in the results from glmm with slight changes in dataset

I am running a GLMM in R using the glmmPQL code from the MASS package. I have been getting a convergence error from the following code. ...
1
vote
0answers
85 views

Mathematical equation for a generalized linear mixed model with interactions

I'd like to know how following R-code for a binomial GLMM (package lme4) translates into a proper mathematical notation form? ...
6
votes
0answers
154 views

Generalized Linear Mixed Models: Diagnostics

I have a random intercept logistic regression (due to repeated measurements) and I would like to do some diagnostics, specifically concerning outliers and influential observations. I looked at ...
1
vote
1answer
221 views

Generalized linear mixed models: model selection

This question/topic came up in a discussion with a colleague and I was looking for some opinions on this: I am modeling some data using a random effects logistic regression, more precisely a random ...
1
vote
0answers
117 views

Random slopes and intercepts using glmmadmb in R - negative binomial

I am new to mixed modeling and using the glmmadmb package. I am interested in using glmmadmb to model count data with a negative binomial distribution, using year and site sampled as my random effect ...
6
votes
1answer
317 views

What are Hommel Hochberg corrections?

I have recently been introduced to to Hommel Hochberg corrections. I am trying to find a simple explanation about what this actually is/does, but am having no luck. Can anyone please give a brief and ...
4
votes
1answer
136 views

How to calculate estimated proportions and their confidence intervals from a mixed model?

I have an experiment with two treatments. It is a split plot experiment, with the structure Block/Treatment1/Treatment2. Each treatment has 2 levels. The dependent variable is presence/absence data ...
5
votes
3answers
577 views

How to assess the fit of a binomial GLMM fitted with lme4 (> 1.0)?

I have a GLMM with a binomial distribution and a logit link function and I have the feeling that an important aspect of the data is not well represented in the model. To test this, I would like to ...
0
votes
0answers
89 views

Correct estimation of arguments for glmmLasso function

I am using glmmLasso for variable selection. In my case, n is slightly less than p and ...
3
votes
0answers
96 views

Plotting fit for binomial lme

I've been asked by a reviewer on a manuscript to provide plots of a model fit for a binomial lme which is specified as follows: ...