[R-sig-ME] extracting p values for main effects of binomial glmm
Doogan, Nathan
doogan.1 at osu.edu
Wed Mar 4 23:09:23 CET 2015
At last check, the summary() of a glmer() object does report p-values. It's output from lmer() that does not. I suppose I could be working with an older version...
-Nate
--
Nathan J. Doogan, Ph.D.
Post-Doctoral Researcher
College of Public Health
The Ohio State University
________________________________________
From: R-sig-mixed-models [r-sig-mixed-models-bounces at r-project.org] on behalf of Megan Kutzer [makutzer at gmail.com]
Sent: Wednesday, March 04, 2015 3:11 PM
To: r-sig-mixed-models at r-project.org
Subject: [R-sig-ME] extracting p values for main effects of binomial glmm
Hi,
I'm fairly new to mixed models and have done a lot of reading without much
success. Unfortunately there is no one at my institution who is really
familiar with them so I thought I would try this list.
I'm running a binomial generalized linear mixed effects model and I need
p-values for the main effects. I know this isn't entirely correct with this
type of model but my supervisor wants the p-values!
The model is:
glmer (Proportion hatched ~ Diet * Infection status * Day + (1|SubjectID) +
(1|Day), family=binomial)
where,
Proportion hatched = cbind(Offspring, Eggs-Offspring)
Diet is a factor with 2 levels
Infection status is a factor with 4 levels
Day is a factor with 3 levels
Using Subject ID number and Day as random effects is supposed to control
for pseudoreplication in the model, although I am not entirely sure that
this is specified in the correct way. I wanted to include experimental
replicate here too but the model failed to converge.
My question is: is there a way to get p-values for the main fixed effects
of Diet, Infection and Day?
If you need more specific model information or the model output I would be
happy to provide it.
Thanks,
Megan
[[alternative HTML version deleted]]
_______________________________________________
R-sig-mixed-models at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
More information about the R-sig-mixed-models
mailing list