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

Ben Bolker bbolker at gmail.com
Wed Mar 4 23:46:26 CET 2015


-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

On 15-03-04 05:09 PM, Doogan, Nathan wrote:
> 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
> 

 Yes (and has more or less for forever), although these are based on
Wald Z-scores.  Thanks to the reminder from Jonathan Baron about
?pvalues ...


> 
> 
> -- 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
> 
> _______________________________________________ 
> R-sig-mixed-models at r-project.org mailing list 
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
> 

-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.11 (GNU/Linux)

iQEcBAEBAgAGBQJU94tCAAoJEOCV5YRblxUHnPcIAKgay6VpUxGahQfi7JMsj43Z
iDPdR2d97ANFLqOXaZJpGktquGKtFbYGmua9XSyz4DSdl5bpX9n2hb+NYT4Secdo
1cVc28noVbP83Od6wHdJgyRWmzJ85onZ+m0t8CK2BtnJO4C/1GMUZwX9H4+F8g6N
ljbIaAD07AOmx14VF8l4XpW9guw3NNEOgKcupZBJ5K7zsLkkdtnLoAJZFCcxpypg
WdGcrcft82nJorqwcEYrjpqF/X9CzFvRIPAhf+cbZbDe06W98coM4RVNyLRXrLNE
2NKOKK2Pg5tHVAW6PnlaOBBuSCFgpSeWTQINDLda/Q3j5OUqmoCeLjpjjTb3A30=
=iWlw
-----END PGP SIGNATURE-----



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