[R-sig-ME] extracting p values for main effects of binomial glmm

Ken Beath ken.beath at mq.edu.au
Wed Mar 4 23:34:22 CET 2015


Did yo mean to have interactions between all 3 as "Diet * Infection status
* Day". With interactions it isn't possible to test for the effect of main
effects.

On 5 March 2015 at 07:11, Megan Kutzer <makutzer at gmail.com> wrote:

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



-- 

*Ken Beath*
Lecturer
Statistics Department
MACQUARIE UNIVERSITY NSW 2109, Australia

Phone: +61 (0)2 9850 8516

Building E4A, room 526
http://stat.mq.edu.au/our_staff/staff_-_alphabetical/staff/beath,_ken/

CRICOS Provider No 00002J
This message is intended for the addressee named and may...{{dropped:9}}



More information about the R-sig-mixed-models mailing list